00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 #include "SundanceDSVWriter.hpp"
00032 #include "PlayaExceptions.hpp"
00033 #include "SundanceOut.hpp"
00034 #include "PlayaTabs.hpp"
00035 #include "Teuchos_XMLObject.hpp"
00036 #include <iomanip>
00037
00038
00039
00040 namespace Sundance
00041 {
00042
00043 using namespace Teuchos;
00044
00045 void DSVWriter::write() const
00046 {
00047 using std::setw;
00048 using std::setprecision;
00049 int width = precision_+7;
00050
00051 std::ofstream os(filename().c_str());
00052 int numNodes = mesh().numCells(0);
00053
00054 std::ios_base::fmtflags oldFlags = os.flags();
00055
00056 for (int i=0; i<numNodes; i++)
00057 {
00058 const Point& x = mesh().nodePosition(i);
00059 for (int d=0; d<x.dim(); d++)
00060 {
00061 if (d != 0) os << delim_;
00062 os << setw(width) << setprecision(precision_) << x[d];
00063 }
00064
00065
00066 for (int j=0; j<pointScalarFields().size(); j++)
00067 {
00068 const RCP<FieldBase>& expr = pointScalarFields()[j];
00069 os << setw(width) << setprecision(precision_);
00070 os << expr->getData(0, i, 0);
00071 }
00072 os << std::endl;
00073 }
00074
00075 os.flags(oldFlags);
00076 }
00077
00078
00079 }