Lines Matching refs:dim
6 static PetscErrorCode private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular(PetscInt dim, PetscInt np[… in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular() argument
14 switch (dim) { in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
25 for (d = 0; d < dim; d++) ds[d] = 2.0 / ((PetscReal)np[d]); in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
27 PetscCall(PetscMalloc1(dim * npoints, &xi)); in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
28 switch (dim) { in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
32 xi[dim * cnt + 0] = -1.0 + 0.5 * ds[d] + ii * ds[0]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
41 xi[dim * cnt + 0] = -1.0 + 0.5 * ds[0] + ii * ds[0]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
42 xi[dim * cnt + 1] = -1.0 + 0.5 * ds[1] + jj * ds[1]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
53 xi[dim * cnt + 0] = -1.0 + 0.5 * ds[0] + ii * ds[0]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
54 xi[dim * cnt + 1] = -1.0 + 0.5 * ds[1] + jj * ds[1]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
55 xi[dim * cnt + 2] = -1.0 + 0.5 * ds[2] + kk * ds[2]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular()
67 static PetscErrorCode private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss(PetscInt dim, PetscInt np_1d… in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss() argument
75 PetscCall(PetscDTGaussTensorQuadrature(dim, 1, np_1d, -1.0, 1.0, &quadrature)); in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss()
77 PetscCall(PetscMalloc1(dim * npoints_q, &xi)); in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss()
79 for (d = 0; d < dim; d++) xi[dim * q + d] = quadrature_xi[dim * q + d]; in private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss()
90 PetscInt dim, npoints_q; in private_DMSwarmInsertPointsUsingCellDM_DA_Q1() local
104 PetscCall(DMGetDimension(dm, &dim)); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
109 PetscCall(private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular(dim, np_dir, &npoints_q, &xi)); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
112 PetscCall(private_DMSwarmCreateCellLocalCoords_DA_Q1_Gauss(dim, npoints, &npoints_q, &xi)); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
123 PetscCall(private_DMSwarmCreateCellLocalCoords_DA_Q1_Regular(dim, np_dir, &npoints_q, &xi)); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
130 PetscCall(PetscMalloc1(dim * npe, &elcoor)); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
135 switch (dim) { in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
137 basis[q][0] = 0.5 * (1.0 - xi[dim * q + 0]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
138 basis[q][1] = 0.5 * (1.0 + xi[dim * q + 0]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
141 basis[q][0] = 0.25 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
142 basis[q][1] = 0.25 * (1.0 + xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
143 basis[q][2] = 0.25 * (1.0 + xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
144 basis[q][3] = 0.25 * (1.0 - xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
148 … basis[q][0] = 0.125 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]) * (1.0 - xi[dim * q + 2]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
149 … basis[q][1] = 0.125 * (1.0 + xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]) * (1.0 - xi[dim * q + 2]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
150 … basis[q][2] = 0.125 * (1.0 + xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]) * (1.0 - xi[dim * q + 2]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
151 … basis[q][3] = 0.125 * (1.0 - xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]) * (1.0 - xi[dim * q + 2]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
152 … basis[q][4] = 0.125 * (1.0 - xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]) * (1.0 + xi[dim * q + 2]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
153 … basis[q][5] = 0.125 * (1.0 + xi[dim * q + 0]) * (1.0 - xi[dim * q + 1]) * (1.0 + xi[dim * q + 2]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
154 … basis[q][6] = 0.125 * (1.0 + xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]) * (1.0 + xi[dim * q + 2]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
155 … basis[q][7] = 0.125 * (1.0 - xi[dim * q + 0]) * (1.0 + xi[dim * q + 1]) * (1.0 + xi[dim * q + 2]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
176 for (d = 0; d < dim; d++) elcoor[dim * k + d] = PetscRealPart(_coor[dim * element[k] + d]); in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
180 for (d = 0; d < dim; d++) swarm_coor[dim * pcnt + d] = 0.0; in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
182 for (d = 0; d < dim; d++) swarm_coor[dim * pcnt + d] += basis[q][k] * elcoor[dim * k + d]; in private_DMSwarmInsertPointsUsingCellDM_DA_Q1()
203 PetscInt dim; in private_DMSwarmInsertPointsUsingCellDM_DA() local
207 PetscCall(DMGetDimension(celldm, &dim)); in private_DMSwarmInsertPointsUsingCellDM_DA()
212 …PetscCheck(dim != 1, PetscObjectComm((PetscObject)dm), PETSC_ERR_SUP, "Support only available for … in private_DMSwarmInsertPointsUsingCellDM_DA()