← Back to team overview

dolfin team mailing list archive

Re: [HG DOLFIN] Fix bug in numpy wrapping of CRS data.

 

On Wednesday 18 February 2009 10:23:52 Garth N. Wells wrote:
> Johan Hake wrote:
> > On Wednesday 18 February 2009 09:38:05 DOLFIN wrote:
> >> One or more new changesets pushed to the primary dolfin repository.
> >> A short summary of the last three changesets is included below.
> >>
> >> changeset:   5728:c01d60f80a759d8126dd67471db65d9b733ff1ee
> >> tag:         tip
> >> user:        "Johan Hake <hake@xxxxxxxxx>"
> >> date:        Wed Feb 18 09:37:59 2009 +0100
> >> files:       dolfin/swig/dolfin_la_post.i test/unit/la/python/test.py
> >> description:
> >> Fix bug in numpy wrapping of CRS data.
> >>   - Numpy type needed to be NPY_ULONG instead of NPY_UINT(!?)
> >
> > I noticed a qurious thing while fixing this bug. The sparsity pattern are
> > not the same at a 64 bit computer and on my laptop, a 32 bit one.
> >
> > Any clue why?
>
> Strange. Can you add some simple code to the sandbox to which reproduces
>   this?

Done.

I have checked some more and it turns out that it only applies for the MTL4 
backend. It removes the nonzero values from the pattern on my laptop.

I also have different versions of MTL4 on the two machines:

  64 bit: svn revision 6624
  32 bit: svn revision 6597

Could be this.

Johan


Follow ups

References