Skip to content

Division by zero in Landau-Zener subroutine #48

@danielhollas

Description

@danielhollas

When running the new LZ_ENE test locally, I noticed this message:

Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG

To find out more, I added the -ffpe-trap=invalid compiler flag to FFLAGS in make.vars.
After recompilation, the LZ_ENE test failed with this message:

Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.

Backtrace for this error:
#0  0x2ac4bbea419f in ???
#1  0x52dab6 in __mod_lz_MOD_lz_hop
	at /home/hollas/ABIN-DEV/landau_zener.F90:173
#2  0x4f48fd in abin_dyn
	at /home/hollas/ABIN-DEV/abin.F90:230
#3  0x4f5cb3 in main
	at /home/hollas/ABIN-DEV/abin.F90:20

We probably need some if statement to guard against dividing by zero in landau_zener.F90:173

Metadata

Metadata

Assignees

Labels

Type

No type
No fields configured for issues without a type.

Projects

No projects

Relationships

None yet

Development

No branches or pull requests

Issue actions