testOneDConduction_double.cpp File Reference

#include <string>
#include <iostream>
#include <fstream>
#include "PC.h"
#include "OneDConduction.h"

Include dependency graph for testOneDConduction_double.cpp:

Go to the source code of this file.

Macros
#define	NCELLS   10

Typedefs
typedef PC< 3, 1 >	PCType

Functions
int	main (int argc, char *argv[])

Macro Definition Documentation

#define NCELLS   10

Definition at line 16 of file testOneDConduction_double.cpp.

Referenced by main().

Typedef Documentation

typedef PC<3,1> PCType

Definition at line 14 of file testOneDConduction_double.cpp.

Function Documentation

int main	(	int	argc,
		char *	argv[]
	)

Definition at line 17 of file testOneDConduction_double.cpp.

References PC< ORDER, DIM >::_data, Array< T >::_data, OneDConduction< T >::getSolution(), NCELLS, NTRIES, OneDConduction< T >::solve(), and stdDev().

{

  double dkConst = 1.;
  PCType pCkConst(dkConst);
  pCkConst._data[1] = 0.1;
  
  //OneDConduction<double> dp(NCELLS,dkConst);
  OneDConduction<PCType> pCp(NCELLS,pCkConst);

  //dp.solve();
  pCp.solve();

  //Array<double>& dt = dp.getSolution();
  Array<PCType>& pCt = pCp.getSolution();

  //cout << "double solution" << endl;
  //for(int i=0; i<NCELLS; i++)
  //  cout << i << " " << dt[i] << endl;
  ofstream mean_file, variance_file;
  mean_file.open("/tmp/mean-mc.dat");
  variance_file.open("/tmp/variance-mc.dat");

  double dx=1.0/NCELLS;
  double x = dx/2.0;
  for(int c=0; c<NCELLS; c++)
  {
      double Tmean = pCt[i]._data[0] << "  "
      double T = gsl_stats_variance_m(Tc[c],1,NTRIES,Tmean);
      mean_file << x << "  " << Tmean << endl;
      variance_file << xc[c] << "  " << stdDev(pCt[i]) << endl;
  }
  
  return 0;
}