Added + and * overload to Misfit class

Now I can do things like:

    goal = Data(x, y) + 0.1*Damping(size)
    estimate = goal.fit()

Added support for default solvers for linear and non-linear cases.
Finished srtomo example and added subclasses for Smoothness and
TotalVariation in 1D and 2D cases.