241 …ts, Vec X, Vec Xdot, PetscReal *Theta, PetscReal *dt, Vec *X0loc, RDNode **x0, Vec *Xloc, RDNode *…  in RDGetLocalArrays()  argument
268   PetscCall(DMDAVecGetArray(rd->da, *X0loc, x0));  in RDGetLocalArrays()
274 static PetscErrorCode RDRestoreLocalArrays(RD rd, Vec *X0loc, RDNode **x0, Vec *Xloc, RDNode **x, V…  in RDRestoreLocalArrays()  argument
277   PetscCall(DMDAVecRestoreArray(rd->da, *X0loc, x0));  in RDRestoreLocalArrays()
314   RDNode       *x, *x0, *xdot, *f;  in RDIFunction_FD()  local
321   PetscCall(RDGetLocalArrays(rd, ts, X, Xdot, &Theta, &dt, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIFunction_FD()
332     rad = (1. - Theta) * RDRadiation(rd, &x0[i], 0) + Theta * RDRadiation(rd, &x[i], 0);  in RDIFunction_FD()
335       RDMaterialEnergy(rd, &x0[i], &Em0, NULL);  in RDIFunction_FD()
351       n.E  = (1. - bcTheta) * x0[0].E + bcTheta * x[0].E;  in RDIFunction_FD()
352       n.T  = (1. - bcTheta) * x0[0].T + bcTheta * x[0].T;  in RDIFunction_FD()
353       nx.E = ((1. - bcTheta) * (x0[1].E - x0[0].E) + bcTheta * (x[1].E - x[0].E)) / hx;  in RDIFunction_FD()
354       nx.T = ((1. - bcTheta) * (x0[1].T - x0[0].T) + bcTheta * (x[1].T - x[0].T)) / hx;  in RDIFunction_FD()
369 …PetscScalar diff = (1. - Theta) * RDDiffusion(rd, hx, x0, i, 0) + Theta * RDDiffusion(rd, hx, x, i…  in RDIFunction_FD()
373   PetscCall(RDRestoreLocalArrays(rd, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIFunction_FD()
382   RDNode       *x, *x0, *xdot;  in RDIJacobian_FD()  local
389   PetscCall(RDGetLocalArrays(rd, ts, X, Xdot, &Theta, &dt, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIJacobian_FD()
400     /*rad = (1.-Theta)* */ RDRadiation(rd, &x0[i], 0); /* + Theta* */  in RDIJacobian_FD()
406       RDMaterialEnergy(rd, &x0[i], &Em0, NULL);  in RDIJacobian_FD()
436       n.E  = (1. - bcTheta) * x0[0].E + bcTheta * x[0].E;  in RDIJacobian_FD()
437       n.T  = (1. - bcTheta) * x0[0].T + bcTheta * x[0].T;  in RDIJacobian_FD()
438       nx.E = ((1. - bcTheta) * (x0[1].E - x0[0].E) + bcTheta * (x[1].E - x[0].E)) / hx;  in RDIJacobian_FD()
439       nx.T = ((1. - bcTheta) * (x0[1].T - x0[0].T) + bcTheta * (x[1].T - x[0].T)) / hx;  in RDIJacobian_FD()
477   PetscCall(RDRestoreLocalArrays(rd, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIJacobian_FD()
630   RDNode       *x, *x0, *xdot, *f;  in RDIFunction_FE()  local
637   PetscCall(RDGetLocalArrays(rd, ts, X, Xdot, &Theta, &dt, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIFunction_FE()
654       RDEvaluate(interp, deriv, q, x0, i, &n0, &n0x);  in RDIFunction_FE()
683       n.E  = (1 - bcTheta) * x0[0].E + bcTheta * x[0].E;  in RDIFunction_FE()
684       n.T  = (1 - bcTheta) * x0[0].T + bcTheta * x[0].T;  in RDIFunction_FE()
702   PetscCall(RDRestoreLocalArrays(rd, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIFunction_FE()
716   RDNode       *x, *x0, *xdot;  in RDIJacobian_FE()  local
724   PetscCall(RDGetLocalArrays(rd, ts, X, Xdot, &Theta, &dt, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIJacobian_FE()
762       n.E  = (1 - Theta) * x0[0].E + Theta * x[0].E;  in RDIJacobian_FE()
763       n.T  = (1 - Theta) * x0[0].T + Theta * x[0].T;  in RDIJacobian_FE()
779   PetscCall(RDRestoreLocalArrays(rd, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIJacobian_FE()