File: Atensor.texi

package info (click to toggle)
maxima 5.21.1-2squeeze
  • links: PTS
  • area: main
  • in suites: squeeze
  • size: 94,928 kB
  • ctags: 43,849
  • sloc: lisp: 298,974; fortran: 14,666; perl: 14,325; tcl: 10,494; sh: 4,052; makefile: 2,975; ansic: 471; awk: 24; sed: 7
file content (265 lines) | stat: -rw-r--r-- 11,006 bytes parent folder | download | duplicates (4)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
@c /Atensor.texi/1.8/Mon Nov 21 00:19:56 2005//
@menu
* Introdu@value{cedilha}@~ao ao Pacote atensor::
* Defini@value{cedilha}@~oes para o Pacote atensor::
@end menu

@node Introdu@value{cedilha}@~ao ao Pacote atensor, Defini@value{cedilha}@~oes para o Pacote atensor, Pacote atensor, Pacote atensor
@section Introdu@value{cedilha}@~ao ao Pacote atensor

@code{atensor} @'e um pacote de manipul@value{cedilha}@~ao de tensores alg@'ebricos.  Para usar @code{atensor},
digite @code{load(atensor)}, seguido por uma chamada @`a fun@value{cedilha}@~ao 
@code{init_atensor}.

A ess@^encia de @code{atensor} @'e um conjunto de regras de simplifica@value{cedilha}@~ao para o operador
de produto (ponto) n@~ao comutativo ("@code{.}").  @code{atensor} reconhece
muitos tipos de @'algebra; as regras de simplifica@value{cedilha}@~ao correspondentes s@~ao activadas quando
a fun@value{cedilha}@~ao @code{init_atensor} @'e chamada.

A compatibilidade de @code{atensor} pode ser demonstrada pela
defini@value{cedilha}@~ao da @'algebra de quaterni@~oes como uma
@'algebra de Clifford Cl(0,2) com dois vectores fundamentais.  As tr@^es
unidades quaterni@'onicas imagin@'arias fundamentais s@~ao ent@~ao os
dois vectores base e seu produto, i.e.:

@example
    i = v     j = v     k = v  .  v
         1         2         1    2
@end example

Embora o pacote @code{atensor} tenha uma defini@value{cedilha}@~ao
interna para a @'algebra dos quaterni@~oes, isso n@~ao foi usado nesse
exemplo, no qual n@'os nos esfor@value{cedilha}amos para construir a
tabela de multiplica@value{cedilha}@~ao dos quaterni@~oes como uma
matriz:

@example

(%i1) load(atensor);
(%o1)       /share/tensor/atensor.mac
(%i2) init_atensor(clifford,0,0,2);
(%o2)                                done
(%i3) atensimp(v[1].v[1]);
(%o3)                                 - 1
(%i4) atensimp((v[1].v[2]).(v[1].v[2]));
(%o4)                                 - 1
(%i5) q:zeromatrix(4,4);
                                [ 0  0  0  0 ]
                                [            ]
                                [ 0  0  0  0 ]
(%o5)                           [            ]
                                [ 0  0  0  0 ]
                                [            ]
                                [ 0  0  0  0 ]
(%i6) q[1,1]:1;
(%o6)                                  1
(%i7) for i thru adim do q[1,i+1]:q[i+1,1]:v[i];
(%o7)                                done
(%i8) q[1,4]:q[4,1]:v[1].v[2];
(%o8)                               v  .  v
                                     1    2
(%i9) for i from 2 thru 4 do for j from 2 thru 4 do
      q[i,j]:atensimp(q[i,1].q[1,j]);
(%o9)                                done
(%i10) q;
                   [    1        v         v      v  .  v  ]
                   [              1         2      1    2 ]
                   [                                      ]
                   [   v         - 1     v  .  v    - v    ]
                   [    1                 1    2      2   ]
(%o10)             [                                      ]
                   [   v      - v  .  v     - 1      v     ]
                   [    2        1    2              1    ]
                   [                                      ]
                   [ v  .  v      v        - v       - 1   ]
                   [  1    2      2          1            ]
@end example

@code{atensor} reconhece como bases vectoriais s@'{@dotless{i}}mbolos indexados, onde o s@'{@dotless{i}}mbolo 
@'e aquele armazenado em @code{asymbol} e o i@'{@dotless{i}}ndice est@'a entre 1 e @code{adim}.
Para s@'{@dotless{i}}mbolos indexado, e somente para s@'{@dotless{i}}mbolos indexados, as formas bilineares
@code{sf}, @code{af}, e @code{av} s@~ao avaliadas.  A avalia@value{cedilha}@~ao
substitui os valores  de @code{aform[i,j]} em lugar de @code{fun(v[i],v[j])}
onde @code{v} representa o valor de @code{asymbol} e @code{fun} @'e
ainda @code{af} ou @code{sf}; ou, isso substitui @code{v[aform[i,j]]}
em lugar de @code{av(v[i],v[j])}.

Desnecess@'ario dizer, as fun@value{cedilha}@~oes @code{sf}, @code{af} e @code{av}
podem ser redefinidas.

Quando o pacote @code{atensor} @'e chamado, os seguintes sinalizadores s@~ao configurados:

@example
dotscrules:true;
dotdistrib:true;
dotexptsimp:false;
@end example

Se quiser experimentar com uma @'algebra n@~ao associativa, pode
tamb@'em considerar a configura@value{cedilha}@~ao de @code{dotassoc}
para @code{false}.  Nesse caso, todavia, @code{atensimp} n@~ao star@'a
sempre habilitado a obter as simplifica@value{cedilha}@~oes desejadas.


@c end concepts atensor
@node Defini@value{cedilha}@~oes para o Pacote atensor,  , Introdu@value{cedilha}@~ao ao Pacote atensor, Pacote atensor

@section Defini@value{cedilha}@~oes para o Pacote atensor

@deffn {Fun@value{cedilha}@~ao} init_atensor (@var{alg_type}, @var{opt_dims})
@deffnx {Fun@value{cedilha}@~ao} init_atensor (@var{alg_type})

Inicializa o pacote @code{atensor} com o tipo especificado de @'algebra.  @var{alg_type}
pode ser um dos seguintes:

@code{universal}: A @'algebra universal tendo regras n@~ao comutativas.

@code{grassmann}: A @'algebra de Grassman @'e definida pela rela@value{cedilha}@~ao de 
comuta@value{cedilha}@~ao @code{u.v+v.u=0}.

@code{clifford}: A @'algebra de Clifford @'e definida pela rela@value{cedilha}@~ao
de comuta@value{cedilha}@~ao @code{u.v+v.u=-2*sf(u,v)} onde @code{sf} @'e a fun@value{cedilha}@~ao
valor-escalar sim@'etrico.  Para essa @'algebra, @var{opt_dims} pode ser acima de tr@^es 
inteiros n@~ao negativos, representando o n@'umero de dimens@~oes positivas,
dimens@~oes degeneradas, e dimens@~oes negativas da @'algebra, respectivamente.  Se
quaisquer valores @var{opt_dims} s@~ao fornecidos, @code{atensor} ir@'a configurar os
valores de @code{adim} e @code{aform} apropriadamente.  Caso contr@'ario,
@code{adim} ir@'a por padr@~ao para 0 e @code{aform} n@~ao ser@'a definida.

@code{symmetric}: A @'algebra sim@'etrica @'e definida pela rela@value{cedilha}@~ao de 
comuta@value{cedilha}@~ao @code{u.v-v.u=0}.

@code{symplectic}: A @'algebra simpl@'etica @'e definida pela rela@value{cedilha}@~ao de 
comuta@value{cedilha}@~ao @code{u.v-v.u=2*af(u,v)} onde @code{af} @'e uma fun@value{cedilha}@~ao valor-escalar 
antisim@'etrica.  Para a @'algebra simpl@'etica, @var{opt_dims} pode
mais de dois inteiros n@~ao negativos, representando a dimens@~ao n@~ao degenerada e
e a dimens@~ao degenerada, respectivamente.  Se quaisquer valores @var{opt_dims} s@~ao
fornecidos, @code{atensor} ir@'a configurar os valores de @code{adim} e @code{aform}
apropriadamente.  Caso contr@'ario, @code{adim} ir@'a por padr@~ao para 0 e @code{aform}
n@~ao ser@'a definida.

@code{lie_envelop}: O inv@'olucro da @'algebra de Lie @'e definido pela 
rela@value{cedilha}@~ao de comuta@value{cedilha}@~ao @code{u.v-v.u=2*av(u,v)} onde @code{av} @'e
uma fun@value{cedilha}@~ao antisim@'etrica.

A fun@value{cedilha}@~ao @code{init_atensor} tamb@'em reconhece muitos tipos pr@'e-definidos de 
@'algebra:

@code{complex} implementa a @'algebra de n@'umeros complexos como a
@'algebra de Clifford Cl(0,1).  A chamada @code{init_atensor(complex)} @'e
equivalente a @code{init_atensor(clifford,0,0,1)}.

@code{quaternion} implementa a @'algebra de quaterni@~oes.  A chamada
@code{init_atensor(quaternion)} @'e equivalente a 
@code{init_atensor(clifford,0,0,2)}.

@code{pauli} implementa a @'algebra de spinores de Pauli como a
@'algebra de Clifford Cl(3,0).  Uma chamada a @code{init_atensor(pauli)}
@'e equivalente a @code{init_atensor(clifford,3)}.

@code{dirac} implementa a @'algebra de spinores de Dirac como a
@'algebra de Clifford Cl(3,1).  Uma chamada a @code{init_atensor(dirac)}
@'e equivalente a @code{init_atensor(clifford,3,0,1)}.

@end deffn


@deffn {Fun@value{cedilha}@~ao} atensimp (@var{expr})

Simplifica a express@~ao alg@'ebrica de tensores @var{expr} conforme as
regras configuradas por uma chamada a @code{init_atensor}.
Simplifica@value{cedilha}@~oes incluem aplica@value{cedilha}@~ao
recursiva de rela@value{cedilha}@~oes comutativas e
resolu@value{cedilha}@~oes de chamadas a @code{sf}, @code{af}, e
@code{av} onde for aplic@'avel.  Uma salvaguarda @'e usada para garantir
que a fun@value{cedilha}@~ao sempre termine, mesmo para express@~oes
complexas.

@end deffn

@deffn {Fun@value{cedilha}@~ao} alg_type

O tipo de @'algebra.  Valores v@'alidos s@'ao @code{universal}, @code{grassmann},
@code{clifford}, @code{symmetric}, @code{symplectic} e @code{lie_envelop}.

@end deffn

@defvr {Vari@'avel} adim

A dimensionalidade da @'algebra.  @code{atensor} usa o valor de @code{adim}
para determinar se um objecto indexado @'e uma base vectorial v@'alida. Veja @code{abasep}.

@end defvr

@defvr {Vari@'avel} aform

Valor por omiss@~ao para as formas bilineares @code{sf}, @code{af}, e
@code{av}.  O padr@~ao @'e a matriz identidade @code{ident(3)}.

@end defvr

@defvr {Vari@'avel} asymbol

O s@'{@dotless{i}}mbolo para bases vectoriais.

@end defvr

@deffn {Fun@value{cedilha}@~ao} sf (@var{u}, @var{v})

@'E uma fun@value{cedilha}@~ao escalar sim@'etrica que @'e usada em
rela@value{cedilha}@~oes comutativas.  A implementa@value{cedilha}@~ao
padr@~ao verifica se ambos os argumentos s@~ao bases vectoriais usando
@code{abasep} e se esse for o caso, substitui o valor correspondente da
matriz @code{aform}.

@end deffn

@deffn {Fun@value{cedilha}@~ao} af (@var{u}, @var{v})

@'E uma fun@value{cedilha}@~ao escalar antisim@'etrica que @'e usada em rela@value{cedilha}@~oes comutativas.
A implementa@value{cedilha}@~ao padr@~ao verifica se ambos os argumentos s@~ao bases vectoriais
usando @code{abasep} e se esse for o caso, substitui o
valor correspondente da matriz @code{aform}.

@end deffn

@deffn {Fun@value{cedilha}@~ao} av (@var{u}, @var{v})

@'E uma fun@value{cedilha}@~ao antisim@'etrica que @'e usada em rela@value{cedilha}@~oes comutativas.
A implementa@value{cedilha}@~ao padr@~ao verifica se ambos os argumentos s@~ao bases vectoriais
usando @code{abasep} e se esse for o caso, substitui o
valor correspondente da matriz @code{aform}.

Por exemplo:

@example
(%i1) load(atensor);
(%o1)       /share/tensor/atensor.mac
(%i2) adim:3;
(%o2)                                  3
(%i3) aform:matrix([0,3,-2],[-3,0,1],[2,-1,0]);
                               [  0    3   - 2 ]
                               [               ]
(%o3)                          [ - 3   0    1  ]
                               [               ]
                               [  2   - 1   0  ]
(%i4) asymbol:x;
(%o4)                                  x
(%i5) av(x[1],x[2]);
(%o5)                                 x
                                       3
@end example

@end deffn


@deffn {Fun@value{cedilha}@~ao} abasep (@var{v})

Verifica se esse argumento @'e uma base vectorial @code{atensor} .  

E ser@'a, se ele for
um s@'{@dotless{i}}mbolo indexado, com o s@'{@dotless{i}}mbolo sendo o mesmo que o valor de
@code{asymbol}, e o @'{@dotless{i}}ndice tiver o mesmo valor num@'erico entre 1
e @code{adim}.

@end deffn