Hi,
The function generate_spline_matrix in fouriertransform.C in the dmft
folder shows the comments as follows:
// A is the matrix whose inverse defines spline_matrix
//
// 6 6
// 1 4 1
// 1 4 1
// A = ...
//
// 1 4 1
// -2 0 0 2
However, the following codes,
dense_matrix A =
4*dt/6.*boost::numeric::ublas::identity_matrix<double>(Np1);
for (int i=1; i<Np1-1; i++) {
A(i,i-…
[View More]1) = dt/6.;
A(i,i+1) = dt/6.;
}
A(0,0) = 1.;
A(0, Np1-1) = 1.;
A(Np1-1, 0) = -2.*dt/6.;
A(Np1-1, 1) = -1.*dt/6.;
A(Np1-1, Np1-2) = 1*dt/6.;
A(Np1-1, Np1-1) = 2*dt/6.;
the red lines will make the matrix looks like
// 6 6
// 1 4 1
// 1 4 1
// A = ...
//
// 1 4 1
// -2 *-1* *1* 2
Is the comment correct or the code correct? Any reference for this
algorithm?
Thank you,
Kuang-Shing Chen
[View Less]
Dear ALPS users,
I have a large list of tasks run from a given job file *parm.in.xml*,
which got terminated after some tasks were completed; my *parm.out.xml*
looks like this after the termination:
===========<SNIP-1>=================
/......//
//......//
// <TASK status="finished">//
// <INPUT file="parm.task62.out.xml"/>//
// </TASK>//
// <TASK status="finished">//
// <INPUT file="parm.task63.out.xml"/>//
// </TASK>//
// <TASK status="…
[View More]new">//
// <INPUT file="parm.task64.out.xml"/>//
// </TASK>//
// <TASK status="new">//
//......//
//....../
============<SNIP-1>================
To restart the job from task64, I execute the command (without changing
anything in the output files)
*worm**--Tmin 10 parm.out.xml*
Could you please confirm that this is correct? I ask because nothing
seems to proceed after a while after executing the above:
============<SNIP-2>===============
/.........//
//.........//
//Loading information about run 1 from file parm.task62.out.run1//
//Loading information about run 1 from file parm.task63.out.run1//
//Task 1 finished.//
//........//
//........//
/
/Task 63 finished.//
//Created run 1 locally//
//Starting task 64.//
//Checking if it is finished: not yet, next check in 10 seconds ( 0%
done).//
//Checking if it is finished: not yet, next check in 203 seconds ( 1%
done).//
//Checking if it is finished: not yet, next check in 174 seconds ( 23%
done)./
=============<SNIP-2>==============
>From thereon the simulation just hangs for many hours (I've tried this
procedure repeatedly), which it should not because my SWEEPS and
THERMALIZATION are not excessively large.
Any ideas on why the re-started simulations don't seem to proceed are
appreciated.
With regards,
Vipin
[View Less]
Hello,
My model need to used complex number. Then I wrote the value of a parameter as "s=0.5*sqrt(-1)" and I got the error message "Cannot evaluate I". What is wrong with me? What can I do to input a complex number? Thank you for your reply.
Best wishes
Tianyou
Dear alps users,
I have been trying to build alps package from source on Windows 7
OS. I have been following the outdated instructions on the page
https://alps.comp-phys.org/trac/wiki/AlpsOnWindows.
I could not find the proper versions of Zlib, SZlib, Lpsolve and SQlite.
There are links on this webpage that read "Check out the source code from
<>" where <> is, for instance,
https://alps.comp-phys.org/svn/alps1/trunk/wistaria/windows/lpsolve for
"Lpsolve".…
[View More] I don't know how to download the source codes using these links.
I tried to download the source codes from a trusted webpage like
sourceforge.net, some of them work while others don't and I keep getting
errors.
Please help me.
Regards,
Niraj Ghimire
Graduate Student
Department of Physics
University of Connecticut
[View Less]
Hello,
In my model, some parameters need to be a complexe number. How can I give a complexe number in the model file. Thank you for your reply.
Best wishes
Tianyou
Hello,
I need to create my own Hamiltonian. I have a term of nearest-neighbor interaction, and I write in the model file
</BONDTERM>
<BONDTERM source="i" target="j">
-t*bdag(i)*b(j)
</BONDTERM>
Here in the field of the target should I write target="i+1" or should I write simply target="j". What is the difference between the two cases ? In the second case, the source i will go through all the sites and the target j will go through alll the sites too, I have then interaction of …
[View More]next-nearest-neighbor terms and so on ?
Thank you for your reply.
Best wishes
Tianyou
[View Less]
Dear ALPS users/developers,
After cmake and ccmake, I am trying to 'make' the ALPS 2.2b1 files on my
computer (running on Linux Mint 13 Maya, which is in turn based on
Ubuntu 12.04); however I get the following error message:
/*<SNIP>*/
Linking CXX static library libspinmc_impl.a
[ 67%] Built target spinmc_impl
Scanning dependencies of target spinmc
[ 68%] Building CXX object
applications/mc/spins/CMakeFiles/spinmc.dir/spinmc.C.o
In file included from
/afs/ictp.it/home/v/vvarma/Desktop/…
[View More]alps-2.2.b1-r7195-src-with-boost/boost/boost/python/detail/wrap_python.hpp:50:0,
from
/afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/boost/boost/python/detail/prefix.hpp:13,
from
/afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/boost/boost/python/args.hpp:8,
from
/afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/boost/boost/python.hpp:11,
from
/afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/src/alps/alea/mcdata.hpp:79,
from
/afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/src/alps/alea.h:34,
from
/afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/src/alps/alea/observableset.h:49,
from
/afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/src/alps/scheduler/montecarlo.h:40,
from
/afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/src/alps/scheduler.h:31,
from
/afs/ictp.it/home/v/vvarma/Desktop/alps-2.2.b1-r7195-src-with-boost/alps/applications/mc/spins/spinmc.C:29:
/usr/include/python2.7/pyconfig.h:1161:0: warning: "_POSIX_C_SOURCE"
redefined [enabled by default]
/usr/include/features.h:164:0: note: this is the location of the
previous definition
/usr/include/python2.7/pyconfig.h:1183:0: warning: "_XOPEN_SOURCE"
redefined [enabled by default]
/usr/include/features.h:166:0: note: this is the location of the
previous definition
make[2]: *** No rule to make target `/usr/lib/libhdf5.so', needed by
`applications/mc/spins/spinmc'. Stop.
make[1]: *** [applications/mc/spins/CMakeFiles/spinmc.dir/all] Error 2
make: *** [all] Error 2
/*<SNIP>*/
I have tried reinstalling hdf5 (libhdf5-serial-dev) and the ALPS
multiple times on this computer but there is always some error related
to the hdf5 files, either during the make process or the cmake
configuration. If you could please help me debug this particular error,
that would be greatly appreciated. A few technical notes are appended at
the end.
Thanks so much.
Regards,
Vipin
======================================
1) Surprisingly I am able to build and install ALPS without any problems
on a laptop running native Ubuntu 12.04.
2) *cmake command*
cmake -D
Boost_ROOT_DIR:PATH=~/Desktop/alps-2.2.b1-r7195-src-with-boost/boost
-DALPS_INSTALL_VISTRAILS_PACKAGES=ON -D
VISTRAILS_APP_DIR:PATH=~/Desktop/vistrails-src-nightly-09cb1cfff309/
-DCMAKE_INSTALL_PREFIX=~/Desktop/alps-2.2.b1-r7195-src-with-boost/opt_alps
~/Desktop/alps-2.2.b1-r7195-src-with-boost/alps
3) *Snippet of CMakeCache.txt*
//HDF5 C++ Wrapper compiler. Used only to detect HDF5 compile
// flags.
HDF5_CXX_COMPILER_EXECUTABLE:FILEPATH=HDF5_CXX_COMPILER_EXECUTABLE-NOTFOUND
//HDF5 Wrapper compiler. Used only to detect HDF5 compile flags.
HDF5_C_COMPILER_EXECUTABLE:FILEPATH=HDF5_C_COMPILER_EXECUTABLE-NOTFOUND
//Path to a file.
HDF5_C_INCLUDE_DIR:PATH=/usr/include
//HDF5 file differencing tool.
HDF5_DIFF_EXECUTABLE:FILEPATH=HDF5_DIFF_EXECUTABLE-NOTFOUND
//Path to a file.
HDF5_DLL:FILEPATH=HDF5_DLL-NOTFOUND
//Path to a file.
HDF5_HL_DLL:FILEPATH=HDF5_HL_DLL-NOTFOUND
//HDF5 library compiled with parallel IO support
HDF5_IS_PARALLEL:BOOL=FALSE
//The HDF5_hdf5 library
HDF5_hdf5_LIBRARY:FILEPATH=/usr/include
//Path to a library.
HDF5_hdf5_LIBRARY_DEBUG:FILEPATH=HDF5_hdf5_LIBRARY_DEBUG-NOTFOUND
//Path to a library.
HDF5_hdf5_LIBRARY_RELEASE:FILEPATH=/usr/lib/libhdf5.so
//The HDF5_hdf5_hl library
HDF5_hdf5_hl_LIBRARY:FILEPATH=/usr/include
//Path to a library.
HDF5_hdf5_hl_LIBRARY_DEBUG:FILEPATH=HDF5_hdf5_hl_LIBRARY_DEBUG-NOTFOUND
//Path to a library.
HDF5_hdf5_hl_LIBRARY_RELEASE:FILEPATH=/usr/lib/libhdf5_hl.so
[View Less]
Dear alps users,
I want to use alps package with my 32 bit version of Python 2.7 on
Windows 7 OS. I have been told that I need to build the package from source
but I have no idea how to do this. Could you please provide me with a clear
step-by-step procedure ?
I would really appreciate your help.
Regards,
Niraj Ghimire
Graduate Student
Department of Physics
University of Connecticut
Dear ALPS users,
I am trying to run a very simple worm simulation with disorder
(site-dependent 'mu') for the Bose-Hubbard model using a Python script.
If I simply translate a working parameter file (obtained from the ALPS
website) into the below Python code, 'random()' is unrecognized (which
is not so for the parameter file where the simulation goes through).
But if I include 'import random' at the beginning of the Python script
then the simulation does go through but a different model gets …
[View More]simulated
(i.e. a random 'mu' but equal for all sites).
I was wondering, should not the below Python script translated directly
from a similar working parameter file run as is for the desired model?
/*<SNIP>*/
import pyalps
import matplotlib.pyplot as plt
import pyalps.plot
#prepare the input parameters
parms = []
for DISORDERSEED in [34275, 49802, 82529]:
parms.append(
{
'LATTICE' : "inhomogeneous square lattice periodic",
'MODEL' : "boson Hubbard",
'T' : 0.1,
'L' : 4 ,
't' : 1.0 ,
'DISORDERSEED' : DISORDERSEED,
'mu' : 5.0*2*(random()-0.5),
'U' : 1.0 ,
'NONLOCAL' : 0 ,
'Nmax' : 2 ,
'THERMALIZATION' : 100,
'SWEEPS' : 5000
}
)
/*<SNIP>*/
Thanks,
Vipin
[View Less]
Hello everyone,
I have installed ALPS in Win7 and I also installed Pythonxy2.7.6.0, but I have not installed Vistrails because I cannot download it.
But I do not know how to run ALPS, can anyone help me ?
Thank you very much!
Tianyou Yi