C the same as matlab, not minimizing though.
Have sparse solving working for normal BA