A puzzling segmentation fault occurs when both "embed" and "dimensions" are used together

[SigmaYot Zeta]

Hi,

I am using the immersed boundary method (IBM) to simulate flow within particles gaps. However, the computational domain is elongated (see Figure 1), so I set the x, y, z ratio to 1:1:10 using dimensions function.

I then read the STL, initialized the distance field, initialized the phi values, and the cs and fs values ​​of the  IBM, and finally counted the number of meshes. The program logs indicate that all of these steps were successful. However, the program then immediately encountered a segmentation fault.

This is very confusing. The source code is attached.

Sincerely,

Siyu Zou

[Stephane Popinet]

Dear Siyu,

I have just pushed a patch which should fix the problem:

https://basilisk.fr/src/?changes=20260618090126

cheers,

Stephane