better handling of rhs/solution vectors (rhs is not destroyed after solveing anymore)
storage type conversion is handled correctly (single precision solvers work!), output print level is configurable per component (default to A:size only, b: size only, x: auto), size check after initializing from file is back on.