On Wed, Oct 3, 2012 at 2:27 AM, Robert Cimrman cimr...@ntc.zcu.cz wrote:
On 10/03/2012 11:06 AM, Ondřej Čertík wrote:
On Wed, Sep 26, 2012 at 9:20 AM, Robert Cimrman cimr...@ntc.zcu.cz wrote:
Hi,
this is to inform you about the latest updates of the code.
- The 'hierarchic' branch has been merged into the master. This branch introduces a hierarchical FE basis into the code and some other tools described below. Look at [1] for an example. It's not finished, but I have decided to merge it already, as the general ideas seem to be work, and I needed some of the developments in the master and vice-versa. The caveats are for now:
- Only 2_4 (rectangles) and 3_8 (bricks) elements are supported.
That's exciting. Is it using the Cartesian product on each element?
Yes, but the order in each direction as well as in all the elements is the same.
What is the maximum polynomial order that it can handle?
The 1D base functions are generated up to order 10. But there is a limitation given by the fact that we assemble all elements in a group at once, so unintegrated element matrices have to fit into memory. In 3D for the 1D order 7 the element matrix has shape 512x512, and the number of quadrature points is also 512 so 8 * 512**3 / 1e9 = 1.073741824, that is 1GB for a single element. This will have to be dealt with.
I see. Looks like one cannot precalculate everything on quadrature points in advance. So I guess one will have to deal with it element by element.
Ondrej