Dear Kyusik, If you open your original mesh file in Paraview (version 5.1 detects it as an AVS UCD file) and you clip the mesh to view its section, then I think you'll see the problem. Although the mesh looks fine from the outside, all of the inner vertices of your torus converge on an annulus (so the hexahedra degenerate into prisms). This is not allowed. You'll need to generate a new mesh with higher quality elements before proceeding.
Regards, J-P On Wednesday, August 31, 2016 at 2:59:02 PM UTC+2, hanks0...@gmail.com wrote: > > Hi > > I'm studying step-20 tutorial example, watching Wolfgang's Lectures > > In step20, I would like to solve that Eq in the 3D mesh that I made by > myself(I attached B_coord3D.inp that is UCD format) > > But, It didn't work out. > > The error message is like this.... > > [ 50%] Built target step-20 > make[3]: Warning: File `step-20' has modification time 81 s in the future > [100%] Run step-20 with Debug configuration > Number of active cells: 36864 > Total number of cells: 42048 > Number of degrees of freedom: 156672 (119808+36864) > 53 CG Schur complement iterations to obtain convergence. > > -------------------------------------------------------- > An error occurred in line <790> of file > </user2/hanks318/dealii/dealii-8.3.0/source/fe/mapping_q1.cc> in function > void dealii::MappingQ1<dim, spacedim>::fill_fe_values(const typename > dealii::Triangulation<dim, spacedim>::cell_iterator&, const > dealii::Quadrature<dim>&, typename dealii::Mapping<dim, > spacedim>::InternalDataBase&, std::vector<dealii::Point<spacedim, double>, > std::allocator<dealii::Point<spacedim, double> > >&, std::vector<double, > std::allocator<double> >&, std::vector<dealii::DerivativeForm<1, spacedim, > dim, double>, std::allocator<dealii::DerivativeForm<1, spacedim, dim, > double> > >&, std::vector<dealii::DerivativeForm<2, dim, spacedim, double>, > std::allocator<dealii::DerivativeForm<2, dim, spacedim, double> > >&, > std::vector<dealii::DerivativeForm<1, spacedim, dim, double>, > std::allocator<dealii::DerivativeForm<1, spacedim, dim, double> > >&, > std::vector<dealii::Point<spacedim, double>, > std::allocator<dealii::Point<spacedim, double> > >&, > dealii::CellSimilarity::Similarity&) const [with int dim = 3, int spacedim > = 3] > The violated condition was: > det > 1e-12*Utilities::fixed_power<dim>(cell->diameter()/ > std::sqrt(double(dim))) > The name and call sequence of the exception was: > (typename > Mapping<dim,spacedim>::ExcDistortedMappedCell(cell->center(), det, point)) > Additional Information: > The image of the mapping applied to cell with center [-1.61075 -0.761538 > 0.649647] is distorted. The cell geometry or the mapping are invalid, > giving a non-positive volume fraction of -1.49556e-07 in quadrature point 5. > > Stacktrace: > ----------- > #0 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > dealii::MappingQ1<3, > 3>::fill_fe_values(dealii::TriaIterator<dealii::CellAccessor<3, 3> > > const&, dealii::Quadrature<3> const&, dealii::Mapping<3, > 3>::InternalDataBase&, std::vector<dealii::Point<3, double>, > std::allocator<dealii::Point<3, double> > >&, std::vector<double, > std::allocator<double> >&, std::vector<dealii::DerivativeForm<1, 3, 3, > double>, std::allocator<dealii::DerivativeForm<1, 3, 3, double> > >&, > std::vector<dealii::DerivativeForm<2, 3, 3, double>, > std::allocator<dealii::DerivativeForm<2, 3, 3, double> > >&, > std::vector<dealii::DerivativeForm<1, 3, 3, double>, > std::allocator<dealii::DerivativeForm<1, 3, 3, double> > >&, > std::vector<dealii::Point<3, double>, std::allocator<dealii::Point<3, > double> > >&, dealii::CellSimilarity::Similarity&) const > #1 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > dealii::FEValues<3, 3>::do_reinit() > #2 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: void > dealii::FEValues<3, 3>::reinit<dealii::DoFHandler<3, 3>, > false>(dealii::TriaIterator<dealii::DoFCellAccessor<dealii::DoFHandler<3, > 3>, false> > const&) > #3 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: void > dealii::hp::FEValues<3, 3>::reinit<dealii::DoFHandler<3, 3>, > false>(dealii::TriaIterator<dealii::DoFCellAccessor<dealii::DoFHandler<3, > 3>, false> >, unsigned int, unsigned int, unsigned int) > #4 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: void > dealii::internal::DataOut::ParallelDataBase<3, > 3>::reinit_all_fe_values<dealii::DoFHandler<3, 3> > >(std::vector<boost::shared_ptr<dealii::internal::DataOut::DataEntryBase<dealii::DoFHandler<3, > > > 3> > >, > std::allocator<boost::shared_ptr<dealii::internal::DataOut::DataEntryBase<dealii::DoFHandler<3, > > 3> > > > >&, dealii::TriaIterator<dealii::CellAccessor<3, 3> > const&, > unsigned int) > #5 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > dealii::DataOut<3, dealii::DoFHandler<3, 3> > >::build_one_patch(std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > 3> >, unsigned int> const*, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&, dealii::DataOut<3, dealii::DoFHandler<3, > 3> >::CurvedCellRegion, std::vector<dealii::DataOutBase::Patch<3, 3>, > std::allocator<dealii::DataOutBase::Patch<3, 3> > >&) > #6 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > boost::_mfi::mf5<void, dealii::DataOut<3, dealii::DoFHandler<3, 3> >, > std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, unsigned int> > const*, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&, dealii::DataOut<3, dealii::DoFHandler<3, > 3> >::CurvedCellRegion, std::vector<dealii::DataOutBase::Patch<3, 3>, > std::allocator<dealii::DataOutBase::Patch<3, 3> > > >&>::operator()(dealii::DataOut<3, dealii::DoFHandler<3, 3> >*, > std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, unsigned int> > const*, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&, dealii::DataOut<3, dealii::DoFHandler<3, > 3> >::CurvedCellRegion, std::vector<dealii::DataOutBase::Patch<3, 3>, > std::allocator<dealii::DataOutBase::Patch<3, 3> > >&) const > #7 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: void > boost::_bi::list6<boost::_bi::value<dealii::DataOut<3, > dealii::DoFHandler<3, 3> >*>, boost::arg<1>, boost::arg<2>, boost::arg<3>, > boost::_bi::value<dealii::DataOut<3, dealii::DoFHandler<3, 3> > >::CurvedCellRegion>, > boost::reference_wrapper<std::vector<dealii::DataOutBase::Patch<3, 3>, > std::allocator<dealii::DataOutBase::Patch<3, 3> > > > > >::operator()<boost::_mfi::mf5<void, dealii::DataOut<3, > dealii::DoFHandler<3, 3> >, > std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, unsigned int> > const*, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&, dealii::DataOut<3, dealii::DoFHandler<3, > 3> >::CurvedCellRegion, std::vector<dealii::DataOutBase::Patch<3, 3>, > std::allocator<dealii::DataOutBase::Patch<3, 3> > >&>, > boost::_bi::list3<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> > >, unsigned int>* const&, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&> >(boost::_bi::type<void>, > boost::_mfi::mf5<void, dealii::DataOut<3, dealii::DoFHandler<3, 3> >, > std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, unsigned int> > const*, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&, dealii::DataOut<3, dealii::DoFHandler<3, > 3> >::CurvedCellRegion, std::vector<dealii::DataOutBase::Patch<3, 3>, > std::allocator<dealii::DataOutBase::Patch<3, 3> > >&>&, > boost::_bi::list3<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> > >, unsigned int>* const&, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&>&, int) > #8 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: void > boost::_bi::bind_t<void, boost::_mfi::mf5<void, dealii::DataOut<3, > dealii::DoFHandler<3, 3> >, > std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, unsigned int> > const*, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&, dealii::DataOut<3, dealii::DoFHandler<3, > 3> >::CurvedCellRegion, std::vector<dealii::DataOutBase::Patch<3, 3>, > std::allocator<dealii::DataOutBase::Patch<3, 3> > >&>, > boost::_bi::list6<boost::_bi::value<dealii::DataOut<3, > dealii::DoFHandler<3, 3> >*>, boost::arg<1>, boost::arg<2>, boost::arg<3>, > boost::_bi::value<dealii::DataOut<3, dealii::DoFHandler<3, 3> > >::CurvedCellRegion>, > boost::reference_wrapper<std::vector<dealii::DataOutBase::Patch<3, 3>, > std::allocator<dealii::DataOutBase::Patch<3, 3> > > > > > >::operator()<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>* const, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> > >(std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, unsigned > int>* const&, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&) > #9 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > boost::detail::function::void_function_obj_invoker3<boost::_bi::bind_t<void, > boost::_mfi::mf5<void, dealii::DataOut<3, dealii::DoFHandler<3, 3> >, > std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, unsigned int> > const*, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&, dealii::DataOut<3, dealii::DoFHandler<3, > 3> >::CurvedCellRegion, std::vector<dealii::DataOutBase::Patch<3, 3>, > std::allocator<dealii::DataOutBase::Patch<3, 3> > >&>, > boost::_bi::list6<boost::_bi::value<dealii::DataOut<3, > dealii::DoFHandler<3, 3> >*>, boost::arg<1>, boost::arg<2>, boost::arg<3>, > boost::_bi::value<dealii::DataOut<3, dealii::DoFHandler<3, 3> > >::CurvedCellRegion>, > boost::reference_wrapper<std::vector<dealii::DataOutBase::Patch<3, 3>, > std::allocator<dealii::DataOutBase::Patch<3, 3> > > > > >, void, > std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, unsigned int>* > const&, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, > 3>&>::invoke(boost::detail::function::function_buffer&, > std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, unsigned int>* > const&, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&) > #10 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > boost::function3<void, > std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, unsigned int>* > const&, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, > 3>&>::operator()(std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> > >, unsigned int>* const&, dealii::internal::DataOut::ParallelData<3, 3>&, > dealii::DataOutBase::Patch<3, 3>&) const > #11 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> > >::operator()(tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > const&) > #12 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > boost::_mfi::mf1<void, > dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> >, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > > const&>::operator()(dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> >&, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > const&) const > #13 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: void > boost::_bi::list2<boost::reference_wrapper<dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> > >, boost::arg<1> > >::operator()<boost::_mfi::mf1<void, > dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> >, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > const&>, > boost::_bi::list1<tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > >&> >(boost::_bi::type<void>, boost::_mfi::mf1<void, > dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> >, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > const&> const&, > boost::_bi::list1<tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > >&>&, int) const > #14 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: void > boost::_bi::bind_t<void, boost::_mfi::mf1<void, > dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> >, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > const&>, > boost::_bi::list2<boost::reference_wrapper<dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> > >, boost::arg<1> > > >::operator()<tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > > >(tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > >&) const > #15 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > tbb::interface6::internal::start_for<tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > >, boost::_bi::bind_t<void, boost::_mfi::mf1<void, > dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> >, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > const&>, > boost::_bi::list2<boost::reference_wrapper<dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> > >, boost::arg<1> > >, > tbb::auto_partitioner > const>::run_body(tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > >&) > #16 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: void > tbb::interface6::internal::partition_type_base<tbb::interface6::internal::auto_partition_type>::execute<tbb::interface6::internal::start_for<tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > >, boost::_bi::bind_t<void, boost::_mfi::mf1<void, > dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> >, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > const&>, > boost::_bi::list2<boost::reference_wrapper<dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> > >, boost::arg<1> > >, > tbb::auto_partitioner const>, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > > >(tbb::interface6::internal::start_for<tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > >, boost::_bi::bind_t<void, boost::_mfi::mf1<void, > dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> >, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > const&>, > boost::_bi::list2<boost::reference_wrapper<dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> > >, boost::arg<1> > >, > tbb::auto_partitioner const>&, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > >&) > #17 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > tbb::interface6::internal::start_for<tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > >, boost::_bi::bind_t<void, boost::_mfi::mf1<void, > dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> >, > tbb::blocked_range<__gnu_cxx::__normal_iterator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>* const*, > std::vector<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*, > std::allocator<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, 3> >, > unsigned int>*> > > > const&>, > boost::_bi::list2<boost::reference_wrapper<dealii::WorkStream::internal::Implementation3::WorkerAndCopier<std::pair<dealii::TriaIterator<dealii::CellAccessor<3, > > 3> >, unsigned int>*, dealii::internal::DataOut::ParallelData<3, 3>, > dealii::DataOutBase::Patch<3, 3> > >, boost::arg<1> > >, > tbb::auto_partitioner const>::execute() > #18 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > tbb::internal::custom_scheduler<tbb::internal::IntelSchedulerTraits>::local_wait_for_all(tbb::task&, > > tbb::task*) > #19 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > tbb::internal::arena::process(tbb::internal::generic_scheduler&) > #20 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > tbb::internal::market::process(rml::job&) > #21 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > tbb::internal::rml::private_worker::run() > #22 /user2/hanks318/dealii/dealii_pre/lib/libdeal_II.g.so.8.3.0: > tbb::internal::rml::private_worker::thread_routine(void*) > -------------------------------------------------------- > > > I can't understand the meaning of distorted image of mapping, non positive > volume etc... > > The below picture is the 3D mesh I used for the calculation(it is refined > twice by triangulation.refine_global(2);) > > > <https://lh3.googleusercontent.com/-mg4mtHj0cPo/V8bL5xBmV3I/AAAAAAAAAAc/AuOBTpAZ5aAnbd3fj4DQWylbPIOh5L2jgCLcB/s1600/mesh_3D.png> > > From this picture, It seems that The Mesh is fine. I don't know what the > problem is. > > And FYI, I used "SphericalManifold<dim> boundary" in this code, even > though the boundary is not sphere. I used this to make the boundary > smoother after every refinement. > > But, When I didn't use this, The problem is not solved. So, I thought it > is not related to this problem > > I also attach the step-20.cc that I revised a little bit.(I commented what > I change in step-20.cc) > > please help me to solve this problem.... > > Thank you > > Regards, > Kyusik > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
