Code implementing a model repair algorithm using biharmonic fields.