Re: [deal.II] Exporting hdf5 file in a loop

2022-02-16 Thread Wolfgang Bangerth
On 2/16/22 10:44, Uclus Heis wrote: It seems this "solution.h5" instead of f_output was a misspelling mistake, I have right f_output in my code. I can see that the problem here is that there are processes without cells, which may be a bug of my implementation as I would like to use all the

Re: [deal.II] Exporting hdf5 file in a loop

2022-02-16 Thread Uclus Heis
Dear Wolfgang, It seems this "solution.h5" instead of f_output was a misspelling mistake, I have right f_output in my code. I can see that the problem here is that there are processes without cells, which may be a bug of my implementation as I would like to use all the processes. How does

Re: [deal.II] Exporting hdf5 file in a loop

2022-02-16 Thread Wolfgang Bangerth
Thank you for the suggestions. I updated the code as you mentioned. I do not fully understand what you mentioned about that I am not writing into the f_output stream. In the code snippet you showed, you are always writing to "solution.h5", not to the f_output stream. Is

Re: [deal.II] Exporting hdf5 file in a loop

2022-02-16 Thread Uclus Heis
Dear Wolfgang, Thank you for the suggestions. I updated the code as you mentioned. I do not fully understand what you mentioned about that I am not writing into the f_output stream. Is data_out.write_hdf5_parallel() not writing the data? However, I still have the error when running more than

Re: [deal.II] Exporting hdf5 file in a loop

2022-02-16 Thread Wolfgang Bangerth
Uclus, the usual style is to create the DataOut object every time you are writing into a file, and then let the variable die at the end of the scope again, rather than keeping it around for the next time you want to create output. You might want to take a look at how time dependent programs

Re: [deal.II] Exporting hdf5 file in a loop

2022-02-16 Thread Uclus Heis
Dear Timo, Thank you for the comments. I am doing now one file per iteration. However, when running with mpirun more tha 1 MPI ranks I get an error in the writing function. My code right now looks : DataOut data_out; data_out.attach_dof_handler(dof_handler); for(int freq_iter=0 ... )

Re: [deal.II] Exporting hdf5 file in a loop

2022-02-15 Thread Timo Heister
What we do typically in this situation is that we write one visualization output file per iteration with a different filename. This is done in many of the examples, especially for time dependent problems. I don't think there is an easy way to append data to an existing file like you suggested.

Re: [deal.II] Exporting hdf5 file in a loop

2022-02-15 Thread Uclus Heis
Thank you very much for the answer. I need to evaluate a high number of frequencies, so in that case I would need to have a large number of vectors to track the results, which is not optimal in my case. Is it any other way to do that? Would be possible to call data_out.add_data_vector() and

Re: [deal.II] Exporting hdf5 file in a loop

2022-02-15 Thread Timo Heister
The call to data_out.add_data_vector() does not copy the contents of the vector but it just keeps track of it until the data is actually written. You will need to store your solutions in different vectors without touching the old ones. On Tue, Feb 15, 2022, 11:05 Uclus Heis wrote: > Good

[deal.II] Exporting hdf5 file in a loop

2022-02-15 Thread Uclus Heis
Good afternoon, I want to store my results in a hdf5 file using a distributed implementation. I am computing different frequencies, so I have a loop in my run() function where I solve for each frequency. When for example computing 5 frequencies, I get 5 results with the same value in my hdf5