Re: [Libmesh-users] question about the numeric vector

2018-01-25 Thread Roy Stogner
On Thu, 25 Jan 2018, Michael Povolotskyi wrote: I added a plain System class, and it worked. Thanks! That was fast! (possibly causality-violatingly so? I could swear I saw this message header in my inbox *before* I hit send on my own previous message...) --- Roy Michael. On 01/25/2018

Re: [Libmesh-users] question about the numeric vector

2018-01-25 Thread Roy Stogner
On Thu, 25 Jan 2018, Michael Povolotskyi wrote: Can I add just a plain System to the EquationSystem? That's actually a good question, and if you'd try it and report back I'd appreciate it. The answer used to be "no", which is why everybody habitually recommends adding ExplicitSystem for stor

Re: [Libmesh-users] question about the numeric vector

2018-01-25 Thread Michael Povolotskyi
I added a plain System class, and it worked. Michael. On 01/25/2018 12:17 PM, Roy Stogner wrote: On Thu, 25 Jan 2018, Michael Povolotskyi wrote: I will add another system and will not solve it. Make sure it's ExplicitSystem (no, wait, maybe just plain System? ExplicitSystem does still all

Re: [Libmesh-users] question about the numeric vector

2018-01-25 Thread Roy Stogner
On Thu, 25 Jan 2018, Michael Povolotskyi wrote: I will add another system and will not solve it. Make sure it's ExplicitSystem (no, wait, maybe just plain System? ExplicitSystem does still allocate an rhs, not just a solution). ImplicitSystem will try to allocate a Jacobian for you whether y

Re: [Libmesh-users] question about the numeric vector

2018-01-25 Thread Roy Stogner
On Wed, 24 Jan 2018, Paul T. Bauman wrote: Sorry for slow response. I have nothing better to offer other than 3 vectors, but maybe others will. You can add another ExplicitSystem to the mesh, giving it a vector-valued variable corresponding to your solution system's scalar-valued variable. T

Re: [Libmesh-users] question about the numeric vector

2018-01-24 Thread Paul T. Bauman
Sorry for slow response. I have nothing better to offer other than 3 vectors, but maybe others will. On Mon, Jan 22, 2018 at 12:40 PM, Michael Povolotskyi wrote: > Dear Libmesh developers and users, > > I need to compute a field on mesh based on a solution. > > To do so I usually add a vector to

[Libmesh-users] question about the numeric vector

2018-01-22 Thread Michael Povolotskyi
Dear Libmesh developers and users, I need to compute a field on mesh based on a solution. To do so I usually add a vector to a system and fill it with data. This works perfectly for a scalar field. What would you recommend to do for a vector field. Of course, I can add 3 vectors, one for each