Lines Matching refs:nx
146 static void RDDiffusionCoefficient(RD rd, PetscBool limit, RDNode *n, RDNode *nx, PetscScalar *D_R,… in RDDiffusionCoefficient() argument
152 denom = 3. * rd->rho * sigma_R + (int)limit * PetscAbsScalar(nx->E) / n->E; in RDDiffusionCoefficient()
153 ddenom.E = -(int)limit * PetscAbsScalar(nx->E) / PetscSqr(n->E); in RDDiffusionCoefficient()
155 dxdenom.E = (int)limit * (PetscRealPart(nx->E) < 0 ? -1. : 1.) / n->E; in RDDiffusionCoefficient()
349 RDNode n, nx; in RDIFunction_FD() local
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()
357 RDDiffusionCoefficient(rd, rd->bclimit, &n, &nx, &D_R, 0, 0); in RDIFunction_FD()
358 f[0].E = hx * (n.E - 2. * D_R * nx.E - rd->Eapplied); in RDIFunction_FD()
433 RDNode n, nx; in RDIJacobian_FD() local
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()
442 RDDiffusionCoefficient(rd, rd->bclimit, &n, &nx, &D, &dD, &dxD); in RDIJacobian_FD()
443 … + 0] = (bcTheta / Theta) * hx * (1. - 2. * D * (-1. / hx) - 2. * nx.E * dD.E + 2. * nx.E * dxD.E … in RDIJacobian_FD()
444 K[0][1 * 2 + 1] = (bcTheta / Theta) * hx * (-2. * nx.E * dD.T); in RDIJacobian_FD()
445 …2 * 2 + 0] = (bcTheta / Theta) * hx * (-2. * D * (1. / hx) - 2. * nx.E * dD.E - 2. * nx.E * dxD.E … in RDIJacobian_FD()
488 …interp[][2], PetscReal deriv[][2], PetscInt q, const RDNode x[], PetscInt i, RDNode *n, RDNode *nx) in RDEvaluate() argument
493 nx->E = 0; in RDEvaluate()
494 nx->T = 0; in RDEvaluate()
498 nx->E += deriv[q][j] * x[i + j].E; in RDEvaluate()
499 nx->T += deriv[q][j] * x[i + j].T; in RDEvaluate()
652 RDNode n, n0, nx, n0x, nt, ntx; in RDIFunction_FE() local
653 RDEvaluate(interp, deriv, q, x, i, &n, &nx); in RDIFunction_FE()
669 RDDiffusionCoefficient(rd, PETSC_TRUE, &n, &nx, &D_R, 0, 0); in RDIFunction_FE()
671 …f[i + j].E += (deriv[q][j] * weight[q] * ((1. - Theta) * D0_R * n0x.E + Theta * D_R * nx.E) + inte… in RDIFunction_FE()
681 RDNode n, nx; in RDIFunction_FE() local
685 nx.E = (x[1].E - x[0].E) / hx; in RDIFunction_FE()
686 nx.T = (x[1].T - x[0].T) / hx; in RDIFunction_FE()
687 RDDiffusionCoefficient(rd, PETSC_TRUE, &n, &nx, &D_R, 0, 0); in RDIFunction_FE()
688 RDDiffusionCoefficient(rd, rd->bclimit, &n, &nx, &D_R_bc, 0, 0); in RDIFunction_FE()
738 RDNode n, nx, nt, ntx, drad, dD_R, dxD_R, dEm; in RDIJacobian_FE() local
739 RDEvaluate(interp, deriv, q, x, i, &n, &nx); in RDIJacobian_FE()
742 RDDiffusionCoefficient(rd, PETSC_TRUE, &n, &nx, &D_R, &dD_R, &dxD_R); in RDIJacobian_FE()
746 … interp[q][k] + deriv[q][j] * weight[q] * ((D_R + dxD_R.E * nx.E) * deriv[q][k] + dD_R.E * nx.E * … in RDIJacobian_FE()
747 … interp[q][k]) + deriv[q][j] * weight[q] * (dxD_R.T * deriv[q][k] + dD_R.T * interp[q][k]) * nx.E); in RDIJacobian_FE()
760 RDNode n, nx; in RDIJacobian_FE() local
764 nx.E = (x[1].E - x[0].E) / hx; in RDIJacobian_FE()
765 nx.T = (x[1].T - x[0].T) / hx; in RDIJacobian_FE()
766 RDDiffusionCoefficient(rd, PETSC_TRUE, &n, &nx, &D_R, 0, 0); in RDIJacobian_FE()
767 RDDiffusionCoefficient(rd, rd->bclimit, &n, &nx, &D_R_bc, 0, 0); in RDIJacobian_FE()
870 RDNode n, nx; in RDTestDifferentiation() local
898 nx.E = 1.; in RDTestDifferentiation()
900 RDDiffusionCoefficient(rd, rd->bclimit, &n, &nx, &D0, &dD, &dxD); in RDTestDifferentiation()
903 nx.E = 1.; in RDTestDifferentiation()
905 RDDiffusionCoefficient(rd, rd->bclimit, &n, &nx, &D, 0, 0); in RDTestDifferentiation()
909 nx.E = 1.; in RDTestDifferentiation()
911 RDDiffusionCoefficient(rd, rd->bclimit, &n, &nx, &D, 0, 0); in RDTestDifferentiation()
915 nx.E = 1. + epsilon; in RDTestDifferentiation()
917 RDDiffusionCoefficient(rd, rd->bclimit, &n, &nx, &D, 0, 0); in RDTestDifferentiation()
921 nx.E = 1.; in RDTestDifferentiation()
923 RDDiffusionCoefficient(rd, rd->bclimit, &n, &nx, &D, 0, 0); in RDTestDifferentiation()