From: Denis D. <dan...@ya...> - 2003-05-05 12:00:01
|
Hi, Does anybody use the libmesh for simulations of time-dependent problems with adaptive mesh refinement? For example a heat transfer equation. When I have refined my mesh I have to transfer the solution vector from the old mesh to the new refined one. Has the library such functionality? Best regards, Denis |
From: Benjamin S. K. <be...@cf...> - 2003-05-06 13:47:17
|
Yes... I just added a project_vector() member last week that projects a vector defined on one mesh to another (refined) mesh. Right now it only works with Lagrange & discintinuous elements since AMR is only functional with these element types. I'm reworking the EquationSystems & System classes right now. They got a little confusing for my taste. I expect to be done with this by the end of the week. When it is done I'll write a new example showing how to use AMR and project between meshes. So, how do you like C++? -Ben Denis Danilov wrote: > Hi, > > Does anybody use the libmesh for simulations of time-dependent problems > with adaptive mesh refinement? For example a heat transfer equation. > When I have refined my mesh I have to transfer the solution vector from > the old mesh to the new refined one. Has the library such functionality? > > Best regards, > Denis > > > ------------------------------------------------------- > This sf.net email is sponsored by:ThinkGeek > Welcome to geek heaven. > http://thinkgeek.com/sf > _______________________________________________ > Libmesh-devel mailing list > Lib...@li... > https://lists.sourceforge.net/lists/listinfo/libmesh-devel |
From: Denis D. <dan...@ya...> - 2003-05-07 14:39:12
Attachments:
probe_1.C
|
Thank you for answer. On Tue, May 06, 2003 at 08:13:45AM -0500, Benjamin S. Kirk wrote: > Yes... I just added a project_vector() member last week that projects a > vector defined on one mesh to another (refined) mesh. Right now it only > works with Lagrange & discintinuous elements since AMR is only > functional with these element types. I am glad to hear that library can this. > > I'm reworking the EquationSystems & System classes right now. They got > a little confusing for my taste. I expect to be done with this by the > end of the week. When it is done I'll write a new example showing how > to use AMR and project between meshes. An example is sometimes best adviser. An example is sometimes best adviser. What about an example showing how to use GeneralSystem::attach_init_function? > > So, how do you like C++? Hmm... hard to say... What can you say about such small testing program and such results? What I need to put in destructor else? refine level: 6 Mesh Information: mesh_dimension()=2 spatial_dimension()=3 n_nodes()=4225 n_elem()=5461 n_active_elem()=4096 n_subdomains()=1 n_processors()=1 processor_id()=0 EquationSystems<General>: n_systems()=1 System "phase-field system" Variables="phi" Finite Element Types="0" Approximation Orders="1" n_dofs()=4225 n_local_dofs()=4225 n_constrained_dofs()=0 Flags: "General" Parameters: "linear solver maximum iterations"=5000 "linear solver tolerance"=1e-12 ---------------------------------------------------------------------------- | Reference count information | ---------------------------------------------------------------------------- | 10SystemBase class reference count information: | Creations: 1 | Destructions: 1 | 4Elem class reference count information: | Creations: 70395 | Destructions: 70395 | 4Node class reference count information: | Creations: 21844 | Destructions: 21844 | t12SparseMatrix1Zd class reference count information: | Creations: 1 | Destructions: 1 | t13NumericVector1Zd class reference count information: | Creations: 5 | Destructions: 5 | t21LinearSolverInterface1Zd class reference count information: | Creations: 1 | Destructions: 1 ---------------------------------------------------------------------------- ERROR: libMesh::close() called prematurely. Library objects still in scope! [0] /home/denis/sci/libmesh-0.3.3-gcc2.95/include/libmesh.h, line 154, compiled on May 7 2003 at 14:25:19 zsh: abort ./probe_1 Best regards, Denis |