dolfin team mailing list archive

[Johan Hake <hake@xxxxxxxxx>]

On Wednesday 18 February 2009 11:00:36 Garth N. Wells wrote:
> Johan Hake wrote:
> > 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.
>
> Are you looking at the sparsity pattern or the assembled matrix? The
> MTL4 wrappers cheat a bit and don't use information from the sparsity
> pattern other than the matrix dimensions.

On the assembled matrix.

Johan

>
> Garth
>
> > 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