almost finished with mosaic convolution
This will not work with multithreads. I need to implement multithread capability in ModelCalcutor class instead of FuncLmdif class. I suspect that 2D interpolation is fast enough that a single thread can do the chi square calculation very fast.