#include "xyz2shc.h"

int XYZ2SHC::InitShc(int given_order){
	if (given_order <= 0){
		cerr << BOLDRED << "error ==> " << RESET << "spherical order must be bigger than zero: " << given_order << endl;
		return -1;
	}
	else{
		shc_order_ = given_order;
		shc_num_ = (shc_order_+1)*(shc_order_+2);
		half_shc_num_ = (shc_order_+1)*(shc_order_+2)/2;
		//初始化shc_value_
		shc_value_.resize(shc_num_,0.0);
		//初始化Anm_和Bnm_
		Anm_ = get_Anm_array(shc_order_+1);
		Bnm_ = get_Bnm_array(shc_order_+1);
		return 0;
	}
}