IO.h

// assemble matrix and vector
#ifndef IO_H
#define IO_H

#include "utils.h"
#include "advection.h"

using namespace std;


// write out the sparse matrix
int write_matrix(cuint P,
				 cuint Q,
				 double** U,
				 char* file_name);

// write out the sparse matrix
int write_vector( cuint P,
				  double* F,
				  char* file_name);

// write out the results
int write_results(  double* U,
					double* V,
					double* W,
					double* P,
					cuint n_dof,
					cuint nx,
					cuint ny,
					cuint nz,
					cdouble xmin,
					cdouble ymin,
					cdouble zmin,
					cdouble hx,
					cdouble hy,
					cdouble hz,
					cuint ts,
					cdouble bcs[][6]
					);

// write out 3d data for debuggin purpose
int write_3d_data( double* U,
				   cuint nx, cuint ny, cuint nz,
				   char* file_name );
#endif //IO_H