Lines Matching refs:vec

92   Vec                vec[6];  in MatNullSpaceCreateRigidBody()  local
110 PetscCall(VecCreate(PetscObjectComm((PetscObject)coords), &vec[0])); in MatNullSpaceCreateRigidBody()
111 PetscCall(VecSetSizes(vec[0], dim * n, dim * N)); in MatNullSpaceCreateRigidBody()
112 PetscCall(VecSetBlockSize(vec[0], dim)); in MatNullSpaceCreateRigidBody()
113 PetscCall(VecSetUp(vec[0])); in MatNullSpaceCreateRigidBody()
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()
148 for (i = 0; i < nmodes; i++) PetscCall(VecRestoreArray(vec[i], &v[i])); in MatNullSpaceCreateRigidBody()
152 PetscCall(VecMDot(vec[i], i, vec, dots)); in MatNullSpaceCreateRigidBody()
154 PetscCall(VecMAXPY(vec[i], i, dots, vec)); in MatNullSpaceCreateRigidBody()
155 PetscCall(VecNormalize(vec[i], NULL)); in MatNullSpaceCreateRigidBody()
157 … PetscCall(MatNullSpaceCreate(PetscObjectComm((PetscObject)coords), PETSC_FALSE, nmodes, vec, sp)); in MatNullSpaceCreateRigidBody()
158 for (i = 0; i < nmodes; i++) PetscCall(VecDestroy(&vec[i])); in MatNullSpaceCreateRigidBody()
341 PetscErrorCode MatNullSpaceRemove(MatNullSpace sp, Vec vec) in MatNullSpaceRemove() argument
349 PetscValidHeaderSpecific(vec, VEC_CLASSID, 2); in MatNullSpaceRemove()
352 PetscCall(VecGetSize(vec, &N)); in MatNullSpaceRemove()
354 PetscCall(VecSum(vec, &sum)); in MatNullSpaceRemove()
356 PetscCall(VecShift(vec, sum)); in MatNullSpaceRemove()
361 PetscCall(VecMDot(vec, sp->n, sp->vecs, sp->alpha)); in MatNullSpaceRemove()
363 PetscCall(VecMAXPY(vec, sp->n, sp->alpha, sp->vecs)); in MatNullSpaceRemove()
366 if (sp->remove) PetscCall((*sp->remove)(sp, vec, sp->rmctx)); in MatNullSpaceRemove()