TODO #2
some arithmetic errors happen, for now check it
Fixes #1, rename cbelt1 -> cbref1 & include performance tests