Home
last modified time | relevance | path

Searched refs:theta_c (Results 1 – 3 of 3) sorted by relevance

/petsc/src/ts/tutorials/phasefield/
H A Dbiharmonic3.c40 PetscReal theta_c; member
73 ctx.theta_c = 1.0; in main()
75 PetscCall(PetscOptionsGetReal(NULL, NULL, "-theta_c", &ctx.theta_c, NULL)); in main()
238 …theta * (-PetscLogScalar(ctx->tol) + PetscLogScalar((1.0 - x[i].u) / 2.0)) + ctx->theta_c * x[i].u; in FormFunction()
239 …theta * (-PetscLogScalar((1.0 + x[i].u) / 2.0) + PetscLogScalar(ctx->tol)) + ctx->theta_c * x[i].u; in FormFunction()
240 …tscLogScalar((1.0 + x[i].u) / 2.0) + PetscLogScalar((1.0 - x[i].u) / 2.0)) + ctx->theta_c * x[i].u; in FormFunction()
253 …[i].u += 2.0 * ctx->theta_c * x[i].u * (x[i + 1].u - x[i - 1].u) * (x[i + 1].u - x[i - 1].u) * .25… in FormFunction()
H A Dbiharmonic2.c40 PetscReal theta_c; member
75 ctx.theta_c = 1.0; in main()
77 PetscCall(PetscOptionsGetReal(NULL, NULL, "-theta_c", &ctx.theta_c, NULL)); in main()
238 …->theta * (-PetscLogReal(ctx->tol) + PetscLogScalar((1.0 - x[i].u) / 2.0)) + ctx->theta_c * x[i].u; in FormFunction()
239 …->theta * (-PetscLogScalar((1.0 + x[i].u) / 2.0) + PetscLogReal(ctx->tol)) + ctx->theta_c * x[i].u; in FormFunction()
240 …tscLogScalar((1.0 + x[i].u) / 2.0) + PetscLogScalar((1.0 - x[i].u) / 2.0)) + ctx->theta_c * x[i].u; in FormFunction()
H A Dbiharmonic.c59 PetscReal theta, theta_c; member
94 ctx.theta_c = 1.0; in main()
96 PetscCall(PetscOptionsGetReal(NULL, NULL, "-theta_c", &ctx.theta_c, NULL)); in main()
205 …PetscReal tol = ctx->tol, theta = ctx->theta, theta_c = ctx->theta_c, a, b; /* a and b are used… in FormFunction() local
272 f[i] += -theta_c * (x[i - 1] + x[i + 1] - 2.0 * x[i]) * sx; in FormFunction()
472 …PetscReal tol = ctx->tol, theta = ctx->theta, theta_c = ctx->theta_c, a, b; /* a and b a… in MyMonitor() local
530 yy[1] = .5 * theta_c * PetscRealPart(1.0 - u[i] * u[i]); in MyMonitor()
595 yy[1] = theta_c * PetscRealPart(-(u[i - 1] + u[i + 1] - 2.0 * u[i])) * sx; in MyMonitor()
697 len = -.5 * theta_c * PetscRealPart(-(u[i - 1] + u[i + 1] - 2.0 * u[i]) * sx / max); in MyMonitor()