Hello,
In my code I have to solve a system of N-2 linear equations K+1 times (with k<10). What this means is that my matrix A (N-2*N-2) will be structurally the same for all K, only with different values. Initially I factored matrix A (with DSS_FACTOR_COMPLEX) every time I needed to solve the system for a different K. But since K<10, I was hoping to allocate the handle ( type(MKL_DSS_HANDLE), allocatable :: handle(:) -- allocate(handle(0:K)) ) and use DSS_FACTOR only once for every K and then use DSS_SOLVE_COMPLEX(handle(k), MKL_DSS_DEFAULTS, RHS, 1, solutionout) to solve my system. It's also important to understand that the code then solves all K+1 systems thousands of times.
Unfortunately, the code exits with:
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
ruth6.exe 0000000110171F6E for__signal_handl Unknown Unknown
libsystem_platfor 00007FFF8AA7F52A _sigtramp Unknown Unknown
ruth6.exe 000000010F97FAD6 mkl_spblas_lp64_a Unknown Unknown
The most interesting here is that it actually solves all K+1 systems a few times (last time solved ~300 times) before exiting with this error.
Please help!!