def METHOD_NAME(self, i, model, positions, index0s, fnorm=1.0):
    """gets the damping matrix for CDAMP2"""
    bi = self.B[i]
    k = bi * array([[1, -1,],
                    [-1, 1]])
    c1, c2 = self.components[i, :]
    n1, n2 = self.node_ids[i, :]
    #print('c1, c2 = %s %s' % (c1, c2))
    delta1 = 0 if c1 in [0, 1, 2, 3] else 3
    delta2 = 0 if c2 in [0, 1, 2, 3] else 3
    #print('delta1, delta2 = %s %s' % (delta1, delta2))
    n_ijv = [
        (n1, 1 + delta1),
        (n2, 1 + delta2),
    ]
    dofs = n_ijv
    return (k, dofs, n_ijv)