aprender-core 0.29.1

use super::*;

// ==================== Additional Coverage Tests ====================

#[test]
fn test_global_max_pool_batched() {
    let x = Tensor::new(&[1.0, 2.0, 5.0, 4.0, 3.0, 6.0, 7.0, 0.0], &[4, 2]);
    let batch = vec![0, 0, 1, 1]; // 2 graphs, 2 nodes each

    let pooled = global_max_pool(&x, Some(&batch));

    assert_eq!(pooled.shape(), &[2, 2]);
    let data = pooled.data();
    // Graph 0: max([1,5], [2,4]) = [5, 4]
    assert!((data[0] - 5.0).abs() < 1e-6);
    assert!((data[1] - 4.0).abs() < 1e-6);
    // Graph 1: max([3,7], [6,0]) = [7, 6]
    assert!((data[2] - 7.0).abs() < 1e-6);
    assert!((data[3] - 6.0).abs() < 1e-6);
}

#[test]
fn test_global_sum_pool_batched() {
    let x = Tensor::new(&[1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0], &[4, 2]);
    let batch = vec![0, 0, 1, 1];

    let pooled = global_sum_pool(&x, Some(&batch));

    assert_eq!(pooled.shape(), &[2, 2]);
    let data = pooled.data();
    // Graph 0: sum([1,3], [2,4]) = [4, 6]
    assert!((data[0] - 4.0).abs() < 1e-6);
    assert!((data[1] - 6.0).abs() < 1e-6);
    // Graph 1: sum([5,7], [6,8]) = [12, 14]
    assert!((data[2] - 12.0).abs() < 1e-6);
    assert!((data[3] - 14.0).abs() < 1e-6);
}

#[test]
fn test_replace_neg_infinity() {
    let mut values = vec![1.0, f32::NEG_INFINITY, 3.0, f32::NEG_INFINITY];
    replace_neg_infinity(&mut values);
    assert!((values[0] - 1.0).abs() < 1e-6);
    assert!((values[1] - 0.0).abs() < 1e-6);
    assert!((values[2] - 3.0).abs() < 1e-6);
    assert!((values[3] - 0.0).abs() < 1e-6);
}

#[test]
fn test_gcn_parameters_mut() {
    let mut gcn = GCNConv::new(4, 8);
    let params_mut = gcn.parameters_mut();
    assert_eq!(params_mut.len(), 2);
}

#[test]
fn test_gat_parameters_mut() {
    let mut gat = GATConv::new(4, 8, 2);
    let params_mut = gat.parameters_mut();
    assert_eq!(params_mut.len(), 4);
}

#[test]
fn test_gin_parameters_mut() {
    let mut gin = GINConv::new(4, 16, 8);
    let params_mut = gin.parameters_mut();
    assert_eq!(params_mut.len(), 4);
}

#[test]
fn test_graphsage_parameters_mut() {
    let mut sage = GraphSAGEConv::new(4, 8);
    let params_mut = sage.parameters_mut();
    assert_eq!(params_mut.len(), 2);
}

#[test]
fn test_edgeconv_parameters_mut() {
    let mut edge_conv = EdgeConv::new(4, 16, 8);
    let params_mut = edge_conv.parameters_mut();
    assert_eq!(params_mut.len(), 4);
}

#[test]
fn test_gat_with_disconnected_node() {
    // Test GAT with a node that has no neighbors (except self-loop)
    let gat = GATConv::new(4, 8, 1);
    let x = Tensor::ones(&[4, 4]);
    // Node 3 is disconnected
    let edges = vec![(0, 1), (1, 2)];

    let out = gat.forward_gnn(&x, &edges);

    assert_eq!(out.shape(), &[4, 8]);
    // All outputs should still be finite
    for &v in out.data() {
        assert!(v.is_finite());
    }
}

#[test]
fn test_gin_with_varying_features() {
    let gin = GINConv::new(2, 8, 4);
    let x = Tensor::new(&[1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0], &[4, 2]);
    let edges = vec![(0, 1), (1, 2), (2, 3)];

    let out = gin.forward_gnn(&x, &edges);

    assert_eq!(out.shape(), &[4, 4]);
    for &v in out.data() {
        assert!(v.is_finite());
    }
}

#[test]
fn test_graphsage_with_isolated_node() {
    let sage = GraphSAGEConv::new(4, 8);
    let x = Tensor::ones(&[5, 4]);
    // Node 4 is isolated
    let edges = vec![(0, 1), (1, 2), (2, 3)];

    let out = sage.forward_gnn(&x, &edges);

    assert_eq!(out.shape(), &[5, 8]);
    // Isolated node should still produce output (zeros for neighbors)
    for &v in out.data() {
        assert!(v.is_finite());
    }
}

#[test]
fn test_graphsage_max_with_single_neighbor() {
    let sage = GraphSAGEConv::new(4, 8).with_aggregation(SAGEAggregation::Max);
    let x = Tensor::new(
        &[
            1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0,
        ],
        &[3, 4],
    );
    // Node 0 only connected to node 1
    let edges = vec![(0, 1)];

    let out = sage.forward_gnn(&x, &edges);

    assert_eq!(out.shape(), &[3, 8]);
    for &v in out.data() {
        assert!(v.is_finite());
    }
}

#[test]
fn test_global_mean_pool_empty_batch() {
    let x = Tensor::new(&[1.0, 2.0, 3.0, 4.0], &[2, 2]);
    // Empty batch (no nodes assigned to any graph) - edge case
    let batch = vec![0, 0];

    let pooled = global_mean_pool(&x, Some(&batch));

    // Should have 1 graph
    assert_eq!(pooled.shape(), &[1, 2]);
}

#[test]
fn test_gcn_debug_format() {
    let gcn = GCNConv::new(4, 8);
    let debug_str = format!("{:?}", gcn);
    assert!(debug_str.contains("GCNConv"));
}

#[test]
fn test_gat_debug_format() {
    let gat = GATConv::new(4, 8, 2);
    let debug_str = format!("{:?}", gat);
    assert!(debug_str.contains("GATConv"));
}

#[test]
fn test_gin_debug_format() {
    let gin = GINConv::new(4, 16, 8);
    let debug_str = format!("{:?}", gin);
    assert!(debug_str.contains("GINConv"));
}

#[test]
fn test_graphsage_debug_format() {
    let sage = GraphSAGEConv::new(4, 8);
    let debug_str = format!("{:?}", sage);
    assert!(debug_str.contains("GraphSAGEConv"));
}

#[test]
fn test_edgeconv_debug_format() {
    let edge_conv = EdgeConv::new(4, 16, 8);
    let debug_str = format!("{:?}", edge_conv);
    assert!(debug_str.contains("EdgeConv"));
}

#[test]
fn test_sage_aggregation_equality() {
    let agg1 = SAGEAggregation::Mean;
    let agg2 = SAGEAggregation::Mean;
    let agg3 = SAGEAggregation::Max;

    assert_eq!(agg1, agg2);
    assert_ne!(agg1, agg3);
}

#[test]
fn test_sage_aggregation_clone() {
    let agg = SAGEAggregation::Sum;
    let cloned = agg;
    assert_eq!(cloned, SAGEAggregation::Sum);
}

#[test]
fn test_sage_aggregation_debug() {
    let agg = SAGEAggregation::Max;
    let debug_str = format!("{:?}", agg);
    assert!(debug_str.contains("Max"));
}

#[test]
#[should_panic(expected = "requires graph structure")]
fn test_gin_forward_panics() {
    let gin = GINConv::new(4, 16, 8);
    let x = Tensor::ones(&[3, 4]);
    let _ = gin.forward(&x);
}

#[test]
#[should_panic(expected = "requires graph structure")]
fn test_gat_forward_panics() {
    let gat = GATConv::new(4, 8, 2);
    let x = Tensor::ones(&[3, 4]);
    let _ = gat.forward(&x);
}

#[test]
fn test_gcn_with_dense_graph() {
    let gcn = GCNConv::new(4, 8);
    let x = Tensor::ones(&[4, 4]);
    // Fully connected graph
    let edges = vec![(0, 1), (0, 2), (0, 3), (1, 2), (1, 3), (2, 3)];

    let out = gcn.forward_gnn(&x, &edges);

    assert_eq!(out.shape(), &[4, 8]);
    for &v in out.data() {
        assert!(v.is_finite());
    }
}

#[test]
fn test_gat_multiple_heads_aggregation() {
    let gat = GATConv::new(4, 4, 4); // 4 heads, 4 features per head
    let x = Tensor::new(
        &[
            1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0,
        ],
        &[3, 4],
    );
    let edges = vec![(0, 1), (1, 2), (2, 0)];

    let out = gat.forward_gnn(&x, &edges);

    // Output should be num_heads * out_features = 4 * 4 = 16
    assert_eq!(out.shape(), &[3, 16]);
    for &v in out.data() {
        assert!(v.is_finite());
    }
}

#[test]
fn test_edgeconv_self_loop_behavior() {
    let edge_conv = EdgeConv::new(2, 8, 4);
    let x = Tensor::new(&[1.0, 2.0, 3.0, 4.0], &[2, 2]);
    // No edges - should create self-loops internally
    let edges: Vec<EdgeIndex> = vec![];

    let out = edge_conv.forward_gnn(&x, &edges);

    assert_eq!(out.shape(), &[2, 4]);
    for &v in out.data() {
        assert!(v.is_finite());
    }
}

#[test]
fn test_accumulate_functions() {
    // Test accumulate_single indirectly through global_sum_pool
    let x = Tensor::new(&[1.0, 2.0, 3.0, 4.0, 5.0, 6.0], &[3, 2]);
    let pooled = global_sum_pool(&x, None);
    let data = pooled.data();
    assert!((data[0] - 9.0).abs() < 1e-6); // 1+3+5
    assert!((data[1] - 12.0).abs() < 1e-6); // 2+4+6
}

#[test]
fn test_global_pool_single_node() {
    let x = Tensor::new(&[5.0, 10.0], &[1, 2]);

    let mean = global_mean_pool(&x, None);
    let sum = global_sum_pool(&x, None);
    let max = global_max_pool(&x, None);

    // With single node, mean = sum = max = original
    assert_eq!(mean.shape(), &[1, 2]);
    assert!((mean.data()[0] - 5.0).abs() < 1e-6);
    assert!((sum.data()[0] - 5.0).abs() < 1e-6);
    assert!((max.data()[0] - 5.0).abs() < 1e-6);
}

#[test]
fn test_gin_with_nonzero_eps() {
    let mut gin = GINConv::new(4, 16, 8);
    gin.set_eps(0.5);

    let x = Tensor::ones(&[3, 4]);
    let edges = vec![(0, 1), (1, 2)];

    let out = gin.forward_gnn(&x, &edges);

    assert_eq!(out.shape(), &[3, 8]);
    // eps affects self-contribution weight
    for &v in out.data() {
        assert!(v.is_finite());
    }
}