zcy
/
openEMS
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

SAR: increase accuracy in intermediate values 

Signed-off-by: Thorsten Liebig <Thorsten.Liebig@gmx.de>
pull/1/head
Thorsten Liebig 2013-01-30 14:16:23 +01:00
parent 944fd3c528
commit 920118ea5f
2 changed files with 24 additions and 24 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

38
tools/sar_calculation.cpp
View File

@ -205,11 +205,11 @@ float*** SAR_Calculation::CalcLocalSAR(float*** SAR)
} }
int SAR_Calculation::FindFittingCubicalMass(unsigned int pos[3], float box_size, unsigned int start[3], unsigned int stop[3], int SAR_Calculation::FindFittingCubicalMass(unsigned int pos[3], float box_size, unsigned int start[3], unsigned int stop[3],
float partial_start[3], float partial_stop[3], float &mass, float &volume, float &bg_ratio, int disabledFace, bool ignoreFaceValid) float partial_start[3], float partial_stop[3], double &mass, double &volume, double &bg_ratio, int disabledFace, bool ignoreFaceValid)
{ {
unsigned int mass_iterations = 0; unsigned int mass_iterations = 0;
float old_mass=0; double old_mass=0;
float old_box_size=0; double old_box_size=0;
bool face_valid; bool face_valid;
bool mass_valid; bool mass_valid;
bool voxel_valid; bool voxel_valid;
@ -265,8 +265,8 @@ int SAR_Calculation::FindFittingCubicalMass(unsigned int pos[3], float box_size,
return -1; return -1;
} }
bool SAR_Calculation::GetCubicalMass(unsigned int pos[3], float box_size, unsigned int start[3], unsigned int stop[3], bool SAR_Calculation::GetCubicalMass(unsigned int pos[3], double box_size, unsigned int start[3], unsigned int stop[3],
float partial_start[3], float partial_stop[3], float &mass, float &volume, float &bg_ratio, int disabledFace) float partial_start[3], float partial_stop[3], double &mass, double &volume, double &bg_ratio, int disabledFace)
{ {
if ((box_size<=0) || isnan(box_size) || isinf(box_size)) if ((box_size<=0) || isnan(box_size) || isinf(box_size))
{ {
@ -343,8 +343,8 @@ bool SAR_Calculation::GetCubicalMass(unsigned int pos[3], float box_size, unsign
mass = 0; mass = 0;
volume = 0; volume = 0;
float bg_volume=0; double bg_volume=0;
float weight[3]; double weight[3];
unsigned int f_pos[3]; unsigned int f_pos[3];
bool face_intersect[6] = {false,false,false,false,false,false}; bool face_intersect[6] = {false,false,false,false,false,false};
for (f_pos[0]=start[0];f_pos[0]<=stop[0];++f_pos[0]) for (f_pos[0]=start[0];f_pos[0]<=stop[0];++f_pos[0])
@ -400,9 +400,9 @@ bool SAR_Calculation::GetCubicalMass(unsigned int pos[3], float box_size, unsign
float SAR_Calculation::CalcCubicalSAR(float*** SAR, unsigned int pos[3], unsigned int start[3], unsigned int stop[3], float partial_start[3], float partial_stop[3], bool assignUsed) float SAR_Calculation::CalcCubicalSAR(float*** SAR, unsigned int pos[3], unsigned int start[3], unsigned int stop[3], float partial_start[3], float partial_stop[3], bool assignUsed)
{ {
float power_mass=0; double power_mass=0;
float mass=0; double mass=0;
float weight[3]; double weight[3];
unsigned int f_pos[3]; unsigned int f_pos[3];
for (f_pos[0]=start[0];f_pos[0]<=stop[0];++f_pos[0]) for (f_pos[0]=start[0];f_pos[0]<=stop[0];++f_pos[0])
{ {
@ -489,13 +489,13 @@ float*** SAR_Calculation::CalcAveragedSAR(float*** SAR)
m_Vx_Used = Create3DArray<bool>(m_numLines); m_Vx_Used = Create3DArray<bool>(m_numLines);
m_Vx_Valid = Create3DArray<bool>(m_numLines); m_Vx_Valid = Create3DArray<bool>(m_numLines);
float voxel_volume; double voxel_volume;
float total_mass; double total_mass;
unsigned int start[3]; unsigned int start[3];
unsigned int stop[3]; unsigned int stop[3];
float partial_start[3]; float partial_start[3];
float partial_stop[3]; float partial_stop[3];
float bg_ratio; double bg_ratio;
int EC=0; int EC=0;
// debug counter // debug counter
@ -551,7 +551,7 @@ float*** SAR_Calculation::CalcAveragedSAR(float*** SAR)
{ {
cerr << "Number of invalid cubes (case 1): " << cnt_case1 << endl; cerr << "Number of invalid cubes (case 1): " << cnt_case1 << endl;
cerr << "Number of invalid cubes (case 2): " << cnt_case2 << endl; cerr << "Number of invalid cubes (case 2): " << cnt_case2 << endl;
cerr << "Number of invalid cubes (falied to converge): " << cnt_NoConvergence << endl; cerr << "Number of invalid cubes (failed to converge): " << cnt_NoConvergence << endl;
} }
// count all used and unused etc. + special handling of unused voxels!! // count all used and unused etc. + special handling of unused voxels!!
@ -570,10 +570,10 @@ float*** SAR_Calculation::CalcAveragedSAR(float*** SAR)
++m_Unused; ++m_Unused;
SAR[pos[0]][pos[1]][pos[2]] = 0; SAR[pos[0]][pos[1]][pos[2]] = 0;
float unused_volumes[6]; double unused_volumes[6];
float unused_SAR[6]; float unused_SAR[6];
float min_Vol=FLT_MAX; double min_Vol=FLT_MAX;
// special handling of unused voxels: // special handling of unused voxels:
for (int n=0;n<6;++n) for (int n=0;n<6;++n)
@ -614,18 +614,18 @@ float*** SAR_Calculation::CalcAveragedSAR(float*** SAR)
return SAR; return SAR;
} }
float SAR_Calculation::CellVolume(unsigned int pos[3]) double SAR_Calculation::CellVolume(unsigned int pos[3])
{ {
if (m_cell_volume) if (m_cell_volume)
return m_cell_volume[pos[0]][pos[1]][pos[2]]; return m_cell_volume[pos[0]][pos[1]][pos[2]];
float volume=1; double volume=1;
for (int n=0;n<3;++n) for (int n=0;n<3;++n)
volume*=m_cellWidth[n][pos[n]]; volume*=m_cellWidth[n][pos[n]];
return volume; return volume;
} }
float SAR_Calculation::CellMass(unsigned int pos[3]) double SAR_Calculation::CellMass(unsigned int pos[3])
{ {
return m_cell_density[pos[0]][pos[1]][pos[2]]*CellVolume(pos); return m_cell_density[pos[0]][pos[1]][pos[2]]*CellVolume(pos);
} }

10
tools/sar_calculation.h
View File

@ -104,16 +104,16 @@ protected:
float*** CalcAveragedSAR(float*** SAR); float*** CalcAveragedSAR(float*** SAR);
int FindFittingCubicalMass(unsigned int pos[3], float box_size, unsigned int start[3], unsigned int stop[3], int FindFittingCubicalMass(unsigned int pos[3], float box_size, unsigned int start[3], unsigned int stop[3],
float partial_start[3], float partial_stop[3], float &mass, float &volume, float &bg_ratio, int disabledFace=-1, bool ignoreFaceValid=false); float partial_start[3], float partial_stop[3], double &mass, double &volume, double &bg_ratio, int disabledFace=-1, bool ignoreFaceValid=false);
bool GetCubicalMass(unsigned int pos[3], float box_size, unsigned int start[3], unsigned int stop[3], bool GetCubicalMass(unsigned int pos[3], double box_size, unsigned int start[3], unsigned int stop[3],
float partial_start[3], float partial_stop[3], float &mass, float &volume, float &bg_ratio, int disabledFace=-1); float partial_start[3], float partial_stop[3], double &mass, double &volume, double &bg_ratio, int disabledFace=-1);
float CalcCubicalSAR(float*** SAR, unsigned int pos[3], unsigned int start[3], unsigned int stop[3], float partial_start[3], float partial_stop[3], bool assignUsed=false); float CalcCubicalSAR(float*** SAR, unsigned int pos[3], unsigned int start[3], unsigned int stop[3], float partial_start[3], float partial_stop[3], bool assignUsed=false);
/****** end SAR averaging and all necessary methods ********/ /****** end SAR averaging and all necessary methods ********/
bool CheckValid(); bool CheckValid();
float CellVolume(unsigned int pos[3]); double CellVolume(unsigned int pos[3]);
float CellMass(unsigned int pos[3]); double CellMass(unsigned int pos[3]);
}; };
#endif // SAR_CALCULATION_H #endif // SAR_CALCULATION_H