You can subscribe to this list here.
2003 
_{Jan}

_{Feb}

_{Mar}

_{Apr}

_{May}

_{Jun}

_{Jul}

_{Aug}

_{Sep}
(2) 
_{Oct}
(2) 
_{Nov}
(27) 
_{Dec}
(31) 

2004 
_{Jan}
(6) 
_{Feb}
(15) 
_{Mar}
(33) 
_{Apr}
(10) 
_{May}
(46) 
_{Jun}
(11) 
_{Jul}
(21) 
_{Aug}
(15) 
_{Sep}
(13) 
_{Oct}
(23) 
_{Nov}
(1) 
_{Dec}
(8) 
2005 
_{Jan}
(27) 
_{Feb}
(57) 
_{Mar}
(86) 
_{Apr}
(23) 
_{May}
(37) 
_{Jun}
(34) 
_{Jul}
(24) 
_{Aug}
(17) 
_{Sep}
(50) 
_{Oct}
(24) 
_{Nov}
(10) 
_{Dec}
(60) 
2006 
_{Jan}
(47) 
_{Feb}
(46) 
_{Mar}
(127) 
_{Apr}
(19) 
_{May}
(26) 
_{Jun}
(62) 
_{Jul}
(47) 
_{Aug}
(51) 
_{Sep}
(61) 
_{Oct}
(42) 
_{Nov}
(50) 
_{Dec}
(33) 
2007 
_{Jan}
(60) 
_{Feb}
(55) 
_{Mar}
(77) 
_{Apr}
(102) 
_{May}
(82) 
_{Jun}
(102) 
_{Jul}
(169) 
_{Aug}
(117) 
_{Sep}
(80) 
_{Oct}
(37) 
_{Nov}
(51) 
_{Dec}
(43) 
2008 
_{Jan}
(71) 
_{Feb}
(94) 
_{Mar}
(98) 
_{Apr}
(125) 
_{May}
(54) 
_{Jun}
(119) 
_{Jul}
(60) 
_{Aug}
(111) 
_{Sep}
(118) 
_{Oct}
(125) 
_{Nov}
(119) 
_{Dec}
(94) 
2009 
_{Jan}
(109) 
_{Feb}
(38) 
_{Mar}
(93) 
_{Apr}
(88) 
_{May}
(29) 
_{Jun}
(57) 
_{Jul}
(53) 
_{Aug}
(48) 
_{Sep}
(68) 
_{Oct}
(151) 
_{Nov}
(23) 
_{Dec}
(35) 
2010 
_{Jan}
(84) 
_{Feb}
(60) 
_{Mar}
(184) 
_{Apr}
(112) 
_{May}
(60) 
_{Jun}
(90) 
_{Jul}
(23) 
_{Aug}
(70) 
_{Sep}
(119) 
_{Oct}
(27) 
_{Nov}
(47) 
_{Dec}
(54) 
2011 
_{Jan}
(22) 
_{Feb}
(19) 
_{Mar}
(92) 
_{Apr}
(93) 
_{May}
(35) 
_{Jun}
(91) 
_{Jul}
(32) 
_{Aug}
(61) 
_{Sep}
(7) 
_{Oct}
(69) 
_{Nov}
(81) 
_{Dec}
(23) 
2012 
_{Jan}
(64) 
_{Feb}
(95) 
_{Mar}
(35) 
_{Apr}
(36) 
_{May}
(63) 
_{Jun}
(98) 
_{Jul}
(70) 
_{Aug}
(171) 
_{Sep}
(149) 
_{Oct}
(64) 
_{Nov}
(67) 
_{Dec}
(126) 
2013 
_{Jan}
(108) 
_{Feb}
(104) 
_{Mar}
(171) 
_{Apr}
(133) 
_{May}
(108) 
_{Jun}
(100) 
_{Jul}
(93) 
_{Aug}
(126) 
_{Sep}
(74) 
_{Oct}
(59) 
_{Nov}
(145) 
_{Dec}
(93) 
2014 
_{Jan}
(38) 
_{Feb}
(45) 
_{Mar}
(26) 
_{Apr}
(41) 
_{May}
(125) 
_{Jun}
(70) 
_{Jul}
(61) 
_{Aug}
(66) 
_{Sep}
(60) 
_{Oct}
(110) 
_{Nov}
(27) 
_{Dec}
(30) 
2015 
_{Jan}
(43) 
_{Feb}
(67) 
_{Mar}
(71) 
_{Apr}
(92) 
_{May}
(39) 
_{Jun}
(15) 
_{Jul}
(3) 
_{Aug}

_{Sep}

_{Oct}

_{Nov}

_{Dec}

From: KIRK, BENJAMIN (JSCEG) (NASA) <benjamin.kirk1@na...>  20051214 19:03:16

You might want to check out SLEPc, an interface on top op PETSc specifically for Eigenvalue problems. Stefan added an interface to it from libMesh some time ago, and he might have some comments on this thread too. Ben _____ From: libmeshusersadmin@... [mailto:libmeshusersadmin@...] On Behalf Of Michael Povolotskyi Sent: Wednesday, December 14, 2005 12:22 PM To: Roy Stogner; libmesh  ML Subject: Re: [Libmeshusers] Eigenvalue problem 1) Yes, I mean Hx = gSx, sorry for the error. 2) Does it mean that if one wants to use Libmesh just for one complex problem, then all the other problems also have to be complex? thank you, Michael. Roy Stogner wrote: On Wed, 14 Dec 2005, Michael Povolotskyi wrote: I have to solve an eigenvalue problem for a quantum physics calculations. My questions are: 1) Is it possible to solve a general eigenvalue problem like Hx = Sx, where H and S are known matrixes. Are you sure you don't mean Hx = gSx, where H and S are known but g is unknown? That's the generalized eigenvalue problem built into SLEPc, which is currently the only eigenvalue solve libMesh interfaces with. Our SLEPc support currently only handles the simple Hx = gx problem, and a new interface would have to be written to get generalized eigenvalues. 2) Is it possible to have complex (selfconjugated) matrixes? If both PETSc and libMesh have been compiled with complex number support on, then the default PetscScalar and libMesh Number types will both be std::complex<double>, and PETSc should handle the linear algebra appropriately. All my code is realvalued, though; I'm not sure who's keeping the libMesh complex number support up to date.  Roy Stogner   Michael Povolotskyi, Ph.D. University of Rome "Tor Vergata" Department of Electronic Engineering Viale Politecnico, 1  00133 Rome  Italy Phone + 39 06 72597781 Fax + 39 06 2020519 http://www.optolab.uniroma2.it/pages/moshe/moshe.html <http://www.optolab.uniroma2.it/pages/moshe/moshe.html>;   This SF.net email is sponsored by: Splunk Inc. Do you grep through log files for problems? Stop! Download the new AJAX search engine that makes searching your log files as easy as surfing the web. DOWNLOAD SPLUNK! http://ads.osdn.com/?ad_id=7637&alloc_id=16865&op=click _______________________________________________ Libmeshusers mailing list Libmeshusers@... https://lists.sourceforge.net/lists/listinfo/libmeshusers 
From: Roy Stogner <roystgnr@ic...>  20051214 18:47:33

On Wed, 14 Dec 2005, John Peterson wrote: > Michael Povolotskyi writes: > > 2) Does it mean that if one wants to use Libmesh just for one complex > > problem, then all the other problems also have to be complex? <br> > > Well, if you have enough disk space, you could have a libmesh_complex > build in a separate directory which is separate and only used for > complex problems. Wouldn't you also need a separate petsc_complex build, and scripts (like the "module load petsc/some_version" we use in CFDlab) to switch around the environment variables PETSc uses? It's not a good solution, but I don't think it'll change any time soon. We could build libMesh with simultaneous Real/Complex support just by adding another template argument to everything which isn't already templated around datatype... but PETSc is C and Fortran, and making that simultaneously compile with multiple datatypes would be pretty tedious.  Roy 
From: John Peterson <peterson@cf...>  20051214 18:38:11

Michael Povolotskyi writes: > 2) Does it mean that if one wants to use Libmesh just for one complex > problem, then all the other problems also have to be complex? <br> Well, if you have enough disk space, you could have a libmesh_complex build in a separate directory which is separate and only used for complex problems. J 
From: Michael Povolotskyi <povolotskyi@in...>  20051214 18:22:29

<!DOCTYPE html PUBLIC "//W3C//DTD HTML 4.01 Transitional//EN"> <html> <head> <meta content="text/html;charset=ISO88591" httpequiv="ContentType"> <title></title> </head> <body bgcolor="#cccccc" text="#000000"> 1) Yes, I mean Hx = gSx, sorry for the error.<br> 2) Does it mean that if one wants to use Libmesh just for one complex problem, then all the other problems also have to be complex? <br> <br> thank you,<br> Michael.<br> <br> <br> <br> Roy Stogner wrote: <blockquote cite="midPine.LNX.4.63.0512141154180.4247@..." type="cite">On Wed, 14 Dec 2005, Michael Povolotskyi wrote: <br> <br> <blockquote type="cite">I have to solve an eigenvalue problem for a quantum physics calculations. <br> My questions are: <br> 1) Is it possible to solve a general eigenvalue problem like Hx = Sx, <br> where H and S are known matrixes. <br> </blockquote> <br> Are you sure you don't mean Hx = gSx, where H and S are known but g is <br> unknown? That's the generalized eigenvalue problem built into SLEPc, <br> which is currently the only eigenvalue solve libMesh interfaces with. <br> Our SLEPc support currently only handles the simple Hx = gx problem, <br> and a new interface would have to be written to get generalized <br> eigenvalues. <br> <br> <blockquote type="cite">2) Is it possible to have complex (selfconjugated) matrixes? <br> </blockquote> <br> If both PETSc and libMesh have been compiled with complex number <br> support on, then the default PetscScalar and libMesh Number types will <br> both be std::complex<double>, and PETSc should handle the linear <br> algebra appropriately. All my code is realvalued, though; I'm not <br> sure who's keeping the libMesh complex number support up to date. <br>  <br> Roy Stogner <br> <br> </blockquote> <br> <br> <pre class="mozsignature" cols="250">  Michael Povolotskyi, Ph.D. University of Rome "Tor Vergata" Department of Electronic Engineering Viale Politecnico, 1  00133 Rome  Italy Phone + 39 06 72597781 Fax + 39 06 2020519 <a class="moztxtlinkfreetext" href="http://www.optolab.uniroma2.it/pages/moshe/moshe.html">http://www.optolab.uniroma2.it/pages/moshe/moshe.html</a>;  </pre> </body> </html> 
From: Roy Stogner <roystgnr@ic...>  20051214 18:10:35

On Wed, 14 Dec 2005, Michael Povolotskyi wrote: > I have to solve an eigenvalue problem for a quantum physics calculations. > My questions are: > 1) Is it possible to solve a general eigenvalue problem like Hx = Sx, > where H and S are known matrixes. Are you sure you don't mean Hx = gSx, where H and S are known but g is unknown? That's the generalized eigenvalue problem built into SLEPc, which is currently the only eigenvalue solve libMesh interfaces with. Our SLEPc support currently only handles the simple Hx = gx problem, and a new interface would have to be written to get generalized eigenvalues. > 2) Is it possible to have complex (selfconjugated) matrixes? If both PETSc and libMesh have been compiled with complex number support on, then the default PetscScalar and libMesh Number types will both be std::complex<double>, and PETSc should handle the linear algebra appropriately. All my code is realvalued, though; I'm not sure who's keeping the libMesh complex number support up to date.  Roy Stogner 
From: Michael Povolotskyi <povolotskyi@in...>  20051214 17:41:00

<!DOCTYPE html PUBLIC "//W3C//DTD HTML 4.01 Transitional//EN"> <html> <head> <meta content=3D"text/html;charset=3DKOI8R" httpequiv=3D"ContentType= "> <title></title> </head> <body bgcolor=3D"#cccccc" text=3D"#000000"> Dear Libmesh developers,<br> I have to solve an eigenvalue problem for a quantum physics calculations.<br> My questions are:<br> 1) Is it possible to solve a general eigenvalue problem like Hx =3D Sx, where H and S are known matrixes.<br> 2) Is it possible to have complex (selfconjugated) matrixes?<br> <br> Thank you,<br> Michael.<br> =9A<br> </body> </html> 
From: John Peterson <peterson@cf...>  20051213 17:59:12

hi, Due to recent vandalism, I removed the ability for anonymous users to create an account (I think!) Furthermore, anonymous users are also not allowed to edit the pages. If you want to sign up to edit the wiki, just send me an email on this list and I'll see if I can add a new user. This is a bit draconian, and not in the spirit of cooperative documentation, but the alternative is me constantly removing ads for Cialis from our wiki. J Ondrej Certik writes: > That's me, unfortunately, I don't have time right now. The program can be found > at: > > http://ondrej.certik.cz/libmesh/heat.tar.bz2 > > and some documentation at: > > http://ondrej.certik.cz/libmesh/fem.ps > > I also have some other programs: > > http://ondrej.certik.cz/libmesh/laplace.tar.bz2 > http://ondrej.certik.cz/libmesh/wave.tar.bz2 > > Could someone (Stefano) please write some info about that to the wiki? > > > Ondrej > > On Tue, Dec 13, 2005 at 08:19:49AM 0600, KIRK, BENJAMIN (JSCEG) (NASA) wrote: > > Yes, in fact there used to be a nice link from the wiki describing just > > that. Unfortunately, that page did not survive the migration to the new > > mysql database. > > > > Could the original author restore the page? > > > > Sorry. > > > > Ben > > > > _____ > > 
From: Ondrej Certik <ondrej@ce...>  20051213 14:28:38

That's me, unfortunately, I don't have time right now. The program can be found at: http://ondrej.certik.cz/libmesh/heat.tar.bz2 and some documentation at: http://ondrej.certik.cz/libmesh/fem.ps I also have some other programs: http://ondrej.certik.cz/libmesh/laplace.tar.bz2 http://ondrej.certik.cz/libmesh/wave.tar.bz2 Could someone (Stefano) please write some info about that to the wiki? Ondrej On Tue, Dec 13, 2005 at 08:19:49AM 0600, KIRK, BENJAMIN (JSCEG) (NASA) wrote: > Yes, in fact there used to be a nice link from the wiki describing just > that. Unfortunately, that page did not survive the migration to the new > mysql database. > > Could the original author restore the page? > > Sorry. > > Ben > > _____ > > From: libmeshusersadmin@... > [mailto:libmeshusersadmin@...] On Behalf Of Stefano > Sent: Tuesday, December 13, 2005 7:47 AM > To: libmeshusers@... > Subject: [Libmeshusers] heat equation 2D > > > Has any of you have an example of employment of libmesh for the solution a > 2D Heat equation problem? > Regards > > Stefano > 
From: KIRK, BENJAMIN (JSCEG) (NASA) <benjamin.kirk1@na...>  20051213 14:19:49

Yes, in fact there used to be a nice link from the wiki describing just that. Unfortunately, that page did not survive the migration to the new mysql database. Could the original author restore the page? Sorry. Ben _____ From: libmeshusersadmin@... [mailto:libmeshusersadmin@...] On Behalf Of Stefano Sent: Tuesday, December 13, 2005 7:47 AM To: libmeshusers@... Subject: [Libmeshusers] heat equation 2D Has any of you have an example of employment of libmesh for the solution a 2D Heat equation problem? Regards Stefano 
From: Stefano <sboschi75personal@ya...>  20051213 13:46:53

Has any of you have an example of employment of libmesh for the solution a 2D Heat equation problem? Regards Stefano 
From: Roy Stogner <roystgnr@ic...>  20051210 00:21:50

On Sat, 10 Dec 2005, David Knezevic wrote: > I've been trying to have a play with libMesh's hierarchic finite elements on > a simple square mesh but I'm not sure what is required to get them working. > What kind of a mesh does one need to generate in order to do, say: > "system.add_variable("u",FIRST,HIERARCHIC)" You'll need second order elements: TRI6 or QUAD9. Other than that there shouldn't be any special requirements  what problems are you having?  Roy 
From: David Knezevic <david.knezevic@ba...>  20051210 00:09:38

Hi, I've been trying to have a play with libMesh's hierarchic finite elements on a simple square mesh but I'm not sure what is required to get them working. What kind of a mesh does one need to generate in order to do, say: "system.add_variable("u",FIRST,HIERARCHIC)" Cheers, Dave 
From: li pan <li76pan@ya...>  20051209 17:19:02

hi, could anybody tell me, how the format of vector and matrix in Newmark system look like? For example: dim = 3; number of nodes is 4; variable u means "displacement" I'm not sure whether u has format (u1_x u2_x u3_x u4_x u1_y u2_y u3_y u4_y u1_z u2_z u3_z u4_z) or (u1_x u1_y u1_z u2_x u2_y u2_z u3_x u3_y u3_z u4_x u4_y u4_z) And the same with stiffness matrix. Whether Ke has format Kuu Kuv Kuw Kvu Kvv Kvw Kwu Kwv Kww (Kuu and other submatrix are 4*4 matrix) or K11 K12 K13 K14 K21 K22 K23 K24 K31 K32 K33 K34 K41 K42 K43 K44 (K11 and other submatrix are 3*3 matrix)? Can anybody help me? Thanx a lot. best pan __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com 
From: li pan <li76pan@ya...>  20051209 13:50:20

the following problem has been solved. I just didn't realize that there was only one variable declared. Sorry. best pan hi everyone, I want to use NewmarkSystem to do integration. I checked ex8. The cube was 3 dimentional. But all the variables, i.e. displacement, velocity, accelereation, stiffness matrix etc, are one dimentional. It seems that only one scalar of the vector was calculated. How can I deal with 3D problem? Does anybody have an idea? best pan __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com 
From: li pan <li76pan@ya...>  20051209 12:59:39

hi everyone, I want to use NewmarkSystem to do integration. I checked ex8. The cube was 3 dimentional. But all the variables, i.e. displacement, velocity, accelereation, stiffness matrix etc, are one dimentional. It seems that only one scalar of the vector was calculated. How can I deal with 3D problem? Does anybody have an idea? best pan __________________________________________________ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com 
From: Roy Stogner <roystgnr@ic...>  20051207 15:55:12

On Wed, 7 Dec 2005, Roy Stogner wrote: > On Wed, 7 Dec 2005, li pan wrote: > >> only one question. should I install any other things >> beside libmesh? > > Actually, that might be the problem. I know at one time we an example > (possibly ex13) that the LASPACK solver failed on. Even for > nonparallel solves we recommend using PETSc. That's supposed to say "we had an example"... but now that I've looked it up, I found that the example that failed with LASPACK was ex8, not ex13. I'm not sure what could be going wrong for you, but when I find some time today I'll try changing the mesh type myself and see if I can replicate the problem with PETSc and the CVS libMesh version.  Roy Stogner 
From: Roy Stogner <roystgnr@ic...>  20051207 15:38:27

On Wed, 7 Dec 2005, li pan wrote: > only one question. should I install any other things > beside libmesh? Actually, that might be the problem. I know at one time we an example (possibly ex13) that the LASPACK solver failed on. Even for nonparallel solves we recommend using PETSc.  Roy 
From: Roy Stogner <roystgnr@ic...>  20051207 14:25:44

On Wed, 7 Dec 2005, li pan wrote: > I read the code in ex13. It used the cubic element. So > I wonder, if I use tetrahedron element, should I > change something in the function "assemble_stokes"? I don't think so. The trick with the velocitypressure formulation of NavierStokes is choosing a stable pair of finite elements for your discretization, but quadratic velocity + linear pressure should work as well as biquadratic velocity + bilinear pressure. If you wanted to mix cubic and tetrahedral elements, there might be a problem  you'd need pyramids and possibly prisms in a hybrid mesh, and last time I checked (which admittedly was a while back) libMesh only supported first order functions on pyramids.  Roy Stogner 
From: Roy Stogner <roystgnr@ic...>  20051207 14:03:47

On Wed, 7 Dec 2005, li pan wrote: > Hi everyone, > I installed libmesh0.5.0. As I tried to start ex13, I > got this output: > EquationSystems > n_systems()=1 > System "NavierStokes" > Type "TransientLinearImplicit" > Variables="u" "v" "p" > Finite Element Types="0" "0" "0" > Approximation Orders="2" "2" "1" > n_dofs()=3803 > n_local_dofs()=3803 > n_constrained_dofs()=0 > n_vectors()=3 > It's different from the online document ex13.php, > which I attach here: > EquationSystems > n_systems()=1 > System "NavierStokes" > Type "TransientImplicit" > Variables="u" "v" "p" > Finite Element Types="0", "12" "0", "12" "0", "12" > > Infinite Element Mapping="0" "0" "0" > Approximation Orders="2", "3" "2", "3" "1", "3" > n_dofs()=3803 > n_local_dofs()=3803 > n_constrained_dofs()=0 > n_vectors()=1 > n_parameters()=2 > Parameters: > "linear solver maximum iterations"=250 > "linear solver tolerance"=0.001 > Can anybody help me? I don't think the differences are anything to worry about. The online documents are generated from the CVS head rather than the last development snapshot, and nearly all the differences are due to interface changes or additional features in the CVS version. The remaining differences appear to be because the version used for the online document was compiled with infinite element support turned on.  Roy 
From: li pan <li76pan@ya...>  20051207 10:29:24

Hi everyone, I installed libmesh0.5.0. As I tried to start ex13, I got this output: EquationSystems n_systems()=1 System "NavierStokes" Type "TransientLinearImplicit" Variables="u" "v" "p" Finite Element Types="0" "0" "0" Approximation Orders="2" "2" "1" n_dofs()=3803 n_local_dofs()=3803 n_constrained_dofs()=0 n_vectors()=3 It's different from the online document ex13.php, which I attach here: EquationSystems n_systems()=1 System "NavierStokes" Type "TransientImplicit" Variables="u" "v" "p" Finite Element Types="0", "12" "0", "12" "0", "12" Infinite Element Mapping="0" "0" "0" Approximation Orders="2", "3" "2", "3" "1", "3" n_dofs()=3803 n_local_dofs()=3803 n_constrained_dofs()=0 n_vectors()=1 n_parameters()=2 Parameters: "linear solver maximum iterations"=250 "linear solver tolerance"=0.001 Can anybody help me? best pan __________________________________________ Yahoo! DSL Something to write home about. Just $16.99/mo. or less. dsl.yahoo.com 
From: li pan <li76pan@ya...>  20051207 10:29:24

Hi everyone, I installed libmesh0.5.0. As I tried to start ex13, I got this output: EquationSystems n_systems()=1 System "NavierStokes" Type "TransientLinearImplicit" Variables="u" "v" "p" Finite Element Types="0" "0" "0" Approximation Orders="2" "2" "1" n_dofs()=3803 n_local_dofs()=3803 n_constrained_dofs()=0 n_vectors()=3 It's different from the online document ex13.php, which I attach here: EquationSystems n_systems()=1 System "NavierStokes" Type "TransientImplicit" Variables="u" "v" "p" Finite Element Types="0", "12" "0", "12" "0", "12" Infinite Element Mapping="0" "0" "0" Approximation Orders="2", "3" "2", "3" "1", "3" n_dofs()=3803 n_local_dofs()=3803 n_constrained_dofs()=0 n_vectors()=1 n_parameters()=2 Parameters: "linear solver maximum iterations"=250 "linear solver tolerance"=0.001 Can anybody help me? best pan __________________________________________ Yahoo! DSL Something to write home about. Just $16.99/mo. or less. dsl.yahoo.com 
From: li pan <li76pan@ya...>  20051207 09:32:12

Hi Roy, I read the code in ex13. It used the cubic element. So I wonder, if I use tetrahedron element, should I change something in the function "assemble_stokes"? best pan On Tue, 6 Dec 2005, li pan wrote: > can anybody tell me, whether I can use libmesh to > solve elastic problem of solid body? Thanx before > hands. It's certainly possible. The only trick I can think of is that libMesh only supports scalar variables, so you'll have to add two or three "u", "v", "w" variables to your system rather than a single vector "u". This isn't hard, though  for similar code, see the NavierStokes equations in example 13 on the website.  Roy __________________________________________ Yahoo! DSL Something to write home about. Just $16.99/mo. or less. dsl.yahoo.com 
From: Roy Stogner <roystgnr@ic...>  20051206 18:54:17

On Tue, 6 Dec 2005, li pan wrote: > can anybody tell me, whether I can use libmesh to > solve elastic problem of solid body? Thanx before > hands. It's certainly possible. The only trick I can think of is that libMesh only supports scalar variables, so you'll have to add two or three "u", "v", "w" variables to your system rather than a single vector "u". This isn't hard, though  for similar code, see the NavierStokes equations in example 13 on the website.  Roy 
From: li pan <li76pan@ya...>  20051206 18:09:06

Hello everyone, can anybody tell me, whether I can use libmesh to solve elastic problem of solid body? Thanx before hands. best pan __________________________________________ Yahoo! DSL Something to write home about. Just $16.99/mo. or less. dsl.yahoo.com 
From: Roy Stogner <roystgnr@ic...>  20051201 20:17:05

On Thu, 1 Dec 2005, Michael Povolotskyi wrote: > I'd like to say that for Lagrange elements one can code a noniterative > algorithm in order to find if a point is inside a polyhedral or not. Does this apply even if the region is not polyhedral or polygonal? It's easy for me to see how to quickly test if a point is within a collection of flat sides, but with secondorder elements it's possible for libMesh users to use curved sides, and with hexes and prisms it's almost impossible for users to avoid curved sides.  Roy 