File: poisson.mac

package info (click to toggle)
maxima 5.27.0-3
  • links: PTS
  • area: main
  • in suites: wheezy
  • size: 120,648 kB
  • sloc: lisp: 322,503; fortran: 14,666; perl: 14,343; tcl: 11,031; sh: 4,146; makefile: 2,047; ansic: 471; awk: 24; sed: 10
file content (22 lines) | stat: -rw-r--r-- 704 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
 
/* written by Furuya <go.maxima@gmail.com> 
#  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.                
*/

/* need to do load("diag") */
/* poisson braket operator */

p_braket(_f,_g):= block([_x,_a,_b,_j,_jj:[],basis ],
	mode_declare([_x,_a,_b],any),
	basis:[],
	for i:1 thru dim do (
	basis:endcons(basis[i],basis)
	),
	_x:coefmatrix([d(_f)],basis),
	_j:matrix([0,1],[-1,0]),
	for i:1 thru dim/2 do (
	_jj:cons(_j,_jj)),
	_b:(_x.diag(_jj)).basis,
	inner(d(_g),_b));