int YAPCAReduce<eltype>::load_map(std::istream &in) {
int header_err = this->read_header(in);
if (header_err !=0 && header_err != 312)
return header_err;
int err_flag = read(in,eigen_vectors);
if (err_flag != 0)
return err_flag;
err_flag = read(in,this->eigen_values);
if (err_flag != 0)
return err_flag;
err_flag = read(in,column_means);
if (err_flag != 0)
return err_flag;
if (!in || in.eof() || eigen_vectors.numel()==0 ||
column_means.numel()==0) {
ya_addwarn(310,10,method_id(),"Improperly formatted map file.");
return 310;
}
return header_err;
}
int ya_string_to_mat<eltype>::convert(const string &value, ya_type &mat) {
int err_flag = 0;
int ascii_flag = 0;
string temp=value;
// Transpose the matrix?
bool transposeme=false;
if (temp[temp.length()-1]=='\'') {
temp[temp.length()-1]=' ';
transposeme=true;
}
// Get the rows
size_t tL=temp.length();
for (size_t i=0; i<tL; i++)
if (temp[i]=='[' || temp[i]==']' || temp[i]==',')
temp[i]=' ';
vector<string> row_str;
ya_get_tokens(';',temp,row_str);
if (row_str.empty()) {
mat.clear();
return err_flag;
}
ya_sizet row_size=row_str.size();
if (ya_whitespace(row_str[row_size-1]))
row_size--;
if (row_size==0) {
mat.clear();
return err_flag;
}
// Get the columns
for (ya_sizet i=0; i<row_size; i++) {
vector<string> col_str;
ya_get_tokens(row_str[i], col_str);
if (col_str.empty()) {
ostringstream owarn;
owarn << "Error parsing matrix string:\n" << value << "\nat row token: "
<< row_str[i];
ya_addwarn(307,11,"YA_BaseT",owarn.str());
err_flag = 307;
if (i==0)
return err_flag;
}
YA_DynT rvec;
size_t cS=col_str.size();
for (size_t j=0; j<cS; j++) {
vector<string> e_tok;
ya_get_tokens(':',col_str[j],e_tok);
if (e_tok.size()>3) {
ostringstream owarn;
owarn << "Error parsing matrix string:\n" << value << "\nat token: "
<< col_str[j];
ya_addwarn(307,11,"YA_BaseT",owarn.str());
err_flag=307;
if (i==0)
return err_flag;
}
eltype jump=1;
ya_sizet high_index=1;
if (e_tok.size()==3) {
high_index=2;
ascii_flag=ya_ascii_to(e_tok[1],jump);
}
if (e_tok.size()>1) {
eltype one,two;
ascii_flag=std::max(ascii_flag,ya_ascii_to(e_tok[0],one));
ascii_flag=std::max(ascii_flag,ya_ascii_to(e_tok[high_index],two));
double intp;
if (modf(static_cast<double>(two-one)/double(jump),&intp) != 0)
two=static_cast<eltype>(static_cast<double>(jump)*intp+one);
if ( (one>two && static_cast<int>(jump)>0) ||
(one<two && static_cast<int>(jump)<0) || jump==0) {
ostringstream owarn;
owarn << "Error parsing matrix string:\n" << value << "\nat token: "
<< col_str[j];
ya_addwarn(307,11,"YA_BaseT",owarn.str());
err_flag=307;
if (i==0)
return err_flag;
}
for (eltype k=one; ; k+=jump) {
rvec.push_back(k);
if (k==two)
break;
}
} else {
eltype one;
ascii_flag=std::max(ascii_flag,ya_ascii_to(e_tok[0],one));
rvec.push_back(one);
}
}
if (i==0)
mat.setup(row_size,rvec.numel());
else
if (rvec.numel()!=mat.cols()) {
ostringstream owarn;
owarn << "Error parsing matrix string:\n" << value << "\nat token: "
<< row_str[i] << "\nAll rows must have the same number of "
<< "columns";
ya_addwarn(307,11,"YA_BaseT",owarn.str());
err_flag=307;
}
mat(i,":")=rvec;
}
if (transposeme)
ip_transpose(mat);
return std::max(ascii_flag,err_flag);
}
