Searched refs:xComp (Results 1 – 1 of 1) sorted by relevance
| /petsc/src/snes/tutorials/ |
| H A D | ex7.c | 328 Vec *xComp, *pComp, xTmp, pTmp; in ComputeFFT() local
332 PetscCall(PetscMalloc2(Nc, &xComp, Nc, &pComp)); in ComputeFFT()
339 PetscCall(VecCreate(PetscObjectComm((PetscObject)x), &xComp[i])); in ComputeFFT()
341 PetscCall(VecSetType(xComp[i], vtype)); in ComputeFFT()
342 PetscCall(VecSetSizes(xComp[i], n / Nc, N / Nc)); in ComputeFFT()
343 PetscCall(VecDuplicate(xComp[i], &pComp[i])); in ComputeFFT()
346 PetscCall(VecStrideGatherAll(x, xComp, INSERT_VALUES)); in ComputeFFT()
348 PetscCall(VecScatterPetscToFFTW(FT, xComp[i], xTmp)); in ComputeFFT()
354 PetscCall(VecDestroy(&xComp[i])); in ComputeFFT()
359 PetscCall(PetscFree2(xComp, pComp)); in ComputeFFT()
|