ok <div>actually reading the documentation is not such a bad idea.....</div><div>building my own postprocessor class for the Simple model now and using the statistics object.</div><div><br></div><div>@time i solved my issue with const GeometryModel::SphericalShell<dim> *geometry = NULL;</div>
<div><div><br></div><div>cheers</div></div><div>Thomas</div><div><br><br><div class="gmail_quote">On Sat, Jan 7, 2012 at 1:28 AM, Wolfgang Bangerth <span dir="ltr"><<a href="mailto:bangerth@math.tamu.edu">bangerth@math.tamu.edu</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
Hi Thomas,<br>
I've taken a look at the changes you made and I think you're making your life much harder than necessary. If I see this right, for every piece of data you want to compute, you open (or keep open) a file and litter the output directory with a bunch of different files. The idea of these postprocessors is different, however: simply put the data you generate into the statistics object every postprocessor function gets as argument, and it will show up in the statistics file in the output directory chosen in the input file. The data you produce will then simply be one column in this file.<br>
<br>
For example, rather than writing code like<br>
if (this->Vrmsout== NULL) this->Vrmsout = fopen("bin/Vrms.dat", "w");<br>
fprintf(Vrmsout,"%e %e\n", this->get_time(), vrms);<br>
and having the Vrmsout handle as a member variable, simply do<br>
statistics.add_value("vrms", vrms);<br>
<br>
If you run this, take a look at bin/statistics (or whichever directory this is for you) and find a column labeled "vrms" that you can then plot just as easily as the data above.<br>
<br>
The reason to do it this way is (i) it's simpler, no need to deal with output files by hand; (ii) it's designed to also work in the case where you checkpoint the program and restart it later, which in your scheme will lead to a loss of the data you've generated so far.<br>
<br>
I think it would be useful to transition the code you wrote in revisions 569 and following to the scheme above. Let me know if you have questions about this and we'll try to help out. I guess I should also write a bit of text in section 7.1.6 of the manual on how this is all supposed to work...<br>
<br>
Cheers<br>
W.<br>
<br>
------------------------------<u></u>------------------------------<u></u>------------<span class="HOEnZb"><font color="#888888"><br>
Wolfgang Bangerth email: <a href="mailto:bangerth@math.tamu.edu" target="_blank">bangerth@math.tamu.edu</a><br>
www: <a href="http://www.math.tamu.edu/~bangerth/" target="_blank">http://www.math.tamu.edu/~<u></u>bangerth/</a><br>
<br>
</font></span></blockquote></div><br></div>