this version of code is going to be used for a test case of calculationg...

        {res[i] = 0.;}

        //number of nonzeroes in upper triangular matrix, always including diagonal

        MKL_INT nnz = ncsr.get_n_of_links() + newN ;//The adjacency matrix is symmetric but the Normallized laplacian is not

        MKL_INT nnz = ncsr.get_n_of_links()/2 + newN ;//The adjacency matrix is symmetric but the Normallized laplacian is not

        double  *a  = (double*)mkl_malloc(nnz*sizeof(double),64);

        MKL_INT *ja = (MKL_INT*)mkl_malloc(nnz*sizeof(MKL_INT),64);

        MKL_INT *ia = (MKL_INT*)mkl_malloc((newN+1)*sizeof(MKL_INT),64);