gistfile1.txt

multi method KnotVector.N0_index(Int $p, $u, KnotBasisDirection $direction = Left)
{
    given $direction
    {
        when Left { UpperBound(@.knots, $u) - $p - 1; }
        when Right { LowerBound(@.knots, $u) - $p - 1; }
    }
}

multi method Nubs.Direction($t)
{
    $t < ([+] self.ParameterRange) / 2 ?? Left !! Right;
}

multi method Nubs.evaluate($t, KnotBasisDirection $direction = self.Direction($t))
{
    my $n0 = $.knot_vector.N0_index($.degree, $t, $direction);
    return [+] ($.knot_vector.N_local($n0, $.degree, $t) 
                >>*<< @.control_points[$n0 .. ($n0 + $.degree)]);
}