Computes the Newton jacobian
INPUT = nodeNb,nonzeroMat, pointMa,flagDir,
diagMass,matMass, diagPoisson,matPoisson,pot, elecTemp,elecDens
- nonzeroMat[nodeNb+1] = index of non null matrix element
- pointMat[ndimMat] = index of non null matrix element
- flagDir[nodeNb] = Dirichlet boundary conditions flag
- diagMass [nodeNb] = diagonal of the Mass matrix
- matMass[ndimMat] = Mass matrix (without the diagonal)
- diagPoisson [nodeNb] = diagonal of the poisson matrix
- matPoisson[ndimMat] = poisson matrix (without the diagonal)
- pot [nodeNb] = potential
- elecTemp = Electronic temperature (eV)
- elecDens = Electronic density (M-3)
- elecTemp2 = Electronic temperature 2nd population [eV]
- elecDens2 = Electronic density 2nd population [m-3]
OUTPUT = diagNewton, matNewton
- diagNewton [nodeNb] = diagonal of the jacobian
- matNewton[ndimMat] = Jacobian matrix (without the diagonal)
Versions and anomalies correction :
Version numberAuthor (name, e-mail) Corrections/Modifications
0.1 F.