Function MultAdd

linalg/multivector.cpp:139–148

  template <class T>
  void MultAdd (const class BaseMatrix & mat, FlatVector<T> s, const MultiVector & x, MultiVector & y)
  {
    if constexpr (std::is_same<T,double>::value)
      {
        mat.MultAdd (s, x, y);
      }
    else
      for (auto i : Range(x.Size()))
        mat.MultAdd (s(i), *x[i], *y[i]);
  }
  
  template void MultAdd (const BaseMatrix & mat, FlatVector<double> s, const MultiVector & x, MultiVector & y);
  template void MultAdd (const BaseMatrix & mat, FlatVector<Complex> s, const MultiVector & x, MultiVector & y);