* check matrix type, chose appropriate function (dune or dumux) * MultiTypeBlock matrix only works sequnetially so far
