On Wed, Mar 24, 2010 at 3:01 PM, (Rebecca) Xuefei YUAN <xy2102 at columbia.edu>wrote:
> Dear Jed and Matt, > This is a genuine bug. It is in DMMGComputeJacobianWithFD() and only occurs when IS_COLORING_GHOSTED is true. So, using IS_COLORING_GLOBAL would probably work here (can't look up the option right now). Barry, what is supposed to happen here? Clearly a local vector is being passed where a global vector is expected (at least part of the time). Matt > Yes, X and F are global and localFIELD is serial. I ran with np=2. > > The error I get is: > > > ************************************************** > 0 SNES Function norm 1.095445115010e+01 > [1]PETSC ERROR: VecScatterBegin() line 1538 in src/vec/vec/utils/vscat.c > Vector wrong size 100 for scatter 60 (scatter forward and vector from != ctx > from size) > [0]PETSC ERROR: VecScatterBegin() line 1538 in src/vec/vec/utils/vscat.c > Vector wrong size 100 for scatter 60 (scatter forward and vector from != ctx > from size) > > > from gdb, I get: > > > Program received signal SIGABRT, Aborted. > [Switching to Thread 0xb7c396b0 (LWP 6301)] > 0xb7f09410 in __kernel_vsyscall () > (gdb) where > #0 0xb7f09410 in __kernel_vsyscall () > #1 0xb7c89085 in raise () from /lib/tls/i686/cmov/libc.so.6 > #2 0xb7c8aa01 in abort () from /lib/tls/i686/cmov/libc.so.6 > #3 0x087d591d in PetscAbortErrorHandler (line=1538, > fun=0x88fe26d "VecScatterBegin", file=0x88fd9f7 "vscat.c", > dir=0x88fd9ff "src/vec/vec/utils/", n=60, p=1, > mess=0xbfd31124 "Vector wrong size 100 for scatter 60 (scatter forward > and v > ector from != ctx from size)", ctx=0x0) at errabort.c:62 > #4 0x0874a88e in PetscError (line=1538, func=0x88fe26d "VecScatterBegin", > file=0x88fd9f7 "vscat.c", dir=0x88fd9ff "src/vec/vec/utils/", n=60, p=1, > mess=0x88fe3a0 "Vector wrong size %D for scatter %D (scatter forward and > vec > tor from != ctx from size)") at err.c:482 > #5 0x086558cc in VecScatterBegin (inctx=0x8a393b0, x=0x8b39cc0, > y=0x8b3b310, > addv=INSERT_VALUES, mode=SCATTER_FORWARD) at vscat.c:1538 > #6 0x0828b54c in DAGlobalToLocalBegin (da=0x8a2d360, g=0x8b39cc0, > mode=INSERT_VALUES, l=0x8b3b310) at dagtol.c:50 > #7 0x080f1b25 in FormFunction (snes=0x8a47840, X=0x8b39cc0, F=0x8adb2f0, > dummg=0x8a2b810) at twgcqt2unffnictv.c:8382 > #8 0x0860c63d in MatFDColoringApply_AIJ (J=0x8a6bbb0, coloring=0x8aa7b40, > x1=0x8b12a60, flag=0xbfd32164, sctx=0x8a47840) at fdmatrix.c:680 > #9 0x0860abaf in MatFDColoringApply (J=0x8a6bbb0, coloring=0x8aa7b40, > x1=0x8b12a60, flag=0xbfd32164, sctx=0x8a47840) at fdmatrix.c:521 > #10 0x08122e45 in SNESDefaultComputeJacobianColor (snes=0x8a47840, > ---Type <return> to continue, or q <return> to quit--- > x1=0x8b12a60, J=0x8a47910, B=0x8a47914, flag=0xbfd32164, ctx=0x8aa7b40) > at snesj2.c:49 > #11 0x0811c7cf in DMMGComputeJacobianWithFD (snes=0x8a47840, x1=0x8a39c20, > J=0x8a47910, B=0x8a47914, flag=0xbfd32164, ctx=0x8a2b810) at > damgsnes.c:365 > #12 0x0811a883 in DMMGComputeJacobian_Multigrid (snes=0x8a47840, > X=0x8a39c20, > J=0x8a47910, B=0x8a47914, flag=0xbfd32164, ptr=0x8a2b6b0) at > damgsnes.c:60 > #13 0x080fc610 in SNESComputeJacobian (snes=0x8a47840, X=0x8a39c20, > A=0x8a47910, B=0x8a47914, flg=0xbfd32164) at snes.c:1188 > #14 0x08124471 in SNESSolve_LS (snes=0x8a47840) at ls.c:189 > #15 0x08103fac in SNESSolve (snes=0x8a47840, b=0x0, x=0x8a39c20) at > snes.c:2242 > #16 0x0811db31 in DMMGSolveSNES (dmmg=0x8a2b6b0, level=0) at damgsnes.c:510 > #17 0x08117029 in DMMGSolve (dmmg=0x8a2b6b0) at damg.c:313 > #18 0x08052ecc in Solve (dmmg=0x8a2b6b0) at twgcqt2unffnictv.c:679 > #19 0x0804def5 in main (argc=Cannot access memory at address 0x189d > ) at twgcqt2unffnictv.c:303 > > > Same things happened to ex25.c from > > > http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/src/snes/examples/tutorials/ex25.c.html > > > with > > ierr = DMMGSetSNES(dmmg,FormFunction,0);CHKERRQ(ierr); > // ierr = > DMMGSetSNESLocal(dmmg,FormFunctionLocal,0,ad_FormFunctionLocal,0);CHKERRQ(ierr); > > > Thanks a lot! > > Rebecca > > > > > Quoting Jed Brown <jed at 59A2.org>: > > On Wed, 24 Mar 2010 15:41:38 -0400, "(Rebecca) Xuefei YUAN" < >> xy2102 at columbia.edu> wrote: >> >>> Hi,all, >>> >>> I have an error from >>> >>> ierr = DAGlobalToLocalBegin(dafield,X,INSERT_VALUES,localFIELD); >>> CHKERRQ(ierr); >>> >>> where >>> >>> ierr = DMMGSetSNES(dmmg, FormFunction,0);CHKERRQ(ierr); is used for >>> set up the SNES. >>> >>> So I check up the vector size of X, F, localFIELD where >>> ierr = DAGetLocalVector(dafield,&localFIELD);CHKERRQ(ierr); >>> PetscInt nlocalFIELD,nX,nF; >>> ierr = VecGetSize(localFIELD,&nlocalFIELD);CHKERRQ(ierr); >>> ierr = VecGetSize(X,&nX);CHKERRQ(ierr); >>> ierr = VecGetSize(F,&nF);CHKERRQ(ierr); >>> ierr = DAGlobalToLocalBegin(dafield,X,INSERT_VALUES,localFIELD); >>> CHKERRQ(ierr); >>> ierr = >>> DAGlobalToLocalEnd(dafield,X,INSERT_VALUES,localFIELD);CHKERRQ(ierr); >>> ierr = >>> DAVecGetArray(dafield,localFIELD,&localfield);CHKERRQ(ierr); >>> >>> >>> (gdb) disp nX >>> 1: nX = 120 >>> (gdb) disp nF >>> 2: nF = 120 >>> (gdb) disp nlocalFIELD >>> 3: nlocalFIELD = 100 >>> >> >> Is this run in parallel? Note that the sizes of X and F are global, >> while localFIELD is serial. What error did you get? >> >> >> Matt, it's clearly FormFunction and not FormFunctionLocal because the >> function prototype has the SNES. >> >> Jed >> >> >> > > > -- > (Rebecca) Xuefei YUAN > Department of Applied Physics and Applied Mathematics > Columbia University > Tel:917-399-8032 > www.columbia.edu/~xy2102 > > -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100324/c44d6064/attachment.htm>
