Re: [Libmesh-users] Trouble with boundary integration

2017-08-26 Thread Mike Marchywka
at each node if you expect it to be uniform rather than oscillating :) From: Renato Poli <rebp...@gmail.com> Sent: Friday, August 25, 2017 7:46 PM To: Mike Marchywka Cc: John Peterson; libmesh-users Subject: Re: [Libmesh-users] Trouble with bo

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread Mike Marchywka
<rebp...@gmail.com> Sent: Friday, August 25, 2017 6:12 PM To: Mike Marchywka Cc: John Peterson; libmesh-users Subject: Re: [Libmesh-users] Trouble with boundary integration Hi Mike, Thanks for the reply. Any hint on how to easily reproduce the shape function calculations by hand (that

Re: [Libmesh-users] Trouble with boundary integration

2017-08-25 Thread Mike Marchywka
I've been trying to learn this myself and encounter a variety of issues like this. If you have an exact solution, why not just print all the pieces- values and normals etc- and see how they differ point by point instead of trying to do an inverse problem of guessing what is wrong with an

Re: [Libmesh-users] "Unknown convergence flag!"

2017-08-16 Thread Mike Marchywka
From: Jed Brown <j...@jedbrown.org> Sent: Wednesday, August 16, 2017 1:39 PM To: Mike Marchywka; John Peterson Cc: libmesh-users@lists.sourceforge.net Subject: Re: [Libmesh-users] "Unknown convergence flag!" Mike Marchywka <marchy...@hotmail.com&

Re: [Libmesh-users] "Unknown convergence flag!"

2017-08-16 Thread Mike Marchywka
ions later. From: Jed Brown <j...@jedbrown.org> Sent: Monday, August 14, 2017 7:39 PM To: Mike Marchywka; John Peterson Cc: libmesh-users@lists.sourceforge.net Subject: Re: [Libmesh-users] "Unknown convergence flag!" Please, please use KSPConvergedReasons[reason]

Re: [Libmesh-users] "Unknown convergence flag!"

2017-08-14 Thread Mike Marchywka
From: John Peterson <jwpeter...@gmail.com> Sent: Monday, August 14, 2017 2:51 PM To: Mike Marchywka Cc: Roy Stogner; libmesh-users@lists.sourceforge.net Subject: Re: [Libmesh-users] "Unknown convergence flag!" On Mon, Aug 14, 2017 a

Re: [Libmesh-users] "Unknown convergence flag!"

2017-08-14 Thread Mike Marchywka
( can not get hotmail to mark or preserve markings on reply text, comments at bottom ) From: John Peterson <jwpeter...@gmail.com> Sent: Monday, August 14, 2017 12:25 PM To: Mike Marchywka Cc: libmesh-users@lists.sourceforge.net Subject: Re: [Libmesh

[Libmesh-users] "Unknown convergence flag!"

2017-08-14 Thread Mike Marchywka
message without actually stepping through it, https://libmesh.github.io/doxygen/classlibMesh_1_1PetscLinearSolver.html 1681 default : 1682 libMesh::err << "Unknown convergence flag!" << std::endl; 1683 return UNKNOWN_FLAG; note new address Mike Marchywka 306 Ch

Re: [Libmesh-users] Meshrefinment snapping nodes to surface of geometry

2017-06-06 Thread Mike Marchywka
One of the first things I tried to do was make a 2D mesh generator and maintain a specific relationship between nodes and (rectangular lol ) boundary. It looked like nodes on the boundary created some issues and I could deal with special elements at the interfaces. I tried to do this using

[Libmesh-users] dealing with spurious oscillations in xvection diffusion equations

2017-01-15 Thread Mike Marchywka
s not seem to be a lot different from anything else you may add to a system. Changing the mesh and comparing points would be a bit more involved I imagine. Thanks. note new address Mike Marchywka 306 Charles Cox Drive Canton, GA 30115 2295 Collinworth Drive Marietta GA 30062. formerly 487 Salem W

Re: [Libmesh-users] TetGen mesh interface's "triangulate_pointset()" stop running

2016-09-30 Thread Mike Marchywka
(hotmail just switched formats and it looks like they are trying to force html on everyone, I'm top posting since I don't have any idea if the reply text is marked etc ) Is there any indication this is a problem with floating point precision as someone pointed out on another recent thread ?

Re: [Libmesh-users] TetGen mesh interface's "triangulate_pointset()" stop running

2016-09-24 Thread Mike Marchywka
> To: libmesh-users@lists.sourceforge.net > From: thiago...@gmail.com > Date: Sat, 24 Sep 2016 17:19:56 +0200 > Subject: [Libmesh-users] TetGen mesh interface's "triangulate_pointset()" > stop running > > Hello! > > I'm having a weird problem with the

Re: [Libmesh-users] DenseSubMatrix and DenseSubVector lack of ctor and operators

2016-08-22 Thread Mike Marchywka
ts.sourceforge.net > > On Mon, Aug 22, 2016 at 1:08 PM Mike Marchywka > <marchy...@hotmail.com<mailto:marchy...@hotmail.com>> wrote: > I tried that and the compiler did not like it because of the way you > designed the classes... > > Compiler didn't like wha

Re: [Libmesh-users] DenseSubMatrix and DenseSubVector lack of ctor and operators

2016-08-22 Thread Mike Marchywka
> > Now: that's not a reason to be sloppy... But libMesh has endured for > ~15 years at this point and it's anything but sloppy. > > Just some advice... > > Derek > On Mon, Aug 22, 2016 at 10:51 AM Mike Marchywka > <marchy...@hotmail.com<mailto:marchy...

Re: [Libmesh-users] DenseSubMatrix and DenseSubVector lack of ctor and operators

2016-08-22 Thread Mike Marchywka
> > > > On Mon, Aug 22, 2016 at 7:53 AM, Mike Marchywka > <marchy...@hotmail.com<mailto:marchy...@hotmail.com>> wrote: > I was trying to move some code around related to the vars typically > called Ke and Fe > but encountered several errors due to missing c

Re: [Libmesh-users] what is the deal with TypeVector size()?

2016-08-21 Thread Mike Marchywka
DIM; i++) > std::cout<<grad_phi[i][qp](i)<<std::endl; > > Derek > > On Sun, Aug 21, 2016 at 11:51 AM Mike Marchywka > <marchy...@hotmail.com<mailto:marchy...@hotmail.com>> wrote: > > > > > > -

Re: [Libmesh-users] what is the deal with TypeVector size()?

2016-08-21 Thread Mike Marchywka
> Date: Sun, 21 Aug 2016 10:17:41 -0500 > From: royst...@ices.utexas.edu > To: marchy...@hotmail.com > CC: libmesh-users@lists.sourceforge.net > Subject: Re: [Libmesh-users] what is the deal with TypeVector size()? > > > O

[Libmesh-users] what is the deal with TypeVector size()?

2016-08-21 Thread Mike Marchywka
dphi vector would be denormal? Haven't thought about memory map much and IIRC the binary dump did not suggest they were any  memory init pattern to flag that problem.  Thanks. note new address   Mike Marchywka 306 Charles Cox Drive Canton, GA 30115  2295 Collinworth  Drive Marietta GA 30062

Re: [Libmesh-users] Where to start learning Libmesh?

2016-08-19 Thread Mike Marchywka
> Date: Fri, 19 Aug 2016 16:35:07 -0500 > From: royst...@ices.utexas.edu > To: rwood...@utexas.edu > CC: libmesh-users@lists.sourceforge.net > Subject: Re: [Libmesh-users] Where to start learning Libmesh? > > > On Wed, 17 Aug 2016, Ryan Woodall wrote: >

Re: [Libmesh-users] Example files won't run a second time

2016-08-17 Thread Mike Marchywka
> AutoPtr(element_type * p = 0) > ^ > ../../../include/libmesh/auto_ptr.h:185:3: note: no known conversion for > argument 1 from ‘libMesh::UniquePtr<libMesh::FEGenericBase> {aka > std::unique_ptr<libMesh::FEGenericBase, > std::default_delete<libMesh::FEGenericBas

Re: [Libmesh-users] Example files won't run a second time

2016-08-16 Thread Mike Marchywka
] > AutoPtr(element_type * p = 0) > ^ > ../../../include/libmesh/auto_ptr.h:185:3: note: no known conversion for > argument 1 from ‘libMesh::UniquePtr<libMesh::FEGenericBase> {aka > std::unique_ptr<libMesh::FEGenericBase, > std::default_delete<libMesh::FEGenericBase&

[Libmesh-users] FW: Example files won't run a second time

2016-08-12 Thread Mike Marchywka
oontheg...@gmail.com; libmesh-users@lists.sourceforge.net >> Subject: Re: [Libmesh-users] Example files won't run a second time >> >> >> On Fri, 12 Aug 2016, Mike Marchywka wrote: >> >>> Again, the compile seems normal but maybe someone could explain th

Re: [Libmesh-users] Example files won't run a second time

2016-08-12 Thread Mike Marchywka
Again, I would not worry about adding an example but rather just try to build and test standalone. I just reinstalled libmesh adding petsc and other things, repeated the process of running "make -n example-opt" and extracting the compile and link lines and it seems to work ok, Silent failures can

Re: [Libmesh-users] FW: libmesh problem: libmesh_nullptr was not declared in scope

2016-08-09 Thread Mike Marchywka
> From: harroontheg...@gmail.com > Date: Tue, 9 Aug 2016 13:19:07 -0700 > Subject: Re: [Libmesh-users] FW: libmesh problem: libmesh_nullptr was > not declared in scope > To: jwpeter...@gmail.com > CC: marchy...@hotmail.com;

Re: [Libmesh-users] FW: libmesh problem: libmesh_nullptr was not declared in scope

2016-08-09 Thread Mike Marchywka
> -L/usr/lib/gcc/x86_64-linux-gnu/4.8 > -Wl,-rpath,/usr/lib/gcc/x86_64-linux-gnu/4.8 > -Wl,-rpath,/usr/lib/x86_64-linux-gnu > CXXLD example-opt > rm: cannot remove 'example-opt': Permission denied > rm: cannot remove '.libs/example-opt': Permission denied > rm: cannot remove

Re: [Libmesh-users] FW: libmesh problem: libmesh_nullptr was not declared in scope

2016-08-09 Thread Mike Marchywka
obably a good idea.  --- end of comments > > > Harry Pearce > > > > On Mon, Aug 8, 2016 at 11:20 AM, Mike Marchywka > <marchy...@hotmail.com<mailto:marchy...@hotmail.com>> wrote: > > > > > >> From:

Re: [Libmesh-users] FW: libmesh problem: libmesh_nullptr was not declared in scope

2016-08-08 Thread Mike Marchywka
rectory of my project under > SUBDIRS but when I tried to run the makefile the compiler said: > > Making all in eigenproblems/eigenproblems_femviblibtests_ex1 > make[1]: Entering directory > `/opt/libmesh-1.0.0-rc1/examples/eigenproblems/eigenproblems_femviblibtests_ex1' > &g

Re: [Libmesh-users] FW: libmesh problem: libmesh_nullptr was not declared in scope

2016-08-05 Thread Mike Marchywka
m, dof_indices); > > fe->reinit (elem); > > Ke.resize (dof_indices.size(), dof_indices.size()); > > Me.resize (dof_indices.size(), dof_indices.size()); > > */ > > // for (unsigned int qp=0; qp<qrule.n_points(); qp++) > > for