Segmentation fault in

I am learning to use Bempp-cl. Wen I ran the code of tutorial example “Electromagnetic scattering from flat screens”, I got segmentation fault (core dumped) at the statement

incident_field_data = incident_field(points)

By looking at the function “incident_field(x)”, I can’t see what the cause of the error is. Is FMM used in this example? I failed to get exafmm-t installed on my Ubuntu 20.04.

Thank you for any helps.

Best regards,

I need to correct myself. By debugging, I found that the code line that causes segmentation fault is

scattered_field_data = -slp_pot * lambda_data

The solver succeeded but calculation of scattered electric field failed. If a small number of elements is used, the program runs fine. The largest number of edges that can let the program succeeds is about 2000, My machine has 32 GB RAM, so 2000 edges should not cause any problems.

Could you anybody help to resolve this issue? I hope bempp could be a useful tool for my project.

Thank you,


Before changing to OpenCl, I could not perform operations with many edges (maybe around thousands) since it would fail the execution (typically at this sort of lines that you indicate).

Then, I installed OpenCL (for windows: OpenCl in Windows 10), and now I can run with many more adges and in faster execution time.



Also, it helps reducing the number of points indicated to compute in : slp_pot(…,NumPoints…)

Hi Marc,

Thank you for the info. I indeed install openCL on my Ubuntu. I’ll check if openCL is enabled.