multequal.c (25da5fcfbc1fce3d855b9e613e14e74389d518e2) multequal.c (09f3b4e5628a00a1eaf17d80982cfbcc515cc9c1)
1#define PETSCMAT_DLL
2
3#include "src/mat/matimpl.h" /*I "petscmat.h" I*/
4
5#undef __FUNCT__
6#define __FUNCT__ "MatMultEqual"
7/*@
8 MatMultEqual - Compares matrix-vector products of two matrices.

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

57 ierr = VecNorm(s1,NORM_INFINITY,&r1);CHKERRQ(ierr);
58 } else {
59 ierr = VecAXPY(s2,none,s1);CHKERRQ(ierr);
60 ierr = VecNorm(s2,NORM_INFINITY,&r1);CHKERRQ(ierr);
61 r1 /= r2;
62 }
63 if (r1 > tol) {
64 *flg = PETSC_FALSE;
1#define PETSCMAT_DLL
2
3#include "src/mat/matimpl.h" /*I "petscmat.h" I*/
4
5#undef __FUNCT__
6#define __FUNCT__ "MatMultEqual"
7/*@
8 MatMultEqual - Compares matrix-vector products of two matrices.

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

57 ierr = VecNorm(s1,NORM_INFINITY,&r1);CHKERRQ(ierr);
58 } else {
59 ierr = VecAXPY(s2,none,s1);CHKERRQ(ierr);
60 ierr = VecNorm(s2,NORM_INFINITY,&r1);CHKERRQ(ierr);
61 r1 /= r2;
62 }
63 if (r1 > tol) {
64 *flg = PETSC_FALSE;
65 ierr = PetscLogInfo((0,"Error: %D-th MatMult() %g\n",k,r1));CHKERRQ(ierr);
65 ierr = PetscVerboseInfo((0,"Error: %D-th MatMult() %g\n",k,r1));CHKERRQ(ierr);
66 break;
67 }
68 }
69 ierr = PetscRandomDestroy(rctx);CHKERRQ(ierr);
70 ierr = VecDestroy(x);CHKERRQ(ierr);
71 ierr = VecDestroy(s1);CHKERRQ(ierr);
72 ierr = VecDestroy(s2);CHKERRQ(ierr);
73 PetscFunctionReturn(0);

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

128 ierr = VecNorm(s1,NORM_INFINITY,&r1);CHKERRQ(ierr);
129 } else {
130 ierr = VecAXPY(s2,none,s1);CHKERRQ(ierr);
131 ierr = VecNorm(s2,NORM_INFINITY,&r1);CHKERRQ(ierr);
132 r1 /= r2;
133 }
134 if (r1 > tol) {
135 *flg = PETSC_FALSE;
66 break;
67 }
68 }
69 ierr = PetscRandomDestroy(rctx);CHKERRQ(ierr);
70 ierr = VecDestroy(x);CHKERRQ(ierr);
71 ierr = VecDestroy(s1);CHKERRQ(ierr);
72 ierr = VecDestroy(s2);CHKERRQ(ierr);
73 PetscFunctionReturn(0);

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

128 ierr = VecNorm(s1,NORM_INFINITY,&r1);CHKERRQ(ierr);
129 } else {
130 ierr = VecAXPY(s2,none,s1);CHKERRQ(ierr);
131 ierr = VecNorm(s2,NORM_INFINITY,&r1);CHKERRQ(ierr);
132 r1 /= r2;
133 }
134 if (r1 > tol) {
135 *flg = PETSC_FALSE;
136 ierr = PetscLogInfo((0,"Error: %d-th MatMultAdd() %g\n",k,r1));CHKERRQ(ierr);
136 ierr = PetscVerboseInfo((0,"Error: %d-th MatMultAdd() %g\n",k,r1));CHKERRQ(ierr);
137 break;
138 }
139 }
140 ierr = PetscRandomDestroy(rctx);CHKERRQ(ierr);
141 ierr = VecDestroy(x);CHKERRQ(ierr);
142 ierr = VecDestroy(y);CHKERRQ(ierr);
143 ierr = VecDestroy(s1);CHKERRQ(ierr);
144 ierr = VecDestroy(s2);CHKERRQ(ierr);

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

198 ierr = VecNorm(s1,NORM_INFINITY,&r1);CHKERRQ(ierr);
199 } else {
200 ierr = VecAXPY(s2,none,s1);CHKERRQ(ierr);
201 ierr = VecNorm(s2,NORM_INFINITY,&r1);CHKERRQ(ierr);
202 r1 /= r2;
203 }
204 if (r1 > tol) {
205 *flg = PETSC_FALSE;
137 break;
138 }
139 }
140 ierr = PetscRandomDestroy(rctx);CHKERRQ(ierr);
141 ierr = VecDestroy(x);CHKERRQ(ierr);
142 ierr = VecDestroy(y);CHKERRQ(ierr);
143 ierr = VecDestroy(s1);CHKERRQ(ierr);
144 ierr = VecDestroy(s2);CHKERRQ(ierr);

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

198 ierr = VecNorm(s1,NORM_INFINITY,&r1);CHKERRQ(ierr);
199 } else {
200 ierr = VecAXPY(s2,none,s1);CHKERRQ(ierr);
201 ierr = VecNorm(s2,NORM_INFINITY,&r1);CHKERRQ(ierr);
202 r1 /= r2;
203 }
204 if (r1 > tol) {
205 *flg = PETSC_FALSE;
206 ierr = PetscLogInfo((0,"Error: %d-th MatMultTranspose() %g\n",k,r1));CHKERRQ(ierr);
206 ierr = PetscVerboseInfo((0,"Error: %d-th MatMultTranspose() %g\n",k,r1));CHKERRQ(ierr);
207 break;
208 }
209 }
210 ierr = PetscRandomDestroy(rctx);CHKERRQ(ierr);
211 ierr = VecDestroy(x);CHKERRQ(ierr);
212 ierr = VecDestroy(s1);CHKERRQ(ierr);
213 ierr = VecDestroy(s2);CHKERRQ(ierr);
214 PetscFunctionReturn(0);

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

269 ierr = VecNorm(s1,NORM_INFINITY,&r1);CHKERRQ(ierr);
270 } else {
271 ierr = VecAXPY(s2,none,s1);CHKERRQ(ierr);
272 ierr = VecNorm(s2,NORM_INFINITY,&r1);CHKERRQ(ierr);
273 r1 /= r2;
274 }
275 if (r1 > tol) {
276 *flg = PETSC_FALSE;
207 break;
208 }
209 }
210 ierr = PetscRandomDestroy(rctx);CHKERRQ(ierr);
211 ierr = VecDestroy(x);CHKERRQ(ierr);
212 ierr = VecDestroy(s1);CHKERRQ(ierr);
213 ierr = VecDestroy(s2);CHKERRQ(ierr);
214 PetscFunctionReturn(0);

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

269 ierr = VecNorm(s1,NORM_INFINITY,&r1);CHKERRQ(ierr);
270 } else {
271 ierr = VecAXPY(s2,none,s1);CHKERRQ(ierr);
272 ierr = VecNorm(s2,NORM_INFINITY,&r1);CHKERRQ(ierr);
273 r1 /= r2;
274 }
275 if (r1 > tol) {
276 *flg = PETSC_FALSE;
277 ierr = PetscLogInfo((0,"Error: %d-th MatMultTransposeAdd() %g\n",k,r1));CHKERRQ(ierr);
277 ierr = PetscVerboseInfo((0,"Error: %d-th MatMultTransposeAdd() %g\n",k,r1));CHKERRQ(ierr);
278 break;
279 }
280 }
281 ierr = PetscRandomDestroy(rctx);CHKERRQ(ierr);
282 ierr = VecDestroy(x);CHKERRQ(ierr);
283 ierr = VecDestroy(y);CHKERRQ(ierr);
284 ierr = VecDestroy(s1);CHKERRQ(ierr);
285 ierr = VecDestroy(s2);CHKERRQ(ierr);
286 PetscFunctionReturn(0);
287}
278 break;
279 }
280 }
281 ierr = PetscRandomDestroy(rctx);CHKERRQ(ierr);
282 ierr = VecDestroy(x);CHKERRQ(ierr);
283 ierr = VecDestroy(y);CHKERRQ(ierr);
284 ierr = VecDestroy(s1);CHKERRQ(ierr);
285 ierr = VecDestroy(s2);CHKERRQ(ierr);
286 PetscFunctionReturn(0);
287}