Segmentation Fault in Elemsd_SetJs on GCC 16.x.x

[shishiousan]

The part of error message is given below.

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x7f77e32249a6 in ???
#1  0x7f77e3223555 in ???
#2  0x7f77dc43e8ef in ???
#3  0x7f77edb58e31 in __elemshapedata_setmethods_MOD_elemsd_setjs
	at /home/shion/Dropbox/easifem/base/src/submodules/ElemshapeData/src/ElemshapeData_SetMethods@Methods.F90:104
#4  0x7f77edb4c787 in __elemshapedata_setmethods_MOD_elemsd_set1
	at /home/shion/Dropbox/easifem/base/src/submodules/ElemshapeData/src/ElemshapeData_SetMethods@Methods.F90:237
#5  0x7f77ef888e4b in __abstract1dsdmixedfem_class_MOD_obj_setelemsdforspace

The error occurs here. (ElemshapeData_SetMethods@Methods.F90:104)

  DO ips = 1, obj%nips
    obj%js(ips) = Det(obj%jacobian(1:obj%nsd, 1:obj%xidim, ips))
  END DO

My source code where the error occur is this.

TYPE(ElemShapeData_), Pointer :: elemsd, linElemsd

CALL ElemsdSet(obj=elemsd, val=xij(1:1, :), N=linElemsd%N(1:2, 1:nips), &
               dNdXi=linElemsd%dNdXi(1:2, 1:1, 1:nips))

This error occurs only when the base and other libraries are compiled with gcc 16.x.x.

I checked the all components of elemsd in my source code as follows

TYPE(ElemShapeData_), Pointer :: elemsd, linElemsd
 
CALL Display(elemsd%nips, "nips:: ")  ! it returns 3 in my case. 
CALL Display(elemsd%jacobian(1:elemsd%nsd, 1:elemsd%xidim,  1:elemsd%nips), "jacobian :: ")
CALL Display(elemsd%js(1:elemsd%nips), "js :: ")

!! will error occur ?
CALL Display(Det(elemsd%jacobian(1:elemsd%nsd, 1:elemsd%xidim, 1)), "det1::")
CALL Display(Det(elemsd%jacobian(1:elemsd%nsd, 1:elemsd%xidim, 2)), "det2::")
CALL Display(Det(elemsd%jacobian(1:elemsd%nsd, 1:elemsd%xidim, 3)), "det3::")

!! CALL ElemsdSet(obj=elemsd, val=xij(1:1, :), N=linElemsd%N(1:2, 1:nips), &
               dNdXi=linElemsd%dNdXi(1:2, 1:1, 1:nips))

In this raw code, no error occurs.

[vickysharma0812]

@shishiousan Can you post here your code which causes the error.

[shishiousan]

@vickysharma0812

MODULE PROCEDURE obj_SetElemsdForSpace
CHARACTER(*), PARAMETER :: myName = "obj_SetElemsdForSpace()"
INTEGER(I4B) :: nips, nns, cellOrder(1), order, cellOrient(1), ii

TYPE(ElemShapeData_), POINTER :: elemsd => NULL(), linElemsd => NULL()
REAL(DFP) :: refElemCoord(1, 2)
LOGICAL(LGT) :: isok

#ifdef DEBUG_VER
CALL e%RaiseInformation(modName//'::'//myName//' - '// &
                        '[START] ')
#endif

#ifdef DEBUG_VER
isok = ivar .LE. defaultOpt%MAX_NUM_FIELD
CALL AssertError1(isok, myName, "ivar must be less than MAX_NUM_FIELD")
#endif

elemsd => obj%elemsdForSpace(ivar)%ptr
linElemsd => obj%linElemsdForSpace(ivar)%ptr

order = obj%spaceOrder(ivar)%ptr%val(spaceElemNum)
nips = Quad_Size(obj%quadForSpace, 2)
nns = obj%totalDOFSpace(ivar)%ptr%val(spaceElemNum)
refElemCoord(1, 1) = math%minus_one
refElemCoord(1, 2) = math%one
cellOrient = math%one_i

CALL ElemsdAllocate(obj=elemsd, nsd=math%one_i, xidim=math%one_i, &
                    nns=nns, nips=nips)

CALL LagrangeElemShapeData(obj=linElemsd, &
                           quad=obj%quadForSpace, &
                           nsd=elemsd%nsd, &
                           xidim=elemsd%xidim, &
                           elemtype=elem%line, &
                           refelemCoord=refelemCoord, &
                           domainName="B", &
                           order=math%one_i)

SELECT CASE (obj%baseInterpolationForSpace)
CASE ("LAGR")

  CALL LagrangeElemShapeData(obj=elemsd, &
                             quad=obj%quadForSpace, &
                             nsd=elemsd%nsd, &
                             xidim=elemsd%xidim, &
                             elemtype=elem%line, &
                             refelemCoord=refelemCoord, &
                             domainName="B", &
                             order=order, &
                             ipType=obj%ipTypeForSpace, &
                             basisType=obj%baseTypeForSpace)

  obj%spaceShapeFuncBndy(ivar, 1:nns, 1) = math%zero
  obj%spaceShapeFuncBndy(ivar, 1, 1) = math%one

  obj%spaceShapeFuncBndy(ivar, 1:nns, 2) = math%zero
  obj%spaceShapeFuncBndy(ivar, 2, 2) = math%one

CASE ("HIER", "HEIR")

  cellOrder = order
  CALL HierarchicalElemShapeData(obj=elemsd, &
                                 quad=obj%quadForSpace, &
                                 nsd=elemsd%nsd, &
                                 xidim=elemsd%xidim, &
                                 elemtype=elem%line, &
                                 refelemCoord=refelemCoord, &
                                 domainName="B", &
                                 cellOrder=cellOrder, &
                                 cellOrient=cellOrient)

  obj%spaceShapeFuncBndy(ivar, 1:nns, 1) = math%zero
  obj%spaceShapeFuncBndy(ivar, 1, 1) = math%one

  obj%spaceShapeFuncBndy(ivar, 1:nns, 2) = math%zero
  obj%spaceShapeFuncBndy(ivar, 2, 2) = math%one

CASE DEFAULT
  CALL e%RaiseError(modName//'::'//myName//' - '// &
                    'no case found for baseInterpolationForSpace')

END SELECT

CALL ElemsdSet(obj=elemsd, val=xij(1:1, :), N=linElemsd%N(1:2, 1:nips), &
               dNdXi=linElemsd%dNdXi(1:2, 1:1, 1:nips))

#ifdef DEBUG_VER
CALL e%RaiseInformation(modName//'::'//myName//' - '// &
                        '[END] ')
#endif

END PROCEDURE obj_SetElemsdForSpace

Error occurs at the bottom

[vickysharma0812]

Did you perform following after updating the gfortran:

easifem clean base classes
easifem install base classes --no-download?

[vickysharma0812]

@shishiousan Also, can you nullify elemsd and linElemsd before exiting obj_SetElemsdForSpace method

[vickysharma0812]

@shishiousan Can you check if the following pointers are associated?

elemsd => obj%elemsdForSpace(ivar)%ptr
linElemsd => obj%linElemsdForSpace(ivar)%ptr

[shishiousan]

@vickysharma0812

Did you perform following after updating the gfortran:

easifem clean base classes easifem install base classes --no-download?

Yes I did.

@shishiousan Also, can you nullify elemsd and linElemsd before exiting obj_SetElemsdForSpace method

I added. But the same error.

@shishiousan Can you check if the following pointers are associated?

elemsd => obj%elemsdForSpace(ivar)%ptr linElemsd => obj%linElemsdForSpace(ivar)%ptr

Both are associated. I checked by the following

CALL Display(ASSOCIATED(elemsd), "associated(elemsd) :: ")   ! return TRUE
CALL Display(ASSOCIATED(linElemsd), "associated(linelemsd) :: ") ! return TRUE

[vickysharma0812]

@shishiousan Does your code compiles if you comment out the ElemsdSet method. If it is true, then please check the bounds related to SetJs method, where following code is used.

  DO ips = 1, obj%nips
    obj%js(ips) = Det(obj%jacobian(1:obj%nsd, 1:obj%xidim, ips))
  END DO

So please check the size of elemsd%js, size(elemsd%jacobian, 3) and elemsd%nips.

[shishiousan]

@vickysharma0812 This error occurred when the part is commented out.

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x7f5beaa249a6 in ???
#1  0x7f5beaa23555 in ???
#2  0x7f5be3c3e8ef in ???
#3  0x7f5bf535692d in __elemshapedata_setmethods_MOD_elemsd_setdndxt
	at /home/shion/Dropbox/easifem/base/src/submodules/ElemshapeData/src/ElemshapeData_SetMethods@Methods.F90:129
#4  0x7f5bf534c796 in __elemshapedata_setmethods_MOD_elemsd_set1
	at /home/shion/Dropbox/easifem/base/src/submodules/ElemshapeData/src/ElemshapeData_SetMethods@Methods.F90:238
#5  0x7f5bf7088e4b in __abstract1dsdmixedfem_class_MOD_obj_setelemsdforspace

This is related to inv of jacobian

DO ips = 1, obj%nips
  CALL Inv(InvA=invJacobian, A=obj%jacobian(1:obj%nsd, 1:obj%nsd, ips)) !!! HERE

  obj%dNdXt(1:obj%nns, 1:obj%nsd, ips) = &
    MATMUL(obj%dNdXi(1:obj%nns, 1:obj%nsd, ips), &
           invJacobian(1:obj%nsd, 1:obj%nsd))
END DO

I already checked the size of everything.

CALL Display(SIZE(elemsd%js), "size of js :: ") ! return 3
CALL Display(SIZE(elemsd%jacobian, 3), "size of jacobian :: ") !return 3
CALL Display(elemsd%nips, "size of nips :: ") ! return 3

[shishiousan]

I could not check these sizes in base routine since the routines are pure subroutines.
But at least, in my code, elemsd has correct sizes.

[vickysharma0812]

@vickysharma0812

MODULE PROCEDURE obj_SetElemsdForSpace
CHARACTER(*), PARAMETER :: myName = "obj_SetElemsdForSpace()"
INTEGER(I4B) :: nips, nns, cellOrder(1), order, cellOrient(1), ii

TYPE(ElemShapeData_), POINTER :: elemsd => NULL(), linElemsd => NULL()
REAL(DFP) :: refElemCoord(1, 2)
LOGICAL(LGT) :: isok

#ifdef DEBUG_VER
CALL e%RaiseInformation(modName//'::'//myName//' - '// &
'[START] ')
#endif

#ifdef DEBUG_VER
isok = ivar .LE. defaultOpt%MAX_NUM_FIELD
CALL AssertError1(isok, myName, "ivar must be less than MAX_NUM_FIELD")
#endif

elemsd => obj%elemsdForSpace(ivar)%ptr
linElemsd => obj%linElemsdForSpace(ivar)%ptr

order = obj%spaceOrder(ivar)%ptr%val(spaceElemNum)
nips = Quad_Size(obj%quadForSpace, 2)
nns = obj%totalDOFSpace(ivar)%ptr%val(spaceElemNum)
refElemCoord(1, 1) = math%minus_one
refElemCoord(1, 2) = math%one
cellOrient = math%one_i

CALL ElemsdAllocate(obj=elemsd, nsd=math%one_i, xidim=math%one_i, &
nns=nns, nips=nips)

CALL LagrangeElemShapeData(obj=linElemsd, &
quad=obj%quadForSpace, &
nsd=elemsd%nsd, &
xidim=elemsd%xidim, &
elemtype=elem%line, &
refelemCoord=refelemCoord, &
domainName="B", &
order=math%one_i)

SELECT CASE (obj%baseInterpolationForSpace)
CASE ("LAGR")

CALL LagrangeElemShapeData(obj=elemsd, &
quad=obj%quadForSpace, &
nsd=elemsd%nsd, &
xidim=elemsd%xidim, &
elemtype=elem%line, &
refelemCoord=refelemCoord, &
domainName="B", &
order=order, &
ipType=obj%ipTypeForSpace, &
basisType=obj%baseTypeForSpace)

obj%spaceShapeFuncBndy(ivar, 1:nns, 1) = math%zero
obj%spaceShapeFuncBndy(ivar, 1, 1) = math%one

obj%spaceShapeFuncBndy(ivar, 1:nns, 2) = math%zero
obj%spaceShapeFuncBndy(ivar, 2, 2) = math%one

CASE ("HIER", "HEIR")

cellOrder = order
CALL HierarchicalElemShapeData(obj=elemsd, &
quad=obj%quadForSpace, &
nsd=elemsd%nsd, &
xidim=elemsd%xidim, &
elemtype=elem%line, &
refelemCoord=refelemCoord, &
domainName="B", &
cellOrder=cellOrder, &
cellOrient=cellOrient)

obj%spaceShapeFuncBndy(ivar, 1:nns, 1) = math%zero
obj%spaceShapeFuncBndy(ivar, 1, 1) = math%one

obj%spaceShapeFuncBndy(ivar, 1:nns, 2) = math%zero
obj%spaceShapeFuncBndy(ivar, 2, 2) = math%one

CASE DEFAULT
CALL e%RaiseError(modName//'::'//myName//' - '// &
'no case found for baseInterpolationForSpace')

END SELECT

CALL ElemsdSet(obj=elemsd, val=xij(1:1, :), N=linElemsd%N(1:2, 1:nips), &
dNdXi=linElemsd%dNdXi(1:2, 1:1, 1:nips))

#ifdef DEBUG_VER
CALL e%RaiseInformation(modName//'::'//myName//' - '// &
'[END] ')
#endif

END PROCEDURE obj_SetElemsdForSpace

Error occurs at the bottom

@shishiousan Which step in the above code causes error