v(i) = i;
norm += std::fabs(1.*i)*std::fabs(1.*i);
}
- v.compress (VectorOperation::add);
+ v.compress (VectorOperation::insert);
// then check the norm
const double eps=typeid(PetscScalar)==typeid(double) ? 1e-14 : 1e-5;
for (unsigned int i=0; i<v.size(); ++i)
if (i%3 == 0)
v(i) = i+1.;
- v.compress (VectorOperation::add);
+ v.compress (VectorOperation::insert);
// then resize with setting to zero
v.reinit (13, false);
y(i) = i+3.;
}
- v.compress (VectorOperation::add);
- w.compress (VectorOperation::add);
- x.compress (VectorOperation::add);
- y.compress (VectorOperation::add);
+ v.compress (VectorOperation::insert);
+ w.compress (VectorOperation::insert);
+ x.compress (VectorOperation::insert);
+ y.compress (VectorOperation::insert);
v.sadd (1.5, 2, w, 3, x, 4, y);
if ((i+2*j+1) % 3 == 0)
m.set (i,j, i*j*.5+.5);
- m.compress (VectorOperation::add);
+ m.compress (VectorOperation::insert);
// compare against the exact value of the
// linfty-norm (max row-sum)
w(i) = i;
}
- m.compress (VectorOperation::add);
- v.compress (VectorOperation::add);
- w.compress (VectorOperation::add);
+ m.compress (VectorOperation::insert);
+ v.compress (VectorOperation::insert);
+ w.compress (VectorOperation::insert);
// w:=Mv
m.Tvmult_add (w,v);
PETScWrappers::Vector f(dim);
PETScWrappers::Vector u(dim);
f = 1.;
- A.compress (VectorOperation::add);
- f.compress (VectorOperation::add);
- u.compress (VectorOperation::add);
+ A.compress (VectorOperation::insert);
+ f.compress (VectorOperation::insert);
PETScWrappers::SolverCG solver(control);
PETScWrappers::PreconditionILU preconditioner(A);