Start on TpetraMatrix - works in serial