Could we please change the PETSc option to allow setting matrix entries that have not been preallocated?

I second this. I've had to make this change locally for quite a while.

It's easy to make the change locally.

Proper fixes for periodic bcs are welcome!

Den Oct 26, 2012 kl. 4:55 PM skrev Garth Wells:

> Your question #212407 on DOLFIN changed:
> https://answers.launchpad.net/dolfin/+question/212407
>
> Status: Open => Answered
>
> Garth Wells proposed the following answer:
> It's easy to make the change locally.

Obviously, but applications may have users that don't compile dolfin themselves.

Hans Petter has a keynote lecture today about computational turbulence (Fenics and cbc.pdesys) at the 25th Nordic seminar on computational mechanics. It is great PR for Fenics, but unfortunate that the applications he talks about don't work with binaries.

>
> Proper fixes for periodic bcs are welcome!

They certainly are. Meanwhile...

Mikael

>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/dolfin/+question/212407/+confirm?answer_id=1
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/dolfin/+question/212407
>
> You received this question notification because you asked the question.

On Fri, Oct 26, 2012 at 03:26:14PM -0000, Mikael Mortensen wrote:
> Question #212407 on DOLFIN changed:
> https://answers.launchpad.net/dolfin/+question/212407
>
> Status: Answered => Open
>
> Mikael Mortensen is still having a problem:
>
> Den Oct 26, 2012 kl. 4:55 PM skrev Garth Wells:
>
> > Your question #212407 on DOLFIN changed:
> > https://answers.launchpad.net/dolfin/+question/212407
> >
> > Status: Open => Answered
> >
> > Garth Wells proposed the following answer:
> > It's easy to make the change locally.
>
> Obviously, but applications may have users that don't compile dolfin
> themselves.
>
> Hans Petter has a keynote lecture today about computational turbulence
> (Fenics and cbc.pdesys) at the 25th Nordic seminar on computational
> mechanics. It is great PR for Fenics, but unfortunate that the
> applications he talks about don't work with binaries.
>
> >
> > Proper fixes for periodic bcs are welcome!
>
> They certainly are. Meanwhile...

Is it possible to call

  MatSetOption(*A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE);

to change the parameter after the matrix has been created?

Then PeriodicBC.cpp could call that function as a temporary fix on the
PETSc Mat pointer.

--
Anders

From my reading of http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatSetOption.html, it seems that Anders' suggestion could be implemented. How would you check that a given GenericMatrix is actually a PETScMatrix?

On Mon, Oct 29, 2012 at 03:35:44PM -0000, Patrick Farrell wrote:
> Question #212407 on DOLFIN changed:
> https://answers.launchpad.net/dolfin/+question/212407
>
> Patrick Farrell posted a new comment:
> >From my reading of http://www.mcs.anl.gov/petsc/petsc-
> current/docs/manualpages/Mat/MatSetOption.html, it seems that Anders'
> suggestion could be implemented. How would you check that a given
> GenericMatrix is actually a PETScMatrix?

if (has_type<Foo>(A))

where Foo is something like PETScMatrix or boost::shared_ptr<PETScMatrix>.

--
Anders

I've proposed a merge of Anders' suggestion at https://code.launchpad.net/~pefarrell/dolfin/periodic-malloc/+merge/132257 .