From: Lorenzo B. <bot...@gm...> - 2010-03-30 09:49:35
|
Hi, I'm sending a debugged (it runs without complaining in dbg mode) version of the dG SIP example and a patch that allows projection of MONOMIALS in case of p-refinement and coarsening. |
From: Roy S. <roy...@ic...> - 2010-03-30 13:34:34
|
On Tue, 30 Mar 2010, Lorenzo Botti wrote: > I'm sending a debugged (it runs without complaining in dbg mode) > version of the dG SIP example and a patch that allows projection of > MONOMIALS in > case of p-refinement and coarsening. Looks like the mailing list blocked the patches? Would you send directly to me? --- Roy |
From: Roy S. <roy...@ic...> - 2010-03-30 17:05:23
|
On Tue, 30 Mar 2010, Lorenzo Botti wrote: > a patch that allows projection of MONOMIALS in case of p-refinement > and coarsening. Would you try my latest commit to system_projection.C? I think it should be enough to get MONOMIAL projection working without setting up a special case code path. --- Roy |
From: Lorenzo B. <bot...@gm...> - 2010-03-30 17:58:58
|
It should work... I'm going to give it a try tomorrow. And it is more consistent with the previous implementation! No wonder, you did the rest of the code ;-). Lorenzo Il giorno 30/mar/2010, alle ore 19.05, Roy Stogner ha scritto: > > On Tue, 30 Mar 2010, Lorenzo Botti wrote: > >> a patch that allows projection of MONOMIALS in case of p-refinement >> and coarsening. > > Would you try my latest commit to system_projection.C? I think it > should be enough to get MONOMIAL projection working without setting > up a special case code path. > --- > Roy |
From: Roy S. <roy...@ic...> - 2010-03-30 18:09:43
|
On Tue, 30 Mar 2010, Lorenzo Botti wrote: > It should work... I'm going to give it a try tomorrow. Thanks! > And it is more consistent with the previous implementation! > No wonder, you did the rest of the code ;-). Yeah, and I'm starting to wonder if I did it right. While I was preparing this patch I think I noticed a bug that would kill the convergence rates for hp refinement on C0 elements in 2D/3D - a situation in which I never did get exponential convergence. No point in attempting a fix now, though, since I won't have time to test it. The app code I'm working on now is LAGRANGE-only and so for us is p-refinement-incompatible; and without that I've got no excuse to spend as much time as it would take to get hp working well by myself. I'd try and conscript Vikram, but he's got a tough enough problem with h refinement alone on some of these non-conforming QoIs. I don't even know if Taylor-Hood elements are LBB-stable under adaptive hp refinement. --- Roy |
From: Lorenzo B. <bot...@gm...> - 2010-03-31 16:41:47
|
With ex21 and some coarsening the error is [0]PETSC ERROR: --------------------- Error Message ------------------------------------ [0]PETSC ERROR: Argument out of range! [0]PETSC ERROR: Negative index -1 at 640 location! [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: Petsc Release Version 2.3.3, Patch 10, Mon Mar 17 20:22:47 CDT 2008 HG revision: 630a491f532632cb3bc52c528f1560e3975eaac6 [0]PETSC ERROR: See docs/changes/index.html for recent updates. [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. [0]PETSC ERROR: See docs/index.html for manual pages. [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: ./ex21-opt on a linux-gnu named digihabia by ihabia Wed Mar 31 18:42:47 2010 [0]PETSC ERROR: Libraries linked from /home/ihabia/src2/petsc-2.3.3-p10//lib/linux-gnu-c-debug [0]PETSC ERROR: Configure run at Thu Mar 27 17:03:05 2008 [0]PETSC ERROR: Configure options --with-cc=gcc --with-fc=gfortran --download-f-blas-lapack=1 --with-mpi-dir=/usr/local/ --with-shared=0 [0]PETSC ERROR: ------------------------------------------------------------------------ [0]PETSC ERROR: VecScatterCheckIndices_Private() line 26 in src/vec/vec/utils/vscat.c [0]PETSC ERROR: VecScatterCreate() line 928 in src/vec/vec/utils/vscat.c [0]PETSC ERROR: User provided function() line 774 in unknowndirectory/src/numerics/petsc_vector.C I'm running in serial. The error appear at line 114 of system_projection.C in the parallel vector implementation old_v.localize(*local_old_vector, projection_list.send_list); compiling with --enable-ghosted the error disappear because we rely on the serial vector implementation Is this correct? Lorenzo ---------- Forwarded message ---------- From: Lorenzo Botti <bot...@gm...> Date: 2010/3/31 Subject: Re: [Libmesh-devel] [Libmesh-users] dg example To: Roy Stogner <roy...@ic...> Cc: Derek Gaston <fri...@gm...>, libmesh-devel < lib...@li...> > It seems that there is a bug enabling coarsening, but I get no errors in dbg mode Lorenzo |
From: Roy S. <roy...@ic...> - 2010-03-31 16:46:46
|
On Wed, 31 Mar 2010, Lorenzo Botti wrote: > With ex21 and some coarsening the error is > > [0]PETSC ERROR: --------------------- Error Message ------------------------------------ > [0]PETSC ERROR: Argument out of range! > [0]PETSC ERROR: Negative index -1 at 640 location! > [0]PETSC ERROR: ------------------------------------------------------------------------ > [0]PETSC ERROR: Petsc Release Version 2.3.3, Patch 10, Mon Mar 17 20:22:47 CDT 2008 HG revision: 630a491f532632cb3bc52c528f1560e3975eaac6 > [0]PETSC ERROR: See docs/changes/index.html for recent updates. > [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting. > [0]PETSC ERROR: See docs/index.html for manual pages. > [0]PETSC ERROR: ------------------------------------------------------------------------ > [0]PETSC ERROR: ./ex21-opt on a linux-gnu named digihabia by ihabia Wed Mar 31 18:42:47 2010 > [0]PETSC ERROR: Libraries linked from /home/ihabia/src2/petsc-2.3.3-p10//lib/linux-gnu-c-debug > [0]PETSC ERROR: Configure run at Thu Mar 27 17:03:05 2008 > [0]PETSC ERROR: Configure options --with-cc=gcc --with-fc=gfortran --download-f-blas-lapack=1 --with-mpi-dir=/usr/local/ --with-shared=0 > [0]PETSC ERROR: ------------------------------------------------------------------------ > [0]PETSC ERROR: VecScatterCheckIndices_Private() line 26 in src/vec/vec/utils/vscat.c > [0]PETSC ERROR: VecScatterCreate() line 928 in src/vec/vec/utils/vscat.c > [0]PETSC ERROR: User provided function() line 774 in unknowndirectory/src/numerics/petsc_vector.C > > I'm running in serial. > The error appear at line 114 of system_projection.C in the parallel vector implementation > > old_v.localize(*local_old_vector, projection_list.send_list); > > compiling with --enable-ghosted the error disappear because we rely on the serial vector implementation > Is this correct? I'm not sure what the problem is at this point. Can you send me an ex21.in that reproduces it? Thanks, --- Roy |
From: Roy S. <roy...@ic...> - 2010-03-31 16:55:19
|
On Wed, 31 Mar 2010, Lorenzo Botti wrote: > the ex21.in in the svn caused the problem on my laptop... What was your exact set of configure options, then? ex21 worked fine in our BuildBot runs (except for with --enable-complex on, and I think I've fixed that). But you may be using a different set of options; I'll try those and see if I can replicate the problem here. --- Roy |