[petsc-dev] coding style

Jed Brown jed at jedbrown.org
Thu Aug 18 11:52:06 CDT 2016

Barry Smith <bsmith at mcs.anl.gov> writes:

>    Seems reasonable. Note for SNES we have something similar in that user can solve either F(x) = b or F(x) - b = 0 depending on if they provide the b explicitly.

I'd just say that the algorithm should still work correctly if the user
doesn't separate everything since that can be inefficient (multiple
communications and passes over memory) and/or unnatural.  How does the
algorithm work if the user has only g(x)?  Would they need to define

  H := \nabla^2 g(x)
  f(x) = g(x) - x'*H*x

so that you can put it back together?
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20160818/53843508/attachment.sig>

More information about the petsc-dev mailing list