Dear All, I use worm with a harmonic trap defined in a standard version as in ALPS2.0 paper with the following addition in <SITETERM> (for 2D vesion)
K#*n(i)*((x-(L-1)/2)^2+(y-(L-1)/2)^2)
works perfectly. I tried dwa getting the error message:
cannot evaluate y
Trying the one dimensional version with just K#*n(i)*(x-(L-1)/2)^2
works for both worm and dwa.
any idea why 2D harmonic trap cannot be used with dwa?
All the best kuba
Jakub Zakrzewski
Zakład Optyki Atomowej, Instytut Fizyki im. Mariana Smoluchowskiego, UJ ul. Lojasiewicza 11, 30-348 Kraków, Poland +48-12-6644555, fax +48126644905 http://chaos.if.uj.edu.pl/~kuba