ListGetElementRealVec yields output of zeroes
Posted: 20 Nov 2023, 19:30
This function is to be found in file Lists.F90.
Its 2nd parameter, ngp, is the number of Gaussian integration points. ngp is a function of the element type. It will thus vary as I loop over all active elements in the mesh. Now when I call this function, setting this 2nd parameter to a constant equal to say the maximum of ngp over all element types then there is no problem and the result will be whatever the handle points to in the SIF. This is the desired behaviour. On the other hand, if I set it to the ngp of the current element being scanned, (the normal way to use the function obviously), then initially everything is OK until the ngp suddenly increases for a certain element , then the result will be all 0's totalling in number to ngp . This will happen again for all elements with that same greater value of ngp.
As this was never reported before, it may date back to the earliest version of ListGetElementRealVec.
It you need more details, I will gladly provide the simplest case that displays this issue.
Its 2nd parameter, ngp, is the number of Gaussian integration points. ngp is a function of the element type. It will thus vary as I loop over all active elements in the mesh. Now when I call this function, setting this 2nd parameter to a constant equal to say the maximum of ngp over all element types then there is no problem and the result will be whatever the handle points to in the SIF. This is the desired behaviour. On the other hand, if I set it to the ngp of the current element being scanned, (the normal way to use the function obviously), then initially everything is OK until the ngp suddenly increases for a certain element , then the result will be all 0's totalling in number to ngp . This will happen again for all elements with that same greater value of ngp.
As this was never reported before, it may date back to the earliest version of ListGetElementRealVec.
It you need more details, I will gladly provide the simplest case that displays this issue.