Thanks again for your prompt and helpful reply.
I just call the corresponding routines and there is no compilation error
I am running a Ubuntu machine with gfortran. I enclosed you the makefile.in
The question is that I am trying to replicate in FEAP the benchmark example of calculating the inverse found in
https://support.nag.com/numeric/nl/nagdoc_26/nagdoc_fl26/pdf/f07/f07ajf.pdffor doing so, I am creating a dummy user element where I create the entrey matrix aaux(4,4)
when I try to compute the inverse via
Call dgetrf(n,n,a,lda,ipiv,info)
Call dgetri(n,a,lda,ipiv,work,lwork,info)
I GET a complete different result in comparison with the becnhmark.
I am not able to see where I made the mistake
Enclosed you have the makefile.in, the I1tension (input) and the user element (elm01.f)