Re: [ESPResSo-devel] p3m kspace force calc, p3m stress tensor

[Tobias Richter]

Hi Axel,

die Komponenten haben vergleichbare Groesse (10^4) sowohl im Realteil der
Kraft (E-Feld) als auch im Imaginaerteil. Die Mpi Meldung, dass der
Imaginaerteil nicht null ist, habe ich auch schon bekommen, als ich hier
und da die Reihenfolge der FFTs veraendert habe, den Fehler bekomme ich
aber nicht, wenn die vorgesehene Reihenfolge an FFTs eingehalten wird (und
trotzdem riesige Zahlen im Imaginaerteil stehen).

Die eingehende Ladungsverteilung kann ich mir aus unerklaerlichen Gruenden
nicht waehrend der Kraftberechnung ausgeben. Ich habe mir daher die
Ruecktrafo angeschaut. Hier sind am "Ende" des Arrays grosse Werte im
Imaginaerteil (25) gegenueber dem Realteil (150). Der Realteil sieht aber
ansonsten so aus wie ich das erwarten wuerde. Zum Testen habe ich eine
einzelne Ladung simuliert.

Kann es sein das Real und Imaginaerteil nicht nebeneinander stehen, also
gerade real, ungerade imaginaer, wie das in der p3m_calc_kspace_forces()
angedeutet ist??

Viele Gruesse,
Tobias

On Thu, April 30, 2015 21:59, Axel Arnold wrote:
> Hallo Tobias,
>
> wie groß relativ zu den eigentlichen Beiträgen ist das denn? Eigentlich
> checkt der Code auch, dass im Imaginärteil nix drin steht, da gibt es
> eine extra Warnung, schau da mal. Und die hat in der Vergangenheit immer
> mal wieder in seltenen Fällen angeschlagen, wir haben aber nie
> verstanden, warum.
>
> Ansonsten kannst Du mal schauen, ob die einlaufende Ladungsverteilung auch
> brav reell ist, sonst wäre das kein Wunder, aber auch kein Problem - die
> FT ist bezüglich der beiden Komponenten unabhängig, man kann sogar im
> Prinzip in Real- und Imaginärteil zwei Ladungsverteilungen gleichzeitig
> berechnen, was man beim Interlacing ausnutzt.
>
> Viele Grüße,
> Axel
>
>
>> Am 30.04.2015 um 15:03 schrieb Tobias Richter
>> <address@hidden <mailto:address@hidden>>:
>>
>> Hello,
>>
>> I am currently trying to repair the contributions to the stress tensor,
>> arising from the longrange electrostatics, so the k-space part.
>> For this I essentially need the forces on the mesh, or more precisely
>> the
>> electric field, since I evaluate the Maxwellian stress tensor.
>>
>> The problem I'm facing now is, that some components in the electric
>> field
>> still have non negligible imaginary components, which should
>> theoretically
>> be zero, since the acting force has to be real valued. I find those
>> nonzero imaginary parts in the p3m_calc_kspace_forces function, before
>> the
>> forces get assigned to the grid, so when the mesh should contain the
>> electric field, which (as I believe) should be strictly real valued.
>>
>> So is this just my misunderstanding or is there something wrong in the
>> calculation?
>>
>> Kind regards,
>> Tobias
>>
>
> ------------------------------------------------
> Dr. Axel Arnold
> ICP, Universität Stuttgart
> Allmandring 3
> 70569 Stuttgart, Germany
> Email: address@hidden <mailto:address@hidden>
> Phone: +49 711 685 67609
>
>
>
> --------------------------------------------
> Axel Arnold
> Martha-Schmidtmann-Str. 7
> 70374 Stuttgart, Germany
> Email: address@hidden
> Phone: +49 173 870 6659
>
>