Re: [Libmesh-devel] VTK output fails in parallel

2011-02-21 Thread Roman Vetter
> EquationSystems::allgather() is actually a second step up in > sophistication from where VTK output is - it gathers geometry data and > degree of freedom numbering in the case of a distributed mesh. The > gathering (and sadly, the decimation of higher-order solutions) of > vector data happens in

Re: [Libmesh-devel] VTK output fails in parallel

2011-02-21 Thread Roy Stogner
On Mon, 21 Feb 2011, Roman Vetter wrote: >> EquationSystems::allgather() is actually a second step up in >> sophistication from where VTK output is - it gathers geometry data and >> degree of freedom numbering in the case of a distributed mesh. The >> gathering (and sadly, the decimation of high

Re: [Libmesh-devel] VTK output fails in parallel

2011-02-21 Thread Roy Stogner
On Mon, 21 Feb 2011, Roy Stogner wrote: > Patch didn't apply cleanly for me (submit as an attachment next time? > my email system must have mangled the whitespace or something) but it > looks simple enough to apply by hand. And I'm glad I looked at it more closely - shouldn't we also be using v

[Libmesh-devel] A bug in Epetra_Matrix::add_matrix ?

2011-02-21 Thread Ming Q
Hi, I don't know it is a bug or not. But in my computer the SumIntoGlobalValues function took the default COLUMN_MAJOR value for the the last parameter. Therefore line 317 in file trilinos_epetra_matrix.C: 00317 _mat->SumIntoGlobalValues(m, (int *)&rows[0], n, (int *)&cols[0], &dm.get_values

Re: [Libmesh-devel] A bug in Epetra_Matrix::add_matrix ?

2011-02-21 Thread Derek Gaston
Ha! Thanks! I think you'll find that tons of stuff is wrong / incomplete with the Trilinos support. I did a lot of that initial wok years ago... and it hasn't been used much since then. Please let us know about anything you find like this and I'll merge it back into trunk. Derek Sent from my