Mark Adams mfad...@lbl.gov writes:
My hash function is:
PetscScalar v = ((PetscScalar)((PETSC_HASH_FACT*idx)%100) -
49.5)/50.0;
ierr = VecSetValues(bb, 1, idx, v,
INSERT_VALUES);CHKERRQ(ierr);
include/petscctable.h:#define PETSC_HASH_FACT 79943
Still grid-aligned
Barry Smith bsm...@mcs.anl.gov writes:
Don't use drand48 directly; use PetscRandom. Matt promised to fix up
PetscRandom to use it by default and work on all systems so all you need to
do is use PetscRandom everywhere.
Yes, that's what I meant. But we need a deterministic seed.
Mark Adams mfad...@lbl.gov writes:
Ah, you're right.
Would you consider this a decent hash?
PetscScalar v = ((PetscScalar)((79943*idx)%6871) - 3435.)/3435. +
((PetscScalar)((7919*idx)%829) - 415.)/415.;
Still bullshit. Plot this and you'll see that there are lots of
patterns and it's
Ah, you're right.
Would you consider this a decent hash?
PetscScalar v = ((PetscScalar)((79943*idx)%6871) - 3435.)/3435. +
((PetscScalar)((7919*idx)%829) - 415.)/415.;
I thought drand48 has been taken off the table because of Windows. If it's
not I can use that.
Mark
On Thu, Aug 27, 2015
On Aug 27, 2015, at 1:58 PM, Jed Brown j...@jedbrown.org wrote:
Mark Adams mfad...@lbl.gov writes:
Ah, you're right.
Would you consider this a decent hash?
PetscScalar v = ((PetscScalar)((79943*idx)%6871) - 3435.)/3435. +
((PetscScalar)((7919*idx)%829) - 415.)/415.;
Still
On Aug 27, 2015, at 2:15 PM, Jed Brown j...@jedbrown.org wrote:
Barry Smith bsm...@mcs.anl.gov writes:
Don't use drand48 directly; use PetscRandom. Matt promised to fix up
PetscRandom to use it by default and work on all systems so all you need to
do is use PetscRandom everywhere.
On Aug 27, 2015, at 4:38 PM, Mark Adams mfad...@lbl.gov wrote:
OK, I can put PetscRandom back into GAMG and add a check for DIVERGE_ITS in
the Cheby solve. Is that necessary and sufficient?
Branch?
Mark
On Thu, Aug 27, 2015 at 4:25 PM, Matthew Knepley knep...@gmail.com wrote:
On Thu, Aug 27, 2015 at 2:17 PM, Barry Smith bsm...@mcs.anl.gov wrote:
On Aug 27, 2015, at 2:15 PM, Jed Brown j...@jedbrown.org wrote:
Barry Smith bsm...@mcs.anl.gov writes:
Don't use drand48 directly; use PetscRandom. Matt promised to fix up
PetscRandom to use it by default and work
OK, I can put PetscRandom back into GAMG and add a check for DIVERGE_ITS in
the Cheby solve. Is that necessary and sufficient?
Mark
On Thu, Aug 27, 2015 at 4:25 PM, Matthew Knepley knep...@gmail.com wrote:
On Thu, Aug 27, 2015 at 2:17 PM, Barry Smith bsm...@mcs.anl.gov wrote:
On Aug 27,
Mark Adams mfad...@lbl.gov writes:
And orthogonal to the high modes of say the Laplacian I assume,
Look at them. They're pretty smooth and grid-aligned (or diagonal)
which is especially bad if you have an anisotropic operator.
signature.asc
Description: PGP signature
Mark Adams mfad...@lbl.gov writes:
OK, I can put PetscRandom back into GAMG and add a check for DIVERGE_ITS in
the Cheby solve. Is that necessary and sufficient?
Yes, but use a deterministic seed (that includes the process rank so all
ranks don't produce the same numbers).
signature.asc
On Aug 27, 2015, at 4:50 PM, Jed Brown j...@jedbrown.org wrote:
Mark Adams mfad...@lbl.gov writes:
OK, I can put PetscRandom back into GAMG and add a check for DIVERGE_ITS in
the Cheby solve. Is that necessary and sufficient?
Yes, but use a deterministic seed (that includes the
Matthew Knepley knep...@gmail.com writes:
I will make them local to the object, and we will all be fine then.
Yes, thanks.
signature.asc
Description: PGP signature
On Aug 27, 2015, at 5:35 PM, Jed Brown j...@jedbrown.org wrote:
Barry Smith bsm...@mcs.anl.gov writes:
This is the default behavior if one just creates a PetscRandom and uses
it, so don't be going and setting seeds inside GAMG or Chebyshev.
Yikes, you are right all three of our
On Wed, Aug 26, 2015 at 3:30 PM, Jed Brown j...@jedbrown.org wrote:
Mark Adams mfad...@lbl.gov writes:
I really think that it does not matter. Just as long as it is noisy.
IMO,
but we have no data, and I'm not sure how we could study it. What is the
chance that your generator will
My hash function is:
PetscScalar v = ((PetscScalar)((PETSC_HASH_FACT*idx)%100) -
49.5)/50.0;
ierr = VecSetValues(bb, 1, idx, v,
INSERT_VALUES);CHKERRQ(ierr);
include/petscctable.h:#define PETSC_HASH_FACT 79943
On Thu, Aug 27, 2015 at 7:53 AM, Mark Adams mfad...@lbl.gov
16 matches
Mail list logo