Method InnerProductD

parallel/parallelvvector.cpp:556–588

    }
                         
    Matrix<> InnerProductD (const MultiVector & v2) const override
    {
      Matrix<double> res(v2.Size(), Size());
      if (res.Height()*res.Width()==0) return res;
      
      auto status1 = (*this)[0] -> GetParallelStatus();
      auto status2 = (v2)[0] -> GetParallelStatus();

      if (auto pv2 = dynamic_cast<const ParallelMultiVector*>(&v2))
        if ( ((status1 == CUMULATED) && (status2 == DISTRIBUTED)) ||
             ((status2 == CUMULATED) && (status1 == DISTRIBUTED)) )
          {
            ArrayMem<double*,32> ptrs1(Size());
            ArrayMem<double*,32> ptrs2(v2.Size());
            MakeSameStatus();
            pv2->MakeSameStatus();
            
            for (int i = 0; i < Size(); i++)
              ptrs1[i] = (*this)[i]->FVDouble().Data();
            for (int i = 0; i < v2.Size(); i++)
              ptrs2[i] = v2[i]->FVDouble().Data();
            
            ngbla::PairwiseInnerProduct((*this)[0]->FVDouble().Size(), ptrs2, ptrs1, res);
            
            paralleldofs->GetCommunicator()
              .AllReduce(FlatArray(res.Height()*res.Width(), res.Data()), NG_MPI_SUM);
            
            return res;
          }

      // fallback
      return MultiVector::InnerProductD(v2);
    }

    
  };