Another possibility is to use the FEMSystem framework, see fem_system
example 3 for example, specifically the ' bool
ElasticitySystem::side_time_derivative (bool request_jacobian, DiffContext
& context) ' function in elasticity_system.C.
This function takes in the unsteady contributions to the res
On Mon, Jun 13, 2016 at 2:43 PM, Ryan Woodall wrote:
> Hi everyone,
>
> I've got a problem where I'm try to implement 2 types of boundary
> conditions on a FEM model. The first type of BC i have is a zero-flux BC,
> which is pretty easy to implement. However, I'm also trying to implement
> time-d
Hi everyone,
I've got a problem where I'm try to implement 2 types of boundary
conditions on a FEM model. The first type of BC i have is a zero-flux BC,
which is pretty easy to implement. However, I'm also trying to implement
time-dependant Neumann BC's, and I'm not sure how to do that. They are
e
On Mon, Jun 13, 2016 at 11:03 AM, Xujun Zhao wrote:
> Hi John,
>
> Thank you for your information. I am using the recent version of libMesh.
>
> Also I finally figured out the culprit last weekend, which is due to hypre
> preconditioner used in one of my fieldsplit variables.
> The BoomerAMG uses
On Mon, Jun 13, 2016 at 11:02 AM, li.luo wrote:
> Thanks.
> I thought of another way may work, that is, to use p2p2 for velocity and
> pressure, then "disregard" those dofs belong to 2nd order elements, and
> set the related rows in matrix to be identities. That I use a larger
> problem, but v a
Hi John,
Thank you for your information. I am using the recent version of libMesh.
Also I finally figured out the culprit last weekend, which is due to hypre
preconditioner used in one of my fieldsplit variables.
The BoomerAMG uses very large memory during the KSPSolve()
-Xujun
On Mon, Jun 13,
On Mon, Jun 13, 2016 at 12:00 AM, wrote:
> Dear developers,
>
> I created a 2nd order finite element, TET or CUBE, however for some
> reasons, I want to identify the 2nd order of dofs in this element.
> For example, the continuity equaiton in N-S equations, if I set 2nd order
> for pressure, but
On Fri, Jun 10, 2016 at 9:21 AM, 张江 wrote:
> Hi,
>
> I am trying to load vtk files using libmesh. But there is always error that
>
> ERROR: In /Users/xxx/VTK-7.0.0/IO/XMLParser/vtkXMLParser.cxx, line 397
> vtkXMLDataParser (0x7f9341c16140): Error parsing XML in stream at line 1,
> column 1, byte
On Thu, Jun 9, 2016 at 3:29 PM, Xujun Zhao wrote:
> Derek,
>
> Excellent analysis! It really helps.
> I further looked at the PETSc log summary when running on 1 CPU, the max
> memory Petsc allocated is about 22.2G (close to what you predicted).
> However, the total memory usage is up to 100G. Th
On Fri, Jun 10, 2016 at 10:01 AM, wrote:
> Dear developers,
>
> I enable the block storage for matrices, however, if I use P2P1 FE type
> for Navier-Stokes equatioins, the block type reduces to non-block type.
>
Hmm.. I don't think this is possible. I just looked at the documentation
for DofM
On Tue, Jun 7, 2016 at 6:19 PM, Jayaraman, Vegnesh
wrote:
> Hi
>
> Can the matrix->zero_rows() function not be used during assembly? I get an
> "Object Not Initialized " error during runtime if I use the above function
> in assembly.
>
zero_rows() can't be called interspersed with the MatSetValu
On Thu, Jun 9, 2016 at 4:47 PM, David Knezevic
wrote:
> On Thu, Jun 9, 2016 at 6:01 PM, Bhalla, Amneet Pal S >
> wrote:
>
> > Hi Folks,
> >
> >
> > I am trying to import a 2d mesh that I created in trelis to libmesh, but
> I
> > am getting error of the type
> >
> > File opened successfully.
> >
On Sat, Jun 11, 2016 at 8:30 PM, Jayaraman, Vegnesh
wrote:
> Hi
>
> I am trying to install Libmesh with DTK and Trilinos
>
> I am getting the following error
>
> /usr/bin/ld:
> /home/vegnesh/codes/trilinos/testbuild/DataTransferKit/packages/Adapters/ClassicDTK/src/libclassic_dtk.a(DTK_CommIndexe
13 matches
Mail list logo
