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