def METHOD_NAME(self, p, bcs):
    '''!
    @brief 计算 p[0], p[1], p[2] 组成的三角形的面积，
        以及每个三角形上 bcs 对应的笛卡尔坐标
    '''
    tm = np.linalg.norm(np.cross(p[1]-p[0], p[2]-p[0]), axis=1)
    point = np.einsum('qi, ifj->qfj', bcs, p)
    return tm, point