00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019 #ifndef GEOM_MAT2_HH
00020 #define GEOM_MAT2_HH
00021
00022 #include <cstddef>
00023 #include <ostream>
00024 #include <cassert>
00025 #include <stdexcept>
00026
00027 #include <boost/operators.hpp>
00028
00029 #include <ost/geom/module_config.hh>
00030
00031 namespace geom {
00032
00033 class DLLEXPORT_OST_GEOM Mat2:
00034 private boost::equality_comparable<Mat2>,
00035 private boost::additive<Mat2>,
00036 private boost::multiplicative2<Mat2, Real>
00037 {
00038 public:
00039 static Mat2 Identity();
00040
00042 Mat2();
00043
00045
00053 Mat2(Real i00, Real i01, Real i10, Real i11);
00054
00056 Mat2(const Mat2& m);
00057
00059 explicit Mat2(const Real[4]);
00060
00062 Mat2& operator=(const Mat2& m);
00063
00065 bool operator==(const Mat2& rhs) const;
00066
00067
00069 const Real& At(std::size_t r, std::size_t c) const
00070 {
00071 if (r>1 || c>1) {
00072 throw std::out_of_range("indices must be smaller than 2");
00073 }
00074 return data_[r][c];
00075 }
00076
00078 Real& At(std::size_t r, std::size_t c)
00079 {
00080 if (r>1 || c>1) {
00081 throw std::out_of_range("indices must be smaller than 2");
00082 }
00083 return data_[r][c];
00084 }
00085
00087 Real& operator()(std::size_t r, std::size_t c)
00088 {
00089 assert(r<2 && c<2);
00090 return data_[r][c];
00091 }
00093 const Real& operator()(std::size_t r, std::size_t c) const
00094 {
00095 assert(r<2 && c<2);
00096 return data_[r][c];
00097 }
00098
00100 Mat2& operator+=(const Mat2& rhs);
00102 Mat2& operator-=(const Mat2& rhs);
00103
00104 Mat2& operator*=(const Real d);
00105 Mat2& operator/=(const Real d);
00106
00107 Mat2& operator*=(const Mat2& m);
00108
00109 Real* Data() {return &data_[0][0];}
00110 const Real* Data() const {return &data_[0][0];}
00111
00112 private:
00113 Real data_[2][2];
00114
00115 void set(Real i00, Real i01, Real i10, Real i11);
00116 };
00117
00118 DLLEXPORT_OST_GEOM std::ostream& operator<<(std::ostream& os, const Mat2& m);
00119
00120 }
00121
00122
00123 #endif