Re: [Libmesh-users] writing/reading data in Reduced Basis

2016-05-31 Thread Lorenzo Zanon
On 5/31/16 5:41 PM, John Peterson wrote: > > > On Tue, May 31, 2016 at 3:23 AM, Lorenzo Zanon > mailto:za...@aices.rwth-aachen.de>> wrote: > > Hello, > > I solved the problem. During the "online" run, I was writing down the > coefficients

Re: [Libmesh-users] writing/reading data in Reduced Basis

2016-05-31 Thread Lorenzo Zanon
; Would I get better results using the new RBDataSerialization? > > > My problem also involves EIM, but the EIM data are read-in also > before the offline RB-Greedy, so I don't think this is relevant. > > > Thanks a lot > > Lorenzo > > > -- > =

Re: [Libmesh-users] crash in RBEIMConstruction

2016-05-10 Thread Lorenzo Zanon
Hello, Calling empty() on the DenseVector is causing the segfault. If I copy the DenseVector into a std::vector, then it's fine. How can I see if there is something wrong with my build/compiler? Everything runs fine apart from that error. The libMesh version is 0.9.5 (according to the first li

[Libmesh-users] crash in RBEIMConstruction

2016-05-10 Thread Lorenzo Zanon
ntil exit from function _ZTVN10__cxxabiv120__si_class_type_infoE, which has no line number information. PetscSignalHandler_Private (sig=0) at /rwthfs/rz/cluster/home/lz474069/Desktop/libmesh_2016/petsc/src/sys/error/signal.c:42 42{ (gdb) -- ======== Lorenzo Zanon Graduate S

[Libmesh-users] Scalar warning

2016-03-19 Thread Lorenzo Zanon
Hello, I just installed the most recent version of libMesh. I'm getting a series of warning when compiling in dbg mode, all of the form: In file included from /home/lz474069/Desktop/libmesh_2016/libmesh/installed/include/libmesh/point.h(24), from /home/lz474069/Desktop/libmes

Re: [Libmesh-users] Scalar warning

2016-03-19 Thread Lorenzo Zanon
>>> ... Thanks! Lorenzo On 3/17/16 3:25 PM, Roy Stogner wrote: > > On Thu, 17 Mar 2016, Lorenzo Zanon wrote: > >> /home/lz474069/Desktop/libmesh_2016/libmesh/installed/include/libmesh/type_vector.h(82): >> >> >> warning #488: template parameter "

Re: [Libmesh-users] RB inclusions

2014-10-31 Thread Lorenzo Zanon
Yes, eim_classes.h, rb_classes.h, and assembly.h all start and end with something like #ifndef __eim_classes_h__ #define __eim_classes_h__ ... #endif There are no other .h files... On 10/31/14 4:51 PM, John Peterson wrote: > On Fri, Oct 31, 2014 at 9:03 AM, Lorenzo Zanon > wrote: &g

[Libmesh-users] RB inclusions

2014-10-31 Thread Lorenzo Zanon
sis_ex4.C:(.text+0xfe50): multiple definition of `SimpleEIMEvaluation::initialize_eim_theta_objects()' rb_classes.x86_64-unknown-linux-gnu.opt.o:rb_classes.C:(.text+0xf790): first defined here for all the SimpleEIM and SimpleRB functions. How could I fix that? Thanks! Lorenzo --

Re: [Libmesh-users] Exodus not working

2014-06-05 Thread Lorenzo Zanon
it becomes compatible with netcdf-C 4.3.2? What about netCDF-Fortran? I don't have to install it myself for using libMesh, do I? Thanks! Lorenzo > > On Mon, May 19, 2014 at 5:35 AM, Lorenzo Zanon > mailto:za...@aices.rwth-aachen.de>> wrote: > > Hi, > >

Re: [Libmesh-users] Exodus not working

2014-05-19 Thread Lorenzo Zanon
-k3 -b -o c0_4.nc $srcdir/c0.cdl FAIL: run_nc4_tests.sh 2 of 2 tests failed Please report to support-net...@unidata.ucar.edu Thanks! Lorenzo On 5/13/14 6:28 PM, John Peterson wrote: > On Tue, May 13, 2014 at 9:40 AM, Lorenzo Zano

Re: [Libmesh-users] Exodus not working

2014-05-13 Thread Lorenzo Zanon
picxx-openmpi-gcc48 CC=mpicc-openmpi-gcc48 FC=mpif90-openmpi-gcc48 F77=mpif90-openmpi-gcc48 \ CXXFLAGS=-I/usr/X11/include LDFLAGS=-L/usr/X11/lib \ Thanks! On 5/13/14 5:29 PM, John Peterson wrote: > > > > On Tue, May 13, 2014 at 9:16 AM, Lorenzo Zanon > mailto:za...@aices.rw

Re: [Libmesh-users] Exodus not working

2014-05-13 Thread Lorenzo Zanon
() #1 0x0061 in ?? () #2 0x000102c19218 in __func__.14673 () I also didn't do anything involving std::complex... On 5/13/14 4:54 PM, John Peterson wrote: > > > > On Tue, May 13, 2014 at 8:30 AM, Lorenzo Zanon > mailto:za...@aices.rwth-aachen.de

[Libmesh-users] Exodus not working

2014-05-13 Thread Lorenzo Zanon
Hi, Besides, since I updated to version 0.9.3, I cannot use Exodus any more. The program crashes at the line e.g. ExodusII_IO(mesh).write_equation_systems ("out_3D.e", equation_systems); Maybe I missed something in the configuration, although I didn't change anything w.r.t. last time. How can

Re: [Libmesh-users] Libmesh with Arpack

2014-05-13 Thread Lorenzo Zanon
few seconds. But shouldn't the solver be the same? Thanks for helping Lorenzo On 4/30/14 9:15 PM, John Peterson wrote: On Wed, Apr 30, 2014 at 12:04 PM, Lorenzo Zanon mailto:za...@aices.rwth-aachen.de>> wrote: Hello, Now I should have slepc correctly installed with ar

[Libmesh-users] Libmesh with Arpack

2014-04-30 Thread Lorenzo Zanon
Hello, Now I should have slepc correctly installed with arpack. But still it looks like I cannot use arpack in libmesh, since in the file slepc_eigen_solver.C the arpack lines are commented out: // case ARPACK: // ierr = EPSSetType (_eps, (char*) EPSARPACK); LIBMESH_CHKERRABORT(ierr); return;

[Libmesh-users] eigensolver

2014-04-10 Thread Lorenzo Zanon
Hello, What do you think is the best SLEPC eigensolver_type for getting the smallest (magnitude) eigenvalue and -mode? The eigenvalues in the spectrum are all real positive. So far I've always used LAPACK, but switching to a 3D problem, that becomes too slow and expensive… Thanks! Lorenzo

Re: [Libmesh-users] installation-petsc problem on iMac

2014-01-22 Thread Lorenzo Zanon
On Jan 20, 2014, at 6:34 PM, John Peterson wrote: > > > > On Mon, Jan 20, 2014 at 9:34 AM, Lorenzo Zanon > wrote: > Hello, > > Our IT team developed a new procedure for installing libmesh on our iMac > machines, so that petsc, slepc, glpk, and eigen are calle

[Libmesh-users] installation-petsc problem on iMac

2014-01-20 Thread Lorenzo Zanon
Hello, Our IT team developed a new procedure for installing libmesh on our iMac machines, so that petsc, slepc, glpk, and eigen are called from the macports that are controlled by them. The procedure - here below - works fine, I think the girl called Vera wrote here to come up with it. There's

Re: [Libmesh-users] Newton iteration nonlinear elasticity slow

2013-12-10 Thread Lorenzo Zanon
0001, absolute=1e-50, divergence=1 > left preconditioning > using PRECONDITIONED norm type for convergence test > PC Object: 1 MPI processes > type: ilu > ILU: out-of-place factorization > 0 levels of fill > tolerance for zero pivo

Re: [Libmesh-users] Newton iteration nonlinear elasticity slow

2013-12-09 Thread Lorenzo Zanon
Hello, I'm sorry it took so long. There were actually a couple of mistakes in my jacobian and residual. Now for the small problem without specifying any options I get: > Running ./example-opt > > Mesh Information: > mesh_dimension()=3 > spatial_dimension()=3 > n_nodes()=66 > n_local

Re: [Libmesh-users] Newton iteration nonlinear elasticity slow

2013-11-25 Thread Lorenzo Zanon
do without LU, but still it's very slow. Thanks Lorenzo On Nov 25, 2013, at 4:01 PM, Kirk, Benjamin (JSC-EG311) wrote: > > On Nov 25, 2013, at 6:23 AM, Lorenzo Zanon > wrote: > >> Hello, >> >> Thanks for the detailed answer! >> >> I've j

Re: [Libmesh-users] Newton iteration nonlinear elasticity slow

2013-11-25 Thread Lorenzo Zanon
terations you're taking and then make sure that both your > libMesh and FEAP implementations are taking the same (or VERY similar) number > of both nonlinear _and_ linear iterations > > > Derek > > > > On Fri, Nov 22, 2013 at 10:08 AM, Lorenzo Zanon > wro

[Libmesh-users] Newton iteration nonlinear elasticity slow

2013-11-22 Thread Lorenzo Zanon
Hello, I have an example based on miscellaneous_ex3.C, where I implemented a nonlinear elastic problem with StVenant stress-strain law on a 3D domain 5x1x1, blocked at x=0 and with an applied load at x=5. The problem is, the computation of the displacement (along x y and z) is really slow (hou

Re: [Libmesh-users] libMesh on Eclipse

2013-09-04 Thread Lorenzo Zanon
On Sep 4, 2013, at 5:31 PM, John Peterson wrote: > > > > On Wed, Sep 4, 2013 at 5:17 AM, Lorenzo Zanon > wrote: > Hello, > > We are having some troubles with using libMesh on iMac workstations. The > problem is, the Macports are regularly updated by our IT team

[Libmesh-users] libMesh on Eclipse

2013-09-04 Thread Lorenzo Zanon
Hello, We are having some troubles with using libMesh on iMac workstations. The problem is, the Macports are regularly updated by our IT team, and this could lead to missing libraries and other issues from time to time. We had a procedure to have libMesh installed on iMac and working on Xcode,

Re: [Libmesh-users] makefile

2013-08-13 Thread Lorenzo Zanon
tting a compile error: ex23aa.C(29): catastrophic error: cannot open source file "libmesh.h" What's the right way to link then? Thanks Lorenzo On Aug 13, 2013, at 1:28 PM, Lorenzo Zanon wrote: > Hi, > > I re-installed libmesh as you indicated. I didn't use --with-m

Re: [Libmesh-users] makefile

2013-08-13 Thread Lorenzo Zanon
lled in my code directory? Thanks Lorenzo On Aug 12, 2013, at 6:03 PM, John Peterson wrote: > > > > On Mon, Aug 12, 2013 at 9:57 AM, Lorenzo Zanon > wrote: > Hello, > > To be more specific, my code now is based on libmesh version 0.7.3. I would > like to run it

Re: [Libmesh-users] makefile

2013-08-12 Thread Lorenzo Zanon
, how do I link this new version of libmesh to my .C and .h files? Thanks for helping Lorenzo On Aug 9, 2013, at 4:17 PM, Lorenzo Zanon wrote: > Hello, > > I just re-installed the newest version of Libmesh. All the examples seem > running, but the Make.common file has not been created

[Libmesh-users] makefile

2013-08-09 Thread Lorenzo Zanon
Hello, I just re-installed the newest version of Libmesh. All the examples seem running, but the Make.common file has not been created in the main directory, so I cannot link my own code to this newly installed version. Does that mean that there's something wrong with the configuration? Thanks

[Libmesh-users] nonlinear solver

2013-07-11 Thread Lorenzo Zanon
=3; equation_systems.get_system("StVen").solve(); But I think here the initial guess is always the 0 function! How can I change it to the solution at each step? Thanks a lot, Lorenzo Lorenzo Zanon, M.Sc. Doctoral Candidate Graduate School AICES

Re: [Libmesh-users] Newton iteration

2013-03-08 Thread Lorenzo Zanon
Thanks for the many answers! Still, according to the PETSc manual "matrix-free methods cannot be used if a direct solver is to be employed", I think this means I cannot use -snes_mf_operator and "-ksp_type preonly -pc_type lu -pc_factor_nonzeros_along_diagonal" at the same time. Do I need some

Re: [Libmesh-users] Newton iteration

2013-03-07 Thread Lorenzo Zanon
ield to > plastic flow, concentration inhomogeneities diffuse to uniformity. Systems > tend to subside to very simple states,independent of their specific history." > Herbert Callen . Thermodynamics and an Introduction to Thermostatics. > > - Original Message - > From: &q

Re: [Libmesh-users] Newton iteration

2013-03-06 Thread Lorenzo Zanon
damp out, internal strains yield to > plastic flow, concentration inhomogeneities diffuse to uniformity. Systems > tend to subside to very simple states,independent of their specific history." > Herbert Callen . Thermodynamics and an Introduction to Thermostatics. > > - Original

[Libmesh-users] Newton iteration

2013-03-06 Thread Lorenzo Zanon
Hello, I'm modifying miscellaneous_ex3.C (2D Laplace-Young Problem Using Nonlinear Solvers, in my current libmesh version). After adapting residual and jacobian to my own problem, although the tolerance is not met, the Newton solver stops after one iteration with an all-zero (incorrect) solutio