[Octopus-devel] possible bug in TDOutput

Nicole Helbig nehelbig at gmail.com
Sat Mar 6 21:11:15 WET 2010


Dear all,

we found a very strange behavior with the TDOutput today. If you run a
1D calculation and specify TDOutput = energy + multipoles + potential
you get an error message saying

Node     0: Assertion "gr%mesh%sb%dim.ne.1" failed in line   481 in
file "states_calc_inc.F90"

Now, I know that "potential" is not an available option for TDOutput,
so legally I shouldn't do what I am doing. However, I checked and the
error message is actually coming from trying to write the angular
momentum which of course doesn't make sense in 1D. I am just wondering
how the calculation ever got to this part of the code in the first
place. I found that angular momentum is classed as Experimental, so I
thought I let the person who is doing that experiment know about this.

Cheers,

Nicole


More information about the Octopus-devel mailing list