[petsc-users] Could be a bug on hypre with petsc

[Sun, Hui]

Thank you Barry. I've just reported the problem to Hypre. 

Hui

________________________________________
From: Barry Smith [bsmith at mcs.anl.gov]
Sent: Monday, April 20, 2015 3:23 PM
To: Sun, Hui
Cc: petsc-users at mcs.anl.gov
Subject: Re: [petsc-users] Could be a bug on hypre with petsc

  It could be a bug or it could just be sloppyness where they send uninitialized data that they never use at the other end. Either way you should submit it to the hypre team so they can decide if they want to fix it.

  Barry

> On Apr 19, 2015, at 6:55 PM, Sun, Hui <hus003 at ucsd.edu> wrote:
>
> I'm basically solving Poisson, and using valgrind to check memory issues. I set pctype to PCHYPRE. It works fine if I use one core. But if I use two cores, it gives me the following error messages about the memory, and then it computes results fine. It seems that it has to do with hypre, for if I set pctype to something else such as PCGAMG or PCJACOBI, it does not output any error message at all, regardless of one core or two cores. By the way, I use DM objects for the matrix and ksp.
>
>
>
> ==59097== Syscall param writev(vector[...]) points to uninitialised byte(s)
>
> ==59097==    at 0x2F19FF2: writev (in /usr/lib/system/libsystem_kernel.dylib)
>
> ==59097==    by 0x2AF470D: MPIDU_Sock_writev (sock_immed.i:610)
>
> ==59097==    by 0x2ADE3E5: MPIDI_CH3_iStartMsgv (ch3_istartmsgv.c:110)
>
> ==59097==    by 0x2ACEEF8: MPIDI_CH3_EagerContigSend (ch3u_eager.c:182)
>
> ==59097==    by 0x2AD556F: MPID_Send (mpid_send.c:143)
>
> ==59097==    by 0x2920057: MPI_Send (send.c:127)
>
> ==59097==    by 0x1A712B1: hypre_MPI_Send (mpistubs.c:988)
>
> ==59097==    by 0x1A6E77B: hypre_DataExchangeList (exchange_data.c:392)
>
> ==59097==    by 0x194D7C1: hypre_NewCommPkgCreate_core (new_commpkg.c:267)
>
> ==59097==    by 0x195108B: hypre_MatvecCommPkgCreate (par_csr_communication.c:478)
>
> ==59097==    by 0x192020B: hypre_BoomerAMGCreateS (par_strength.c:190)
>
> ==59097==    by 0x189D6AD: hypre_BoomerAMGSetup (par_amg_setup.c:883)
>
> ==59097==    by 0x187BD7A: HYPRE_BoomerAMGSetup (HYPRE_parcsr_amg.c:52)
>
> ==59097==    by 0x11EFCE6: PCSetUp_HYPRE (hypre.c:157)
>
> ==59097==    by 0x12766D4: PCSetUp (precon.c:918)
>
> ==59097==    by 0x1385670: KSPSetUp (itfunc.c:330)
>
> ==59097==    by 0x138798B: KSPSolve (itfunc.c:542)
>
> ==59097==    by 0x10003C854: fourPoissonSolve (FourPoisson3D.h:383)
>
> ==59097==    by 0x1000DEA42: main (test.c:473)
>
> ==59097==  Address 0x100a295d8 is 8 bytes inside a block of size 28 alloc'd
>
> ==59097==    at 0xC258: malloc (vg_replace_malloc.c:300)
>
> ==59097==    by 0x1A6F501: hypre_MAlloc (hypre_memory.c:62)
>
> ==59097==    by 0x1A6DDBA: hypre_DataExchangeList (exchange_data.c:212)
>
> ==59097==    by 0x194D7C1: hypre_NewCommPkgCreate_core (new_commpkg.c:267)
>
> ==59097==    by 0x195108B: hypre_MatvecCommPkgCreate (par_csr_communication.c:478)
>
> ==59097==    by 0x192020B: hypre_BoomerAMGCreateS (par_strength.c:190)
>
> ==59097==    by 0x189D6AD: hypre_BoomerAMGSetup (par_amg_setup.c:883)
>
> ==59097==    by 0x187BD7A: HYPRE_BoomerAMGSetup (HYPRE_parcsr_amg.c:52)
>
> ==59097==    by 0x11EFCE6: PCSetUp_HYPRE (hypre.c:157)
>
> ==59097==    by 0x12766D4: PCSetUp (precon.c:918)
>
> ==59097==    by 0x1385670: KSPSetUp (itfunc.c:330)
>
> ==59097==    by 0x138798B: KSPSolve (itfunc.c:542)
>
> ==59097==    by 0x10003C854: fourPoissonSolve (FourPoisson3D.h:383)
>
> ==59097==    by 0x1000DEA42: main (test.c:473)
>
>
>
> Best,
>
> Hui
>