Lines Matching refs:lambda
113 PetscScalar *u, *lambda, *w, *fu, *fw, *flambda, d, h; in FormFunction() local
127 PetscCall(DMDAVecGetArray(user->da1, vlambda, &lambda)); in FormFunction()
134 fw[0] = -2. * d * lambda[0]; in FormFunction()
139 if (i == 0) flambda[0] = h * u[0] + 2. * d * lambda[0] - d * lambda[1]; in FormFunction()
140 else if (i == 1) flambda[1] = 2. * h * u[1] + 2. * d * lambda[1] - d * lambda[2]; in FormFunction()
141 else if (i == N - 1) flambda[N - 1] = h * u[N - 1] + 2. * d * lambda[N - 1] - d * lambda[N - 2]; in FormFunction()
142 …else if (i == N - 2) flambda[N - 2] = 2. * h * u[N - 2] + 2. * d * lambda[N - 2] - d * lambda[N - … in FormFunction()
143 else flambda[i] = 2. * h * u[i] - d * (lambda[i + 1] - 2.0 * lambda[i] + lambda[i - 1]); in FormFunction()
157 PetscCall(DMDAVecRestoreArray(user->da1, vlambda, &lambda)); in FormFunction()
169 Vec w, u, lambda, U, F; in Monitor() local
173 PetscCall(DMCompositeGetAccess(user->packer, U, &w, &u, &lambda)); in Monitor()
175 PetscCall(VecView(lambda, user->lambda_viewer)); in Monitor()
176 PetscCall(DMCompositeRestoreAccess(user->packer, U, &w, &u, &lambda)); in Monitor()
179 PetscCall(DMCompositeGetAccess(user->packer, F, &w, &u, &lambda)); in Monitor()
181 PetscCall(VecView(lambda, user->flambda_viewer)); in Monitor()
182 PetscCall(DMCompositeRestoreAccess(user->packer, F, &w, &u, &lambda)); in Monitor()