Thanks for the feedback. I will try out the suggestions. Is it more more common to construct the global matrix (in sparse format) than have a linear operator that iterates over all the elements for each operator evaluation? I can see that storing the global matrix reduces the time per iteration, but that way also increases the memory usage. Greetings Arne Rekdal