Re: [Libmesh-users] Adaptative Mesh Refinement

2018-07-15 Thread Derek Gaston
MOOSE has had this for a long time - you can see how we do it here: https://github.com/idaholab/moose/blob/devel/framework/src/problems/FEProblemBase.C#L4651 Basically: 1. Indicate/Mark 2. Adapt 3. Reinit and re-project 4. Repeat Derek On Sun, Jul 15, 2018 at 1:55 PM Vinicius C. Reis wrot

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-07-15 Thread Vinicius C. Reis
Hi Roy, my problem is pretty simple so the workaround you've suggested will work perfectly fine. For now, since I'm testing only simple geometries it might suffice to do all refinement before solving, so no post solve projection would even be necessary. If I have the chance to implement any complex

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-29 Thread Roy Stogner
On Thu, 28 Jun 2018, Roy Stogner wrote: On Tue, 12 Jun 2018, Vinicius C. Reis wrote: Hi John, I pasted a somewhat smaller version of my code below. By playing with the refinement steps limit and the commented out init() and reinit() lines you should be able to reproduce the two exceptions

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-28 Thread Roy Stogner
On Tue, 12 Jun 2018, Vinicius C. Reis wrote: Hi John, I pasted a somewhat smaller version of my code below. By playing with the refinement steps limit and the commented out init() and reinit() lines you should be able to reproduce the two exceptions I got. Got it reproduced. The exceptions

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-20 Thread Roy Stogner
On Tue, 12 Jun 2018, Vinicius C. Reis wrote: Hi John, I pasted a somewhat smaller version of my code below. By playing with the refinement steps limit and the commented out init() and reinit() lines you should be able to reproduce the two exceptions I got. Could you email me a copy of the m

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-11 Thread Vinicius C. Reis
Hi John, I pasted a somewhat smaller version of my code below. By playing with the refinement steps limit and the commented out init() and reinit() lines you should be able to reproduce the two exceptions I got. I hope this clarifies. Thanks for your time! //---

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-11 Thread John Peterson
On Sun, Jun 10, 2018 at 10:02 PM, Vinicius C. Reis wrote: > Hi again, it sounds possible that the issue is the hanging nodes. if I only > do one refinement step no error occurs. Yes, but presumably there are still hanging nodes even after one level of refinement is performed, unless you are ref

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-10 Thread Vinicius C. Reis
Hi again, it sounds possible that the issue is the hanging nodes. if I only do one refinement step no error occurs. Adding another step, same exception. Since I am hard flagging the elements to be refined, quite a few hanging nodes appear already in the first refinement step. I tried calling init()

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-05 Thread John Peterson
On Tue, Jun 5, 2018 at 1:04 PM, Vinicius C. Reis wrote: > Sorry about that, I should have imagined... Well, exception and stack trace > from a fresh run in the text dump below. I managed to reproduce it in a > quite smaller code, what would be the best way to share, paste it into an > email or th

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-05 Thread Vinicius C. Reis
Sorry about that, I should have imagined... Well, exception and stack trace from a fresh run in the text dump below. I managed to reproduce it in a quite smaller code, what would be the best way to share, paste it into an email or through git repository? Thanks, VCR --

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-05 Thread Roy Stogner
On Mon, 4 Jun 2018, Vinicius C. Reis wrote: The exception is thrown from within EquationSystems.init(). I'll do some cleanup and get back with a somewhat more concise reproduction. I've attached two text files with the exception and the stack trace. I'm afraid the attachments got stripped by

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-04 Thread Vinicius C. Reis
The exception is thrown from within EquationSystems.init(). I'll do some cleanup and get back with a somewhat more concise reproduction. I've attached two text files with the exception and the stack trace. Thanks for the help, VCR 2018-06-04 15:20 GMT-03:00 John Peterson : > > > On Mon, Jun 4,

Re: [Libmesh-users] Adaptative Mesh Refinement

2018-06-04 Thread John Peterson
On Mon, Jun 4, 2018 at 11:44 AM, Vinicius C. Reis wrote: > Hi, > > I'm trying to perform a h-refinement and a subdomain > evaluation/reevaluation for each refinement iteration, prior to > initializing an EquationSystems object. Currently when call init() I get a > out_of_range exception. I suspec

[Libmesh-users] Adaptative Mesh Refinement

2018-06-04 Thread Vinicius C. Reis
Hi, I'm trying to perform a h-refinement and a subdomain evaluation/reevaluation for each refinement iteration, prior to initializing an EquationSystems object. Currently when call init() I get a out_of_range exception. I suspect it is thrown from within its node or the element loop. Everything wo