positive_real_power_explicit.cpp

来自「矩阵运算源码最新版本」· C++ 代码 · 共 57 行

#include <iostream>
#include <cmath>

#include <boost/numeric/linear_algebra/operators.hpp>
#include <boost/numeric/linear_algebra/concepts.hpp>
#include <boost/numeric/linear_algebra/identity.hpp>
#include <boost/numeric/linear_algebra/inverse.hpp>
#include <boost/numeric/linear_algebra/is_invertible.hpp>

#include <libs/numeric/linear_algebra/test/power_explicit.hpp>
#include <libs/numeric/linear_algebra/test/positive_real.hpp>
#include <libs/numeric/linear_algebra/test/positive_real_power.hpp>

// User defined data types and operators

using mtl::positive_real;
  
// We assume that 0 and infinity is not to guarantee invertibility

# ifdef __GXX_CONCEPTS__
  namespace algebra { 
      concept_map SemiGroup< semigroup_mult, positive_real > {};
      concept_map Monoid< monoid_mult, positive_real > {};
      concept_map Group< group_mult, positive_real > {};
  }
# endif


 

int main(int, char* []) 
{
    using mtl::power;
    using math::mult;

    positive_real          value(1.1), zero(0.0);

    compute_power(value, 777, mult<positive_real>(), "Magma");
    std::cout << '\n';
 
    compute_power(value, 777, semigroup_mult(), "SemiGroup");
    std::cout << '\n';

    compute_power(value, 777, monoid_mult(), "Monoid");
    compute_power(value, -777, monoid_mult(), "Monoid");
    std::cout << '\n';
 
    compute_power(value, 777, group_mult(), "Group");
    compute_power(value, -777, group_mult(), "Group");
    compute_power(zero, 777, group_mult(), "Group");
    compute_power(zero, -777, group_mult(), "Group");
    std::cout << '\n';
   
    return 0;
}