template <class DESTN,class SRC>
class MKLSparseRepresentation
{
	private:
		MKL_INT *rowindices;
		DESTN *values;
		MKL_INT *colptrs;
	public :
		MKLSparseRepresentation(SRC &A)
		{
			int numRows=0;
			rowindices=(MKL_INT *)malloc(sizeof(MKL_INT)*A.n_nonzero);
			values=(DESTN *)malloc(sizeof(DESTN)*A.n_nonzero);
			#pragma omp for
			for(int i=0; i<A.n_nonzero;i++)
			{
				rowindices[i]=A.row_indices[i];
				values[i]=A.values[i];
			}
			colptrs=(MKL_INT *)malloc(sizeof(MKL_INT)*A.n_cols+1);
			#pragma omp for
			for(int i=0;i<A.n_cols+1;i++)
			{
				colptrs[i]=A.col_ptrs[i];
			}

		}
		DESTN* getValues()
		{
			return values;
		}
		MKL_INT* getRowIndices()
		{
			return rowindices;
		}
		MKL_INT* getcolptrs()
		{
			return colptrs;
		}
		~MKLSparseRepresentation()
		{
			free(rowindices);
			free(colptrs);
			free(values);
		}

};