Lines Matching refs:i
93 PetscInt n, N, dim, nmodes, i, j; in MatNullSpaceCreateRigidBody() local
114 for (i = 1; i < nmodes; i++) PetscCall(VecDuplicate(vec[0], &vec[i])); in MatNullSpaceCreateRigidBody()
115 for (i = 0; i < nmodes; i++) PetscCall(VecGetArray(vec[i], &v[i])); in MatNullSpaceCreateRigidBody()
117 for (i = 0; i < n; i++) { in MatNullSpaceCreateRigidBody()
119 v[0][i * 2 + 0] = sN; in MatNullSpaceCreateRigidBody()
120 v[0][i * 2 + 1] = 0.; in MatNullSpaceCreateRigidBody()
121 v[1][i * 2 + 0] = 0.; in MatNullSpaceCreateRigidBody()
122 v[1][i * 2 + 1] = sN; in MatNullSpaceCreateRigidBody()
124 v[2][i * 2 + 0] = -x[i * 2 + 1]; in MatNullSpaceCreateRigidBody()
125 v[2][i * 2 + 1] = x[i * 2 + 0]; in MatNullSpaceCreateRigidBody()
127 v[0][i * 3 + 0] = sN; in MatNullSpaceCreateRigidBody()
128 v[0][i * 3 + 1] = 0.; in MatNullSpaceCreateRigidBody()
129 v[0][i * 3 + 2] = 0.; in MatNullSpaceCreateRigidBody()
130 v[1][i * 3 + 0] = 0.; in MatNullSpaceCreateRigidBody()
131 v[1][i * 3 + 1] = sN; in MatNullSpaceCreateRigidBody()
132 v[1][i * 3 + 2] = 0.; in MatNullSpaceCreateRigidBody()
133 v[2][i * 3 + 0] = 0.; in MatNullSpaceCreateRigidBody()
134 v[2][i * 3 + 1] = 0.; in MatNullSpaceCreateRigidBody()
135 v[2][i * 3 + 2] = sN; in MatNullSpaceCreateRigidBody()
137 v[3][i * 3 + 0] = x[i * 3 + 1]; in MatNullSpaceCreateRigidBody()
138 v[3][i * 3 + 1] = -x[i * 3 + 0]; in MatNullSpaceCreateRigidBody()
139 v[3][i * 3 + 2] = 0.; in MatNullSpaceCreateRigidBody()
140 v[4][i * 3 + 0] = 0.; in MatNullSpaceCreateRigidBody()
141 v[4][i * 3 + 1] = -x[i * 3 + 2]; in MatNullSpaceCreateRigidBody()
142 v[4][i * 3 + 2] = x[i * 3 + 1]; in MatNullSpaceCreateRigidBody()
143 v[5][i * 3 + 0] = x[i * 3 + 2]; in MatNullSpaceCreateRigidBody()
144 v[5][i * 3 + 1] = 0.; in MatNullSpaceCreateRigidBody()
145 v[5][i * 3 + 2] = -x[i * 3 + 0]; in MatNullSpaceCreateRigidBody()
148 for (i = 0; i < nmodes; i++) PetscCall(VecRestoreArray(vec[i], &v[i])); in MatNullSpaceCreateRigidBody()
150 for (i = dim; i < nmodes; i++) { in MatNullSpaceCreateRigidBody()
152 PetscCall(VecMDot(vec[i], i, vec, dots)); in MatNullSpaceCreateRigidBody()
153 for (j = 0; j < i; j++) dots[j] *= -1.; in MatNullSpaceCreateRigidBody()
154 PetscCall(VecMAXPY(vec[i], i, dots, vec)); in MatNullSpaceCreateRigidBody()
155 PetscCall(VecNormalize(vec[i], NULL)); in MatNullSpaceCreateRigidBody()
158 for (i = 0; i < nmodes; i++) PetscCall(VecDestroy(&vec[i])); in MatNullSpaceCreateRigidBody()
189 PetscInt i; in MatNullSpaceView() local
196 for (i = 0; i < sp->n; i++) PetscCall(VecView(sp->vecs[i], viewer)); in MatNullSpaceView()
233 PetscInt i; in MatNullSpaceCreate() local
238 for (i = 0; i < n; i++) PetscValidHeaderSpecific(vecs[i], VEC_CLASSID, 4); in MatNullSpaceCreate()
241 for (i = 0; i < n; i++) { in MatNullSpaceCreate()
243 PetscCall(VecLockReadPush(vecs[i])); in MatNullSpaceCreate()
248 for (i = 0; i < n; i++) { in MatNullSpaceCreate()
250 PetscCall(VecNorm(vecs[i], NORM_2, &norm)); in MatNullSpaceCreate()
251 …tscObjectComm((PetscObject)vecs[i]), PETSC_ERR_ARG_WRONG, "Vector %" PetscInt_FMT " must have 2-no… in MatNullSpaceCreate()
254 for (i = 0; i < n; i++) { in MatNullSpaceCreate()
256 PetscCall(VecSum(vecs[i], &sum)); in MatNullSpaceCreate()
257 …((PetscObject)vecs[i]), PETSC_ERR_ARG_WRONG, "Vector %" PetscInt_FMT " must be orthogonal to const… in MatNullSpaceCreate()
261 for (i = 0; i < n - 1; i++) { in MatNullSpaceCreate()
263 PetscCall(VecMDot(vecs[i], n - i - 1, vecs + i + 1, dots)); in MatNullSpaceCreate()
264 for (j = 0; j < n - i - 1; j++) { in MatNullSpaceCreate()
265 …cObject)vecs[i]), PETSC_ERR_ARG_WRONG, "Vector %" PetscInt_FMT " must be orthogonal to vector %" P… in MatNullSpaceCreate()
286 for (i = 0; i < n; i++) { in MatNullSpaceCreate()
287 PetscCall(PetscObjectReference((PetscObject)vecs[i])); in MatNullSpaceCreate()
288 sp->vecs[i] = vecs[i]; in MatNullSpaceCreate()
310 PetscInt i; in MatNullSpaceDestroy() local
320 for (i = 0; i < (*sp)->n; i++) PetscCall(VecLockReadPop((*sp)->vecs[i])); in MatNullSpaceDestroy()
344 PetscInt i, N; in MatNullSpaceRemove() local
362 for (i = 0; i < sp->n; i++) sp->alpha[i] = -sp->alpha[i]; in MatNullSpaceRemove()