Add fast trap_burden mode with lightweight bulk/original-basis computation
#371
+203
−46
trap_burden mode with lightweight bulk/original-basis computation
#371