Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2014-06-23 Thread Roy Stogner
On Mon, 23 Jun 2014, Onur Solmaz wrote: > First, I wanted to ask whether there has been any development on giving > Dirichlet boundary conditions node-wise, on irregular meshes. Yes. If you assign a DirichletBoundary to a particular boundary id, and you give a node that boundary id, then that d

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2014-06-23 Thread Onur Solmaz
It has been two years since I last replied to this thread, and again I have the same problem. First, I wanted to ask whether there has been any development on giving Dirichlet boundary conditions node-wise, on irregular meshes. I have been doing it using MeshData which will be deprecated (actually

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-08 Thread Onur Solmaz
It seems that the problem was due to unused nodes which didn't belong to any element. mesh.prepare_for_use(false); Solves the problem in some way, since it jettisons those nodes. I'm going to need to solve the problem of matching those soln values back to my original nodes. Thanks Onur Onur 201

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-07 Thread Onur Solmaz
Deprecating MeshData would require the methods of UNVIO etc. to be rewritten. I'm glad I'm not the one to figure out the classes and methods going to be involved when transferring from a file to an explicitsystem. As for me, I am using my own function to transfer from my data structures to the mes

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-07 Thread Roy Stogner
On Tue, 7 Aug 2012, Onur Solmaz wrote: > I'm using Tet4 elements. > How do I compile libmesh in debug mode? I couldn't figure out. Using a Bourne-type shell, "export METHOD=dbg" before you compile libMesh or run any examples; or just replace "make" with "METHOD=dbg make". > I'm using the most r

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-07 Thread Onur Solmaz
I'm using Tet4 elements. How do I compile libmesh in debug mode? I couldn't figure out. I'm using the most recent svn build by the way. Onur 2012/8/7 Roy Stogner > > On Tue, 7 Aug 2012, Onur Solmaz wrote: > > It would be great if there were an example on that ExplicitSystem. >> > > *Definitely

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-07 Thread Roy Stogner
On Tue, 7 Aug 2012, Onur Solmaz wrote: > It would be great if there were an example on that ExplicitSystem. *Definitely*. Trouble is finding someone with both the knowhow and the time to sit back and write it up. > I also had another weird problem. Would there be a reason that the > number of

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-07 Thread Onur Solmaz
It would be great if there were an example on that ExplicitSystem. I also had another weird problem. Would there be a reason that the number of dofs won't be equal to the number of nodes? It is weird because my program works perfectly fine with my first input file (n_dofs=n_nodes), but it crashed

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-07 Thread Roy Stogner
On Tue, 7 Aug 2012, Onur Solmaz wrote: Thanks for the detailed response! No trouble. Try to hit "Reply to All" rather than just "Reply" to libmesh messages, though - I like making sure these threads get archived so Google can help others with similar issues in the future. Calling Lib

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-07 Thread Roy Stogner
On Tue, 7 Aug 2012, Onur Solmaz wrote: Calling LibmeshInit in my main function gives me segfaults for obvious reasons. I'm afraid not obvious to me - all our examples and nearly all my application codes create a LibMeshInit object in main(). And I can't allocate from the heap: LibMeshInit *

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-06 Thread Roy Stogner
On Tue, 7 Aug 2012, Onur Solmaz wrote: I finally did it, however I started to have some other problems. In a short list: 1. I am trying to call libmesh from another function than main(). This works, but I'm not sure if this is the best way to do so. I use fake arguments for LibMeshInit()   i

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-04 Thread Roy Stogner
On Sat, 4 Aug 2012, Onur Solmaz wrote: FYI the quadrature rule defined for the faces does not get used in example 4 (used in example 3).   QGauss qface(dim-1, FIFTH);   fe_face->attach_quadrature_rule (&qface); These are the only lines, and removing these does not affect the result. We used

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-03 Thread Roy Stogner
On Fri, 3 Aug 2012, John Peterson wrote: > You are mistaken. The quadrature points on the face of the tet will > not coincide with the nodes of the tet unless you are using > trapezoidal quadrature. > > But I don't think the DirichletBoundary class of Introduction ex4 does > quadrature. It sets

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-03 Thread John Peterson
On Fri, Aug 3, 2012 at 1:13 PM, Onur Solmaz wrote: > That is what I tried at first, since the quadrature points are the vertices > of the tetrahedra (If I'm not mistaken). You are mistaken. The quadrature points on the face of the tet will not coincide with the nodes of the tet unless you are us

Re: [Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-03 Thread John Peterson
On Fri, Aug 3, 2012 at 2:10 AM, Onur Solmaz wrote: > > The problem is that the example (I think) uses the exact solution function > to impose boundary conditions in the element loop. I couldn't figure out a > workaround, so I am asking here, how can I impose nodal boundary conditions > in the asse

[Libmesh-users] Introduction Example 4 with Nodal BC's

2012-08-03 Thread Onur Solmaz
Hi I am having trouble writing a modified version of the assembly function of Introduction Ex 4. Example 4 demonstrates the method to solve the Poisson equation. My domain consists of a 3D unstructured mesh with tetrahedral elements. I haven't experienced any problems converting my data structures