Lines Matching refs:u
45 const PetscScalar *u, *udot; in IFunction() local
48 PetscCall(VecGetArrayRead(U, &u)); in IFunction()
51 f[0] = udot[0] - k1(ctx, t) * u[2] + ctx->k2 * u[0]; in IFunction()
52 f[1] = udot[1] - k1(ctx, t) * u[2] + ctx->k3 * u[1] * u[3] - ctx->sigma2; in IFunction()
53 f[2] = udot[2] - ctx->k3 * u[1] * u[3] + k1(ctx, t) * u[2]; in IFunction()
54 f[3] = udot[3] - ctx->k2 * u[0] + ctx->k3 * u[1] * u[3]; in IFunction()
55 PetscCall(VecRestoreArrayRead(U, &u)); in IFunction()
65 const PetscScalar *u, *udot; in IJacobian() local
68 PetscCall(VecGetArrayRead(U, &u)); in IJacobian()
75 J[1][1] = a + ctx->k3 * u[3]; in IJacobian()
77 J[1][3] = ctx->k3 * u[1]; in IJacobian()
79 J[2][1] = -ctx->k3 * u[3]; in IJacobian()
81 J[2][3] = -ctx->k3 * u[1]; in IJacobian()
83 J[3][1] = ctx->k3 * u[3]; in IJacobian()
85 J[3][3] = a + ctx->k3 * u[1]; in IJacobian()
87 PetscCall(VecRestoreArrayRead(U, &u)); in IJacobian()
115 PetscScalar *u; in main() local
141 PetscCall(VecGetArrayWrite(ctx.initialsolution, &u)); in main()
142 u[0] = 0.0; in main()
143 u[1] = 1.3e8; in main()
144 u[2] = 5.0e11; in main()
145 u[3] = 8.0e11; in main()
146 PetscCall(VecRestoreArrayWrite(ctx.initialsolution, &u)); in main()