* constructor to disallow the
* creation of such an object.
*/
- Number values[(dim!=0) ? (dim) : 1];
+ Number values[(dim!=0) ? (dim) : (dim+1)];
/**
* Help function for unroll. If
values[1] = p.values[1];
break;
case 3:
+
values[0] = p.values[0];
values[1] = p.values[1];
values[2] = p.values[2];
unsigned int d2,
unsigned int d3)
{
- int nodes[8];
+ int nodes[1<<dim];
nodes[GeometryInfo<dim>::dx_to_deal[0]] = start;
nodes[GeometryInfo<dim>::dx_to_deal[1]] = start+d1;
if (dim>=2)
unsigned int d2,
unsigned int d3)
{
- int nodes[8];
+ int nodes[1<<dim];
nodes[GeometryInfo<dim>::ucd_to_deal[0]] = start;
nodes[GeometryInfo<dim>::ucd_to_deal[1]] = start+d1;
if (dim>=2)
(void)comm;
AssertThrow(false, ExcMessage ("HDF5 support is disabled."));
#else
- int world_size;
hid_t h5_file_id, plist_id;
hid_t node_dataspace, node_dataset, node_file_dataspace, node_memory_dataspace;
hid_t cell_dataspace, cell_dataset, cell_file_dataspace, cell_memory_dataspace;
// If HDF5 is not parallel and we're using multiple processes, abort
#ifndef H5_HAVE_PARALLEL
#ifdef DEAL_II_COMPILER_SUPPORTS_MPI
+ int world_size;
MPI_Comm_size(comm, &world_size);
AssertThrow (world_size <= 1,
ExcMessage ("Serial HDF5 output on multiple processes is not yet supported."));