create matrix C
new estimation.py passes
check estimation
dev and opt results difference fixed
resolve conflict between branch dev and opt