Lines Matching refs:daf
39 static PetscErrorCode DMCreateInterpolation_DA_1D_Q1(DM dac, DM daf, Mat *A) in DMCreateInterpolation_DA_1D_Q1() argument
54 …PetscCall(DMDAGetInfo(daf, NULL, &mx, NULL, NULL, NULL, NULL, NULL, &dof, NULL, NULL, NULL, NULL, … in DMCreateInterpolation_DA_1D_Q1()
63 PetscCall(DMDAGetCorners(daf, &i_start, NULL, NULL, &m_f, NULL, NULL)); in DMCreateInterpolation_DA_1D_Q1()
64 PetscCall(DMDAGetGhostCorners(daf, &i_start_ghost, NULL, NULL, &m_ghost, NULL, NULL)); in DMCreateInterpolation_DA_1D_Q1()
65 PetscCall(DMGetLocalToGlobalMapping(daf, <og_f)); in DMCreateInterpolation_DA_1D_Q1()
170 static PetscErrorCode DMCreateInterpolation_DA_1D_Q0(DM dac, DM daf, Mat *A) in DMCreateInterpolation_DA_1D_Q0() argument
185 …PetscCall(DMDAGetInfo(daf, NULL, &mx, NULL, NULL, NULL, NULL, NULL, &dof, NULL, NULL, NULL, NULL, … in DMCreateInterpolation_DA_1D_Q0()
196 PetscCall(DMDAGetCorners(daf, &i_start, NULL, NULL, &m_f, NULL, NULL)); in DMCreateInterpolation_DA_1D_Q0()
197 PetscCall(DMDAGetGhostCorners(daf, &i_start_ghost, NULL, NULL, &m_ghost, NULL, NULL)); in DMCreateInterpolation_DA_1D_Q0()
198 PetscCall(DMGetLocalToGlobalMapping(daf, <og_f)); in DMCreateInterpolation_DA_1D_Q0()
256 static PetscErrorCode DMCreateInterpolation_DA_2D_Q1(DM dac, DM daf, Mat *A) in DMCreateInterpolation_DA_2D_Q1() argument
272 …PetscCall(DMDAGetInfo(daf, NULL, &mx, &my, NULL, NULL, NULL, NULL, &dof, NULL, NULL, NULL, NULL, N… in DMCreateInterpolation_DA_2D_Q1()
292 PetscCall(DMDAGetCorners(daf, &i_start, &j_start, NULL, &m_f, &n_f, NULL)); in DMCreateInterpolation_DA_2D_Q1()
293 …PetscCall(DMDAGetGhostCorners(daf, &i_start_ghost, &j_start_ghost, NULL, &m_ghost, &n_ghost, NULL)… in DMCreateInterpolation_DA_2D_Q1()
294 PetscCall(DMGetLocalToGlobalMapping(daf, <og_f)); in DMCreateInterpolation_DA_2D_Q1()
313 PetscCallMPI(MPI_Comm_size(PetscObjectComm((PetscObject)daf), &size_f)); in DMCreateInterpolation_DA_2D_Q1()
314 PetscCallMPI(MPI_Comm_rank(PetscObjectComm((PetscObject)daf), &rank_f)); in DMCreateInterpolation_DA_2D_Q1()
318 …MatPreallocateBegin(PetscObjectComm((PetscObject)daf), m_f * n_f, col_scale * m_c * n_c, dnz, onz); in DMCreateInterpolation_DA_2D_Q1()
348 PetscCall(MatCreate(PetscObjectComm((PetscObject)daf), &mat)); in DMCreateInterpolation_DA_2D_Q1()
484 static PetscErrorCode DMCreateInterpolation_DA_2D_Q0(DM dac, DM daf, Mat *A) in DMCreateInterpolation_DA_2D_Q0() argument
500 …PetscCall(DMDAGetInfo(daf, NULL, &mx, &my, NULL, NULL, NULL, NULL, &dof, NULL, NULL, NULL, NULL, N… in DMCreateInterpolation_DA_2D_Q0()
505 …PetscCheck(ratioi * Mx == mx, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_WRONG, "Fine grid p… in DMCreateInterpolation_DA_2D_Q0()
506 …PetscCheck(ratioj * My == my, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_WRONG, "Fine grid p… in DMCreateInterpolation_DA_2D_Q0()
507 …PetscCheck(ratioi == 2, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_WRONG, "Coarsening factor… in DMCreateInterpolation_DA_2D_Q0()
508 …PetscCheck(ratioj == 2, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_WRONG, "Coarsening factor… in DMCreateInterpolation_DA_2D_Q0()
510 PetscCall(DMDAGetCorners(daf, &i_start, &j_start, NULL, &m_f, &n_f, NULL)); in DMCreateInterpolation_DA_2D_Q0()
511 …PetscCall(DMDAGetGhostCorners(daf, &i_start_ghost, &j_start_ghost, NULL, &m_ghost, &n_ghost, NULL)… in DMCreateInterpolation_DA_2D_Q0()
512 PetscCall(DMGetLocalToGlobalMapping(daf, <og_f)); in DMCreateInterpolation_DA_2D_Q0()
531 PetscCallMPI(MPI_Comm_size(PetscObjectComm((PetscObject)daf), &size_f)); in DMCreateInterpolation_DA_2D_Q0()
532 PetscCallMPI(MPI_Comm_rank(PetscObjectComm((PetscObject)daf), &rank_f)); in DMCreateInterpolation_DA_2D_Q0()
536 …MatPreallocateBegin(PetscObjectComm((PetscObject)daf), m_f * n_f, col_scale * m_c * n_c, dnz, onz); in DMCreateInterpolation_DA_2D_Q0()
560 PetscCall(MatCreate(PetscObjectComm((PetscObject)daf), &mat)); in DMCreateInterpolation_DA_2D_Q0()
605 static PetscErrorCode DMCreateInterpolation_DA_3D_Q0(DM dac, DM daf, Mat *A) in DMCreateInterpolation_DA_3D_Q0() argument
624 …PetscCall(DMDAGetInfo(daf, NULL, &mx, &my, &mz, NULL, NULL, NULL, &dof, NULL, NULL, NULL, NULL, NU… in DMCreateInterpolation_DA_3D_Q0()
628 …PetscCheck(ratioi * Mx == mx, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_WRONG, "Fine grid p… in DMCreateInterpolation_DA_3D_Q0()
629 …PetscCheck(ratioj * My == my, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_WRONG, "Fine grid p… in DMCreateInterpolation_DA_3D_Q0()
630 …PetscCheck(ratiol * Mz == mz, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_WRONG, "Fine grid p… in DMCreateInterpolation_DA_3D_Q0()
631 …PetscCheck(ratioi == 2 || ratioi == 1, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_WRONG, "Co… in DMCreateInterpolation_DA_3D_Q0()
632 …PetscCheck(ratioj == 2 || ratioj == 1, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_WRONG, "Co… in DMCreateInterpolation_DA_3D_Q0()
633 …PetscCheck(ratiol == 2 || ratiol == 1, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_WRONG, "Co… in DMCreateInterpolation_DA_3D_Q0()
635 PetscCall(DMDAGetCorners(daf, &i_start, &j_start, &l_start, &m_f, &n_f, &p_f)); in DMCreateInterpolation_DA_3D_Q0()
636 …PetscCall(DMDAGetGhostCorners(daf, &i_start_ghost, &j_start_ghost, &l_start_ghost, &m_ghost, &n_gh… in DMCreateInterpolation_DA_3D_Q0()
637 PetscCall(DMGetLocalToGlobalMapping(daf, <og_f)); in DMCreateInterpolation_DA_3D_Q0()
656 PetscCallMPI(MPI_Comm_size(PetscObjectComm((PetscObject)daf), &size_f)); in DMCreateInterpolation_DA_3D_Q0()
657 PetscCallMPI(MPI_Comm_rank(PetscObjectComm((PetscObject)daf), &rank_f)); in DMCreateInterpolation_DA_3D_Q0()
661 …MatPreallocateBegin(PetscObjectComm((PetscObject)daf), m_f * n_f * p_f, col_scale * m_c * n_c * p_… in DMCreateInterpolation_DA_3D_Q0()
689 PetscCall(MatCreate(PetscObjectComm((PetscObject)daf), &mat)); in DMCreateInterpolation_DA_3D_Q0()
734 static PetscErrorCode DMCreateInterpolation_DA_3D_Q1(DM dac, DM daf, Mat *A) in DMCreateInterpolation_DA_3D_Q1() argument
751 …PetscCall(DMDAGetInfo(daf, NULL, &mx, &my, &mz, NULL, NULL, NULL, &dof, NULL, NULL, NULL, NULL, NU… in DMCreateInterpolation_DA_3D_Q1()
786 PetscCall(DMDAGetCorners(daf, &i_start, &j_start, &l_start, &m_f, &n_f, &p_f)); in DMCreateInterpolation_DA_3D_Q1()
787 …PetscCall(DMDAGetGhostCorners(daf, &i_start_ghost, &j_start_ghost, &l_start_ghost, &m_ghost, &n_gh… in DMCreateInterpolation_DA_3D_Q1()
788 PetscCall(DMGetLocalToGlobalMapping(daf, <og_f)); in DMCreateInterpolation_DA_3D_Q1()
1007 PetscErrorCode DMCreateInterpolation_DA(DM dac, DM daf, Mat *A, Vec *scale) in DMCreateInterpolation_DA() argument
1016 PetscValidHeaderSpecific(daf, DM_CLASSID, 2); in DMCreateInterpolation_DA()
1021 …PetscCall(DMDAGetInfo(daf, &dimf, &Mf, &Nf, &Pf, &mf, &nf, &pf, &doff, &sf, &bxf, &byf, &bzf, &stf… in DMCreateInterpolation_DA()
1022 …PetscCheck(dimc == dimf, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "Dimensions of D… in DMCreateInterpolation_DA()
1023 …PetscCheck(dofc == doff, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "DOF of DMDA do … in DMCreateInterpolation_DA()
1024 …PetscCheck(sc == sf, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "Stencil width of DM… in DMCreateInterpolation_DA()
1025 …PetscCheck(bxc == bxf && byc == byf && bzc == bzf, PetscObjectComm((PetscObject)daf), PETSC_ERR_AR… in DMCreateInterpolation_DA()
1026 …PetscCheck(stc == stf, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "Stencil type diff… in DMCreateInterpolation_DA()
1033 PetscCall(DMCreateInterpolation_DA_1D_Q1(dac, daf, A)); in DMCreateInterpolation_DA()
1035 PetscCall(DMCreateInterpolation_DA_2D_Q1(dac, daf, A)); in DMCreateInterpolation_DA()
1037 PetscCall(DMCreateInterpolation_DA_3D_Q1(dac, daf, A)); in DMCreateInterpolation_DA()
1038 …} else SETERRQ(PetscObjectComm((PetscObject)daf), PETSC_ERR_SUP, "No support for this DMDA dimensi… in DMCreateInterpolation_DA()
1041 PetscCall(DMCreateInterpolation_DA_1D_Q0(dac, daf, A)); in DMCreateInterpolation_DA()
1043 PetscCall(DMCreateInterpolation_DA_2D_Q0(dac, daf, A)); in DMCreateInterpolation_DA()
1045 PetscCall(DMCreateInterpolation_DA_3D_Q0(dac, daf, A)); in DMCreateInterpolation_DA()
1046 …} else SETERRQ(PetscObjectComm((PetscObject)daf), PETSC_ERR_SUP, "No support for this DMDA dimensi… in DMCreateInterpolation_DA()
1048 if (scale) PetscCall(DMCreateInterpolationScale(dac, daf, *A, scale)); in DMCreateInterpolation_DA()
1052 static PetscErrorCode DMCreateInjection_DA_1D(DM dac, DM daf, VecScatter *inject) in DMCreateInjection_DA_1D() argument
1067 …PetscCall(DMDAGetInfo(daf, NULL, &mx, NULL, NULL, NULL, NULL, NULL, &dof, NULL, NULL, NULL, NULL, … in DMCreateInjection_DA_1D()
1076 PetscCall(DMDAGetCorners(daf, &i_start, NULL, NULL, &m_f, NULL, NULL)); in DMCreateInjection_DA_1D()
1077 PetscCall(DMDAGetGhostCorners(daf, &i_start_ghost, NULL, NULL, &m_ghost, NULL, NULL)); in DMCreateInjection_DA_1D()
1078 PetscCall(DMGetLocalToGlobalMapping(daf, <og_f)); in DMCreateInjection_DA_1D()
1098 …PetscCall(ISCreateBlock(PetscObjectComm((PetscObject)daf), dof, nc, cols, PETSC_OWN_POINTER, &isf)… in DMCreateInjection_DA_1D()
1100 PetscCall(DMGetGlobalVector(daf, &vecf)); in DMCreateInjection_DA_1D()
1103 PetscCall(DMRestoreGlobalVector(daf, &vecf)); in DMCreateInjection_DA_1D()
1108 static PetscErrorCode DMCreateInjection_DA_2D(DM dac, DM daf, VecScatter *inject) in DMCreateInjection_DA_2D() argument
1123 …PetscCall(DMDAGetInfo(daf, NULL, &mx, &my, NULL, NULL, NULL, NULL, &dof, NULL, NULL, NULL, NULL, N… in DMCreateInjection_DA_2D()
1139 PetscCall(DMDAGetCorners(daf, &i_start, &j_start, NULL, &m_f, &n_f, NULL)); in DMCreateInjection_DA_2D()
1140 …PetscCall(DMDAGetGhostCorners(daf, &i_start_ghost, &j_start_ghost, NULL, &m_ghost, &n_ghost, NULL)… in DMCreateInjection_DA_2D()
1141 PetscCall(DMGetLocalToGlobalMapping(daf, <og_f)); in DMCreateInjection_DA_2D()
1164 …PetscCall(ISCreateBlock(PetscObjectComm((PetscObject)daf), dof, nc, cols, PETSC_OWN_POINTER, &isf)… in DMCreateInjection_DA_2D()
1166 PetscCall(DMGetGlobalVector(daf, &vecf)); in DMCreateInjection_DA_2D()
1169 PetscCall(DMRestoreGlobalVector(daf, &vecf)); in DMCreateInjection_DA_2D()
1174 static PetscErrorCode DMCreateInjection_DA_3D(DM dac, DM daf, VecScatter *inject) in DMCreateInjection_DA_3D() argument
1193 …PetscCall(DMDAGetInfo(daf, NULL, &mx, &my, &mz, NULL, NULL, NULL, &dof, NULL, NULL, NULL, NULL, NU… in DMCreateInjection_DA_3D()
1217 PetscCall(DMDAGetCorners(daf, &i_start, &j_start, &k_start, &m_f, &n_f, &p_f)); in DMCreateInjection_DA_3D()
1218 …PetscCall(DMDAGetGhostCorners(daf, &i_start_ghost, &j_start_ghost, &k_start_ghost, &m_ghost, &n_gh… in DMCreateInjection_DA_3D()
1219 PetscCall(DMGetLocalToGlobalMapping(daf, <og_f)); in DMCreateInjection_DA_3D()
1254 …PetscCall(ISCreateBlock(PetscObjectComm((PetscObject)daf), dof, nc, cols, PETSC_OWN_POINTER, &isf)… in DMCreateInjection_DA_3D()
1256 PetscCall(DMGetGlobalVector(daf, &vecf)); in DMCreateInjection_DA_3D()
1259 PetscCall(DMRestoreGlobalVector(daf, &vecf)); in DMCreateInjection_DA_3D()
1264 PetscErrorCode DMCreateInjection_DA(DM dac, DM daf, Mat *mat) in DMCreateInjection_DA() argument
1273 PetscValidHeaderSpecific(daf, DM_CLASSID, 2); in DMCreateInjection_DA()
1277 …PetscCall(DMDAGetInfo(daf, &dimf, &Mf, &Nf, &Pf, &mf, &nf, &pf, &doff, &sf, &bxf, &byf, &bzf, &stf… in DMCreateInjection_DA()
1278 …PetscCheck(dimc == dimf, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "Dimensions of D… in DMCreateInjection_DA()
1279 …PetscCheck(dofc == doff, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "DOF of DMDA do … in DMCreateInjection_DA()
1280 …PetscCheck(sc == sf, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "Stencil width of DM… in DMCreateInjection_DA()
1281 …PetscCheck(bxc == bxf && byc == byf && bzc == bzf, PetscObjectComm((PetscObject)daf), PETSC_ERR_AR… in DMCreateInjection_DA()
1282 …PetscCheck(stc == stf, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "Stencil type diff… in DMCreateInjection_DA()
1288 PetscCall(DMCreateInjection_DA_1D(dac, daf, &inject)); in DMCreateInjection_DA()
1290 PetscCall(DMCreateInjection_DA_2D(dac, daf, &inject)); in DMCreateInjection_DA()
1292 PetscCall(DMCreateInjection_DA_3D(dac, daf, &inject)); in DMCreateInjection_DA()
1305 PetscErrorCode DMCreateAggregates(DM dac, DM daf, Mat *mat) in DMCreateAggregates() argument
1307 return DMDACreateAggregates(dac, daf, mat); in DMCreateAggregates()
1332 PetscErrorCode DMDACreateAggregates(DM dac, DM daf, Mat *rest) in DMDACreateAggregates() argument
1356 PetscValidHeaderSpecificType(daf, DM_CLASSID, 2, DMDA); in DMDACreateAggregates()
1360 …PetscCall(DMDAGetInfo(daf, &dimf, &Mf, &Nf, &Pf, &mf, &nf, &pf, &doff, &sf, &bxf, &byf, &bzf, &stf… in DMDACreateAggregates()
1361 …PetscCheck(dimc == dimf, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "Dimensions of D… in DMDACreateAggregates()
1362 …PetscCheck(dofc == doff, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "DOF of DMDA do … in DMDACreateAggregates()
1363 …PetscCheck(sc == sf, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "Stencil width of DM… in DMDACreateAggregates()
1364 …PetscCheck(bxc == bxf && byc == byf && bzc == bzf, PetscObjectComm((PetscObject)daf), PETSC_ERR_AR… in DMDACreateAggregates()
1365 …PetscCheck(stc == stf, PetscObjectComm((PetscObject)daf), PETSC_ERR_ARG_INCOMP, "Stencil type diff… in DMDACreateAggregates()
1376 PetscCall(DMDAGetCorners(daf, &i_start, &j_start, &l_start, &m_f, &n_f, &p_f)); in DMDACreateAggregates()
1377 …PetscCall(DMDAGetGhostCorners(daf, &i_start_ghost, &j_start_ghost, &l_start_ghost, &m_ghost, &n_gh… in DMDACreateAggregates()
1379 PetscCall(DMGetLocalToGlobalMapping(daf, <ogmf)); in DMDACreateAggregates()
1399 …PetscCall(MatCreateAIJ(PetscObjectComm((PetscObject)daf), m_c * n_c * p_c * dofc, m_f * n_f * p_f … in DMDACreateAggregates()