uikit0/mfem

Automatically exported from code.google.com/p/mfem

mfem's Issues

How to build mfem in windows 7 64bit + VS2010/VS2012?

Dear all,
Can I use the MFEM in windows 7 64bit + VS2010/VS2012? How to build it?

Thanks,
Tang Laoya

Original issue reported on code.google.com by [email protected] on 22 Dec 2013 at 11:43

Testing mfem-dev hookup

Test

Original issue reported on code.google.com by tzanio on 2 Oct 2013 at 8:52

How to compute gradient of displacement and other functions in ex2.cpp

I am new to MFEM. In example ex2.cpp I would like to compute gradient of 
displacement field u and then stress sigma(u), strain eps(u), strain energy U 
and elementwise strain energy Ue. The functional form of these quantities as a 
function of u is given by

eps(u) = (grad(u) + grad(u)^T)/2

sigma(u) = 2*mu*eps + lambda*trace(eps) I

U = lmbda/2.0*(trace(eps))^2 + mu*trace(eps^2) 

Ue= \int_(Omega)(U/vol *v dOmega) , where vol is cell volume and v are trial 
functions 

Can you provide an example or illustration for this?
Also, how is it possible to apply a point load at specific location on free end 
of the cantilever beam (ex2.cpp)?

Original issue reported on code.google.com by [email protected] on 17 Mar 2014 at 10:36

mfem data directory

How do you generate .mesh file (e.g. beam-tet.mesh) in the data directory of 
the MFEM code? Thanks.

Original issue reported on code.google.com by [email protected] on 13 May 2013 at 9:27

Mesh file format for MFEM

I see that NETGEN with OCC can generate mesh which could be exported in 
different file formats. Can you kindly tell me which file format would be 
useful or do I need to add any extra line like MFEM mesh v1.0.... etc?

Original issue reported on code.google.com by [email protected] on 29 Jul 2013 at 6:04

setting HYPRE_DIR

Shouldn't HYPRE_DIR be set before MPICC ?

# Parallel compiler
MPICC      = mpicxx
MPIOPTS    = $(CCOPTS) -I$(HYPRE_DIR)/include
MPIDEBUG   = $(DEBUG_OPTS) -I$(HYPRE_DIR)/include

# The HYPRE library (needed to build the parallel version)
HYPRE_DIR  = ../../hypre-2.8.0b/src/hypre

Regards,
Jian

Original issue reported on code.google.com by [email protected] on 16 Feb 2012 at 8:14

Ex1 help

For EX1, actually it is Delta u = const which is Poisson solver and not Laplace 
solver as in head comment. Now suppose, if I consider that I have a solid rod 
and I want to set outer boundary temperature is say 1000.0 degrees. I want to 
see distribution radially. 

To do this, shall I change:

ConstantCoefficient Temp(1000.0);

and change b->AddDomainIntegrator in LinearForm and a->AddDomainIntegrator in 
Bilinear to incorporate Temp instead of one?

Thanks and regs

Original issue reported on code.google.com by [email protected] on 30 Jul 2013 at 11:37

Report for Dev release

In OpenSuse 12.3, I tried to run this dev series while I got following error 
which may help you to fix the development issues.
Best,
Abhijit

make
cd lib; g++  -O3 -I/usr/X11R6/include  -DGLVIS_USE_LIBPNG -DGLVIS_USE_FREETYPE 
-I/usr/X11R6/include/freetype2 -I/usr/include/freetype2 -DGLVIS_MULTISAMPLE=4 
-DGLVIS_MS_LINEWIDTH=1.4 -I../../mfem -c aux_gl.cpp
cd lib; g++  -O3 -I/usr/X11R6/include  -DGLVIS_USE_LIBPNG -DGLVIS_USE_FREETYPE 
-I/usr/X11R6/include/freetype2 -I/usr/include/freetype2 -DGLVIS_MULTISAMPLE=4 
-DGLVIS_MS_LINEWIDTH=1.4 -I../../mfem -c aux_vis.cpp
cd lib; gcc  -O3 -I/usr/X11R6/include -I../../mfem -c gl2ps.c
cd lib; g++  -O3 -I/usr/X11R6/include  -DGLVIS_USE_LIBPNG -DGLVIS_USE_FREETYPE 
-I/usr/X11R6/include/freetype2 -I/usr/include/freetype2 -DGLVIS_MULTISAMPLE=4 
-DGLVIS_MS_LINEWIDTH=1.4 -I../../mfem -c material.cpp
cd lib; g++  -O3 -I/usr/X11R6/include  -DGLVIS_USE_LIBPNG -DGLVIS_USE_FREETYPE 
-I/usr/X11R6/include/freetype2 -I/usr/include/freetype2 -DGLVIS_MULTISAMPLE=4 
-DGLVIS_MS_LINEWIDTH=1.4 -I../../mfem -c openglvis.cpp
cd lib; g++  -O3 -I/usr/X11R6/include  -DGLVIS_USE_LIBPNG -DGLVIS_USE_FREETYPE 
-I/usr/X11R6/include/freetype2 -I/usr/include/freetype2 -DGLVIS_MULTISAMPLE=4 
-DGLVIS_MS_LINEWIDTH=1.4 -I../../mfem -c threads.cpp
threads.cpp: In static member function ‘static void* 
communication_thread::execute(void*)’:
threads.cpp:560:73: error: invalid user-defined conversion from 
‘std::basic_istream<char>’ to ‘int’ [-fpermissive]
In file included from /usr/include/c++/4.7/ios:45:0,
                 from /usr/include/c++/4.7/ostream:40,
                 from /usr/include/c++/4.7/iostream:40,
                 from ../../mfem/fem/../general/array.hpp:15,
                 from ../../mfem/fem/fem.hpp:15,
                 from ../../mfem/mfem.hpp:13,
                 from vsdata.hpp:16,
                 from visual.hpp:20,
                 from threads.cpp:16:
/usr/include/c++/4.7/bits/basic_ios.h:113:7: note: candidate is: 
std::basic_ios<_CharT, _Traits>::operator void*() const [with _CharT = char; 
_Traits = std::char_traits<char>] <near match>
/usr/include/c++/4.7/bits/basic_ios.h:113:7: note:   no known conversion for 
implicit ‘this’ parameter from ‘void*’ to ‘int’
threads.cpp:580:72: error: invalid user-defined conversion from ‘std::istream 
{aka std::basic_istream<char>}’ to ‘int’ [-fpermissive]
In file included from /usr/include/c++/4.7/ios:45:0,
                 from /usr/include/c++/4.7/ostream:40,
                 from /usr/include/c++/4.7/iostream:40,
                 from ../../mfem/fem/../general/array.hpp:15,
                 from ../../mfem/fem/fem.hpp:15,
                 from ../../mfem/mfem.hpp:13,
                 from vsdata.hpp:16,
                 from visual.hpp:20,
                 from threads.cpp:16:
/usr/include/c++/4.7/bits/basic_ios.h:113:7: note: candidate is: 
std::basic_ios<_CharT, _Traits>::operator void*() const [with _CharT = char; 
_Traits = std::char_traits<char>] <near match>
/usr/include/c++/4.7/bits/basic_ios.h:113:7: note:   no known conversion for 
implicit ‘this’ parameter from ‘void*’ to ‘int’
make: *** [lib/threads.o] Error 1

Original issue reported on code.google.com by [email protected] on 24 Jul 2013 at 4:22

building library for using LAPACK

Hi,
I need to use methods for DenseMatrix Class, then I have to build (serial) 
library by compiling with LAPACK.

I defined the variable USE_LAPACK       = YES in the makefile, but this setting 
was ignored by the compiler, indeed I see the following compilation line when I 
run make:
cd linalg; g++ -O3  -DMFEM_USE_MEMALLOC      -c densemat.cpp


Finally, my code using DenseMatrixEigensystem doesn't work, and the following 
message is obtained:

DenseMatrixEigensystem::Eval(): Compiled without LAPACK

Could you help me?
Thanks a lot, Pasqua D'Ambra

Original issue reported on code.google.com by [email protected] on 27 Jan 2014 at 9:55

Attachments:

• makefile

wrapper to the hypre eigensolver

Adding a wrapper to the hypre eigensolver and a corresponding modification of 
ex3/ex3p, please? 

Andrew Knyazev

Original issue reported on code.google.com by [email protected] on 9 Aug 2012 at 9:56

ex3.cpp crash when disable/change the mesh refinement

Dear all,

I am trying to use mfem for my problem. I focus on the electromagnetic 
calculation and use ex3.cpp for my test. For debugging and learning the code of 
mfem, I disabled the refinement from the example:

...
   // 2. Refine the mesh to increase the resolution. In this example we do
   //    'ref_levels' of uniform refinement. We choose 'ref_levels' to be the
   //    largest number that gives a final mesh with no more than 50,000
   //    elements.
   if(0){   /// disable refinement for test
      int ref_levels =
         (int)floor(log(50000./mesh->GetNE())/log(2.)/mesh->Dimension());
      for (int l = 0; l < ref_levels; l++)
         mesh->UniformRefinement();
   }
   mesh->ReorientTetMesh();
...


The program crashed. Furthermore, the program still crash even I set a smaller 
number refinement number:

      int ref_levels =
         (int)floor(log(30000./mesh->GetNE())/log(2.)/mesh->Dimension());


Could anyone help me to take a look at it?


Thanks,
Tang Laoya 

Original issue reported on code.google.com by [email protected] on 8 Dec 2014 at 2:18