[aspect-devel] Implementing another solution variable

[Jonathan Perry-Houts]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

To implement this Keller melt transport model, I'm adding another
(optional) solution variable in Aspect that will eventually be solved
along with velocity and pressure in the stokes step. Does it matter in
which order the FESystem components are organized? ie, should I put
porosity before temperature, and then adjust the introspection indices
of temperature and compositions and probably a dozen other hard coded
spots in in the code, or can I tack the new solution component on to
the end and use blocks to access that chunk of the matrix?

Also, since I'll be modifying the system matrix, I'd like to be able
to use the direct solver (which is still a work in progress at this
time, although it looks like Timo's been making good headway on it).
Anyway, that means this new variable will need to be able to be
(optionally) part of the same "block" as velocity and pressure. Is it
possible for example to make block 0 contain components 0-dim and
dim+2+n_compositional_fields?

Thanks!

- -Jonathan
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1
Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/

iQEcBAEBAgAGBQJTekDwAAoJEGe6xJ1FYRpRowoH/REDm5BpnqCrnzgkb5Z2OJ26
3eGOBJJKuqj3KcAb8oB3BUsvJsiQ4zugm5/hIoY61hPQ551izGuE65ezZ8Eejtli
bRikkOKnzxFVmEP1/mGCY7OfhB/BFRKwB526SKE1Zdi5TvGtYbGVit9zRuCCSp+g
Kf8ELnaTbzsOhM4ePOYQjXryNMkhcWT1cFo3o7SkyCh7pDzIxNWS0GDPdbOl/xr+
cw+9HLCSSLx91H+pkDKwUJBSzspWIoxT4VtyvaP6ujC/TPrulBO8sNyWmah7vxJd
JitEkNrsf/Lg6Cwzj9JDDwTgmHg0tg4eAOveRGtdYysitckHdofHH3bDASwHRw0=
=Exul
-----END PGP SIGNATURE-----