// Math functions implementation
    
    u64 BitVecHammingDistance(BitVec *bv1, BitVec *bv2) {
        ValidateBitVec(bv1);
        ValidateBitVec(bv2);

    // Test BitVecHammingDistance basic functionality
    bool test_bitvec_hamming_distance_basic(void) {
        WriteFmt("Testing BitVecHammingDistance basic functionality\n");
    
        BitVec bv1    = BitVecInit();

        BitVecPush(&bv2, true);
    
        u64 distance = BitVecHammingDistance(&bv1, &bv2);
        result       = result && (distance == 0);

        BitVecPush(&bv2, false);
    
        distance = BitVecHammingDistance(&bv1, &bv2);
        result   = result && (distance == 3);

        BitVecPush(&bv2, true); // Same
    
        distance = BitVecHammingDistance(&bv1, &bv2);
        result   = result && (distance == 1);

    // Test BitVecHammingDistance edge cases
    bool test_bitvec_hamming_distance_edge_cases(void) {
        WriteFmt("Testing BitVecHammingDistance edge cases\n");
    
        BitVec bv1    = BitVecInit();

    
        // Test empty bitvectors
        u64 distance = BitVecHammingDistance(&bv1, &bv2);
        result       = result && (distance == 0);

        BitVecPush(&bv2, true);
    
        distance = BitVecHammingDistance(&bv1, &bv2);
        result   = result && (distance == 1); // 1 length difference

        // Test one empty, one non-empty
        BitVecClear(&bv2);
        distance = BitVecHammingDistance(&bv1, &bv2);
        result   = result && (distance == 2); // Length of bv1

    
        // Test that all functions complete without crashing
        u64    hamming     = BitVecHammingDistance(&bv1, &bv2);
        double jaccard     = BitVecJaccardSimilarity(&bv1, &bv2);
        double cosine      = BitVecCosineSimilarity(&bv1, &bv2);

    
    bool test_bitvec_hamming_distance_null_bv1(void) {
        WriteFmt("Testing BitVecHammingDistance(NULL, bv2) - should fatal\n");
        BitVec bv2 = BitVecInit();
        BitVecPush(&bv2, true);

        BitVec bv2 = BitVecInit();
        BitVecPush(&bv2, true);
        BitVecHammingDistance(NULL, &bv2);
        BitVecDeinit(&bv2);
        return true;

    
    bool test_bitvec_hamming_distance_null_bv2(void) {
        WriteFmt("Testing BitVecHammingDistance(bv1, NULL) - should fatal\n");
        BitVec bv1 = BitVecInit();
        BitVecPush(&bv1, true);

        BitVec bv1 = BitVecInit();
        BitVecPush(&bv1, true);
        BitVecHammingDistance(&bv1, NULL);
        BitVecDeinit(&bv1);
        return true;