Lars Ruthotto

Hi Hari, 1. storing the dphi in the mesh structure is a good idea I think. There is a way in MATLAB to have dependent fields (i.e., that are computed...

Perfect. Thank you both for your help! On Mon, Jul 16, 2018 at 3:47 PM, James Herring wrote: > Checked out Hari's branch today. I have some updates to hyperElasticFEM.m...

I cannot reproduce this issue. It seems like there is a problem with the MATLAB path. Could you verify that the kernel folder and its subfolders are in the path?...

OK, this error message also suggests that there is an issue with your matlab path. If FAIR is your current folder you could try >> addpath(genpath(pwd)) and then call getCellCenteredGrid...

Can you try again after pulling the latest version? For some reason, it seems that this file got deleted. I've added it back.

@thraen: You might want to look into the function `sor` and other iterative solvers that are optimized for sparse matrices in https://github.com/lruthotto/KrylovMethods.jl. I use the `sor` a lot as a...

I'd be glad to help work out a common API and a convergence of the packages in terms of efficiency. I'm not sure how much time I'll have to work...

@zimoun: We chose to provide the preconditioner as a function handle that computes the inverse, for example, `M(x) = M\x` to allow as much flexibility to the user as possible....

I agree that we should distinguish between left and right preconditioners. How about having two variables, for example, `Mr` and `Ml` doing this? Regarding the abstraction issue: Can you give...

No, you need the inverse for left- and right-preconditioning. See, e.g., Chapter 9 in Saad's book: www-users.cs.umn.edu/~saad/IterMethBook_2ndEd.pdf Regarding the inner product: Think, also about when you need that **inside** the...