}
    
    double BitVecEntropy(BitVec *bv) {
    ValidateBitVec(bv);

    // Test BitVecEntropy basic functionality
    bool test_bitvec_entropy_basic(void) {
    printf("Testing BitVecEntropy basic functionality\n");
    
    BitVec bv     = BitVecInit();

    BitVecPush(&bv, false);
    
    double entropy = BitVecEntropy(&bv);
    result         = result && (fabs(entropy - 1.0) < 0.001);

    BitVecPush(&bv, true);
    
    entropy = BitVecEntropy(&bv);
    result  = result && (entropy == 0.0);

    // Test BitVecEntropy edge cases
    bool test_bitvec_entropy_edge_cases(void) {
    printf("Testing BitVecEntropy edge cases\n");
    
    BitVec bv     = BitVecInit();

    
    // Test empty bitvector
    double entropy = BitVecEntropy(&bv);
    result         = result && (entropy == 0.0);

    // Test single bit
    BitVecPush(&bv, true);
    entropy = BitVecEntropy(&bv);
    result  = result && (entropy == 0.0);

    u64    dot_prod    = BitVecDotProduct(&bv1, &bv2);
    double correlation = BitVecCorrelation(&bv1, &bv2);
    double entropy1    = BitVecEntropy(&bv1);
    int    align_score = BitVecAlignmentScore(&bv1, &bv2, 1, -1);
    u64    best_align  = BitVecBestAlignment(&bv1, &bv2);

    
    bool test_bitvec_entropy_null(void) {
    printf("Testing BitVecEntropy(NULL) - should fatal\n");
    BitVecEntropy(NULL);
    return true;

    bool test_bitvec_entropy_null(void) {
    printf("Testing BitVecEntropy(NULL) - should fatal\n");
    BitVecEntropy(NULL);
    return true;
    }