Exemplo n.º 1
0
 virtual void apply(const NumericVector<Real> & x, NumericVector<Real> & y) override
 {
   y.pointwise_mult(_diag_inverse, x);
 }
Exemplo n.º 2
0
void TensorShellMatrix<T>::get_diagonal (NumericVector<T>& dest) const
{
  dest.pointwise_mult(_v,_w);
}