Oh, Thank you very much!!!!! The problem is solved!!!
you are right. My mesh was the problem. I was stupid!!! But, Could I have the last question....? I have to show my result to my boss(variable p). And, To show the inside of torus, It seems that I have to use contour plot using Visit.(In the Wolfgang lecture, he used contour plot to see the inside of the calculated result in 3D problem) But, When I used Contour plot in Visit for variable p, It is said that " The Contour plot of variable "p" yielded no data." Is there anyway for me to show the inside of the calculated result to my boss? Please help me.... 2016년 8월 31일 수요일 오후 11시 45분 8초 UTC+9, Jean-Paul Pelteret 님의 말: > > Hi Kyusik, > > You welcome :-) No, deal.II only supports quads/hexes > <https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#can-i-useimplement-trianglestetrahedra-in-dealii> > so > you cannot use prisms. However, since meshing an object with a circular > cross section with a structured mesh is not particularly difficult > <http://www.asergeev.com/files/qatar/openfoam/cases/courant/mesh.jpg>, I > don't think that this should be particularly challenging for you to fix. > > I must admit that I haven't looked at your application of the > SphericalManifold - you should probably deal with one problem at a time, > especially since your mesh quality issue is such a fundamental one and > might have a significant role to play here. If you fix your mesh and > continue to run into issues that you cannot resolve yourself, then perhaps > you could post again (with the updated mesh) and I'll (or someone else) > take another look at it. > > Regards, > J-P > > On Wednesday, August 31, 2016 at 4:24:13 PM UTC+2, hank...@gmail.com > <javascript:> wrote: >> >> First of all, Thank you very much for your answering my Question. >> >> Could I have one more question...? >> >> Then, Is it okay If I use Prism cell instead of the Hexahedron cell in >> UCD format?? Can the deal.II treat the Prism cell?? >> >> And, Do you also agree with that "SphericalManifold<dim> boundary" is not >> related to this problem?? >> >> Thank you. >> >> Regard >> Kyusik >> >> 2016년 8월 31일 수요일 오후 10시 36분 4초 UTC+9, Jean-Paul Pelteret 님의 말: >>> >>> 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, hank...@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.
