223   PetscScalar       *Xarr, *z, *lambdai;  in TaoPDIPMInitializeSolution()  local
228   PetscCall(VecGetArrayWrite(pdipm->X, &Xarr));  in TaoPDIPMInitializeSolution()
232   PetscCall(PetscArraycpy(Xarr, xarr, pdipm->nx));  in TaoPDIPMInitializeSolution()
258   PetscCall(VecRestoreArrayWrite(pdipm->X, &Xarr));  in TaoPDIPMInitializeSolution()
280   const PetscScalar *Xarr, *aa;  in TaoSNESJacobian_PDIPM()  local
290   PetscCall(VecGetArrayRead(X, &Xarr));  in TaoSNESJacobian_PDIPM()
299       vals[1] = Xarr[pdipm->off_lambdai + i] / Xarr[pdipm->off_z + i];  in TaoSNESJacobian_PDIPM()
307       vals[0] = Xarr[pdipm->off_z + i];  in TaoSNESJacobian_PDIPM()
308       vals[1] = Xarr[pdipm->off_lambdai + i];  in TaoSNESJacobian_PDIPM()
362   PetscCall(VecPlaceArray(pdipm->x, Xarr));  in TaoSNESJacobian_PDIPM()
415   PetscCall(VecRestoreArrayRead(X, &Xarr));  in TaoSNESJacobian_PDIPM()
446   const PetscScalar *Xarr, *carr, *zarr, *larr;  in TaoSNESFunction_PDIPM()  local
451   PetscCall(VecGetArrayRead(X, &Xarr));  in TaoSNESFunction_PDIPM()
456   PetscCall(VecPlaceArray(x, Xarr));  in TaoSNESFunction_PDIPM()
469       PetscCall(VecPlaceArray(tao->DI, Xarr + pdipm->off_lambdai));  in TaoSNESFunction_PDIPM()
475     PetscCall(VecPlaceArray(pdipm->lambdai_xb, Xarr + pdipm->off_lambdai + pdipm->nh));  in TaoSNESFunction_PDIPM()
489       PetscCall(VecPlaceArray(tao->DE, Xarr + pdipm->off_lambdae));  in TaoSNESFunction_PDIPM()
495       PetscCall(VecPlaceArray(pdipm->lambdae_xfixed, Xarr + pdipm->off_lambdae + pdipm->ng));  in TaoSNESFunction_PDIPM()
514       larr = Xarr + pdipm->off_lambdai;  in TaoSNESFunction_PDIPM()
515       zarr = Xarr + pdipm->off_z;  in TaoSNESFunction_PDIPM()
525       larr = Xarr + pdipm->off_lambdai;  in TaoSNESFunction_PDIPM()
526       zarr = Xarr + pdipm->off_z;  in TaoSNESFunction_PDIPM()
535   PetscCall(VecRestoreArrayRead(X, &Xarr));  in TaoSNESFunction_PDIPM()
552   const PetscScalar *Xarr = NULL;  in TaoSNESFunction_PDIPM_residual()  local
557   PetscCall(VecGetArrayRead(X, &Xarr));  in TaoSNESFunction_PDIPM_residual()
571         for (i = 0; i < pdipm->nci; i++) tmparr[i] *= Xarr[pdipm->off_z + i];  in TaoSNESFunction_PDIPM_residual()
579         for (i = 0; i < pdipm->nci; i++) tmparr[i] /= Xarr[pdipm->off_z + i];  in TaoSNESFunction_PDIPM_residual()
605   PetscCall(VecRestoreArrayRead(X, &Xarr));  in TaoSNESFunction_PDIPM_residual()
708   PetscScalar       *Xarr, *z, *lambdai, dot, *taosolarr;  in SNESLineSearch_PDIPM()  local
718   PetscCall(VecGetArrayWrite(X, &Xarr));  in SNESLineSearch_PDIPM()
720   z  = Xarr + pdipm->off_z;  in SNESLineSearch_PDIPM()
726   lambdai  = Xarr + pdipm->off_lambdai;  in SNESLineSearch_PDIPM()
736   PetscCall(VecRestoreArrayWrite(X, &Xarr));  in SNESLineSearch_PDIPM()
745   PetscCall(VecGetArrayWrite(X, &Xarr));  in SNESLineSearch_PDIPM()
747   for (i = 0; i < pdipm->nx; i++) Xarr[i] -= alpha_p * dXarr[i];  in SNESLineSearch_PDIPM()
748 …for (i = 0; i < pdipm->nce; i++) Xarr[i + pdipm->off_lambdae] -= alpha_d * dXarr[i + pdipm->off_la…  in SNESLineSearch_PDIPM()
751     Xarr[i + pdipm->off_lambdai] -= alpha_d * dXarr[i + pdipm->off_lambdai];  in SNESLineSearch_PDIPM()
752     Xarr[i + pdipm->off_z] -= alpha_p * dXarr[i + pdipm->off_z];  in SNESLineSearch_PDIPM()
755   PetscCall(PetscArraycpy(taosolarr, Xarr, pdipm->nx));  in SNESLineSearch_PDIPM()
758   PetscCall(VecRestoreArrayWrite(X, &Xarr));  in SNESLineSearch_PDIPM()
846   const PetscScalar *aa, *Xarr;  in TaoSetup_PDIPM()  local
915   PetscCall(VecGetArrayRead(pdipm->X, &Xarr));  in TaoSetup_PDIPM()
917   if (pdipm->Nx) PetscCall(VecCreateMPIWithArray(comm, 1, pdipm->nx, pdipm->Nx, Xarr, &pdipm->x));  in TaoSetup_PDIPM()
920 …if (pdipm->Nce) PetscCall(VecCreateMPIWithArray(comm, 1, pdipm->nce, pdipm->Nce, Xarr + pdipm->off…  in TaoSetup_PDIPM()
924 …PetscCall(VecCreateMPIWithArray(comm, 1, pdipm->ng, pdipm->Ng, Xarr + pdipm->off_lambdae, &tao->DE…  in TaoSetup_PDIPM()
933 …PetscCall(VecCreateMPIWithArray(comm, 1, pdipm->nci, pdipm->Nci, Xarr + pdipm->off_lambdai, &pdipm…  in TaoSetup_PDIPM()
936 … PetscCall(VecCreateMPIWithArray(comm, 1, pdipm->nci, pdipm->Nci, Xarr + pdipm->off_z, &pdipm->z));  in TaoSetup_PDIPM()
940 …if (pdipm->Nh) PetscCall(VecCreateMPIWithArray(comm, 1, pdipm->nh, pdipm->Nh, Xarr + pdipm->off_la…  in TaoSetup_PDIPM()
945   PetscCall(VecRestoreArrayRead(pdipm->X, &Xarr));  in TaoSetup_PDIPM()