adding check for jacobian columns equal to # layout varbs

        # convert to int matrix

        J = sy.Matrix(np.matrix(J).astype(int))

    # ensure jacobian columns are all accounted for

    if len(bool_measured) != J.shape[1]:

        jac_layout_err_msg = f"Jacobian columns: {J.shape[1]} doesn' match layout columns: {len(bool_measured)}"

        raise ValueError(jac_layout_err_msg)

    bool_measured = np.array(bool_measured).astype(bool)

    Jm = J[:, list(np.where(bool_measured)[0])]