527   User               user_ptr = (User)ctx;  in FormFunctionGradient()  local
528   TS                 ts       = user_ptr->ts;  in FormFunctionGradient()
534   user_ptr->mu = y_ptr[0];  in FormFunctionGradient()
541   PetscCall(VecGetArray(user_ptr->U, &x_ptr));  in FormFunctionGradient()
543 …x_ptr[1] = -2.0 / 3.0 + 10.0 / (81.0 * user_ptr->mu) - 292.0 / (2187.0 * user_ptr->mu * user_ptr->…  in FormFunctionGradient()
544   PetscCall(VecRestoreArray(user_ptr->U, &x_ptr));  in FormFunctionGradient()
546   PetscCall(TSSolve(ts, user_ptr->U));  in FormFunctionGradient()
548   PetscCall(VecGetArrayRead(user_ptr->U, &y_ptr));  in FormFunctionGradient()
549 …*f = (y_ptr[0] - user_ptr->ob[0]) * (y_ptr[0] - user_ptr->ob[0]) + (y_ptr[1] - user_ptr->ob[1]) * …  in FormFunctionGradient()
552   PetscCall(VecGetArray(user_ptr->Lambda[0], &x_ptr));  in FormFunctionGradient()
553   x_ptr[0] = 2. * (y_ptr[0] - user_ptr->ob[0]);  in FormFunctionGradient()
554   x_ptr[1] = 2. * (y_ptr[1] - user_ptr->ob[1]);  in FormFunctionGradient()
555   PetscCall(VecRestoreArrayRead(user_ptr->U, &y_ptr));  in FormFunctionGradient()
556   PetscCall(VecRestoreArray(user_ptr->Lambda[0], &x_ptr));  in FormFunctionGradient()
558   PetscCall(VecGetArray(user_ptr->Mup[0], &x_ptr));  in FormFunctionGradient()
560   PetscCall(VecRestoreArray(user_ptr->Mup[0], &x_ptr));  in FormFunctionGradient()
561   PetscCall(TSSetCostGradients(ts, 1, user_ptr->Lambda, user_ptr->Mup));  in FormFunctionGradient()
565   PetscCall(VecGetArray(user_ptr->Mup[0], &x_ptr));  in FormFunctionGradient()
566   PetscCall(VecGetArrayRead(user_ptr->Lambda[0], &y_ptr));  in FormFunctionGradient()
568 …= y_ptr[1] * (-10.0 / (81.0 * user_ptr->mu * user_ptr->mu) + 2.0 * 292.0 / (2187.0 * user_ptr->mu …  in FormFunctionGradient()
569   PetscCall(VecRestoreArray(user_ptr->Mup[0], &x_ptr));  in FormFunctionGradient()
570   PetscCall(VecRestoreArrayRead(user_ptr->Lambda[0], &y_ptr));  in FormFunctionGradient()
577   User           user_ptr = (User)ctx;  in FormHessian()  local
583   PetscCall(Adjoint2(P, harr, user_ptr));  in FormHessian()