[CIG-SEISMO] Floating-point exception - erroneous arithmetic operation
komatitsch at lma.cnrs-mrs.fr
Wed May 27 07:40:46 PDT 2015
It seems you simply forgot to also divide the time step by two to remain
stable when you increased mesh resolution by two:
*** Max stability for P wave velocity (must be below about 0.50 or so)
On 05/27/2015 03:36 PM, Moritz Wehler wrote:
> Hey there,
> I have some issues with increasing the number of elements along the
> z-axis with specfem2D.
> As soon as i double the number of elements from 111 (in total along the
> z-axis) to 222 (in total), I get the following error message:
> "Program received signal SIGFPE: Floating-point exception - erroneous
> arithmetic operation.
> Backtrace for this error:
> #0 0x7F62A8D0D7D7
> #1 0x7F62A8D0DDDE
> #2 0x7F62A8449D3F
> #3 0x433DBE in compute_forces_viscoelastic_
> #4 0x4DF021 in MAIN__ at specfem2D.F90:?
> Gleitkomma-Ausnahme (Speicherabzug geschrieben)"
> I should add that everything worked perfectly fine BEFORE I increased
> the element number. The corresponding par_file and interface_file are
> Hope I didn't forget anything important and that someone can help me.
> Cheers, Mo
> CIG-SEISMO mailing list
> CIG-SEISMO at geodynamics.org
CNRS Research Director (DR CNRS), Laboratory of Mechanics and Acoustics,
UPR 7051, Marseille, France http://komatitsch.free.fr
More information about the CIG-SEISMO