On Aug 6, 2008, at 9:52 PM, Jefferson Thomas wrote: > The problem is that the function fields::get_field works to slow when > used for copying > the field from the entire volume into a new array, which would be > needed for VTK output. > > So is there maybe a way to copy the entire field array into a new > double*** array > or something like that ?
You don't want to do it using the get_field function. You would want to do it using the loop_in_chunks function, just as for output_hdf5. One problem with the VTK file format is that I'm not aware of any parallel I/O libraries for this format (although I could be missing something), so you would have to implement this yourself. You say that calling output_hdf5, then h5tovtk, then deleting the .h5 file, is too slow. Do you have any benchmarks to illustrate this? I would have thought that the time to run h5tovtk would be fine compared to the time for output_hdf5. I disagree that writing an HDF5 file and then converting it is "definitely not elegant". The alternative is for every program to implement support for every major file format, which is definitely not elegant. Writing to a single standard interchange file format and then using converter programs for other formats and programs follows the the spirit of the Unix philosophy. Steven _______________________________________________ meep-discuss mailing list [email protected] http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/meep-discuss
