synaps/linalg/MATRIX3D.h

Go to the documentation of this file.

/*********************************************************************
*      This file is part of the source code of SYNAPS library        *
*   Author(s): L. Deschamps                                          *
*              G. Gatellier                                          *
*              B. Mourrain, GALAAD, INRIA                            *
**********************************************************************/
#ifndef synaps_linalg_TENSOR3D_H
#define synaps_linalg_TENSOR3D_H

/*********************************************************************/

#ifdef HAVE_CONFIG_H
#include <synaps/init.h>
#endif 

#include <iostream>
#include <assert.h>

#include <synaps/linalg/MethodName.h>

__BEGIN_NAMESPACE_SYNAPS

/*********************************************************************/
namespace TENSOR3D {
 
  typedef unsigned size_type;
  
  template < class R >
  inline std::ostream & print( std::ostream & os, const R & m ) {
    typedef typename R::size_type  size_type;
    
    for ( size_type k = 0 ; k < m.nbhgt() ; k++ ) {
      for ( size_type i = 0 ; i < m.nbrow(); ++i ){
        os << "[ ";
        for ( size_type j = 0 ; j < m.nbcol() ; j++ )
          os<< m( i, j, k ) << " ";
        os << "]"<< std::endl;
      }
      os << std::endl;
    }
    return os;
  }

  template < class R >
  inline std::ostream & printsign( std::ostream & os, const R & m ) {
    typedef typename R::size_type  size_type;
    
    for ( size_type k = 0 ; k < m.nbhgt() ; k++ ) {
      for ( size_type i = 0 ; i < m.nbrow(); ++i ){     
        for( size_type u=0;u<k;u++) os<<" ";
        os << "[";
        for ( size_type j = 0 ; j < m.nbcol() ; j++ )
          if(m(i,j,k)<0) os<< "-"; else if (m(i,j,k)>0) os<< "+"; else os << "0"; 
        os << "]"<< std::endl;
      }
      os << std::endl;
    }
    return os;
  }

  
  template< class R >  inline
  std::istream & read( std::istream & is, R & M ){
    typedef typename R::size_type size_type;
    size_type m, n, o;
    is >> m >> n >> o;
    M.rep().resize( m, n );
 
    for ( size_type k = 0; k < o ; k++ )
      for ( size_type i = 0 ; i < m ; i++ )
        for (size_type j = 0 ; j < n ; j++ ) {
          is >> M(i,j);
        }
    return( is );
  }

  template< class R, class C > inline
  void init( R & a, const C & c) { 
    typedef typename R::size_type size_type;

    for ( size_type i = 0 ; i < a.nbrow() ; i++ ) 
      for ( size_type j = 0 ; j < a.nbcol(); j++ )
        for ( size_type k = 0 ; k < a.nbhgt() ; k++)
          a( i, j, k ) = c;
  }

  template < class R1, class R2 >
  inline void assign( R1 & r, const R2 & m ) {
    typedef typename R1::size_type  size_type;
    
    for ( size_type i = 0 ; i < m.nbrow(); ++i )
      for (size_type j = 0 ; j < m.nbcol(); j++ )
        for( size_type k = 0 ; k < m.nbhgt() ; k++ )
          r( i, j, k ) = m( i, j, k );
  }
  
  
} // TENSOR3D

/*********************************************************************/

__END_NAMESPACE_SYNAPS

#endif // synaps_linalg_TENSOR3D_H