Weights are sometimes <0

seismology / mc_kernel

Calculate seismic sensitivity kernels using Axisem

GNU General Public License v3.0

24 stars 8 forks source link

Weights are sometimes <0 #30

Closed sstaehler closed 9 years ago

sstaehler commented 9 years ago

Happens in @kasra-hosseini's huge mesh after 24h runtime, so it is a bit difficult to debug. More of a memo to self to look into it.

ERROR: weights is smaller zero! Check whether point is outside of element
    weight: -0.49999999953433383     
    dx    :   40809.243542678654     
    dy    :  -32979.621791672747     
    dz    :  -69273.540243368130

sstaehler commented 9 years ago

Was caused by integer overflow in halton.f90, which caused negative random numbers, which caused points outside of tetrahedron to be generated, which caused negative weights