Skip to content

Commit

Permalink
Courant Snyder (canonization) done
Browse files Browse the repository at this point in the history
  • Loading branch information
@mattsignorelli
mattsignorelli committed Jun 25, 2024
1 parent 79a0983 commit 46cbfb4
Show file tree
Hide file tree
Showing 6 changed files with 845 additions and 26 deletions.
41 changes: 41 additions & 0 deletions fpp-ptc-sandbox/code/Makefile
View file
Original file line number Diff line number Diff line change
Expand Up @@ -389,6 +389,19 @@ z_radiation_matt_fake_maps-exe/fast:
$(MAKE) $(MAKESILENT) -f CMakeFiles/z_radiation_matt_fake_maps-exe.dir/build.make CMakeFiles/z_radiation_matt_fake_maps-exe.dir/build
.PHONY : z_radiation_matt_fake_maps-exe/fast

#=============================================================================
# Target rules for targets named z_canonize_matt_fake_maps-exe

# Build rule for target.
z_canonize_matt_fake_maps-exe: cmake_check_build_system
$(MAKE) $(MAKESILENT) -f CMakeFiles/Makefile2 z_canonize_matt_fake_maps-exe
.PHONY : z_canonize_matt_fake_maps-exe

# fast build rule for target.
z_canonize_matt_fake_maps-exe/fast:
$(MAKE) $(MAKESILENT) -f CMakeFiles/z_canonize_matt_fake_maps-exe.dir/build.make CMakeFiles/z_canonize_matt_fake_maps-exe.dir/build
.PHONY : z_canonize_matt_fake_maps-exe/fast

normal_simple.o: normal_simple.f90.o
.PHONY : normal_simple.o

Expand All @@ -413,6 +426,30 @@ normal_simple.f90.s:
$(MAKE) $(MAKESILENT) -f CMakeFiles/normal_simple-exe.dir/build.make CMakeFiles/normal_simple-exe.dir/normal_simple.f90.s
.PHONY : normal_simple.f90.s

z_canonize_matt_fake_maps.o: z_canonize_matt_fake_maps.f90.o
.PHONY : z_canonize_matt_fake_maps.o

# target to build an object file
z_canonize_matt_fake_maps.f90.o:
$(MAKE) $(MAKESILENT) -f CMakeFiles/z_canonize_matt_fake_maps-exe.dir/build.make CMakeFiles/z_canonize_matt_fake_maps-exe.dir/z_canonize_matt_fake_maps.f90.o
.PHONY : z_canonize_matt_fake_maps.f90.o

z_canonize_matt_fake_maps.i: z_canonize_matt_fake_maps.f90.i
.PHONY : z_canonize_matt_fake_maps.i

# target to preprocess a source file
z_canonize_matt_fake_maps.f90.i:
$(MAKE) $(MAKESILENT) -f CMakeFiles/z_canonize_matt_fake_maps-exe.dir/build.make CMakeFiles/z_canonize_matt_fake_maps-exe.dir/z_canonize_matt_fake_maps.f90.i
.PHONY : z_canonize_matt_fake_maps.f90.i

z_canonize_matt_fake_maps.s: z_canonize_matt_fake_maps.f90.s
.PHONY : z_canonize_matt_fake_maps.s

# target to generate assembly for a file
z_canonize_matt_fake_maps.f90.s:
$(MAKE) $(MAKESILENT) -f CMakeFiles/z_canonize_matt_fake_maps-exe.dir/build.make CMakeFiles/z_canonize_matt_fake_maps-exe.dir/z_canonize_matt_fake_maps.f90.s
.PHONY : z_canonize_matt_fake_maps.f90.s

z_coast_matt_fake_maps.o: z_coast_matt_fake_maps.f90.o
.PHONY : z_coast_matt_fake_maps.o

Expand Down Expand Up @@ -902,6 +939,7 @@ help:
@echo "... edit_cache"
@echo "... rebuild_cache"
@echo "... normal_simple-exe"
@echo "... z_canonize_matt_fake_maps-exe"
@echo "... z_coast_matt_fake_maps-exe"
@echo "... z_cut-exe"
@echo "... z_cut_sub-exe"
Expand All @@ -925,6 +963,9 @@ help:
@echo "... normal_simple.o"
@echo "... normal_simple.i"
@echo "... normal_simple.s"
@echo "... z_canonize_matt_fake_maps.o"
@echo "... z_canonize_matt_fake_maps.i"
@echo "... z_canonize_matt_fake_maps.s"
@echo "... z_coast_matt_fake_maps.o"
@echo "... z_coast_matt_fake_maps.i"
@echo "... z_coast_matt_fake_maps.s"
Expand Down
11 changes: 7 additions & 4 deletions fpp-ptc-sandbox/code/z_canonize_matt_fake_maps.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -159,7 +159,8 @@ program example
do i=1,c_%nd2
m%v(i)=m%v(i) !*decrement(i)
enddo
! call print(m)
call print(m)
stop

lielib_print(4)=0
call c_normal(m,normal,dospin=putspin,phase=phase,nu_spin=nu_spin)
Expand All @@ -168,7 +169,7 @@ program example
call c_fast_canonise_clean_julia(normal%atot,a1,dospin=putspin)

call print(a1)

stop

normal%atot=a1
write(6,*) normal%damping(1:3)
Expand Down Expand Up @@ -1259,6 +1260,8 @@ subroutine c_fast_canonise_clean_julia(u,u_c,phase,damping,q_cs,q_as,q_orb,q_rot
type(c_damap) uct
integer ndt,ndptb,ndpt,ndct

! nd = number dimensions
! ndt = number harmonic oscillators = 4 if coastings

ndt=c_%nd2t/2
ndptb=c_%ndptb
Expand All @@ -1284,8 +1287,8 @@ subroutine c_fast_canonise_clean_julia(u,u_c,phase,damping,q_cs,q_as,q_orb,q_rot
b0=0
id=0
do i=1,nd
b0(2*i-1,2*i-1)=1
b0(2*i,2*i)=1
b0(2*i-1,2*i-1)=1 ! I
b0(2*i,2*i)=1 !
s(2*i-1,2*i)=1
s(2*i,2*i-1)=-1
enddo
Expand Down
124 changes: 124 additions & 0 deletions src/canonize/ac.map
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,124 @@
6 Dimensional TPSA/DA map

1, NO = 3, NV = 6, INA = 30
*********************************************

I COEFFICIENT ORDER EXPONENTS
NO = 3 NV = 6
1 1.000312298776967 0.000000000000000 1 0 0 0 0 0
1 -0.4830800052522120E-16 0.000000000000000 0 1 0 0 0 0
1 0.1078435759919717E-01 0.000000000000000 0 0 1 0 0 0
1 -0.1979756005464457E-02 0.000000000000000 0 0 0 1 0 0
1 0.5630792227601156E-02 0.000000000000000 0 0 0 0 0 1
-5 0.000000000000000 0.000000000000000 0 0 0 0 0 0

1, NO = 3, NV = 6, INA = 29
*********************************************

I COEFFICIENT ORDER EXPONENTS
NO = 3 NV = 6
1 0.1306080210285989E-02 0.000000000000000 1 0 0 0 0 0
1 0.9995715517477228 0.000000000000000 0 1 0 0 0 0
1 0.4917593339866650E-02 0.000000000000000 0 0 1 0 0 0
1 0.9879832270829701E-02 0.000000000000000 0 0 0 1 0 0
1 0.4464348946484611E-02 0.000000000000000 0 0 0 0 0 1
-5 0.000000000000000 0.000000000000000 0 0 0 0 0 0

1, NO = 3, NV = 6, INA = 28
*********************************************

I COEFFICIENT ORDER EXPONENTS
NO = 3 NV = 6
1 -0.9852382729493981E-02 0.000000000000000 1 0 0 0 0 0
1 -0.1972277583719225E-02 0.000000000000000 0 1 0 0 0 0
1 0.9965336729127154 0.000000000000000 0 0 1 0 0 0
1 -0.3263630881646084E-16 0.000000000000000 0 0 0 1 0 0
1 0.1239780734022430E-01 0.000000000000000 0 0 0 0 0 1
-5 0.000000000000000 0.000000000000000 0 0 0 0 0 0

1, NO = 3, NV = 6, INA = 27
*********************************************

I COEFFICIENT ORDER EXPONENTS
NO = 3 NV = 6
1 0.4807753199428283E-02 0.000000000000000 1 0 0 0 0 0
1 -0.1084012449016447E-01 0.000000000000000 0 1 0 0 0 0
1 0.1156632979614590E-01 0.000000000000000 0 0 1 0 0 0
1 1.003361696547912 0.000000000000000 0 0 0 1 0 0
1 0.9355219167887742E-02 0.000000000000000 0 0 0 0 0 1
-5 0.000000000000000 0.000000000000000 0 0 0 0 0 0

1, NO = 3, NV = 6, INA = 26
*********************************************

I COEFFICIENT ORDER EXPONENTS
NO = 3 NV = 6
1 -0.4306612093237661E-02 0.000000000000000 1 0 0 0 0 0
1 0.5512437038594799E-02 0.000000000000000 0 1 0 0 0 0
1 -0.9199848978962361E-02 0.000000000000000 0 0 1 0 0 0
1 0.1250395461075949E-01 0.000000000000000 0 0 0 1 0 0
1 1.000000000000000 0.000000000000000 0 0 0 0 1 0
-5 0.000000000000000 0.000000000000000 0 0 0 0 0 0

1, NO = 3, NV = 6, INA = 25
*********************************************

I COEFFICIENT ORDER EXPONENTS
NO = 3 NV = 6
1 1.000000000000000 0.000000000000000 0 0 0 0 0 1
-1 0.000000000000000 0.000000000000000 0 0 0 0 0 0
No Spin Matrix
Quaternion
c_quaternion

1, NO = 3, NV = 6, INA = 15
*********************************************

I COEFFICIENT ORDER EXPONENTS
NO = 3 NV = 6
0 0.9044197722704104 0.000000000000000 0 0 0 0 0 0
1 -1.230805330670363 0.000000000000000 1 0 0 0 0 0
1 0.7593604560872026 0.000000000000000 0 1 0 0 0 0
1 -0.2005031974645732 0.000000000000000 0 0 1 0 0 0
1 0.5873101678547868E-01 0.000000000000000 0 0 0 1 0 0
1 0.1062946080338559 0.000000000000000 0 0 0 0 0 1
-6 0.000000000000000 0.000000000000000 0 0 0 0 0 0

1, NO = 3, NV = 6, INA = 232
*********************************************

I COEFFICIENT ORDER EXPONENTS
NO = 3 NV = 6
0 0.3632814512983805 0.000000000000000 0 0 0 0 0 0
1 1.206286862714071 0.000000000000000 1 0 0 0 0 0
1 -2.545393061216957 0.000000000000000 0 1 0 0 0 0
1 -0.7194219382032653 0.000000000000000 0 0 1 0 0 0
1 0.2321269775156017 0.000000000000000 0 0 0 1 0 0
1 -0.3591574559148977 0.000000000000000 0 0 0 0 0 1
-6 0.000000000000000 0.000000000000000 0 0 0 0 0 0

1, NO = 3, NV = 6, INA = 233
*********************************************

I COEFFICIENT ORDER EXPONENTS
NO = 3 NV = 6
1 -0.1392966253184698 0.000000000000000 1 0 0 0 0 0
1 2.688418743115568 0.000000000000000 0 1 0 0 0 0
1 -3.100864700305056 0.000000000000000 0 0 1 0 0 0
1 2.355391437568686 0.000000000000000 0 0 0 1 0 0
1 0.1387778780781446E-16 0.000000000000000 0 0 0 0 0 1
-5 0.000000000000000 0.000000000000000 0 0 0 0 0 0

1, NO = 3, NV = 6, INA = 234
*********************************************

I COEFFICIENT ORDER EXPONENTS
NO = 3 NV = 6
0 -0.2237218421810473 0.000000000000000 0 0 0 0 0 0
1 -3.016884842915607 0.000000000000000 1 0 0 0 0 0
1 -1.063434273197377 0.000000000000000 0 1 0 0 0 0
1 -1.978759416966205 0.000000000000000 0 0 1 0 0 0
1 0.6143562772582531 0.000000000000000 0 0 0 1 0 0
1 -0.1534955027326110 0.000000000000000 0 0 0 0 0 1
-6 0.000000000000000 0.000000000000000 0 0 0 0 0 0
No Stochastic Radiation
Loading

0 comments on commit 46cbfb4

Please sign in to comment.