topic I found the issue, I had an in IntelĀ® oneAPI Math Kernel Library
https://community.intel.com/t5/InteloneAPIMathKernelLibrary/SIGSEGVwithPoissonSolver/mp/1121488#M24994
<P>I found the issue, I had an error in my makefile. I was not adding the flags when compiling the file that had the calls to the Helmoltz functions. </P>
<P> </P>
Mon, 09 May 2016 19:42:55 GMT
Camilo_P_
20160509T19:42:55Z

SIGSEGV with Poisson Solver
https://community.intel.com/t5/InteloneAPIMathKernelLibrary/SIGSEGVwithPoissonSolver/mp/1121487#M24993
<P>I have been running and debugging a code with the Intel visual fortran composer that uses the MKL poisson solver. It runs without issues in Windows. But when I run it in a UNIX system I get a SIGSEGV error. I am compiling it in 64 bits and I have checked that that is happening. I have tried: heaparrays, ulimit s unlimited and none of them have fixed the problem. I used traceback and it points to the line where I do </P>
<PRE class="brush:fortran;">call D_COMMIT_HELMHOLTZ_3D (f, bd_ax, bd_bx, bd_ay, bd_by, bd_az , bd_bz, xhandle, yhandle, ipar, dpar, stat)</PRE>
<P>The subroutine where this call is:</P>
<PRE class="brush:fortran;"> subroutine solvePoisson3DVorticityPerMKL(f, bd_ax, bd_bx , bd_ay, bd_by, bd_az, bd_bz, nx, ny, nz, ax, bx, ay, by, az, bz, i)
double precision dpar(5*(nx+ny)/2+9)
integer ipar(128), i, stat, nx, ny, nz
double precision :: bd_ax(:,:), bd_bx(:,:) , bd_ay(:,:), bd_by(:,:), bd_az(:,:), bd_bz(:,:)
double precision :: f(:,:,:), finish, start
character(6) BCtype
double precision q, ax, bx, ay, by, az, bz
type(DFTI_DESCRIPTOR), pointer :: xhandle, yhandle
SELECT CASE (i) ! select component
CASE (1)
BCtype = 'PPDDPP'
#ifndef SHEAR_FLOW
bd_ay =0
bd_by =0
#endif
case (2)
BCtype = 'PPNNPP'
case (3)
BCtype = 'PPDDPP'
#ifndef SHEAR_FLOW
print *,"Unidirectional flow"
bd_ay = 0.1*(bxax)/2 ! for unidirectional flow
bd_by =0.1*(bxax)/2
#endif
end select
q=0.0
ipar =0
call D_INIT_HELMHOLTZ_3D(ax, bx, ay, by, az, bz, nx, ny, nz, BCtype , q, ipar, dpar, stat)
ipar(3)=0 ! dissable warnings
call D_COMMIT_HELMHOLTZ_3D (f, bd_ax, bd_bx, bd_ay, bd_by, bd_az , bd_bz, xhandle, yhandle, ipar, dpar, stat)
call D_HELMHOLTZ_3D(f, bd_ax, bd_bx, bd_ay, bd_by, bd_az , bd_bz, xhandle, yhandle, ipar, dpar, stat)
call free_Helmholtz_3D(xhandle, yhandle, ipar, stat)
end subroutine solvePoisson3DVorticityPerMKL</PRE>
<P>The error :</P>
<P>forrtl: severe (174): SIGSEGV, segmentation fault occurred<BR />
Image PC Routine Line Source <BR />
CoupledSystem 0000000000519B15 Unknown Unknown Unknown<BR />
CoupledSystem 00000000005178D7 Unknown Unknown Unknown<BR />
CoupledSystem 00000000004CBB94 Unknown Unknown Unknown<BR />
CoupledSystem 00000000004CB9A6 Unknown Unknown Unknown<BR />
CoupledSystem 000000000047F4B6 Unknown Unknown Unknown<BR />
CoupledSystem 0000000000482AD0 Unknown Unknown Unknown<BR />
libpthread.so.0 00002B3F0728F130 Unknown Unknown Unknown<BR />
libmkl_def.so 00002B3F095743E4 Unknown Unknown Unknown<BR />
libmkl_core.so 00002B3F05D48E59 Unknown Unknown Unknown<BR />
libmkl_intel_ilp6 00002B3F0475F7CF Unknown Unknown Unknown<BR />
CoupledSystem 0000000000447A77 vorticitysolverut 380 FluidUtilities.f90<BR />
CoupledSystem 00000000004579CC vorticitysolverut 2596 FluidUtilities.f90<BR />
CoupledSystem 000000000045648F vorticitysolverut 2567 FluidUtilities.f90<BR />
CoupledSystem 0000000000403865 MAIN__ 135 main.f90<BR />
CoupledSystem 00000000004030DE Unknown Unknown Unknown<BR />
libc.so.6 00002B3F079D1AF5 Unknown Unknown Unknown<BR />
CoupledSystem 0000000000402FE9 Unknown Unknown Unknown</P>
<P><SPAN style="fontsize: 1em; lineheight: 1.5;">I have allocated the vectors previously. Again the code works without an issue on windows. Thanks. </SPAN></P>
<P> </P>
Mon, 09 May 2016 18:08:00 GMT
https://community.intel.com/t5/InteloneAPIMathKernelLibrary/SIGSEGVwithPoissonSolver/mp/1121487#M24993
Camilo_P_
20160509T18:08:00Z

I found the issue, I had an
https://community.intel.com/t5/InteloneAPIMathKernelLibrary/SIGSEGVwithPoissonSolver/mp/1121488#M24994
<P>I found the issue, I had an error in my makefile. I was not adding the flags when compiling the file that had the calls to the Helmoltz functions. </P>
<P> </P>
Mon, 09 May 2016 19:42:55 GMT
https://community.intel.com/t5/InteloneAPIMathKernelLibrary/SIGSEGVwithPoissonSolver/mp/1121488#M24994
Camilo_P_
20160509T19:42:55Z

Hi, I'm facing the same
https://community.intel.com/t5/InteloneAPIMathKernelLibrary/SIGSEGVwithPoissonSolver/mp/1121489#M24995
<P>Hi, I'm facing the same problem..</P><P>Can you explain how did you solve the problem and which flag have you added?</P>
Mon, 15 Jun 2020 11:43:46 GMT
https://community.intel.com/t5/InteloneAPIMathKernelLibrary/SIGSEGVwithPoissonSolver/mp/1121489#M24995
cavaiola__mattia
20200615T11:43:46Z