add matrix.h and coresponding operations, estimate.h and corresponding fusion operations