write_series.cc
f31c5570ab70master

write_series.cc
View Options

			#include "write_series.hh"

			WriteSeries::WriteSeries(Series& series, unsigned int N, unsigned int precision): DumperSeries(series)
			{
			maxiter = N;
			}

			void WriteSeries::setSeparator(std::string separator){
			if (separator == "comma"){
			sep = ", ";
			format = ".csv";
			}
			else if (separator == "pipe"){
			sep = "\|";
			format = ".psv";
			}

			}

			void WriteSeries::dump(std::ostream & os){

			std::string file_name = "series_output";
			file_name.append(format);

			std::ofstream my_file;
			my_file.open(file_name, std::ios::out \| std::ios::trunc);

			double computed = 0.0;
			double analytical = 0.0;
			double error = 0.0;

			for (unsigned long i = 0.0; i <= maxiter; ++i){

			if (std::isnan(series.getAnalyticPrediction())){
			computed = series.compute(i);
			my_file << std::setprecision(dump_precision) << i << "\t" << computed << std::endl;
			}else {
			computed = series.compute(i);
			analytical = series.getAnalyticPrediction();
			error = 100.0*(analytical - computed)/analytical;
			my_file << std::setprecision(dump_precision) << i << "\t" << computed << "\t" << analytical << "\t" << error <<std::endl;
			}

			}

			my_file.close();

			}