arrays.hh

Go to the documentation of this file.

 
#ifndef arrays_hh
#define arrays_hh
 
#include "toolbox/array.hh"
 
namespace concepts {
 
  // ****************************************************** ArrayDeterminant **
 
  template <class F, int dim>
  class ArrayDeterminant : public Array<F> {
  public:
    ArrayDeterminant(const concepts::Array<concepts::Mapping<F,dim> >* array) :
      Array<F>(array ? array->size() : 0) {
      if (array) {
  F* foo = data_;
  const concepts::Mapping<F,dim>* matrix = 
    (const concepts::Mapping<F,dim>*)*array;
  for(uint i = array->size(); i--;)
    *foo++ = matrix++->determinant();
      }
    }
  };
 
  // **************************************************** ArrayMatrixInverse **
 
  template <class F, int dim>
  class ArrayMatrixInverse : public Array<concepts::Mapping<F,dim> > {
  public:
    ArrayMatrixInverse
    (const concepts::Array<concepts::Mapping<F,dim> >* array) :
      Array<concepts::Mapping<F,dim> >(array ? array->size() : 0) {
      if (array) {
  concepts::Mapping<F,dim>* foo = data_;
  const concepts::Mapping<F,dim>* matrix = 
    (const concepts::Mapping<F,dim>*)*array;
  for(uint i = array->size(); i--;)
    *foo++ = matrix++->inverse();
      }
    }
  };
 
  // ************************************************** ArrayMatrixTranspose **
 
  template <class F, int dim>
  class ArrayMatrixTranspose : public Array<concepts::Mapping<F, dim> > {
  public:
    ArrayMatrixTranspose
    (const concepts::Array<concepts::Mapping<F,dim> >* array) :
      Array<concepts::Mapping<F,dim> >(array ? array->size() : 0) {
      if (array) {
  concepts::Mapping<F,dim>* foo = data_;
  const concepts::Mapping<F,dim>* matrix = 
    (const concepts::Mapping<F,dim>*)*array;
  for(uint i = array->size(); i--;)
    *foo++ = matrix++->transpose();
      }
    }
  };
 
  // ******************************************************* ArrayReciprocal **
 
  template <class F>
  class ArrayReciprocal : public Array<F> {
  public:
    ArrayReciprocal(const concepts::Array<F>* array) :
      Array<F>(array ? array->size() : 0) {
      if (array) {
  F* foo = data_; const F* a = *array;
  for(uint i = array->size(); i--;)
    *foo++ = 1.0 / *a++;
      }
    }
  };
 
  // ******************************************************* ArrayGramMatrix **
 
  template <class F, int dim>
  class ArrayGramMatrix : public Array<concepts::Mapping<F,dim> > {
  public:
    ArrayGramMatrix(const concepts::Array<concepts::Mapping<F,dim> >& array) :
      Array<concepts::Mapping<F,dim> >(array.size()) {
      concepts::Mapping<F,dim>* foo = data_, *matrix = *array;
      for(uint i = array.size(); i--;)
  *foo++ = *matrix * matrix++->transpose();
    }
  };
 
} // namespace concepts
 
#endif