00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00056 #ifndef __vtkTensor_h
00057 #define __vtkTensor_h
00058
00059 #include "vtkObject.h"
00060
00061 class VTK_EXPORT vtkTensor : public vtkObject
00062 {
00063 public:
00064 static vtkTensor *New();
00065 vtkTypeMacro(vtkTensor,vtkObject);
00066
00068 void Initialize();
00069
00071 float GetComponent(int i, int j) {return this->T[i+3*j];};
00072
00074 void SetComponent(int i, int j, float v) {if (i > 2 || j > 2) {vtkErrorMacro("trying to set tensor component i or j > 2: i = " << i << ", j = " << j); return;}; this->T[i+3*j] = v;};
00075
00077 void AddComponent(int i, int j, float v) { if (i > 2 || j > 2) {vtkErrorMacro("trying to add tensor component i or j > 2: i = " << i << ", j = " << j); return;}; this->T[i+3*j] += v;};
00078
00081 float *GetColumn(int j) { if (j > 2) {vtkErrorMacro("trying to get tensor column j > 2: j = " << j); return NULL;}; return this->T + 3*j;};
00082
00084 void DeepCopy(vtkTensor *t);
00085
00086 #ifndef VTK_REMOVE_LEGACY_CODE
00087
00088 void DeepCopy(vtkTensor &t);
00089 #endif
00090
00092 operator float*() {return this->T;};
00093
00095 float *T;
00096
00097 protected:
00098 vtkTensor();
00099 ~vtkTensor() {};
00100 vtkTensor(const vtkTensor&) {};
00101 void operator=(const vtkTensor&) {};
00102
00103 float Storage[9];
00104 };
00105
00106 inline void vtkTensor::Initialize()
00107 {
00108 for (int j=0; j<3; j++)
00109 {
00110 for (int i=0; i<3; i++)
00111 {
00112 this->T[i+j*3] = 0.0;
00113 }
00114 }
00115 }
00116
00117 inline void vtkTensor::DeepCopy(vtkTensor *t)
00118 {
00119 for (int j=0; j < 3; j++)
00120 {
00121 for (int i=0; i < 3; i++)
00122 {
00123 this->T[i+3*j] = t->T[i+3*j];
00124 }
00125 }
00126 }
00127
00128 #endif