#include "alglibinternal.h"
#include "linalg.h"
#include "ap.h"
#include <fstream>
#include <iostream>
#include <string>
#include <stdlib.h>
#include <stdio.h>
#include "stdafx.h"

using namespace alglib;
using namespace std;


int main(){

   ifstream archivo;
   archivo.open("salida.txt");
   int n,m,NM,N,M;
   string valor,ns,ms,NMS;
   getline(archivo,valor,' ');
   N=atoi(valor.c_str());
   getline(archivo,valor,' ');
   M=atoi(valor.c_str());
   getline(archivo,NMS);
   NM=atoi(NMS.c_str());
   
   double matrix[N][M];
        double *array;
        real_1d_array b;
        array=new double(N*M);

   for(int i=0;i<N;i++){
      for(int j=0;j<M;j++){
         matrix[i][j]=0;
      }
   }
   
   while(!archivo.eof()){
      getline(archivo,ns,' ');
      getline(archivo,ms,' ');
      n=atoi(ns.c_str());
      m=atoi(ms.c_str());
      getline(archivo,valor,'\n');
      if(archivo.fail())
         break;
      matrix[n-1][m-1]=atof(valor.c_str());
      
   }
   archivo.close();
   
   for(int i=0,l=0;i<N;i++,l++){
      for(int j=0;j<M;j++){
         array[l]=matrix[i][j];
         l++;
      }
   }

   real_1d_array X;
        real_2d_array A; //<<<<<------ [color=#FF0000]PROBLEM[/color]
   ae_int_t max_fila(N);
   ae_int_t max_columna(M);
   ae_int_t info;
   double r1,rinf;
   //densesolverreport& rep; [color=#0000FF]spdmatrixcholesky need this data type but it doesn't work[/color]
   A.setlength(N,N);
   A.setcontent(max_fila,max_columna,array);
   X.setlength(N);
   for(int i=0; i<N;i++){
      array[i]=0;   
   }
   array[7]=1;
   b.setcontent(N,array);
   
   for(int i=1,l=0;i<10;i++,l++){
      for(int j=1;j<10;j++){
         cout<< A(i,j)<<"\t";
      }
   }
   //spdmatrixcholesky(A, max_fila, false,b,info,rep,X);[color=#0000FF] Is this the right way to use this function[/color]?
        delete [] array;
   return 0;


}

g++ -g ap.cpp  bcsstk01.cpp