these are returned as an extra array from the low
level - the high level interface remains the same

Additionally a new routine ccalc is introduced that makes calling
from C easier.