---

CTBdot.icc

Go to the documentation of this file.

/*----------------------------------------------------------------------------*/
/* C Tool Box: Designed and implemented by:                                   */
/*    Walter F.J. Mueller   Gesellschaft fuer Schwerionenforschung (GSI)      */
/*                          Planckstrasse 1, D-64291 Darmstadt, Germany       */
/*                  Email:  W.F.J.Mueller@gsi.de                              */
/*                  WWW:    http://www-kp3.gsi.de/www/kp3/people/mueller.html */
/*----------------------------------------------------------------------------*/

#include "CTBtraits.hxx"
#include "CTBnum.hxx"
#include "CTBexceptionIndexRange.hxx"
#include "CTBexceptionConformity.hxx"

//------------------------------------------+-----------------------------------
template <class T>
inline T CTBdot(const CTBvector<T>& v1, const CTBvector<T>& v2)
{
  const T* p_v1 = v1;
  const T* p_v2 = v2;
  CTBtraits<T>::Sum_t f_sum = 0.;

#ifdef CTB__ConformityCheck
  CTBint n = v1.Size();
  if (n != v2.Size())
    throw CTBexceptionConformity(n, v2.Size(),
                                 "CTBdot(CTBvector&,CTBvector&)");
#else
  CTBint n = CTBmin(v1.Size(), v2.Size());
#endif

  for (CTBint i = 0; i < n; i++) f_sum += *p_v1++ * *p_v2++;
  return f_sum;
}

//------------------------------------------+-----------------------------------
template <class T>
inline T CTBdot(const CTBvector<T>& v1, const CTBvector<T>& v2,
                CTBint i_beg, CTBint i_end)
{
  const T* p_v1 = &v1(i_beg);
  const T* p_v2 = &v2(i_beg);
  CTBtraits<T>::Sum_t f_sum = 0.;

#ifdef CTB__IndexCheck
  CTBint n = CTBmin(v1.Size(), v2.Size());
  if (i_end >= i_beg && (i_end < 0 || i_end >= n)) 
    throw CTBexceptionIndexRange(i_end,n,
                                 "CTBdot(CTBvector&,CTBvector&,CTBint,CTBint)");
#endif

  for (CTBint i = i_beg; i <= i_end; i++) f_sum += *p_v1++ * *p_v2++;
  return f_sum;
}

//------------------------------------------+-----------------------------------
template <class T>
inline T CTBdot(const CTBmatrix<T>& m, CTBint i_row, const CTBvector<T>& v)
{
  const T* p_m = &m(i_row,0);
  const T* p_v = v;
  CTBtraits<T>::Sum_t f_sum = 0.;

#ifdef CTB__ConformityCheck
  CTBint n = m.NColumn();
  if (n != v.Size())
    throw CTBexceptionConformity(n, v.Size(),
                                 "CTBdot(CTBmatrix&,CTBint,CTBvector&)");
#else
  CTBint n = CTBmin(m.NColumn(), v.Size());
#endif

  for (CTBint i = 0; i < n; i++) f_sum += *p_m++ * *p_v++;
  return f_sum;
}

//------------------------------------------+-----------------------------------
template <class T>
inline T CTBdot(const CTBmatrix<T>& m, CTBint i_row, const CTBvector<T>& v,
                CTBint i_beg, CTBint i_end)
{
  const T* p_m = &m(i_row,i_beg);
  const T* p_v = &v(i_beg);
  CTBtraits<T>::Sum_t f_sum = 0.;

#ifdef CTB__IndexCheck
  CTBint n = CTBmin(m.NColumn(), v.Size());
  if (i_end >= i_beg && (i_end < 0 || i_end >= n))
    throw CTBexceptionIndexRange(i_end,n,
            "CTBdot(CTBmatrix&,CTBint,CTBvector&,CTBint,CTBint)");
#endif

  for (CTBint i = i_beg; i <= i_end; i++) f_sum += *p_m++ * *p_v++;
  return f_sum;
}

//------------------------------------------+-----------------------------------
template <class T>
inline T CTBdot(CTBint i_col, const CTBmatrix<T>& m, const CTBvector<T>& v)
{
  const T* p_m  = &m(0,i_col);
  const T* p_v  = v;
  CTBint   i_nc = m.NColumn();
  CTBtraits<T>::Sum_t f_sum = 0.;

#ifdef CTB__ConformityCheck
  CTBint n = m.NRow();
  if (n != v.Size())
    throw CTBexceptionConformity(n, v.Size(),
                                 "CTBdot(CTBint,CTBmatrix&,CTBvector&)");
#else
  CTBint n = CTBmin(m.NRow(), v.Size());
#endif

  for (CTBint i = 0; i < n; i++) {
    f_sum += *p_m * *p_v++;
    p_m   += i_nc;
  }
  
  return f_sum;
}

//------------------------------------------+-----------------------------------
template <class T>
inline T CTBdot(CTBint i_col, const CTBmatrix<T>& m, const CTBvector<T>& v,
                CTBint i_beg, CTBint i_end)
{
  const T* p_m  = &m(i_beg,i_col);
  const T* p_v  = &v(i_beg);
  CTBint   i_nc = m.NColumn();
  CTBtraits<T>::Sum_t f_sum = 0.;

#ifdef CTB__IndexCheck
  CTBint n = CTBmin(m.NRow(), v.Size());
  if (i_end >= i_beg && (i_end < 0 || i_end >= n))
    throw CTBexceptionIndexRange(i_end,n,
            "CTBdot(CTBint,CTBmatrix&,CTBvector&,CTBint,CTBint)");
#endif

  for (CTBint i = i_beg; i <= i_end; i++) {
    f_sum += *p_m * *p_v++;
    p_m   += i_nc;
  }
  
  return f_sum;
}

//------------------------------------------+-----------------------------------
template <class T>
inline T CTBdot(const CTBmatrix<T>& m1, CTBint i_row, 
                CTBint i_col, const CTBmatrix<T>& m2)
{
  const T* p_m1 = &m1(i_row,0);
  const T* p_m2 = &m2(0,i_col);
  CTBint i_nc2 = m2.NColumn();
  CTBtraits<T>::Sum_t f_sum = 0.;

#ifdef CTB__ConformityCheck
  CTBint n = m1.NColumn();
  if (n != m2.NRow()) 
    throw CTBexceptionConformity(n, m2.NRow(),
                                 "CTBdot(CTBmatrix&,CTBint,CTBint,CTBmatrix&)");
#else
  CTBint n = CTBmin(m1.NColumn(), m2.NRow());
#endif

  for (CTBint i = 0; i < n; i++) {
    f_sum += *p_m1++ * *p_m2;
    p_m2  += i_nc2;
  }
  
  return f_sum;
}

//------------------------------------------+-----------------------------------
template <class T>
inline T CTBdot(const CTBmatrix<T>& m1, CTBint i_row, 
                CTBint i_col, const CTBmatrix<T>& m2,
                CTBint i_beg, CTBint i_end)
{
  const T* p_m1 = &m1(i_row,i_beg);
  const T* p_m2 = &m2(i_beg,i_col);
  CTBint i_nc2 = m2.NColumn();
  CTBtraits<T>::Sum_t f_sum = 0.;

#ifdef CTB__IndexCheck
  CTBint n = CTBmin(m1.NColumn(),m2.NRow() );
  if (i_end >= i_beg && (i_end < 0 || i_end >= n))
    throw CTBexceptionIndexRange(i_end,n,
            "CTBdot(CTBmatrix&,CTBint,CTBint,CTBmatrix&,CTBint,CTBint)");
#endif

  for (CTBint i = i_beg; i <= i_end; i++) {
    f_sum += *p_m1++ * *p_m2;
    p_m2  += i_nc2;
  }
  
  return f_sum;
}

---