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}
(46) 
_{Aug}
(63) 
_{Sep}
(84) 
_{Oct}
(82) 
_{Nov}
(69) 
_{Dec}
(45) 
2016 
_{Jan}
(92) 
_{Feb}
(91) 
_{Mar}
(148) 
_{Apr}
(43) 
_{May}
(58) 
_{Jun}
(117) 
_{Jul}
(92) 
_{Aug}
(136) 
_{Sep}

_{Oct}

_{Nov}

_{Dec}

S  M  T  W  T  F  S 




1
(5) 
2
(2) 
3
(1) 
4
(5) 
5
(2) 
6
(13) 
7
(5) 
8
(9) 
9
(6) 
10
(5) 
11
(1) 
12
(1) 
13
(9) 
14
(1) 
15
(3) 
16
(4) 
17
(2) 
18

19

20
(2) 
21

22
(1) 
23
(11) 
24

25

26
(2) 
27
(7) 
28
(12) 
29
(9) 
30
(7) 
31
(2) 

From: Roy Stogner <roystgnr@ic...>  20060301 20:25:33

On Wed, 1 Mar 2006, David Xu wrote: > I'd like to know when you expect to have a working 3D hermite > polynomial in place? Well, the bug hunting hasn't been high on my priority list since the 3D version of my method has a couple harder obstacles to overcome first. I can give it more attention if I know there's other people waiting to use the code, but it won't be fixed any time in the next week. > I'm trying to model 3D timeindependent schrodinger equation and > some references tell me using hermite polynomial with C1 continuity > has clear advantage over Lagrange polynomial. I'm not sure the advantages will be clear unless your weak formulation depends on integrals on a W^{2,p} space. C1 elements can give you similar accuracy with fewer DoFs (I think ~8 per element instead of ~27 when comparing C0 and C1 cubic hexes), I suppose. > What's the highest order of Lagrange polynomial supported? Only linears and quadratics are supported with a Lagrange basis on most elements, and AFAIK only pseudolinears are supported on pyramid elements. We've also got a set of hierarchic elements you might be interesed in. They're not very complete, but they do include cubic 3D hexes, without any rectilinear mesh restrictions. >> No. The higher order elements give you the chance to use higher order >> mapping functions, and give libMesh more places to put degrees of >> freedom (because edge and face DoFs are usually necessary for higher >> order shape functions), but if you just convert a QUAD4 mesh into a >> QUAD9 mesh while specifying bilinear shape functions it shouldn't >> change your results beyond whatever different floatingpoint errors >> do to the mapping functions. > > So, is the number of nodes in an element unrelated to the order of > Lagrange polynomial in libmesh? Not unrelated  if you want quadratic Lagrange elements, then you need second order geometric elements like QUAD9. However, if you want linear elements, then you can still use second order geometric elements (to better resolve a curved boundary, for instance). > How do I sepcify the order of gauss quadrature? When you construct a QGauss object, one of the parameters is the degree of polynomials you want to exactly integrate. You can also call FEType::default_quadrature_rule() if you have an FEType corresponding to an element of order p and you want a quadrature rule that exactly integrates polynomials of degree 2p+1. That's the safest way to go if you want compatibility with future more exotic elements, too. Macroelement basis functions require Gaussian rules on each of their subelements for exact integration, for example.  Roy 
From: David Xu <dxu@my...>  20060301 19:57:15

Roy, Thanks for your quick and thorough answers. I have some followup questions= : > > > 1. Does libmesh support Hermite polynomial shape function? > > Yes, but as far as I know only the 2D case (the BognerFoxSchmidt > rectangle) is well tested. I just caught a bug in the 3D Hermite a > little while ago, and I think there's still a bug remaining. I > haven't tested the 1D case at all except insofar as the 1D code gets > used by the 2D element. > > > If so, how do I specify it in the my code? > > By overwriting the default argument to ::add_variable(), like so: > > system.add_variable("u", THIRD, HERMITE); > > Note that only cubics (THIRD) are supported by the Hermite elements, > and your mesh had better be rectilinear if you want them to work > correctly. I'd like to know when you expect to have a working 3D hermite polynomial in place? I'm trying to model 3D timeindependent schrodinger equation and some references tell me using hermite polynomial with C1 continuity has clear advantage over Lagrange polynomial. What's the highest order of Lagrange polynomial supported? > > I didn't find any indication in the examples. > > That's my fault. I don't want to commit an updated example 15 until > I'm certain the 3D case is working correctly. > > > I assume they all use Lagrange polynomial as default? > > Yes. > > 2. How do I specify the polynomial order (p) in the code? > > In the second argument to add_variable(). > > > 3. In this line equation_systems("Poisson").add_variable("u", SECOND); > > Is "SECOND" the order of polynomial shape function or the order of > > gauss quadrature? > > The polynomial order. > > > 4. Do the element types (QUAD4, QUAD9, etc.) already include > > polynomial order information implicitly? > > No. The higher order elements give you the chance to use higher order > mapping functions, and give libMesh more places to put degrees of > freedom (because edge and face DoFs are usually necessary for higher > order shape functions), but if you just convert a QUAD4 mesh into a > QUAD9 mesh while specifying bilinear shape functions it shouldn't > change your results beyond whatever different floatingpoint errors > do to the mapping functions. So, is the number of nodes in an element unrelated to the order of Lagrange polynomial in libmesh? How do I sepcify the order of gauss quadrature? Thanks, David 
From: Roy Stogner <roystgnr@ic...>  20060301 19:17:39

On Wed, 1 Mar 2006, David Xu wrote: > 1. Does libmesh support Hermite polynomial shape function? Yes, but as far as I know only the 2D case (the BognerFoxSchmidt rectangle) is well tested. I just caught a bug in the 3D Hermite a little while ago, and I think there's still a bug remaining. I haven't tested the 1D case at all except insofar as the 1D code gets used by the 2D element. > If so, how do I specify it in the my code? By overwriting the default argument to ::add_variable(), like so: system.add_variable("u", THIRD, HERMITE); Note that only cubics (THIRD) are supported by the Hermite elements, and your mesh had better be rectilinear if you want them to work correctly. > I didn't find any indication in the examples. That's my fault. I don't want to commit an updated example 15 until I'm certain the 3D case is working correctly. > I assume they all use Lagrange polynomial as default? Yes. > 2. How do I specify the polynomial order (p) in the code? In the second argument to add_variable(). > 3. In this line equation_systems("Poisson").add_variable("u", SECOND); > Is "SECOND" the order of polynomial shape function or the order of > gauss quadrature? The polynomial order. > 4. Do the element types (QUAD4, QUAD9, etc.) already include > polynomial order information implicitly? No. The higher order elements give you the chance to use higher order mapping functions, and give libMesh more places to put degrees of freedom (because edge and face DoFs are usually necessary for higher order shape functions), but if you just convert a QUAD4 mesh into a QUAD9 mesh while specifying bilinear shape functions it shouldn't change your results beyond whatever different floatingpoint errors do to the mapping functions.  Roy Stogner 
From: David Xu <dxu@my...>  20060301 18:59:50

Hi Steffen, I'm looking forward to your example 17 and implementation of generalized eigen value solver with libmesh. Thanks, David On 2/27/06, Steffen Petersen <steffen.petersen@...> wrote: > I have just committed an example 17 including a generalized > eigenvalue problem. Currently, the largest eigenvalues > in the spectrum are computed. The code can be easily changed > in order to compute the smallest values (or use the slepc comman > line options). However, the Lanczos solver slepc offers with this > option appears to be somewhat unstable when used to compute > the smallest values (perhaps this is due to the zero eigenmode in > this example). An option to work around this problem is to > configure and install slepc with arpack and parpack (the arpack > solver works fine for example 17). > > Steffen > > > > David Xu schrieb: > > Hi Steffen and Ondrej, > > > > I read your previous message that you've extended the slepc interface > > to support generalized eigen problems. I'd like to get the details > > from you. I'm new to libmesh and I'm wondering if you can post some > > new tutorial examples on how to solve eigen problems using libmesh > > (eg. solve 1D or 2D simple harmonic oscillator)? That will a great > > help. > > > > Ondrej, I read about your message about your implementation of libmesh > > and pysparse's jacobidavidson solver. That sounds really good. Could > > you please give me the details too? I really appreciate it. > > > > David > > 
From: David Xu <dxu@my...>  20060301 18:54:28

Hi everyone, My questions are: 1. Does libmesh support Hermite polynomial shape function? If so, how do I specify it in the my code? I didn't find any indication in the examples. I assume they all use Lagrange polynomial as default? 2. How do I specify the polynomial order (p) in the code? 3. In this line equation_systems("Poisson").add_variable("u", SECOND); Is "SECOND" the order of polynomial shape function or the order of gauss quadrature? 4. Do the element types (QUAD4, QUAD9, etc.) already include polynomial order information implicitly? Thanks for your time, David 