string tmpstr = JFile::readFile("..\\testLinAlg\\vb.csv");
   alglib::real_2d_array mac;
   int n,m;
   double *b = JStr::toNumArr<double>(tmpstr,n,newline);
   double **a = new double*[n];
   vector<string> al = JFile::readFileLines("..\\testLinAlg\\ma.csv",true);
   if (al.size()!=n) throw JException("Files do not match.");
   for (int r=0;r<n;r++) {
      al[r] = al[r].substr(0,al[r].length()-2);
      a[r] = JStr::toNumArr<double>(al[r],m,",");
   }

   //4 check
   double *d = new double[n*m];
   memset(d,0,n*m*sizeof(double));
   mac.setcontent(n,m,d);

   //pop data
   sparsematrix ma;
    sparsecreate(n, m, ma);
   for (int r=0;r<n;r++){
      for (int c=0;c<m;c++){
         if (a[r][c]!=0) {
            mac[r][c] = a[r][c];
            sparseset(ma, r, c, a[r][c]);
         }
      }
   }
   sparseconverttocrs(ma);
   real_1d_array rhs;
   //rhs.setlength(n);
   //for (int r=0;r<n;r++)rhs[r] = b[r];
   rhs.setcontent(n,b);

   //check with full solver
   alglib::real_1d_array check;
   alglib::ae_int_t      info;
   alglib::densesolverlsreport repc;
   alglib::rmatrixsolvels(mac,n,m,rhs,0.0,info,repc,check);
   cout << "velAdj:" << endl; showMatrix(check,10); //last value: 0.000263801


   //solve
    linlsqrstate s;
    linlsqrreport rep;
    real_1d_array x;
    linlsqrcreate(n, m, s);
   linlsqrsetcond(s, 1.0e-10,1.0e-10,n+m+50);
    linlsqrsolvesparse(s, ma, rhs);
    linlsqrresults(s, x, rep);

   printf("%d\n", int(rep.terminationtype)); // EXPECTED: 4
   //printf("%s\n", x.tostring(6).c_str());
   double *res = x.getcontent();
    cout << JArray::join<double>(res,min(m,10),newline) << endl; //last value: 0.000123301