STK::hidden::MultPointArray< Lhs, Rhs, Result > Struct Template Reference

This structure regroup the products between a point and different kind of array. More...

#include <STK_ArrayByVectorProduct.h>

Public Types
typedef Result::Type	Type

Static Public Member Functions
template<class SubLhs , class SubRhs >
static void	mult (SubLhs const &l, SubRhs const &r, Result &res)
	Compute the product res = l*r with l a point (a row-vector) and r a matrix using the standard formula .
static void	run (ExprBase< Lhs > const &l, ExprBase< Rhs > const &r, Result &res)
	Compute the product res = l*r with l a point (a row-vector) and r a matrix using column decomposition.

Detailed Description

template<typename Lhs, typename Rhs, typename Result>
struct STK::hidden::MultPointArray< Lhs, Rhs, Result >

This structure regroup the products between a point and different kind of array.

Definition at line 126 of file STK_ArrayByVectorProduct.h.

Member Typedef Documentation

Type

template<typename Lhs , typename Rhs , typename Result >

typedef Result::Type STK::hidden::MultPointArray< Lhs, Rhs, Result >::Type

Definition at line 128 of file STK_ArrayByVectorProduct.h.

Member Function Documentation

mult()

template<typename Lhs , typename Rhs , typename Result >

template<class SubLhs , class SubRhs >

static void STK::hidden::MultPointArray< Lhs, Rhs, Result >::mult ( SubLhs const &  l, SubRhs const &  r, Result &  res  )

inlinestatic

Compute the product res = l*r with l a point (a row-vector) and r a matrix using the standard formula .

Definition at line 133 of file STK_ArrayByVectorProduct.h.

  {
    // loop over the column of the right hand side
    for (int j=r.beginCols(); j< r.endCols(); ++j)
    {
      Type sum(0); // compute dot product
      for (int i=l.begin(); i<l.end(); ++i) { sum += l.elt(i) * r.elt(i,j); }
      res.elt(j) += sum;
    }
  }

References STK::sum().

Referenced by STK::hidden::MultPointArray< Lhs, Rhs, Result >::run().

run()

template<typename Lhs , typename Rhs , typename Result >

static void STK::hidden::MultPointArray< Lhs, Rhs, Result >::run ( ExprBase< Lhs > const &  l, ExprBase< Rhs > const &  r, Result &  res  )

inlinestatic

Compute the product res = l*r with l a point (a row-vector) and r a matrix using column decomposition.

Definition at line 146 of file STK_ArrayByVectorProduct.h.

  {
    pointByArrayRowRange rowRange(r.beginRows(), pointByArrayRowSize_);
    for(;rowRange.end()<r.endRows(); rowRange += pointByArrayRowSize_)
    {
      PointByArrayColRange colRange(r.beginCols(), pointByArrayColSize_); // create fixed size col range
      for(; colRange.end() < r.endCols(); colRange += pointByArrayColSize_)
      {
        mult(l.sub(rowRange), r.sub(rowRange, colRange), res);
      }
      Range lastColRange(colRange.begin(), r.lastIdxCols(), 0);
      mult(l.sub(rowRange), r.sub(rowRange, lastColRange ), res);
    }
    // end
    Range lastRowRange(rowRange.begin(), r.lastIdxRows(), 0);
    PointByArrayColRange colRange(r.beginCols(), pointByArrayColSize_); // create fixed size col range
    for(; colRange.end() < r.endCols(); colRange += pointByArrayColSize_)
    { mult(l.sub(lastRowRange), r.sub(lastRowRange, colRange), res);}
    // last term
    Range lastColRange(colRange.begin(), r.lastIdxCols(), 0);
    mult(l.sub(lastRowRange), r.sub(lastRowRange, lastColRange ), res);
  }

References STK::hidden::MultPointArray< Lhs, Rhs, Result >::mult(), STK::hidden::pointByArrayColSize_, and STK::hidden::pointByArrayRowSize_.

---

The documentation for this struct was generated from the following file:

• STK_ArrayByVectorProduct.h