Wouldn't that give me a vector instead of a (sparse) matrix which I need
afterwards?
Am Mittwoch, 21. November 2018 15:57:08 UTC+1 schrieb Wolfgang Bangerth:
>
> On 11/21/18 5:39 AM, 'Maxi Miller' via deal.II User Group wrote:
> > Hmm, but in that case I have an addition (update of the sparse
> jacobian),
> > which I do not know how to handle (yet)
>
> Well, if you also have a previous Jacobian matrix, say J, then you would
> replace...
>
> > void BroydenOperator::vmult (Vector<double> &dst,
> > const Vector<double> &src) const
> > {
> > dst = 0;
>
> ...this line by
> J.vmult (dst, src);
> and keep the rest...
>
> > for (unsigned int k=0; k<n_vectors; ++k)
> > {
> > const double xk_dot_src = x[k] * src;
> > dst.add (xk_dot_src, y[k]);
> > }
> > }
> ...as is. This corresponds to the operator
>
> B = J + sum_k y_k x_k^T
>
> Best
> W.
>
> --
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: bang...@colostate.edu
> <javascript:>
> www: http://www.math.colostate.edu/~bangerth/
>
>
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
