Segfault on call to tgsyl [ solver for Sylvester equation]

[Dharma]

Hello I have written a code for solving a quadratic matrix equation. I get a segmentation fault when i run the code. Can you please help me fix this. I ran the code using intel debugger, I get the following information

Program received signal SIGSEGV
mkl_lapack_ctgsyl ()

I am attaching the code and a make file to test it. QuadMatEqnSolve.F90
Makefile

Thanks
Reddy

[barragan_villanueva_]

Reddy,

In your Makefile there is inconsistency related to using LP64 MKL libraries (see MKL_LIBS_Static) but ILP64 Fortran-mod files (see -I$(MKL_INC)/intel64/ilp64)

But I can still see SegFault:
Program received signal SIGSEGV, Segmentation fault.
0x000000000041a564 in mkl_lapack_ctgsyl ()
(gdb) bt
#0 0x000000000041a564 in mkl_lapack_ctgsyl ()
#1 0x0000000000412968 in ctgsyl_ ()
#2 0x00000000004126e1 in ctgsyl_f95_ ()
#3 0x000000000040e29d in solve_mp_solvequadmateqn_ ()
#4 0x0000000000410ad7 in testsolvequadmateqn () at QuadMatEqnSolve.F90:119
#5 0x0000000000409cbc in main ()

So, it needs more inverstigation to comply with MKL functions interface

[Dharma]

Hello Victor,
The code works but does not converge to a solution for, ijob=1,2,3,4 but when i use ijob=0, i get a segfault. Can you suggest me a way to debug.

Thanks
Reddy