ALF issueshttps://git.physik.uni-wuerzburg.de/groups/ALF/-/issues2020-03-24T19:03:16Zhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/117Update_doc2020-03-24T19:03:16ZFakher F. AssaadUpdate_docHi Jefferson, in the ALF 1.2 and next versions the \hat{Z}_k could be the third Pauli matrix, in which case the field s takes the from \pm 1. That is: \hat{Z}_k | s > = \pm 1 | s >. This would correspond to the vertex type 1. B...Hi Jefferson, in the ALF 1.2 and next versions the \hat{Z}_k could be the third Pauli matrix, in which case the field s takes the from \pm 1. That is: \hat{Z}_k | s > = \pm 1 | s >. This would correspond to the vertex type 1. But now \hat{Z}_k could be for example the position operator. In this case, \hat{Z}_k | Phi > = \phi| \phi > with \phi a real number. In this case, the vertex type is set to three. One should update the doc to take this into account. The only thing to change is to make it very clear to the user that s_{i,\tau} can be a \pm 1 variable or a continuous real field. Thanks Fakher.Jefferson Stafusa E. PortelaJefferson Stafusa E. Portelahttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/116Histogram observables2019-10-03T23:05:37ZJonas SchwabHistogram observablesImplementing a new observable type, tailored to histograms. Also, since one histogram bin should capture the whole space-time configuration, perhaps implement some structure in main.F90 to conveniently measure histograms once per sweep.Implementing a new observable type, tailored to histograms. Also, since one histogram bin should capture the whole space-time configuration, perhaps implement some structure in main.F90 to conveniently measure histograms once per sweep.Jonas SchwabJonas Schwabhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/115Sequential_Sweeps2019-10-02T15:00:26ZJefferson Stafusa E. PortelaSequential_SweepsImplement a switch for deciding whether the program carries out sequential single spin-flip sweeps.
Such a switch was mentioned in earlier versions of the documentation under the name "Sequential_Sweeps".Implement a switch for deciding whether the program carries out sequential single spin-flip sweeps.
Such a switch was mentioned in earlier versions of the documentation under the name "Sequential_Sweeps".https://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/114block diagonal interaction terms2019-09-29T10:50:34ZJohannes Hofmannblock diagonal interaction termsSome models have a block diagonal structure in the interaction terms, similar to the already exploited N_fl structure. However, the non-interacting part may not be block diagonal such that N_fl cannot be used.
This branch aims to impleme...Some models have a block diagonal structure in the interaction terms, similar to the already exploited N_fl structure. However, the non-interacting part may not be block diagonal such that N_fl cannot be used.
This branch aims to implement this structure!
I would propose to redefine the operator type to contain a list of the "dense" blocks of the operators. I.e., P(:) --> P(:,b); O(:,:) --> O(:,:,b) and so on. b should be the index of the block and appear as the last index for memory efficiency as with n_fl. The rest of the code should then be updated accordingly.
It might also be convenient, to use the number of blocks as an optional argument when the operator is allocated, ideally in a way that the existing model without that structure can stay as they are. This way, the more experienced user can make use of this option whereas a "new" user doesn't have to care about this in the beginning.Emilie HuffmanEmilie Huffmanhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/113Cholesky vs maxent2022-05-30T13:12:54ZFlorian GothCholesky vs maxentyou have to touch xqmc1 and Vhelp.you have to touch xqmc1 and Vhelp.Florian GothFlorian Gothhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/112testsuite compilation failure2019-09-24T09:46:31ZJefferson Stafusa E. Portelatestsuite compilation failure### Problem
When the main code has been compiled using, e.g., the Intel compiler, then running **make** in `testsuite/tests` as described in README.md will fail, because it uses gfortran -- which won't recognize the intel-compiled modul...### Problem
When the main code has been compiled using, e.g., the Intel compiler, then running **make** in `testsuite/tests` as described in README.md will fail, because it uses gfortran -- which won't recognize the intel-compiled modules.
**Details**:
Error is `f951: Fatal Error: Reading module 'mymats' at line 1 column 2: Unexpected EOF`,
which disappears when switching from gfortran to the compiler set the configureHPC in `tests/CMakeCache.txt` (CMAKE_Fortran_COMPILER:FILEPATH)
---
### Possible solution
To use the environment variable ``$ALF_FC`` (set by configureHPC.sh) instead.https://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/106Use scientific number format in analysis of Green functions2019-08-31T18:28:14ZMartin HohenadlerUse scientific number format in analysis of Green functionsThe current format (F16.8) in cov_tau.f90 can be insufficient, leading to zero average values and/or errors.The current format (F16.8) in cov_tau.f90 can be insufficient, leading to zero average values and/or errors.https://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/103Suggestion: increment N_auto in powers of two to speed up autocorrelation ana...2019-05-21T12:02:40ZMartin HohenadlerSuggestion: increment N_auto in powers of two to speed up autocorrelation analysisThe documentation points out that the autocorrelation analysis is expensive. However, there is very little use in analyzing and including every bin length between 1 and 500 (in the example of Fig.1). It is common in the literature to inc...The documentation points out that the autocorrelation analysis is expensive. However, there is very little use in analyzing and including every bin length between 1 and 500 (in the example of Fig.1). It is common in the literature to increase N_auto by a factor of two instead.https://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/96Reconsider current use of confin and confout files2019-02-24T15:45:08ZMartin HohenadlerReconsider current use of confin and confout filesFrom my perspective, the potential pitfalls (incorrect error bars, waste of CPU time) of the current way the configuration files are handled outweigh any advantages I can think of.
On the present occasion, problems arose for the followi...From my perspective, the potential pitfalls (incorrect error bars, waste of CPU time) of the current way the configuration files are handled outweigh any advantages I can think of.
On the present occasion, problems arose for the following nontrivial reason: although I used the example job file that includes a call to out_to_in.sh, things failed because the script file did not have the right permissions (as can easily happen when copying files to a new location). The problem remained completely hidden because the output of out_to_in.sh is piped to /dev/null.https://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/94automatically generate zslgemm and zslhemm with arbitary N2019-02-07T22:01:21ZJonas Schwabautomatically generate zslgemm and zslhemm with arbitary NCreate Libraries/Modules/Mat_subroutines.f90 with script during compilation with arbitrary high N.Create Libraries/Modules/Mat_subroutines.f90 with script during compilation with arbitrary high N.Jonas SchwabJonas Schwabhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/93Analysis routines extremely slow on JURECA2018-11-15T11:07:05ZMartin HohenadlerAnalysis routines extremely slow on JURECARoughly since the update of the software stack, even the analysis of equal-time correlations with less than 100 bins of data takes much too long. There seems to be a fundamental issue reflected in an overload of the CPU's. For example,
...Roughly since the update of the software stack, even the analysis of equal-time correlations with less than 100 bins of data takes much too long. There seems to be a fundamental issue reflected in an overload of the CPU's. For example,
15846 hwb03b 20 0 5708480 26612 7648 R 7987 0.0 34:35.29 cov_eq.out
suggesting 8000% of CPU usage.Johannes HofmannJohannes Hofmannhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/91Test the compact WY decomposition in real life2018-10-01T15:55:30ZFlorian GothTest the compact WY decomposition in real lifeMake a branch from Johannes' branch and introduce xGEQRT into udv_state in the case of the logscale where permutations are disabled anyway.Make a branch from Johannes' branch and introduce xGEQRT into udv_state in the case of the logscale where permutations are disabled anyway.Florian GothFlorian Gothhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/90PLASMAAlf2020-09-04T12:08:18ZFlorian GothPLASMAAlfPlasma,
https://bitbucket.org/icl/plasma
seems to have a superior QR decomposition that is highly parralizable.
Let's see where this leadsPlasma,
https://bitbucket.org/icl/plasma
seems to have a superior QR decomposition that is highly parralizable.
Let's see where this leadsFlorian GothFlorian Gothhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/89Single Precision in ALF2020-07-23T18:04:23ZFlorian GothSingle Precision in ALFnobody knows what happens when we compile ALF for using single precision. Let's find out!nobody knows what happens when we compile ALF for using single precision. Let's find out!Florian GothFlorian Gothhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/87Implement checksums ofr MPI communication2018-05-09T10:02:58ZJohannes HofmannImplement checksums ofr MPI communicationIntroduce checksums for MPI communication to test the correctness of the transferred data, starting with tempering moves and observable printout.Introduce checksums for MPI communication to test the correctness of the transferred data, starting with tempering moves and observable printout.Johannes HofmannJohannes Hofmannhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/86ALF-for-Summit2018-05-02T15:21:24ZFlorian GothALF-for-SummitWe intend to port ALF to the Summit HPC System with the help of the MAGMA library.We intend to port ALF to the Summit HPC System with the help of the MAGMA library.Johannes HofmannJohannes Hofmannhttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/79./analysis.sh: line 20: 28106 CPU time limit exceeded $ANNAL/cov_eq.out2017-12-04T19:12:52ZMartin Hohenadler./analysis.sh: line 20: 28106 CPU time limit exceeded $ANNAL/cov_eq.outSince pulling the latest version yesterday and compiling with standard flags (O3), the analysis has become extremely slow on JURECA. The QMC code itself runs normally. The numbers of bins was only 120 and we are talking about equal-time ...Since pulling the latest version yesterday and compiling with standard flags (O3), the analysis has become extremely slow on JURECA. The QMC code itself runs normally. The numbers of bins was only 120 and we are talking about equal-time analysis.https://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/77Examples/Hubbard_SU2_Ising_Square2/parameters2017-10-20T16:34:52ZMartin HohenadlerExamples/Hubbard_SU2_Ising_Square2/parametersreplace &VAR_Hubbard by &VAR_Hub_Isingreplace &VAR_Hubbard by &VAR_Hub_Isinghttps://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/76Automatically set values for certain parameters2017-10-02T07:13:10ZMartin HohenadlerAutomatically set values for certain parametersWould it be possible to enable automatically setting variables such as LOBS_EN for given beta, dtau? In my codes, I usually allow for negative values of such flags to signal automatic values. This avoids accidentally using incorrect values.Would it be possible to enable automatically setting variables such as LOBS_EN for given beta, dtau? In my codes, I usually allow for negative values of such flags to signal automatic values. This avoids accidentally using incorrect values.https://git.physik.uni-wuerzburg.de/ALF/ALF/-/issues/74Optimization for hopping operators2017-09-12T15:59:41ZJohannes HofmannOptimization for hopping operatorsIt might be useful to specialize hopping operator as many of us simulate Ising fields coupled to bond hopping or hopping operators squared in Kondo type couplings.It might be useful to specialize hopping operator as many of us simulate Ising fields coupled to bond hopping or hopping operators squared in Kondo type couplings.