Skip to content

Commit 48731ef

Browse files
author
David Richie
committed
initial commit
0 parents  commit 48731ef

Some content is hidden

Large Commits have some content hidden by default. Use the searchbox below for content that may be hidden.

127 files changed

+46479
-0
lines changed

‎A.out

+102
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,102 @@
1+
Number_of_Trajectories 1
2+
0 0 0 0 0
3+
0.1 -0.19369 -1.13157 0.5231 0.0864218
4+
0.2 -0.836144 0.560665 -0.279843 1.06058
5+
0.3 -0.789374 1.30411 0.459902 1.10694
6+
0.4 -0.747222 1.28108 0.454987 1.17905
7+
0.5 -0.797061 1.33631 0.422966 0.883422
8+
0.6 -0.605418 1.42889 0.563409 0.596592
9+
0.7 -0.326183 1.47847 0.62412 0.281208
10+
0.8 -0.250129 1.4667 0.634468 0.213052
11+
0.9 -0.143535 1.38651 0.513788 0.0886517
12+
1 0.101217 1.3843 0.493163 -0.0965062
13+
1.1 0.0564498 1.40633 0.557509 -0.0645235
14+
1.2 0.200553 1.32969 0.456104 -0.157255
15+
1.3 0.393466 1.26735 0.369889 -0.278685
16+
1.4 0.208679 1.28488 0.422147 -0.149755
17+
1.5 0.324517 1.13053 0.287081 -0.185473
18+
1.6 0.358975 1.17147 0.319782 -0.217105
19+
1.7 0.345095 1.17138 0.323789 -0.200728
20+
1.8 0.4948 1.11153 0.261008 -0.282383
21+
1.9 0.613448 0.998443 0.162683 -0.311769
22+
2 0.686766 0.828151 0.0542823 -0.278904
23+
2.1 0.777057 0.703962 -0.0280322 -0.270001
24+
2.2 0.826076 0.522435 -0.10054 -0.205208
25+
2.3 0.79549 0.417098 -0.109644 -0.154267
26+
2.4 0.743275 0.324541 -0.102263 -0.110348
27+
2.5 0.759219 0.20468 -0.123108 -0.0698916
28+
2.6 0.765835 0.143096 -0.131018 -0.0498692
29+
2.7 0.779562 0.0490327 -0.141901 -0.0164498
30+
2.8 0.736146 -0.0322608 -0.126021 0.0123483
31+
2.9 0.681473 -0.170421 -0.0986821 0.0575459
32+
3 0.616146 -0.23531 -0.0719384 0.0701622
33+
3.1 0.569975 -0.270678 -0.0556615 0.0736033
34+
3.2 0.50906 -0.310053 -0.0349716 0.0742686
35+
3.3 0.462054 -0.32576 -0.0228282 0.0700858
36+
3.4 0.379148 -0.340772 -0.00469007 0.0589274
37+
3.5 0.323144 -0.360333 0.00728077 0.0527592
38+
3.6 0.263873 -0.386448 0.0200163 0.0457099
39+
3.7 0.199466 -0.382639 0.0254314 0.0334115
40+
3.8 0.15146 -0.382424 0.0288566 0.0248891
41+
3.9 0.112465 -0.377133 0.0300069 0.0179568
42+
4 0.0651347 -0.366975 0.029821 0.00949105
43+
4.1 0.0202948 -0.358957 0.0290825 0.00185911
44+
4.2 -0.0169693 -0.336555 0.0252107 -0.00381661
45+
4.3 -0.046084 -0.315917 0.0215816 -0.00756612
46+
4.4 -0.0680313 -0.290444 0.0174459 -0.009595
47+
4.5 -0.0835026 -0.271298 0.0145713 -0.0106967
48+
4.6 -0.105644 -0.244465 0.0104303 -0.0119368
49+
4.7 -0.120593 -0.222794 0.007431 -0.0122762
50+
4.8 -0.136518 -0.197471 0.00411002 -0.0121924
51+
4.9 -0.146557 -0.177999 0.00189579 -0.0117575
52+
5 -0.153735 -0.151502 -0.000491958 -0.0103839
53+
5.1 -0.155445 -0.12678 -0.00205291 -0.00870058
54+
5.2 -0.159634 -0.103581 -0.00349839 -0.00723356
55+
5.3 -0.164622 -0.0826125 -0.00471726 -0.00589259
56+
5.4 -0.163335 -0.0594173 -0.00528974 -0.00411226
57+
5.5 -0.16269 -0.0397021 -0.00564141 -0.00264513
58+
5.6 -0.157688 -0.0220505 -0.00547364 -0.00132611
59+
5.7 -0.152028 -0.00545748 -0.00514195 -0.000156538
60+
5.8 -0.144773 0.00905015 -0.004616 0.000773624
61+
5.9 -0.13502 0.0219761 -0.00388513 0.00147515
62+
6 -0.127563 0.0342834 -0.00327822 0.00208751
63+
6.1 -0.117795 0.0449576 -0.00253869 0.00247117
64+
6.2 -0.107208 0.0536845 -0.00180794 0.0026449
65+
6.3 -0.0962027 0.0616277 -0.00109848 0.00269291
66+
6.4 -0.0848532 0.0673477 -0.000480844 0.00256815
67+
6.5 -0.0740349 0.0717409 2.9093E-05 0.00236653
68+
6.6 -0.0630095 0.0754241 0.000477724 0.00209693
69+
6.7 -0.0521176 0.0763979 0.000769511 0.00173787
70+
6.8 -0.0426332 0.0773455 0.000988495 0.00142659
71+
6.9 -0.0334597 0.0774842 0.00113623 0.00110836
72+
7 -0.0246595 0.0764776 0.00120241 0.00079194
73+
7.1 -0.0160819 0.0752268 0.00122755 0.000488619
74+
7.2 -0.00801827 0.0730326 0.00118758 0.000210712
75+
7.3 5.32563E-06 0.0698926 0.00108929 -4.94386E-05
76+
7.4 0.00605975 0.0661682 0.000961204 -0.000220891
77+
7.5 0.0121745 0.0618265 0.000804568 -0.000371401
78+
7.6 0.0175518 0.0572607 0.000642202 -0.00047721
79+
7.7 0.0220926 0.0527936 0.000489768 -0.00054321
80+
7.8 0.025939 0.0478445 0.00033588 -0.000569726
81+
7.9 0.0291931 0.042806 0.000193544 -0.000567342
82+
8 0.0319077 0.0378254 6.73649E-05 -0.000543014
83+
8.1 0.0337313 0.0328089 -3.62097E-05 -0.000493467
84+
8.2 0.0349714 0.0279092 -0.000118777 -0.000431352
85+
8.3 0.0356989 0.02312 -0.000181662 -0.000361149
86+
8.4 0.0360879 0.0184161 -0.000228537 -0.000286962
87+
8.5 0.0357984 0.0138666 -0.000253005 -0.000210322
88+
8.6 0.0351247 0.00969943 -0.00026093 -0.000140343
89+
8.7 0.0341555 0.00583368 -0.000256457 -7.73506E-05
90+
8.8 0.0329202 0.00218185 -0.000241925 -2.09048E-05
91+
8.9 0.0313287 -0.000953313 -0.00021794 2.3214E-05
92+
9 0.029538 -0.00390455 -0.000188671 6.00773E-05
93+
9.1 0.0276144 -0.00651781 -0.000156812 8.75463E-05
94+
9.2 0.0255347 -0.008852 -0.000123198 0.000106661
95+
9.3 0.0233195 -0.0108119 -8.99343E-05 0.000116779
96+
9.4 0.0210597 -0.0124141 -5.91482E-05 0.000119504
97+
9.5 0.0187451 -0.0137677 -3.07583E-05 0.0001167
98+
9.6 0.0164075 -0.0148168 -6.05143E-06 0.000108845
99+
9.7 0.0141266 -0.015564 1.40231E-05 9.75366E-05
100+
9.8 0.0118704 -0.0160745 3.0099E-05 8.37975E-05
101+
9.9 0.00971182 -0.0163624 4.19222E-05 6.90133E-05
102+
10 0.00766058 -0.0163833 4.92987E-05 5.37773E-05

‎ACML-EULA.pdf

22.1 KB
Binary file not shown.

‎CHANGES

+87
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,87 @@
1+
Changes in version 1.3.5:
2+
3+
1. Removed obsolete compiler options from Complex.cc.
4+
5+
2. Fixed a minor bug in Complex.cc.
6+
7+
3. Fixed a minor bug in State.cc.
8+
9+
4. Fixed a minor bug in ACG.cc
10+
11+
-----------------------------------------------------------------------------
12+
13+
Changes in version 1.3.4:
14+
15+
1.) Fixed minor bug in Makefile.
16+
17+
2.) Improved instructions in README.
18+
19+
-----------------------------------------------------------------------------
20+
21+
Additions in version 1.3.3:
22+
23+
1.) New example program `onespin.cc'; improved comments in the existing
24+
example programs.
25+
26+
2.) The directory `doc' contains, in LaTeX as well as postscript
27+
format, a detailed description of the code, which has been published
28+
in Comp. Phys. Comm. 102 (1997) 210-228.
29+
30+
-----------------------------------------------------------------------------
31+
32+
Major differences between versions 1.2 and 1.3:
33+
34+
1.) The test program `testprog.cc' now uses a simple
35+
platform-independent random-number generator. The results produced by
36+
`testprog.cc' should now be the same for all operating systems and C++
37+
dialects.
38+
39+
2.) The directory `example_output' contains the output of all example
40+
programs. The output was obtained using Gnu g++ Version 2.7.2 on a
41+
DEC Alpha under Digital UNIX. Different compilers or operating systems
42+
may produce different results except for `testprog.cc' (see above).
43+
44+
3.) We added the example program `template.cc'.
45+
46+
4.) The method `Trajectory::plotExp' now allows one to specify a file
47+
into which the final state is saved.
48+
49+
5.) The method `Trajectory::sumExp' now computes the imaginary as well
50+
as the real parts of operator expectation values.
51+
52+
6.) `Makefile' is now be compatible with standard (nonGnu) versions of
53+
`make'.
54+
55+
-----------------------------------------------------------------------------
56+
57+
Major differences between versions 1.1 and 1.2:
58+
59+
1.) The quantum jump and orthogonal jump unravelings have been added
60+
via the "stepper" classes `AdaptiveJump' and `AdaptiveOrthoJump'.
61+
The example program `moving.cc' illustrates how to use `AdaptiveJump'.
62+
63+
2.) A bug involving a wrong nested loop in `sumExp' has been fixed.
64+
65+
-----------------------------------------------------------------------------
66+
67+
Here is a list of major differences between versions 1.0 and 1.1:
68+
69+
1.) The third argument of the function `plotExp' is now a list of file
70+
names, not a list of file pointers as before. The old version of
71+
`plotExp' has been renamed to `plotExp_obsolete'. The example programs
72+
`simple.cc', `spins.cc', and `moving.cc' show how to use `plotExp'.
73+
74+
2.) A new function named `sumExp' permits to compute the mean of
75+
expectation values over many trajectories. The example program
76+
`sums.cc' illustrates how to use 'sumExp'.
77+
78+
3.) The program now compiles and runs under VMS using DEC C++.
79+
80+
4.) A quantum jump stepper has been added, but detailed comments will
81+
not be available before version 1.2.
82+
83+
5.) The files `src/AtomOp.cc' and `include/AtomOp.h' define the new
84+
primary operator class `TransitionOperator'.
85+
86+
6.) A test program has been added to enable users to check correct
87+
installation. See the file `TEST' for details.

‎EXAMPLES

+20
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,20 @@
1+
onespin.cc: One spin, one or more trajectories, very detailed comments.
2+
3+
template.cc: Uses many features, may serve as a template.
4+
5+
testprog.cc: Test program; identical to template.cc except that it uses a
6+
simple portable random-number generator.
7+
8+
simple.cc: Simple harmonic oscillator, one trajectory.
9+
10+
moving.cc: Simple harmonic oscillator, using the moving basis, one
11+
trajectory.
12+
Different unravelings can be chosen by uncommenting or
13+
commenting out the appropriate lines.
14+
15+
spins.cc: Two interacting spins, time-dependent Hamiltonian, one
16+
trajectory.
17+
18+
sums.cc: Simple harmonic oscillator, sum over trajectories.
19+
20+

‎Einstein/Makefile

+165
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,165 @@
1+
# Makefile for QSD code.
2+
3+
QSD_DIR = .
4+
5+
CXX = g++
6+
INCLUDE = /opt/compiler/intel/ict/3.2/mkl/10.1/include
7+
FLAGS = -g -O -Wno-deprecated -I$(INCLUDE)
8+
FC = ifort -c
9+
DFLAGS = -D__INTEL -D__LAPACK -D__FFTW3
10+
##-D__FFTSG -D__FFTW2 -D__parallel -D__BLACS -D__SCALAPACK
11+
FCFLAGS = -O2 -ftz -IPF-fp-relaxed -heap-arrays 64 -I$(INCLUDE)
12+
13+
INC = $(QSD_DIR)/include
14+
SRC = $(QSD_DIR)/src
15+
16+
LOADLIBES = -lm
17+
LDFLAGS = -L/opt/compiler/intel/ict/3.2/mkl/10.1/lib/em64t -lmkl_scalapack
18+
LDFLAGS = -L$(MKLPATH) -lmkl_scalapack -lmkl_solver_ilp64_sequential -Wl,--start-group -lmkl_intel_ilp64 -lmkl_sequential -lmkl_core -Wl,--end-group -lpthread -L/opt/compiler/intel/compiler/11.0/lib/intel64/ -lifcore $(LOADLIBES)
19+
FFTFLAGS = -L/usr/cta/unsupported/fftw/3.1.2-intel/lib -lfftw3
20+
##FFTFLAGS = -L$(HOME)/lib -lfftw3
21+
CXXFLAGS = $(DCC) -I$(INC) -I$(HOME)/include $(FLAGS)
22+
CPPFLAGS = -L$(MKLPATH) -lmkl_lapack -lmkl -lguide
23+
24+
COMPILE = $(CXX) -c $(CPPFLAGS) $(CXXFLAGS)
25+
FCOMPILE = $(FC) $(FCFLAGS) $(DFLAGS)
26+
LINK = icpc $(CPPFLAGS) $(CXXFLAGS) $(LDFLAGS) $(FFTFLAGS)
27+
28+
RAN_OBJ = Random.o Normal.o Uniform.o RNG.o ACG.o MLCG.o CmplxRan.o
29+
RAN_HEAD = $(INC)/Random.h $(INC)/Normal.h $(INC)/Uniform.h $(INC)/RNG.h $(INC)/ACG.h $(INC)/MLCG.h $(INC)/CmplxRan.h
30+
31+
OBJ1 = Operator.o PrimOp.o State.o Complex.o
32+
HEAD1 = $(INC)/Operator.h $(INC)/PrimOp.h $(INC)/State.h $(INC)/Complex.h $(INC)/Mesh.h
33+
OBJ2 = FieldOp.o SpinOp.o AtomOp.o
34+
OBJ3 = Traject.o $(RAN_OBJ)
35+
OBJ4 = mesh.o coeff.o poisson.o
36+
OBJ5 = Mesh.o Coeff.o Epot.o Poisson.o
37+
OBJ = $(OBJ1) $(OBJ2) $(OBJ3) $(OBJ4) $(OBJ5)
38+
#-------------------------------------------------------------------
39+
# Add driver routines here.
40+
41+
ALL = onespin spins simple moving sums testprog template lineob damped qcascade sechar qd
42+
43+
qd: qd.cc $(OBJ)
44+
$(LINK) -o qd qd.cc $(OBJ) $(LOADLIBES)
45+
46+
##qdf: qdf.f90 $(OBJ)
47+
## $(LINK) -o qdf qdf.f90 $(OBJ) $(LOADLIBES)
48+
49+
onespin: onespin.cc $(OBJ)
50+
$(LINK) -o onespin onespin.cc $(OBJ) $(LOADLIBES)
51+
52+
spins: spins.cc $(OBJ)
53+
$(LINK) -o spins spins.cc $(OBJ) $(LOADLIBES)
54+
55+
simple: simple.cc $(OBJ)
56+
$(LINK) -o simple simple.cc $(OBJ) $(LOADLIBES)
57+
58+
lineob: lineob.cc $(OBJ)
59+
$(LINK) -o lineob lineob.cc $(OBJ) $(LOADLIBES)
60+
61+
qcascade: qcascade.cc $(OBJ)
62+
$(LINK) -o qcascade qcascade.cc $(OBJ) $(LOADLIBES)
63+
64+
sechar: sechar.cc $(OBJ)
65+
$(LINK) -o sechar sechar.cc $(OBJ) $(LOADLIBES)
66+
67+
damped: damped.cc $(OBJ)
68+
$(LINK) -o damped damped.cc $(OBJ) $(LOADLIBES)
69+
70+
moving: moving.cc $(OBJ)
71+
$(LINK) -o moving moving.cc $(OBJ) $(LOADLIBES)
72+
73+
sums: sums.cc $(OBJ)
74+
$(LINK) -o sums sums.cc $(OBJ) $(LOADLIBES)
75+
76+
testprog: testprog.cc $(OBJ)
77+
$(LINK) -o testprog testprog.cc $(OBJ) $(LOADLIBES)
78+
79+
template: template.cc $(OBJ)
80+
$(LINK) -o template template.cc $(OBJ) $(LOADLIBES)
81+
82+
#-------------------------------------------------------------------
83+
84+
all: $(ALL)
85+
86+
cleanexe:
87+
rm $(ALL)
88+
89+
cleanrand:
90+
rm $(RAN_OBJ)
91+
92+
clean:
93+
rm *.o
94+
95+
Traject.o: $(SRC)/Traject.cc $(INC)/Traject.h $(HEAD1) $(RAN_HEAD)
96+
$(COMPILE) $(SRC)/Traject.cc
97+
98+
State.o: $(SRC)/State.cc $(HEAD1)
99+
$(COMPILE) $(SRC)/State.cc
100+
101+
Operator.o: $(SRC)/Operator.cc $(HEAD1)
102+
$(COMPILE) $(SRC)/Operator.cc
103+
104+
PrimOp.o: $(SRC)/PrimOp.cc $(HEAD1)
105+
$(COMPILE) $(SRC)/PrimOp.cc
106+
107+
FieldOp.o: $(SRC)/FieldOp.cc $(INC)/FieldOp.h $(HEAD1)
108+
$(COMPILE) $(SRC)/FieldOp.cc
109+
110+
SpinOp.o: $(SRC)/SpinOp.cc $(INC)/SpinOp.h $(HEAD1)
111+
$(COMPILE) $(SRC)/SpinOp.cc
112+
113+
AtomOp.o: $(SRC)/AtomOp.cc $(INC)/AtomOp.h $(HEAD1)
114+
$(COMPILE) $(SRC)/AtomOp.cc
115+
116+
Complex.o: $(SRC)/Complex.cc $(INC)/Complex.h
117+
$(COMPILE) $(SRC)/Complex.cc
118+
119+
Mesh.o: $(SRC)/Mesh.cc $(INC)/Mesh.h
120+
$(COMPILE) $(SRC)/Mesh.cc
121+
122+
Epot.o: $(SRC)/Epot.cc $(INC)/Epot.h
123+
$(COMPILE) $(SRC)/Epot.cc
124+
125+
Coeff.o: $(SRC)/Coeff.cc
126+
$(COMPILE) $(SRC)/Coeff.cc
127+
128+
Poisson.o: $(SRC)/Poisson.cc
129+
$(COMPILE) $(SRC)/Poisson.cc
130+
131+
Random.o: $(SRC)/Random.cc $(RAN_HEAD)
132+
$(COMPILE) $(SRC)/Random.cc
133+
134+
Normal.o: $(SRC)/Normal.cc $(RAN_HEAD)
135+
$(COMPILE) $(SRC)/Normal.cc
136+
137+
Uniform.o: $(SRC)/Uniform.cc $(RAN_HEAD)
138+
$(COMPILE) $(SRC)/Uniform.cc
139+
140+
RNG.o: $(SRC)/RNG.cc $(RAN_HEAD)
141+
$(COMPILE) $(SRC)/RNG.cc
142+
143+
ACG.o: $(SRC)/ACG.cc $(RAN_HEAD)
144+
$(COMPILE) $(SRC)/ACG.cc
145+
146+
MLCG.o: $(SRC)/MLCG.cc $(RAN_HEAD)
147+
$(COMPILE) $(SRC)/MLCG.cc
148+
149+
CmplxRan.o: $(SRC)/CmplxRan.cc $(RAN_HEAD) $(INC)/Complex.h
150+
$(COMPILE) $(SRC)/CmplxRan.cc
151+
152+
coeff.o: $(SRC)/coeff.f90
153+
$(FCOMPILE) $(SRC)/coeff.f90
154+
155+
fftw3.o: $(SRC)/fftw3.f90
156+
$(FCOMPILE) $(SRC)/fftw3.f90
157+
158+
cf.o: $(SRC)/cf.f90
159+
$(FCOMPILE) $(SRC)/cf.f90
160+
161+
poisson.o: $(SRC)/poisson.f90
162+
$(FCOMPILE) $(SRC)/poisson.f90
163+
164+
mesh.o : $(SRC)/mesh.f90
165+
$(FCOMPILE) $(SRC)/mesh.f90

0 commit comments

Comments
 (0)