Lines Matching refs:u
98 const PetscScalar *u, *udot; in IFunction() local
102 PetscCall(VecGetArrayRead(U, &u)); in IFunction()
105 f[0] = udot[0] + ctx->k * u[0] * u[1]; in IFunction()
106 f[1] = udot[1] + ctx->k * u[0] * u[1]; in IFunction()
107 f[2] = udot[2] - ctx->k * u[0] * u[1]; in IFunction()
108 PetscCall(VecRestoreArrayRead(U, &u)); in IFunction()
121 const PetscScalar *u, *udot; in IJacobian() local
124 PetscCall(VecGetArrayRead(U, &u)); in IJacobian()
126 J[0][0] = a + ctx->k * u[1]; in IJacobian()
127 J[0][1] = ctx->k * u[0]; in IJacobian()
129 J[1][0] = ctx->k * u[1]; in IJacobian()
130 J[1][1] = a + ctx->k * u[0]; in IJacobian()
132 J[2][0] = -ctx->k * u[1]; in IJacobian()
133 J[2][1] = -ctx->k * u[0]; in IJacobian()
136 PetscCall(VecRestoreArrayRead(U, &u)); in IJacobian()
154 PetscScalar *u, d0, q; in Solution() local
158 PetscCall(VecGetArrayWrite(U, &u)); in Solution()
162 u[0] = uinit[0] / (1.0 + uinit[1] * q); in Solution()
163 u[1] = u[0] - d0; in Solution()
164 u[2] = uinit[1] + uinit[2] - u[1]; in Solution()
165 PetscCall(VecRestoreArrayWrite(U, &u)); in Solution()
178 PetscScalar *u; in main() local
208 PetscCall(VecGetArrayWrite(ctx.initialsolution, &u)); in main()
209 u[0] = 1; in main()
210 u[1] = .7; in main()
211 u[2] = 0; in main()
212 PetscCall(VecRestoreArrayWrite(ctx.initialsolution, &u)); in main()