basics/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