/home/hauberg/Dokumenter/Capture/humim-tracker-0.1/src/OpenTissue/OpenTissue/core/math/big/big_is_symmetric.h

Go to the documentation of this file.

#ifndef OPENTISSUE_CORE_MATH_BIG_BIG_IS_SYMMETRIC_H
#define OPENTISSUE_CORE_MATH_BIG_BIG_IS_SYMMETRIC_H
//
// OpenTissue Template Library
// - A generic toolbox for physics-based modeling and simulation.
// Copyright (C) 2008 Department of Computer Science, University of Copenhagen.
//
// OTTL is licensed under zlib: http://opensource.org/licenses/zlib-license.php
//
#include <OpenTissue/configuration.h>

#include <boost/cast.hpp>             // needed for boost::numeric_cast

namespace OpenTissue
{
  namespace math
  {
    namespace big
    {

      template<typename matrix_type>
      inline bool is_symmetric(  matrix_type const & A  )
      {
        using std::fabs;

        typedef typename matrix_type::value_type                     value_type;
        typedef typename matrix_type::size_type                      size_type;

        size_type  const & m       = A.size1();
        size_type  const & n       = A.size2();

        assert( m>0         || !"is_symmetric(): m was out of range");
        assert( n>0         || !"is_symmetric(): n was out of range");
        assert( m==n        || !"is_symmetric(): m and n was not equal");

        value_type const precision = ::boost::numeric_cast<value_type>(10e-6);

        for ( size_type i = 0; i < n; ++i )
          for ( size_type j = i+1; j < n; ++j )
          {
            value_type tmp = A(i,j) - A(j,i);
            if( fabs(tmp) >  precision )
              return false;
          }
          return true;
      }


    } // end of namespace big
  } // end of namespace math
} // end of namespace OpenTissue

// OPENTISSUE_CORE_MATH_BIG_BIG_IS_SYMMETRIC_H
#endif