Dear Jakub,

We will look into it

Matthias

On 23 Jul 2015, at 10:00, Jakub Zakrzewski jakub.zakrzewski@gmail.com wrote:

Dear all,

As another example of using dwa-algorithm I tried to run 3D case in trap. The relevant part of "model.xml" is:

<SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> <PARAMETER name="K#" default="K"/> -mu#*n(i)+U#*n(i)*(n(i)-1)/2+K#*n(i)*((x-(L-1)/2)^2+(y-(L-1)/2)^2+(z-(L-1)/2)^2) </SITETERM>

which works with both worm and dwa. But suppose I want a QUARTIC trap and not harmonic so I modify to

<SITETERM site="i"> <PARAMETER name="mu#" default="mu"/> <PARAMETER name="U#" default="U"/> <PARAMETER name="K#" default="K"/> -mu#*n(i)+U#*n(i)*(n(i)-1)/2+K#*n(i)*((x-(L-1)/2)^2+(y-(L-1)/2)^2+(z-(L-1)/2)^2)^2 </SITETERM>

Note that the only difference is a power 2 on squared radius. Worm manages fine. Dwa gives error message "Cannot evaluate y"

What I am doing wrong? Is the failure of dwa due to an additional power 2 (since "dwa" in polish means "two")? I know the answer - dwa fails for any power different than 1....

Help appreciated

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 http://chaos.if.uj.edu.pl/~kuba