ComputeMatrixElements
Si EvanescentJ n'est pas activé, jmin et jmax valent 0 et ns.
Du coup les j+1 et j-1 de:
real gradertm = (enrad[ljp]-enrad[ljm])/(Rmed[i]*sinPhiMed[h]*(thetaMed[j+1]-thetaMed[j-1]));
posent problème.
Je modifie en tenant compte que c'est cyclique ?