fix matrix_toarray bug, crashed for empty rows
add ConstrainedOperator, fixes solve w/constraints
replace axpy by update
replace add_left/right_precon by set_precon