Dear all,
sorry but It's not clear for me. What we have seen is that simulations with particles >1000 gave wrong results. In fact only one particle survives after the first steps, even with very similar initial conditions. With a particle number <1000 the behavior is what we have expected. So, following your emails, I have different questions:
- Do you mean that with "macroparticles" <1000, the SC is off and when you move to simulations with a particle number >1000 the SC is automatically on? Is there any way to track a high number of particles with SC off? It could be interesting to analyze the SC effects
- If there were any wrong units somewhere we would have obtained some wrong results even with a low number of particles, wouldn't we?
Cheers
Conchi
El sáb, 11-02-2017 a las 14:58 +0000, Daniel Winklehner escribió:
I remember having similar errors when using the wrong units somewhere (geometry, particle distribution, …).
Could it be that simple? Can you send me the link to the Dropbox folder with the input files?
Best,
Daniel
From: opal-request AT lists.psi.ch [mailto:opal-request AT lists.psi.ch] On Behalf Of Christof Metzger-Kraus
Sent: Saturday, February 11, 2017 03:50
To: Adelmann Andreas (PSI) <andreas.adelmann AT psi.ch>
Cc: Oliver Amoros, Concepcion <concepcion.oliver AT ciemat.es>; Calvo Portela, Pedro <Pedro.Calvo AT ciemat.es>; opal <Opal AT lists.psi.ch>
Subject: Re: [Opal] Error OPAL
We now know why we have a different behavior when run with 999 particles than with 1000. But the main problem isn't solved: why does OPAL crash when SC is on? It looks to me like the mesh isn't adapted to the particle distribution when computing SC. boundp missing or not at the correct location?
cheers,
christof
On Sat, Feb 11, 2017 at 8:38 AM, Adelmann Andreas (PSI) <andreas.adelmann AT psi.ch> wrote:
Hi is the problem solved/understood?
Andreas
> On 10 Feb 2017, at 19:12, Oliver Amoros, Concepcion <concepcion.oliver AT ciemat.es> wrote:
>
> Hello,
> We use matlab by extracting the informatíon in the .h5 file.
> Conchi
> ________________________________________
> De: opal-request AT lists.psi.ch [opal-request AT lists.psi.ch] en nombre de Adelmann Andreas (PSI) [andreas.adelmann AT psi.ch]
> Enviado el: jueves, 09 de febrero de 2017 12:05
> Para: Calvo Portela, Pedro
> Cc: opal
> Asunto: Re: [Opal] Error OPAL
>
> How do you plot this?
>
> A
> ------
> Dr. sc. math. Andreas (Andy) Adelmann
> Staff Scientist
> Paul Scherrer Institut WBBA/219 CH-5232 Villigen PSI
> Phone Office: xx41 56 310 42 33 Fax: xx41 56 310 31 91
> Phone Home: xx41 62 891 91 44
> -------------------------------------------------------
> Friday: ETH CAB F37 or HPK G 28 +41 44 632 75 22
> ============================================
> The more exotic, the more abstract the knowledge,
> the more profound will be its consequences.
> Leon Lederman
> ============================================
>
> On 9 Feb 2017, at 11:47, Calvo Portela, Pedro <Pedro.Calvo AT ciemat.es<mailto:Pedro.Calvo AT ciemat.es>> wrote:
>
> Dear Andreas
>
> I've include the modification at the Field Solver, and now the error doesn't occur. But the tracking of the particles are absolutely unexpected. I attach two plots of the orbits with 1000 and 999 particles. You could see that the results are complete different.
>
> Best regards
> Pedro
>
> ________________________________
> De: Adelmann Andreas (PSI) [andreas.adelmann AT psi.ch<mailto:andreas.adelmann AT psi.ch>]
> Enviado el: miércoles, 08 de febrero de 2017 21:15
> Para: Calvo Portela, Pedro
> Asunto: Re: [Opal] Error OPAL
>
> Hi Pedro, try to change:
>
>
> FSTYPE=FFT, MX=32, MY=32, MT=64
>
> to
>
> FSTYPE=FFT, MX=16, MY=16, MT=32, ….
>
> You only have 1000 particles and apps. 65k grid cells.
>
> How many cores do you use?
>
> Andreas
>
>
>
> On 31 Jan 2017, at 08:58, Calvo Portela, Pedro <Pedro.Calvo AT ciemat.es<mailto:Pedro.Calvo AT ciemat.es>> wrote:
>
> Dear Andreas.
>
> I'm using the version 1.4.0 of OPAL.
> I share with you the input files in a dropbox folder.
>
> Cheers
> Pedro
>
> ________________________________
> De: Adelmann Andreas (PSI) [andreas.adelmann AT psi.ch<mailto:andreas.adelmann AT psi.ch>]
> Enviado el: viernes, 27 de enero de 2017 23:35
> Para: Calvo Portela, Pedro
> Cc: opal AT lists.psi.ch<mailto:opal AT lists.psi.ch>
> Asunto: Re: [Opal] Error OPAL
>
> Hi Pedro
>
> can you make all the input files available so we can reproduce the problem.
> Which version of OPAL are you using?
> Cheers Andreas
>
> On 27 Jan 2017, at 11:03, Calvo Portela, Pedro <Pedro.Calvo AT ciemat.es<mailto:Pedro.Calvo AT ciemat.es>> wrote:
>
> Dear OPAL users, I write you because we have some problems about the simulation of a big number of particles.
> We simulate without problems our machine with the correct electric and magnetic fields, including the geometry in the simulation. But when we increase the number of particles above 1000, we obtain the following error in the output:
>
> Error> Interpolator::getFieldIter: attempt to access non-local index{[-2147483648:-2147483648:1],[-2147483648:-2147483648:1],[-2147483648:-2147483648:1]} on node 0
> Error> Dumping local owned and allocated domains:
> Error> 0: owned = {[0:31:1],[0:31:1],[0:31:1]}, allocated = {[-1:32:1],[-1:32:1],[-1:32:1]}
> Error> Error occurred for BareField with layout = Domain = {[0:31:1],[0:31:1],[0:31:1]}
> Error> FieldLayoutUsers = 3
> Error> Total number of vnodes = 1
> Error> Local Vnodes = 1
> Error> vnode 0: Node = 0 ; vnode_m = -1 ; Domain = {[0:31:1],[0:31:1],[0:31:1]}
> Error> Remote Vnodes = 0
> Error>
> Error> Calling abort ...
>
> This error don't occur when the number of particles is lower than 1000.
> We have tried to solve this error changing the fieldsolver in the code for a 1000 particles, but the results are completly diferent if we compare with the simulation with a number of particles of 999. I
>
> I will appreciate any suggestion to solve this problem.
> Best Regards
> Pedro Calvo
>
>
> ----------------------------
> Confidencialidad:
> Este mensaje y sus ficheros adjuntos se dirige exclusivamente a su destinatario y puede contener información privilegiada o
> confidencial. Si no es vd. el destinatario indicado, queda notificado de que la utilización, divulgación y/o copia sin autorización está
> prohibida en virtud de la legislación vigente. Si ha recibido este mensaje por error, le rogamos que nos lo comunique
> inmediatamente respondiendo al mensaje y proceda a su destrucción.
>
> Disclaimer:
> This message and its attached files is intended exclusively for its recipients and may contain confidential information. If you received
> this e-mail in error you are hereby notified that any dissemination, copy or disclosure of this communication is strictly prohibited and
> may be unlawful. In this case, please notify us by a reply and delete this email and its contents immediately.
> ----------------------------
>
>
> <orbits_999part_161632.png><orbits_1000part_161632.png>
>
|
