> Should I go for the vector valued approach when coding up the function F? I > then want to implement Newton's method with either finite differences, or > (preferably) automatic differentiation derivatives (thinking to follow the > tutorial and use Sacado). Again, my system has 3 discrete equations; how > could I cast this in the vector-valued framework given in the tutorial?
In order to build the system of equations you need to cast this in a single bilinear form. To see how to do this, you may want to take a look at my paper in SISC, see number 23 here: https://www.math.tamu.edu/~bangerth/publications.html#x-reviewed The paper also has a discussion on how to solve the linear problems. On the other hand, if your system is small enough, you could also use SparseDirectUMFPACK to solve it. Best W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
