Thanks. I'll try to improve "my code" 2016-01-15 14:56 GMT-05:00 Jed Brown <[email protected]>:
> Matthew Knepley <[email protected]> writes: > > The way I read this, you are taking about 23 iterates/solve, and most of > > your work is residual computation which should > > be highly parallelizable/vectorizable. This seems great to me. > > This in the sense that it's up to you to determine whether your > matrix-free residual and preconditioning code is fast. This profile > merely says that almost all of the run-time is in *your code*. If your > code is fast, then this is good performance. If you can use a different > algorithm to converge in fewer iterations, or a different representation > to apply the operator faster, then you could do better. >
