update()

void fl::WaveSimulation1D_Real::update

(

)

Definition at line 108 of file wave_simulation_real.cpp.

                                   {
    i16 *curr = (whichGrid == 0) ? grid1.data() : grid2.data();
    i16 *next = (whichGrid == 0) ? grid2.data() : grid1.data();
 
    // Update boundaries with a Neumann (zero-gradient) condition:
    curr[0] = curr[1];
    curr[length + 1] = curr[length];
 
    // Compute dampening factor as an integer value: 2^(mDampenening)
    i32 dampening_factor = 1 << mDampenening;
 
    i32 mCourantSq32 = static_cast<i32>(mCourantSq);
    // Iterate over each inner cell.
    for (fl::size i = 1; i < length + 1; i++) {
        // Compute the 1D Laplacian:
        // lap = curr[i+1] - 2 * curr[i] + curr[i-1]
        i32 lap =
            (i32)curr[i + 1] - ((i32)curr[i] << 1) + curr[i - 1];
 
        // Multiply the Laplacian by the simulation speed using Q15 arithmetic:
        i32 term = (mCourantSq32 * lap) >> 15;
 
        // Compute the new value:
        // f = -next[i] + 2 * curr[i] + term
        i32 f = -(i32)next[i] + ((i32)curr[i] << 1) + term;
 
        // Apply damping:
        f = f - (f / dampening_factor);
 
        // Clamp the result to the Q15 range [-32768, 32767].
        if (f > 32767)
            f = 32767;
        else if (f < -32768)
            f = -32768;
 
        next[i] = (i16)f;
    }
 
    if (mHalfDuplex) {
        // Set the negative values to zero.
        for (fl::size i = 1; i < length + 1; i++) {
            if (next[i] < 0) {
                next[i] = 0;
            }
        }
    }
 
    // Toggle the active grid.
    whichGrid ^= 1;
}

References grid1, grid2, length, mCourantSq, mDampenening, mHalfDuplex, and whichGrid.