April 2015 - Comp-phys-alps-users

Inconsistent comment in generate_spline_matrix in the DMFT folder
by Kuangshing Chen 22 Jul '15

22 Jul '15

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-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

5 25

Warning: No parameter set has been found. No XML files will be generated
by Alexander Herzog 02 May '15

02 May '15

2 1

Cannot calculate boundary crossing if neither vertex is in the original cell
by Alexander Herzog 02 May '15

02 May '15

2 11

building alps from source
by Tadeusz Wasiutyński 30 Apr '15

30 Apr '15

while cmake looks OK I meet problem during make command: ********************************************************************************************************************************* [ 20%] Building CXX object src/boost/CMakeFiles/boost.dir/home/twasiutynsk/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/once.cpp.o In file included from /home/twasiutynsk/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic.hpp:12:0, from /home/twasiutynsk/alps-2.2.b3-r7462-src-with-boost/boost/boost/thread/pthread/once_atomic.hpp:20, from /home/twasiutynsk/alps-2.2.b3-r7462-src-with-boost/boost/boost/thread/once.hpp:20, from /home/twasiutynsk/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/./once_atomic.cpp:9, from /home/twasiutynsk/alps-2.2.b3-r7462-src-with-boost/boost/libs/thread/src/pthread/once.cpp:8: /home/twasiutynsk/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic/atomic.hpp:202:16: error: ‘uintptr_t’ was not declared in this scope typedef atomic<uintptr_t> atomic_uintptr_t; ^ /home/twasiutynsk/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic/atomic.hpp:202:25: error: template argument 1 is invalid typedef atomic<uintptr_t> atomic_uintptr_t; ^ /home/twasiutynsk/alps-2.2.b3-r7462-src-with-boost/boost/boost/atomic/atomic.hpp:202:43: error: invalid type in declaration before ‘;’ token typedef atomic<uintptr_t> atomic_uintptr_t; ********************************************************************************************************************************** is it connected with my compiler or alps source? In the past on the other machine compilation went smooth. Regards -- Tadeusz Wasiutyński

2 1

model not implemented in file factory.C
by Alexander Herzog 18 Apr '15

18 Apr '15

2 7

Problems with ALPS fulldiag measurement code and suggested fix
by Pikulski Marek 14 Apr '15

14 Apr '15

Dear ALPS Users and Developers, I have been observing some unexpected behavior when using translation symmetry with the ALPS fulldiag application. The detailed description (PDF), as well as a simple example and the output in question (.zip), are attached. More recently, I have discovered another issue: If I consider a single spin-1/2 in the Sz basis subjected to magnetic field pointing along the y-direction, the eigenvectors and energies are calculated correctly but the expectation value of e.g. the identity-operator is zero. I have tracked down the problem to a surplus complex conjugation in the calculation of the scalar product of complex vectors and the attached patch seems to fix it. Moreover, the originally encountered problem (cf. first paragraph) seems to be no more than a different manifestation of the same bug. I hope this information will be useful. Best regards, Marek Pikulski Laboratory for Solid State Physics ETH Zurich

2 2