Go to the documentation of this file.00001
00002 #ifndef vnl_linear_operators_3_h_
00003 #define vnl_linear_operators_3_h_
00004 #ifdef VCL_NEEDS_PRAGMA_INTERFACE
00005 #pragma interface
00006 #endif
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023 #include <vnl/vnl_double_3.h>
00024 #include <vnl/vnl_double_3x3.h>
00025
00026
00027
00028 inline
00029 vnl_double_3 operator* (const vnl_double_3x3& A, const vnl_double_3& x)
00030 {
00031 const double* a = A.data_block();
00032 double r0 = a[0] * x[0] + a[1] * x[1] + a[2] * x[2];
00033 double r1 = a[3] * x[0] + a[4] * x[1] + a[5] * x[2];
00034 double r2 = a[6] * x[0] + a[7] * x[1] + a[8] * x[2];
00035 return vnl_double_3(r0, r1, r2);
00036 }
00037
00038
00039
00040 inline
00041 vnl_double_3 operator+ (const vnl_double_3& a, const vnl_double_3& b)
00042 {
00043 double r0 = a[0] + b[0];
00044 double r1 = a[1] + b[1];
00045 double r2 = a[2] + b[2];
00046 return vnl_double_3(r0, r1, r2);
00047 }
00048
00049 #endif // vnl_linear_operators_3_h_