def METHOD_NAME(x1,v1,x2,v2,u=0.5):
    """Subdivides a hermite curve into two hermite curves."""
    xm = hermite_eval(x1,v1,x2,v2,u)
    vm = hermite_deriv(x1,v1,x2,v2,u)
    return [(x1,vectorops.mul(v1,u),xm,vectorops.mul(vm,u)),(xm,vectorops.mul(vm,1.0-u),x2,vectorops.mul(v2,1.0-u))]