add_sawtooth_wave Subroutine

public subroutine add_sawtooth_wave(tape, track, t1, t2, f, Amp, envelope)

Adds on the track a sawtooth wave with an ADSR envelope:

Arguments

Type IntentOptional Attributes Name
type(tape_recorder), intent(inout) :: tape
integer, intent(in) :: track
real(kind=wp), intent(in) :: t1
real(kind=wp), intent(in) :: t2
real(kind=wp), intent(in) :: f
real(kind=wp), intent(in) :: Amp
type(ADSR_envelope), intent(in), optional :: envelope

Calls

proc~~add_sawtooth_wave~~CallsGraph proc~add_sawtooth_wave add_sawtooth_wave proc~adsr_level ADSR_envelope%ADSR_level proc~add_sawtooth_wave->proc~adsr_level

Called by

proc~~add_sawtooth_wave~~CalledByGraph proc~add_sawtooth_wave add_sawtooth_wave program~all_signals all_signals program~all_signals->proc~add_sawtooth_wave

Source Code

    subroutine add_sawtooth_wave(tape, track, t1, t2, f, Amp, envelope)
        type(tape_recorder), intent(inout) :: tape
        integer, intent(in) :: track
        real(wp), intent(in) :: t1, t2, f, Amp
        type(ADSR_envelope), optional, intent(in) :: envelope
        ! Period in seconds:
        real(wp) :: tau
        ! Time in seconds:
        real(wp) :: t
        real(wp) :: signal
        ! ADSR Envelope value:
        real(wp) :: env
        integer  :: i

        env = 1._wp     ! Default value if no envelope is passed
        tau = 1.0_wp / f

        do concurrent(i = nint(t1*RATE) : nint(t2*RATE)-1)
            t = (i - nint(t1*RATE)) * dt

            if (present(envelope)) env = envelope%get_level(t1+t, t1, t2)

            ! We substract 0.5 for the signal to be centered on 0:
            signal = 2 * (((t/tau) - floor(t/tau)) - 0.5_wp) * Amp * env

            tape%left(track,  i) = tape%left(track,  i) + signal
            tape%right(track, i) = tape%right(track, i) + signal
        end do
    end subroutine add_sawtooth_wave