From: John P. <pet...@cf...> - 2007-07-08 21:18:08
|
Ondrej Certik writes: > > DistributedVector (whose existence I'd forgotten about - thanks, John) > > derives from PetscVector directly too, doesn't it? Would it be a > > DistributedVector derives from NumericVector, at least in 0.6.0-rc2. Yes, AFAIK it has always derived from the abstract NumericVector base. > > partial fix if NumericVector::build() created a DistributedVector by > > default instead of just exiting with an error? > > There are actually two places in the system that exit with the solver > error and I tried to return NULL, which obviously segfaulted. :) But > maybe returning the DistributedVector would do the job. I hope it works. This is also a good opportunity (for anyone out there who wants to) to actually *write* a LibMeshSparseMatrix/DistributedMatrix class and LibMeshLinearSolver. Serial is fine at first, we can work on parallel later ;) -John |