20 lines
625 B
C++
20 lines
625 B
C++
|
#include "xyz2shc.h"
|
||
|
|
|
||
|
|
int XYZ2SHC::OutShc(char* filename){
|
||
|
|
int start_index,flow_index;
|
||
|
|
|
||
|
|
ofstream outfile;
|
||
|
|
if (open_outfile(outfile,filename)) return -1;
|
||
|
|
|
||
|
|
for (int i = 0; i < shc_order_+1; i++){
|
||
|
|
//按列优先排列将球谐系数输出 由列号计算数组中的起始位置
|
||
|
|
start_index = int(0.5*i*(2*shc_order_+3-i));
|
||
|
|
for (int j = 0; j < shc_order_+1; j++){
|
||
|
|
//根据行号计算偏移量
|
||
|
|
flow_index = j - i;
|
||
|
|
outfile << j << " " << i << setprecision(16) << " " << shc_value_[start_index+flow_index] << " " << shc_value_[start_index+flow_index+half_shc_num_] << endl;
|
||
|
|
}
|
||
|
|
}
|
||
|
|
outfile.close();
|
||
|
|
return 0;
|
||
|
|
}
|