9 for (
unsigned int i=0;i<3;i++)
mA[i]=0.;
27 for (
unsigned int i=0;i<3;i++){
33 for (
unsigned int i=0;i<3;i++){
44 for (
unsigned int i=0;i<3;i++){
64 cout <<
"DMatrix3x1:" <<endl;
65 cout <<
" | 0 |" <<endl;
66 cout <<
"----------------------" <<endl;
67 for (
unsigned int i=0;i<3;i++){
68 cout <<
" "<<i<<
" |" << setw(11)<<setprecision(4) <<
mA[i] << endl;
90 for (
unsigned int j=0;j<2;j++){
91 mA->v[j]=_mm_setzero_pd();
109 __m128d GetV(
int pair)
const{
136 return DMatrix3x1(_mm_add_pd(GetV(0),m2.GetV(0)),
137 _mm_add_pd(GetV(1),m2.GetV(1)));
140 mA->v[0]=_mm_add_pd(GetV(0),m2.GetV(0));
141 mA->v[1]=_mm_add_pd(GetV(1),m2.GetV(1));
147 return DMatrix3x1(_mm_sub_pd(GetV(0),m2.GetV(0)),
148 _mm_sub_pd(GetV(1),m2.GetV(1)));
152 return mA->d[0]*
mA->d[0]+
mA->d[1]*
mA->d[1]+
mA->d[2]*
mA->d[2];
156 cout <<
"DMatrix3x1:" <<endl;
157 cout <<
" | 0 |" <<endl;
158 cout <<
"----------------------" <<endl;
159 for (
unsigned int i=0;i<3;i++){
160 cout <<
" "<<i<<
" |" << setw(11)<<setprecision(4) <<
mA->d[i] << endl;
177 scale=_mm_set1_pd(c);
179 _mm_mul_pd(scale,M.GetV(1)));
DMatrix3x1 operator+(const DMatrix3x1 &m2) const
DMatrix2x1 operator*(const double c, const DMatrix2x1 &M)
#define ALIGNED_16_BLOCK_PTR(TYPE, NUM, PTR)
#define ALIGNED_16_BLOCK_WITH_PTR(TYPE, NUM, PTR)
double & operator()(int row)
DMatrix3x1 & operator+=(const DMatrix3x1 &m2)
DMatrix3x1 & operator=(const DMatrix3x1 &m2)
DMatrix3x1(double a1, double a2, double a3)
DMatrix3x1(const DMatrix3x1 &m2)
#define ALIGNED_16_BLOCK(TYPE, NUM, PTR)
DMatrix3x1 operator-(const DMatrix3x1 &m2) const