-
-
Notifications
You must be signed in to change notification settings - Fork 5
/
Copy pathclifford_sta.dem
59 lines (42 loc) · 1.46 KB
/
clifford_sta.dem
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
/* Copyright (C) 2015 Dimiter Prodanov
*
* This program is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public License as
* published by the Free Software Foundation; either version 2 of
* the License, or (at your option) any later version.
*
* This program is distributed in the hope that it will be
* useful, but WITHOUT ANY WARRANTY; without even the implied
* warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
* PURPOSE. See the GNU General Public License for more details.
*
*/
if get('clifford,'version)=false then load("clifford");
("Clifford implements Clifford algebra for Maxima.")$
("Space-time algebra")$
clifford(g,1,3);
("here are all elements of the algebra")$
elements();
("anticommutativity")$
g[2].g[1];
("Electric field vector")$
EE:cvect(E,[[],x,y,z]);
("Magnetic field vector")$
BB:cvect(B,[[],x,y,z]);
("Faraday multi-vector")$
F: EE + %iv.BB;
("Faraday multi-vector full simplification")$
F: F,expand, dotsimpc;
("pseudo-Euclidaen vector segment")$
mtable1([g[1],g[2],g[3],g[4]]);
("Quaternion segment")$
mtable1([g[1] . g[3] . g[4], g[1] . g[2] . g[4],g[1] . g[2] . g[3]]);
("space-time split")$
F1:F.g[1],expand,dotsimpc;
("Quaternion segment")$
mtable1([ g[3] . g[4], g[2] . g[4], g[2] . g[3]]);
("Euclidaen vector segment")$
mtable1([ g[1] . g[2], g[1] . g[3], g[1] . g[4]]);
("Full algebra")$
mtable2();
/* End of demo -- comment line needed by MAXIMA to resume demo menu */