Dear all,
I would like to test my simulation on the C1 elements. But when I try to
project the initial values into the system, I got an error. Could any one tell
me how to do the initial values when using a C1 elements.
Thank you in advance,
/ Ming Q.
Here is my code:
> //_________________________________________________________________
> Number init_solution (const Point &p,
> const Parameters& par,
> const std::string& sys_name,
> const std::string& unk_name)
> {
> Real init_value=0.;
> Real x = p(0);
> Real y = p(1);
> if (unk_name.compare("val1")==0)
> {
> init_value = sin(x)*cos(y);
> }
> return init_value;
> }
> //_________________________________________________________________
> void sys1_init (EquationSystems& es, const std::string& sys_name)
> {
> libmesh_assert ( sys_name == "system1" );
>
> TransientLinearImplicitSystem&
> system=es.get_system<TransientLinearImplicitSystem> ("system1");
>
> system.project_solution(init_solution, NULL, es.parameters);
> }
Then here is the error I got when libMesh doing the projection:
> #0 0x00007fff886a7616 in __kill ()
> #1 0x00007fff88747cca in abort ()
> #2 0x00000001056e7325 in __gnu_cxx::__verbose_terminate_handler ()
> #3 0x000000010567bdba in __cxxabiv1::__terminate ()
> #4 0x000000010567be13 in std::terminate ()
> #5 0x000000010567be96 in __cxa_throw ()
> #6 0x000000010065cc9e in libMesh::System::ProjectSolution::operator()
> (this=0x7fff5fbfe1f0, range=@0x7fff5fbfe1b0) at
> src/systems/system_projection.C:896
> #7 0x0000000100664d3b in
> libMesh::Threads::parallel_for<libMesh::StoredRange<libMesh::MeshBase::const_element_iterator,
> libMesh::Elem const*>, libMesh::System::ProjectSolution>
> (range=@0x7fff5fbfe1b0, body=@0x7fff5fbfe1f0) at threads.h:270
> #8 0x000000010066348f in libMesh::System::project_vector (this=0x106376e40,
> fptr=0x100007f6b <init_solution(libMesh::Point const&, libMesh::Parameters
> const&, std::string const&, std::string const&)>, gptr=0,
> parameters=@0x7fff5fbfe6c8, new_vector=@0x10639bd00) at
> src/systems/system_projection.C:282
> #9 0x00000001006637cc in libMesh::System::project_solution
> (this=0x106376e40, fptr=0x100007f6b <init_solution(libMesh::Point const&,
> libMesh::Parameters const&, std::string const&, std::string const&)>, gptr=0,
> parameters=@0x7fff5fbfe6c8) at src/systems/system_projection.C:250
> #10 0x0000000100007f44 in sys1_init (es=@0x7fff5fbfe6c0,
> sys_name=@0x106376ed8) at source/addinc.cc:56
> #11 0x000000010063c9d1 in libMesh::System::user_initialization
> (this=0x106376e40) at src/systems/system.C:1505
> #12 0x000000010063cc52 in libMesh::System::init (this=0x106376e40) at
> src/systems/system.C:204
> #13 0x00000001005d552d in libMesh::EquationSystems::init
> (this=0x7fff5fbfe6c0) at src/systems/equation_systems.C:119
------------------------------------------------------------------------------
Special Offer-- Download ArcSight Logger for FREE (a $49 USD value)!
Finally, a world-class log management solution at an even better price-free!
Download using promo code Free_Logger_4_Dev2Dev. Offer expires
February 28th, so secure your free ArcSight Logger TODAY!
http://p.sf.net/sfu/arcsight-sfd2d
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
