matnull.c (ea709b57ddb8cc304f11df2e9a39d6a8fdb28b53) matnull.c (ef66eb6987ddfdf4e414d6b820cbc8d8d7d17bc2)
1/*$Id: matnull.c,v 1.38 2001/08/06 21:15:07 bsmith Exp balay $*/
1/*$Id: matnull.c,v 1.39 2001/08/07 03:02:41 balay Exp bsmith $*/
2/*
3 Routines to project vectors out of null spaces.
4*/
5
6#include "src/mat/matimpl.h" /*I "petscmat.h" I*/
7#include "petscsys.h"
8
9#undef __FUNCT__

--- 84 unchanged lines hidden (view full) ---

94
95.keywords: PC, null space, remove
96
97.seealso: MatNullSpaceCreate(), MatNullSpaceDestroy()
98@*/
99int MatNullSpaceRemove(MatNullSpace sp,Vec vec,Vec *out)
100{
101 PetscScalar sum;
2/*
3 Routines to project vectors out of null spaces.
4*/
5
6#include "src/mat/matimpl.h" /*I "petscmat.h" I*/
7#include "petscsys.h"
8
9#undef __FUNCT__

--- 84 unchanged lines hidden (view full) ---

94
95.keywords: PC, null space, remove
96
97.seealso: MatNullSpaceCreate(), MatNullSpaceDestroy()
98@*/
99int MatNullSpaceRemove(MatNullSpace sp,Vec vec,Vec *out)
100{
101 PetscScalar sum;
102 int j,n = sp->n,N,ierr;
103 Vec l = vec;
102 int j,n = sp->n,N,ierr;
103 Vec l = vec;
104
105 PetscFunctionBegin;
106 if (out) {
107 if (!sp->vec) {
108 ierr = VecDuplicate(vec,&sp->vec);CHKERRQ(ierr);
109 }
110 *out = sp->vec;
111 ierr = VecCopy(vec,*out);CHKERRQ(ierr);

--- 88 unchanged lines hidden ---
104
105 PetscFunctionBegin;
106 if (out) {
107 if (!sp->vec) {
108 ierr = VecDuplicate(vec,&sp->vec);CHKERRQ(ierr);
109 }
110 *out = sp->vec;
111 ierr = VecCopy(vec,*out);CHKERRQ(ierr);

--- 88 unchanged lines hidden ---