I have just tried, and the petsc is now also working. Thanks again. Now I will see if I can make some sense of the rest of it :)
Evan On Wed, Feb 11, 2009 at 9:41 PM, Evan Lezar <[email protected]> wrote: > I will pull the latest source and try again. > > Thanks > Evan > > > On Wed, Feb 11, 2009 at 9:23 PM, <[email protected]> wrote: > >> > On Mon, Feb 9, 2009 at 1:05 PM, Anders Logg <[email protected]> wrote: >> > >> >> On Mon, Feb 09, 2009 at 11:58:28AM +0100, Kent Andre wrote: >> >> > >> >> > Is it the code in DofMap that is not understandable ? >> >> > >> >> > The idea behind the code is that the dof_map provides a numbering >> >> > of the unknowns ie. the global dofs are D = [0, 1, .....N-1]. >> >> > In the restricted function space we need a corresponding numbering >> >> > R = [0, 1, .... M-1] where M < N. >> >> > >> >> > The code in restriction is such that it runs through all dofs in the >> >> > restricted mesh (as it is now by runing through the whole mesh except >> >> > for the ones marked by the meshfunction). I run through all the dofs >> >> > and for each dof I check the dof in the global mesh, If I already >> have >> >> > seen this dof I use the same dof_number as last time (stored in >> >> > R[dofs[k]]). If the dof is new I update the dof_counter and add >> >> > the new dof to R. By doing this I get a new dof_map where each >> >> > node is numbered in an increasing sequence. >> >> > >> >> > Check also the demo under dolfin/demo/function/restriction. You can >> >> eg. >> >> > probably begin by changing demo. >> >> > >> >> > >> >> > Kent >> >> >> >> The demo looks a little unfinished. Could you move it to >> >> >> >> dolfin/demo/function/restriction/python/demo.py >> >> >> >> then add the copyright headers etc like in the other demos, clean it >> >> up a little, add a C++ version in the cpp directory (or a copy of the >> >> README file from other demos where the C++ version is missing). You >> >> can also move it to the sandbox until it's finished. >> >> >> >> -- >> >> Anders >> >> >> >> >> > Regarding the demo. >> > >> > I have just tried to run it, and it fails with Petsc reporting a >> > Segmentation Violation. >> > >> > Evan >> > >> >> I fixed a bug yesterday, do you have the newest dev version ? >> (I have only tested the trilinos and uBlas version.) >> >> Kent >> >> >> >> >
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
