Hello,
I am trying to add CPML layers while using the internal mesher. Although Stacey absorbing boundary conditions worked fine, I noticed some reflections at the boundaries. That is why I wanted to try CPML.
However, I am encountering an error when using CPML during the database generation stage:
=================================
Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
Backtrace for this error:
STOP there is an alpha error in mesh of CPML-layer z
[orfoz132:2260164:0:2260164] Caught signal 8 (Floating point exception: floating-point invalid operation)
==== backtrace (tid:2260164) ====
0 /arf/sw/lib/ucx/lib/libucs.so.0(ucs_handle_error+0x294) [0x14dda50a8b44]
1 /arf/sw/lib/ucx/lib/libucs.so.0(+0x35cff) [0x14dda50a8cff]
2 /arf/sw/lib/ucx/lib/libucs.so.0(+0x3610a) [0x14dda50a910a]
3 /lib64/libm.so.6(+0x649c0) [0x14dda610c9c0]
4 /lib64/libm.so.6(pow+0x18) [0x14dda60dd788]
5 ./bin/xgenerate_databases() [0x4715e9]
6 ./bin/xgenerate_databases() [0x429a71]
7 ./bin/xgenerate_databases() [0x49466d]
8 ./bin/xgenerate_databases() [0x42bc09]
9 ./bin/xgenerate_databases() [0x403717]
10 /lib64/libc.so.6(+0x3feb0) [0x14dda5e3feb0]
11 /lib64/libc.so.6(__libc_start_main+0x80) [0x14dda5e3ff60]
12 ./bin/xgenerate_databases() [0x403745]
The mesh is covering an area of 655 x 547 km. Elements used in the XI and ETA direction are equal to 320. I have 60 elements in the z direction for a depth of 100 km. I use 4 doubling layers at 15, 25, 30, and 40. I try adding 20 PML layers with a thickness of 10.d0
Since the error indicates an issue with the CPML layer in the z-direction, I attempted to change the element size along the z-axis and varied the number of extra CPML layers. However, the problem persists.
I would greatly appreciate any guidance or suggestions on how to fix this issue.
Hello,
I am trying to add CPML layers while using the internal mesher. Although Stacey absorbing boundary conditions worked fine, I noticed some reflections at the boundaries. That is why I wanted to try CPML.
However, I am encountering an error when using CPML during the database generation stage:
=================================
Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
Backtrace for this error:
STOP there is an alpha error in mesh of CPML-layer z
[orfoz132:2260164:0:2260164] Caught signal 8 (Floating point exception: floating-point invalid operation)
==== backtrace (tid:2260164) ====
0 /arf/sw/lib/ucx/lib/libucs.so.0(ucs_handle_error+0x294) [0x14dda50a8b44]
1 /arf/sw/lib/ucx/lib/libucs.so.0(+0x35cff) [0x14dda50a8cff]
2 /arf/sw/lib/ucx/lib/libucs.so.0(+0x3610a) [0x14dda50a910a]
3 /lib64/libm.so.6(+0x649c0) [0x14dda610c9c0]
4 /lib64/libm.so.6(pow+0x18) [0x14dda60dd788]
5 ./bin/xgenerate_databases() [0x4715e9]
6 ./bin/xgenerate_databases() [0x429a71]
7 ./bin/xgenerate_databases() [0x49466d]
8 ./bin/xgenerate_databases() [0x42bc09]
9 ./bin/xgenerate_databases() [0x403717]
10 /lib64/libc.so.6(+0x3feb0) [0x14dda5e3feb0]
11 /lib64/libc.so.6(__libc_start_main+0x80) [0x14dda5e3ff60]
12 ./bin/xgenerate_databases() [0x403745]
The mesh is covering an area of 655 x 547 km. Elements used in the XI and ETA direction are equal to 320. I have 60 elements in the z direction for a depth of 100 km. I use 4 doubling layers at 15, 25, 30, and 40. I try adding 20 PML layers with a thickness of 10.d0
Since the error indicates an issue with the CPML layer in the z-direction, I attempted to change the element size along the z-axis and varied the number of extra CPML layers. However, the problem persists.
I would greatly appreciate any guidance or suggestions on how to fix this issue.