one of the last commits, as I have switched over to using the LTPDA repositories.
implementing com-derivatino in code currently gives incorrect answers;
  I do not know whether the code or the math is wrong yet