[deal.II] Boundary support points

2017-06-09 Thread Ignacio Tomas
Hello everybody. I am trying to generate a map (or vectors, whatever is more convenient) that gives me the ''dof coordinate'' (''support point'' in the language of dealii) but only for the boundary dofs. Ideally I am looking for something of the style std::map< types::global_dof_index, Point<

Re: [deal.II] Re: making a package of PETSc vectors

2017-03-24 Thread Ignacio Tomas
I wanted to do something like (in reality I would do this with a for loop) std::vector< std::shared_ptr > VectorPacket ; VectorPacket.resize(3) ; VectorPacket[0].reset(new PETScWrappers::MPI::Vector) ; VectorPacket[1].reset(new PETScWrappers::MPI::Vector) ; VectorPacket[2].reset(new

Re: [deal.II] Re: making a package of PETSc vectors

2017-03-24 Thread Ignacio Tomas
Thanks, I will definitely try it. Ignacio. On Fri, Mar 24, 2017 at 1:50 AM, Jean-Paul Pelteret wrote: > Dear Itomas, > > I would recommend going the route suggested in this recent post >

Re: [deal.II] error when trying to initialize PETSc preconditioners

2016-08-15 Thread Ignacio Tomas
edu> wrote: > Thanks, this should fix it: https://github.com/dealii/dealii/pull/2963 > > On Sat, Aug 13, 2016 at 5:28 PM, Ignacio Tomas <nacho...@gmail.com> wrote: > > This is the whole error message > > > > > > > > > > > > terminate calle

Re: [deal.II] error when trying to initialize PETSc preconditioners

2016-08-13 Thread Ignacio Tomas
This is the whole error message terminate called after throwing an instance of 'dealii::StandardExceptions::ExcInvalidState' what(): An error occurred in line <68> of file in function void

Re: [deal.II] error when trying to initialize PETSc preconditioners

2016-08-12 Thread Ignacio Tomas
ner that works with parallel matrices like > PreconditionBlockJacobi. > > Best, > Timo > > > On Fri, Aug 12, 2016 at 9:46 PM, Ignacio Tomas <nach...@gmail.com > > wrote: > > Hi, everybody. I am just trying to initialize either ICC, Jacobi or > > ParaSa

[deal.II] error when trying to initialize PETSc preconditioners

2016-08-12 Thread Ignacio Tomas
Hi, everybody. I am just trying to initialize either ICC, Jacobi or ParaSalis preconditioner from PETSc, using a PETSc matrix in parallel. So far I am running just one mpi process. After the call to ''.initialize'' I get: [0]PETSC ERROR: No support for this operation for this object type