In fact, I get the same error (IncompatibleDimensions) when I try to
process step-18's *.d2 files using step-19:
bash-3.00$ run ../step-19/step-19 solution-0001.0000.d2 -x gmv test.gmv
-bash: run: command not found
-bash-3.00$ gdb ../step-19/step-19 solution-0001.0000.d2 -x gmv test.gmv
Excess command line arguments ignored. (test.gmv)
GNU gdb Red Hat Linux (6.3.0.0-1.143.el4rh)
Copyright 2004 Free Software Foundation, Inc.
GDB is free software, covered by the GNU General Public License, and you are
welcome to change it and/or distribute copies of it under certain conditions.
Type "show copying" to see the conditions.
There is absolutely no warranty for GDB. Type "show warranty" for details.
This GDB was configured as "x86_64-redhat-linux-gnu"...Using host
libthread_db library "/lib64/tls/libthread_db.so.1".
"/home/deal.ii_7.0.1/deal.II/examples/step-18/solution-0001.0000.d2"
is not a core dump: File format not recognized
gmv: No such file or directory.
(gdb) run solution-0001.0000.d2 -x gmv test.gmv
Starting program: /home/deal.ii_7.0.1/deal.II/examples/step-19/step-19
solution-0001.0000.d2 -x gmv test.gmv
[Thread debugging using libthread_db enabled]
[New Thread 183032137568 (LWP 29735)]
test.gmv
Program received signal SIGSEGV, Segmentation fault.
[Switching to Thread 183032137568 (LWP 29735)]
0x000000346870885c in _Unwind_RaiseException () from /lib64/libgcc_s.so.1
(gdb) bt
#0 0x000000346870885c in _Unwind_RaiseException () from /lib64/libgcc_s.so.1
#1 0x0000003468baf296 in __cxa_throw () from /usr/lib64/libstdc++.so.6
#2 0x0000002a9990cb5c in
dealii::deal_II_exceptions::internals::issue_error_throw<dealii::DataOutReader<3,
3>::ExcIncompatibleDimensions> (file=0x2a9aadbbc0
"/bevo2/deal.ii_7.0.1/deal.II/source/base/data_out_base.cc",
line=5704, function=0x2a9aafa100 "void dealii::DataOutReader<dim,
spacedim>::read(std::basic_istream<char, std::char_traits<char>> &)
[with dim = 3, spacedim = 3]",
cond=0x2a9aaf9e80 "(dimension_info.first == dim) &&
(dimension_info.second == spacedim)", exc_name=0x2a9aaf9e20
"ExcIncompatibleDimensions (dimension_info.first, dim,
dimension_info.second, spacedim)", e=
{<dealii::ExceptionBase> = {<exception> = {_vptr.exception =
0x2a9bcbf3a0}, file = 0x2a9aadbbc0
"/bevo2/deal.ii_7.0.1/deal.II/source/base/data_out_base.cc", line =
5704, function = 0x2a9aafa100 "void dealii::DataOutReader<dim,
spacedim>::read(std::basic_istream<char, std::char_traits<char>> &)
[with dim = 3, spacedim = 3]", cond = 0x2a9aaf9e80
"(dimension_info.first == dim) && (dimension_info.second ==
spacedim)", exc = 0x2a9aaf9e20 "ExcIncompatibleDimensions
(dimension_info.first, dim, dimension_info.second, spacedim)",
stacktrace = 0x7ba1e0, n_stacktrace_frames = 8}, arg1 = 4279959, arg2
= 3, arg3 = 0, arg4 = 3})
at /bevo2/deal.ii_7.0.1/deal.II/include/deal.II/base/exceptions.h:303
#3 0x0000002a999cee3a in dealii::DataOutReader<3, 3>::read
(this=0x7fbfffe768, in=@0x7fbfffe298) at
/bevo2/deal.ii_7.0.1/deal.II/source/base/data_out_base.cc:5704
#4 0x0000000000418cc0 in Step19::do_convert<3, 3> () at step-19.cc:511
#5 0x0000000000411301 in Step19::convert () at step-19.cc:619
#6 0x00000000004114d7 in main (argc=5, argv=0x7fbfffefc8) at step-19.cc:652
On Sun, Apr 15, 2012 at 2:35 PM, pleramorphyllactic
<[email protected]> wrote:
> Also, here's the backtrace for where step-19 is dying on the *.d2 file
> (I'm using version 7.0.1 btw), .... some sort of dimension mismatch
> for the fesystem ....
>
> #0 0x000000346870885c in _Unwind_RaiseException () from /lib64/libgcc_s.so.1
> #1 0x0000003468baf296 in __cxa_throw () from /usr/lib64/libstdc++.so.6
> #2 0x0000002a9990cab4 in
> dealii::deal_II_exceptions::internals::issue_error_throw<dealii::DataOutReader<2,
> 2>::ExcIncompatibleDimensions> (file=0x2a9aadbbc0
> "/bevo2/running/deal.ii_7.0.1/deal.II/source/base/data_out_base.cc",
> line=5704, function=0x2a9aafa060 "void dealii::DataOutReader<dim,
> spacedim>::read(std::basic_istream<char, std::char_traits<char>> &)
> [with dim = 2, spacedim = 2]",
> cond=0x2a9aaf9e80 "(dimension_info.first == dim) &&
> (dimension_info.second == spacedim)", exc_name=0x2a9aaf9e20
> "ExcIncompatibleDimensions (dimension_info.first, dim,
> dimension_info.second, spacedim)", e=
> {<dealii::ExceptionBase> = {<exception> = {_vptr.exception =
> 0x2a9bcbef60}, file = 0x2a9aadbbc0
> "/bevo2/running/deal.ii_7.0.1/deal.II/source/base/data_out_base.cc",
> line = 5704, function = 0x2a9aafa060 "void dealii::DataOutReader<dim,
> spacedim>::read(std::basic_istream<char, std::char_traits<char>> &)
> [with dim = 2, spacedim = 2]", cond = 0x2a9aaf9e80
> "(dimension_info.first == dim) && (dimension_info.second ==
> spacedim)", exc = 0x2a9aaf9e20 "ExcIncompatibleDimensions
> (dimension_info.first, dim, dimension_info.second, spacedim)",
> stacktrace = 0x5a5390, n_stacktrace_frames = 8}, arg1 = 4279959, arg2
> = 2, arg3 = 0, arg4 = 2})
> at
> /bevo2/running/deal.ii_7.0.1/deal.II/include/deal.II/base/exceptions.h:303
> #3 0x0000002a999cdec2 in dealii::DataOutReader<2, 2>::read
> (this=0x7fbfffe5d8, in=@0x7fbfffe108) at
> /bevo2/running/deal.ii_7.0.1/deal.II/source/base/data_out_base.cc:5704
> #4 0x0000000000417fd7 in Step19::do_convert<2, 2> () at step-19.cc:509
> #5 0x00000000004111fd in Step19::convert () at step-19.cc:604
> #6 0x00000000004114d7 in main (argc=5, argv=0x7fbfffee28) at step-19.cc:650
>
>
> On Sat, Apr 14, 2012 at 8:42 PM, pleramorphyllactic
> <[email protected]> wrote:
>> Hi guys,
>>
>> Quick question, I hope someone can help with. I'm using fesystem with
>> two components: one is a scalar component, the other a vector
>> component. In serial I'm having no trouble outputing data doing the
>> following type of thing:
>>
>>
>> DataOut<dim,DoFHandler<dim> > data_out;
>> data_out0.attach_dof_handler (*(dof_handler[k]));
>> data_out0.add_data_vector (subdomain_solution[k],
>> hd<dim>::component_names (k),
>> DataOut<dim, DoFHandler<dim> >::type_dof_data,
>> hd<dim>::component_interpretation ());
>>
>>
>> std::ostringstream filename;
>> filename << "solution0-"
>> << cycle
>> << ".vtk";
>>
>> std::ofstream output (filename.str().c_str());
>> data_out.write_vtk (output);
>>
>> where the component interpretation works fine, and here's the pseudo flow:
>>
>> if(k>=0){
>> std::vector<std::string> names;
>> std::ostringstream parser;
>> parser << "alpha_" << (k);
>> names.push_back(parser.str());
>> for(unsigned int j=0;j<dim;j++){
>> std::ostringstream parser;
>> parser << "sigma_" << (k) << "_<" << (j) << ">";
>> names.push_back(parser.str());
>> }
>>
>> std::vector<DataComponentInterpretation::DataComponentInterpretation>
>> data_component_interpretation;
>>
>> data_component_interpretation.push_back(DataComponentInterpretation::component_is_scalar);
>> for(unsigned int j=0;j<dim;j++){
>>
>> data_component_interpretation.push_back(DataComponentInterpretation::component_is_part_of_vector);
>> }
>> return data_component_interpretation;
>>
>> The problem is, when I'm running over many processors using petsc and
>> mpi, this output gives processor specific data, and isn't really
>> useful. I've tried the style used in step-40 and step-18, but in the
>> former case I get an error on the vtu/pvtu files where paraview can't
>> load them (, and in the ladder step-19 kicks an exception and can't
>> convert. I'm guessing it might have something to do with the fact
>> that the fesystem has one scalar component and one vector component?
>> My preference would be to have something like used in step-18. Here's
>> what I've tried:
>>
>> for(unsigned int k=0; k<alphadim; k++){
>>
>>
>> if(k==0){
>>
>>
>> FilteredDataOut<dim> data_out0(this_mpi_process);
>>
>> data_out0.attach_dof_handler (*(dof_handler[k]));
>> data_out0.add_data_vector (subdomain_solution[k],
>> rmhd<dim>::component_names (k),
>> DataOut<dim, DoFHandler<dim>
>> >::type_dof_data,
>> rmhd<dim>::component_interpretation ());
>>
>> std::vector<unsigned int> partition_int
>> (triangulation.n_active_cells());
>> GridTools::get_subdomain_association (triangulation, partition_int);
>> const Vector<double> partitioning(partition_int.begin(),
>> partition_int.end());
>> data_out0.add_data_vector (partitioning, "partitioning");
>>
>> data_out0.build_patches ();
>>
>> std::ostringstream filename;
>> filename << "solution-";
>> filename << std::setfill('0');
>> filename.setf(std::ios::fixed, std::ios::floatfield);
>> filename << std::setw(9) << std::setprecision(4) << cycle;
>>
>> if (n_mpi_processes != 1)
>> {
>> AssertThrow (n_mpi_processes < 1000, ExcNotImplemented());
>>
>> filename << '-';
>> filename << std::setfill('0');
>> filename << std::setw(3) << this_mpi_process;
>> }
>>
>> filename << data_out0.default_suffix(DataOut<dim>::deal_II_intermediate);
>>
>> std::ofstream output (filename.str().c_str());
>> data_out0.write_deal_II_intermediate (output);
>> }
>>
>> The *.d2 files print out, and look reasonable to me when I open them,
>> but ../../../step-19/step-19 solution-000000000.d1 -x gmv test1.gmv
>> aborts without converting.
>>
>> An error occurred in line <580> of file <step-19.cc> in function
>> void Step19::convert()
>> The violated condition was:
>> input
>> The name and call sequence of the exception was:
>> ExcIO()
>> Additional Information:
>> (none)
>>
>> I'm not sure at present how to fix this.
>>
>> Many thanks,
>> Evan
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
