On 10/19/11 5:02 AM, Lorenzo Alessio Botti wrote:
Looking at the implementation of FE::side_map, it seems like it should
be possible to avoid calling inverse_map() in the reinit(elem,side,...), even when
there are user-supplied points, by doing something like:
// Find where the integratio
On Wed, 19 Oct 2011, Roy Stogner wrote:
> so why the hell haven't I added TBB to FEMSystem yet?"
FEMSystem is now multithreaded (for assembly, postprocessing, qoi
assembly, and qoi derivative assembly) as of r4878. It's passing my
basic initial checks; we'll hit it with some stricter regression
On Wed, 19 Oct 2011, Paul T. Bauman wrote:
> I have up and running tests using FEMSystem for a variety of problems that
> can also be used as guinea pigs. Let me know how I
> can help here. If it's trivial, might be a nice stepping stone into TBB for
> me... ;-)
Yeah, I'm an idiot. When I sa
On Wed, Oct 19, 2011 at 2:27 PM, Roy Stogner wrote:
>
> On Wed, 19 Oct 2011, David Fuentes wrote:
>
> > I probably need my version of libmesh, but
> > Is there an example that uses threading for Matrix assembly for FemSystem
> ?
>
> Not at the moment. I made the DiffContext change to make the
> F
On Wed, 19 Oct 2011, David Fuentes wrote:
> I probably need my version of libmesh, but
> Is there an example that uses threading for Matrix assembly for FemSystem ?
Not at the moment. I made the DiffContext change to make the
FEMSystem API thread-safe as soon as I noticed that issue, but I
have
Hi,
I probably need my version of libmesh, but
Is there an example that uses threading for Matrix assembly for FemSystem ?
thanks,
David
--
The demand for IT networking professionals continues to grow, and the
demand for
On Wed, Oct 19, 2011 at 9:05 AM, Matteo Semplice
wrote:
>
> Well yes, but there are a couple of quirks:
>
> 1) The mesh must be initialized setting its dimension, like
> Mesh mesh(2);
> since the Mesh constructor defauts to 1 and GmshIO::read_mesh chooses dim as
> const unsigned int dim =
On 19/10/2011 14:05, David Knezevic wrote:
> On 10/19/2011 07:41 AM, Matteo Semplice wrote:
>
>> Using gmsh was indeed my first attempt, but something does not work here.
>> Consider the following geometry saved in squarewithhole.geo:
>> -
>> Point(1) = {0, 0, 0};
>> Point(2) = {0, 1, 0
On 10/19/2011 07:41 AM, Matteo Semplice wrote:
> Using gmsh was indeed my first attempt, but something does not work here.
> Consider the following geometry saved in squarewithhole.geo:
> -
> Point(1) = {0, 0, 0};
> Point(2) = {0, 1, 0};
> Point(3) = {1, 1, 0};
> Point(4) = {1, 0, 0};
> Poi
Hi Matteo,
I usually use gmsh. You can set gmsh "physicals" on boundary sides, and
they will be read in by libMesh as boundary IDs.
David
On 10/19/2011 06:48 AM, robert wrote:
> Hello,
>
> I had success using the exodusII file format, which I had generated with
> cubit[1]. In addition to gene
Hello,
I had success using the exodusII file format, which I had generated with
cubit[1]. In addition to generating the geometry and meshing it, you
have to create "sidesets" containing the surfaces on which you want to
apply your BC.
Robert
[1] http://cubit.sandia.gov/
On Wed, Oct 19, 2011 at
Dear all,
I am need to solve a PDE in which different kinds of boundry
conditions are applied in different portions of the boundary. I can do
that if I generate the mesh with libmesh and using
MeshTools::Modification::change_boundary_id. However I will need to deal
with complicated geometries
>
> Looking at the implementation of FE::side_map, it seems like it should
> be possible to avoid calling inverse_map() in the reinit(elem,side,...), even
> when there are user-supplied points, by doing something like:
>
> // Find where the integration points are located on the
> // full eleme
13 matches
Mail list logo
