geometric/internal.py
def d_cross_ab(a, b, da, db):
    """
    Given two vectors a, b and their derivatives w/r.t. a parameter, return the derivative
    of the cross product
    """
    answer = np.zeros((da.shape[0], 3), dtype=float)
    for i in range(da.shape[0]):
        answer[i] = np.cross(a, db[i]) + np.cross(da[i], b)
    return answer

    Given two vectors a, b and their derivatives w/r.t. a parameter, return the derivative
    of the cross product
    
Given two vectors a, b and their derivatives w/r.t. a parameter, return the derivative
of the cross product
Given two vectors a b and their derivatives w / r . t . a parameter return the derivative of the cross product
def d_cross_ab(a, b, da, db):
    answer = np.zeros((da.shape[0], 3), dtype=float)
    for i in range(da.shape[0]):
        answer[i] = np.cross(a, db[i]) + np.cross(da[i], b)
    return answer
def d_cross_ab ( a , b , da , db ) : answer = np . zeros ( ( da . shape [ 0 ] , 3 ) , dtype = float ) for i in range ( da . shape [ 0 ] ) : answer [ i ] = np . cross ( a , db [ i ] ) + np . cross ( da [ i ] , b ) return answer
Given two vectors a, b and their derivatives w/r.t.
Given two vectors a b and their derivatives w / r . t .
"""
    Given two vectors a, b and their derivatives w/r.t. a parameter, return the derivative
    of the cross product
    """
a
b
da
db