RsBundle  Diff

                    .iter()
                    .map(|&cidx| {
                        subdatas
                            .iter()
                            .enumerate()
                            .map(|(fidx, sub)| {
                                let primal = state.aggregated_primal(fidx);
                                sub.evaluate_constraint(&primal, cidx)
                                sub.evaluate_constraint(primal, cidx)
                            })
                            .sum::<Real>()
                    })
                    .enumerate()
                    .filter_map(|(i, sg)| if sg < 1e-3 { Some(i) } else { None })
                    .collect::<Vec<_>>()
            };