organ-pipe-waveguide 0.0.1__tar.gz

organ_pipe_waveguide-0.0.1/LICENSE

@@ -0,0 +1,7 @@
+Copyright (c) 2026 William Brady Call
+
+Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the "Software"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

organ_pipe_waveguide-0.0.1/PKG-INFO

@@ -0,0 +1,39 @@
+Metadata-Version: 2.4
+Name: organ_pipe_waveguide
+Version: 0.0.1
+Summary: A package for physically modelling the sounds of a pipe organ.
+Author-email: William Brady Call <wbradycall@gmail.com>
+License-Expression: MIT
+Project-URL: Homepage, https://github.com/wbradycall/organ_pipe_waveguide
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy>=2.1.3
+Requires-Dist: scipy>=1.17.1
+Dynamic: license-file
+
+# Organ Pipe Waveguide
+
+This project is to be used for physically modeling organ pipes. However, it's not 100% realistic just yet and it needs some tweaking. It also only includes flue pipes for now but in later versions, it will include reed pipes.
+
+To import the flue pipe model, you need to import the 'flue_pipe' class which one can import via "from organ_pipe_waveguide.FluePipe import flue_pipe". Here are the default parameters for the 'flue_pipe' object-oriented class and the methods:
+
+flue_pipe(self,note=0,default_diam=0.1555,rho=1.225,c=343,gain1=0.05,gain2=0.8,k=1.4)
+
+<!-- This method is for generating the turbulence. -->
+jet_noise(self,time,fs=44100,seed=41)
+
+<!-- This method is for creating the IIR filter function and returns a tuple of the list of beta values and the list of alpha values. -->
+create_IIR(self,fs=44100)
+
+<!-- This method is for applying an IIR filter. -->
+def IIR(self,input_array,filter_array,index,beta,alpha)
+
+<!-- This method is for adding an envelope to generated noise so that one can use it to resonate. -->
+create_input(self,time,noise_data,amp=0.5,attack=0.1)
+
+<!-- This method is for applying a non-linear jet-flue interaction simulator. -->
+nonlinearity(self,input_array)
+
+<!-- This method is for simulating the acoustics of the organ pipe and writing it to a .wav file (make to sure add '.wav' or ".wav" to the end of your filename). -->
+simulate(self,filename,duration,init_amp=0.5,fs=44100,seed=41)

organ_pipe_waveguide-0.0.1/README.md

@@ -0,0 +1,25 @@
+# Organ Pipe Waveguide
+
+This project is to be used for physically modeling organ pipes. However, it's not 100% realistic just yet and it needs some tweaking. It also only includes flue pipes for now but in later versions, it will include reed pipes.
+
+To import the flue pipe model, you need to import the 'flue_pipe' class which one can import via "from organ_pipe_waveguide.FluePipe import flue_pipe". Here are the default parameters for the 'flue_pipe' object-oriented class and the methods:
+
+flue_pipe(self,note=0,default_diam=0.1555,rho=1.225,c=343,gain1=0.05,gain2=0.8,k=1.4)
+
+<!-- This method is for generating the turbulence. -->
+jet_noise(self,time,fs=44100,seed=41)
+
+<!-- This method is for creating the IIR filter function and returns a tuple of the list of beta values and the list of alpha values. -->
+create_IIR(self,fs=44100)
+
+<!-- This method is for applying an IIR filter. -->
+def IIR(self,input_array,filter_array,index,beta,alpha)
+
+<!-- This method is for adding an envelope to generated noise so that one can use it to resonate. -->
+create_input(self,time,noise_data,amp=0.5,attack=0.1)
+
+<!-- This method is for applying a non-linear jet-flue interaction simulator. -->
+nonlinearity(self,input_array)
+
+<!-- This method is for simulating the acoustics of the organ pipe and writing it to a .wav file (make to sure add '.wav' or ".wav" to the end of your filename). -->
+simulate(self,filename,duration,init_amp=0.5,fs=44100,seed=41)

organ_pipe_waveguide-0.0.1/pyproject.toml

@@ -0,0 +1,18 @@
+[build-system]
+requires = ["setuptools >= 77.0.3"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "organ_pipe_waveguide"
+version = "0.0.1"
+dependencies = ["numpy>=2.1.3",
+                "scipy>=1.17.1",]
+authors = [{ name="William Brady Call", email="wbradycall@gmail.com"},]
+description = "A package for physically modelling the sounds of a pipe organ."
+readme = "README.md"
+requires-python = ">=3.9"
+license = "MIT"
+license-files = ["LICENSE*"]
+
+[project.urls]
+Homepage = "https://github.com/wbradycall/organ_pipe_waveguide"

organ_pipe_waveguide-0.0.1/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

organ_pipe_waveguide-0.0.1/src/organ_pipe_waveguide/FluePipe.py

@@ -0,0 +1,130 @@
+import numpy as np
+from numpy.fft import fft,ifft,fftfreq
+from scipy import signal
+import scipy.io.wavfile as wavfile
+
+
+class flue_pipe: # Organ flue pipe simulator
+    def __init__(self,note=0,default_diam=0.1555,rho=1.225,c=343,gain1=0.05,gain2=0.8,k=1.4):
+        self.note = note
+        self.f0 = 55 * pow(2,1/4) * pow(2,note/12)
+        self.r = (default_diam / 2) * pow(2,-note/16)
+        self.rho = rho
+        self.c = c
+        self.gain1 = gain1
+        self.gain2 = gain2
+        self.k = k
+
+    def __repr__(self):
+        return f"FluePipe(freq:{self.f0},radius:{self.r},rho:{self.rho},wavespeed:{self.c},gain1:{self.gain1},gain2:{self.gain2},nonlin_coeff:{self.k})"
+
+    def __str__(self):
+        return f"FluePipe(freq:{self.f0},radius:{self.r},wavespeed:{self.c})"
+
+    def jet_noise(self,time,fs=44100,seed=41):
+        N = len(time)
+        np.random.seed(seed)
+        noise_data_derivative3 = np.random.normal(0, 1, N)
+
+        cutoff = 50 * pow(2, self.note / 12)
+        frequencies = fftfreq(N, 1 / fs)
+        omegas = 2 * np.pi * frequencies
+
+        FFT = fft(noise_data_derivative3)
+        newFFT = np.zeros(N, dtype=complex)
+
+        for i in range(N):
+            if np.abs(frequencies[i]) <= cutoff:
+                pass
+            else:
+                newFFT[i] = FFT[i] / (1j * np.abs(omegas[i])) ** 3
+
+        y = ifft(newFFT).real
+        y /= np.max(np.abs(y))
+        return y
+
+    def create_IIR(self,fs=44100):
+        cutoff_freq = 1.84 * self.c / (2 * np.pi * self.r)
+        cutoff_om = 2 * np.pi * cutoff_freq
+        num = [cutoff_om]
+        den = [1, cutoff_om]
+
+        beta, alpha = signal.bilinear(num, den, fs=fs)
+        return beta, alpha
+
+    def IIR(self,input_array,filter_array,index,beta,alpha):
+        b0 = beta[0]
+        b1 = beta[1]
+
+        a0 = alpha[0]
+        a1 = alpha[1]
+
+        return (b0 * input_array[index] + b1 * input_array[index - 1] - a1 * filter_array[index - 1]) / a0
+
+    def create_input(self,time,noise_data,amp=0.5,attack=0.1):
+        N = len(time)
+        duration = np.max(time)
+
+        envelope = np.zeros(N)
+        release1 = (4 / 5) * duration
+        release2 = (22 / 25) * duration
+
+        for i in range(N):
+            if time[i] < attack:
+                envelope[i] = time[i] / attack
+            elif attack <= time[i] and time[i] < release1:
+                envelope[i] = 1
+            elif release1 <= time[i] and time[i] < release2:
+                envelope[i] = np.exp(-5 * (time[i] - release1))
+
+        return amp * noise_data * envelope
+
+    def nonlinearity(self,input_array):
+        return (self.k * input_array) + 0.005 * (self.k * input_array)**2 - (self.k * input_array)**3
+
+    def simulate(self,filename,duration,init_amp=0.5,fs=44100,seed=41):
+        N = int(duration * fs)
+        T = 1 / self.f0
+        delay = T * fs
+        start_resonating = int(delay)
+        beta, alpha = self.create_IIR(fs=fs)
+
+        time = np.linspace(0,duration,N,endpoint=False)
+        noise_data = self.jet_noise(time,fs=fs,seed=seed)
+        input_array = self.create_input(time,noise_data,amp=init_amp)
+        jet_input = np.zeros(N)
+        filter_array = np.zeros(N)
+        output_array = np.zeros(N)
+
+        for i in range(start_resonating):
+            jet_input[i] = self.nonlinearity(self.gain1 * input_array[i])
+            filter_array[i] = self.IIR(jet_input,filter_array,i,beta,alpha)
+            output_array[i] = filter_array[i]
+            print(i)
+
+        for i in range(start_resonating,N):
+            k = i-delay # Fractional delay
+            k_floor = int(np.floor(k))
+            k_ceil = int(np.ceil(k))
+            fraction = k - k_floor
+
+            output_state = (1 - fraction) * output_array[k_floor] + fraction * output_array[k_ceil]
+
+            jet_input[i] = self.nonlinearity(self.gain1 * input_array[i] + self.gain2 * output_state) # Open end
+            filter_array[i] = self.IIR(jet_input,filter_array,i,beta,alpha)
+            output_array[i] = filter_array[i]
+            print(i)
+
+        print("Finished resonating!\n")
+        print("Converting to stereo and writing to .wav file...\n")
+
+        output_array /= np.max(np.abs(output_array))
+        output_array2 = np.zeros(N)
+
+        for i in range(N - 441):
+            output_array2[i + 441] = 0.8 * output_array[i]
+
+        stereo_output = np.vstack((output_array, output_array2)).T
+        wavfile.write(filename, fs, stereo_output)
+
+        print("Finished!")

organ_pipe_waveguide-0.0.1/src/organ_pipe_waveguide/__init__.py

@@ -0,0 +1 @@
+__version__ = "0.0.1"

organ_pipe_waveguide-0.0.1/src/organ_pipe_waveguide.egg-info/PKG-INFO

@@ -0,0 +1,39 @@
+Metadata-Version: 2.4
+Name: organ_pipe_waveguide
+Version: 0.0.1
+Summary: A package for physically modelling the sounds of a pipe organ.
+Author-email: William Brady Call <wbradycall@gmail.com>
+License-Expression: MIT
+Project-URL: Homepage, https://github.com/wbradycall/organ_pipe_waveguide
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy>=2.1.3
+Requires-Dist: scipy>=1.17.1
+Dynamic: license-file
+
+# Organ Pipe Waveguide
+
+This project is to be used for physically modeling organ pipes. However, it's not 100% realistic just yet and it needs some tweaking. It also only includes flue pipes for now but in later versions, it will include reed pipes.
+
+To import the flue pipe model, you need to import the 'flue_pipe' class which one can import via "from organ_pipe_waveguide.FluePipe import flue_pipe". Here are the default parameters for the 'flue_pipe' object-oriented class and the methods:
+
+flue_pipe(self,note=0,default_diam=0.1555,rho=1.225,c=343,gain1=0.05,gain2=0.8,k=1.4)
+
+<!-- This method is for generating the turbulence. -->
+jet_noise(self,time,fs=44100,seed=41)
+
+<!-- This method is for creating the IIR filter function and returns a tuple of the list of beta values and the list of alpha values. -->
+create_IIR(self,fs=44100)
+
+<!-- This method is for applying an IIR filter. -->
+def IIR(self,input_array,filter_array,index,beta,alpha)
+
+<!-- This method is for adding an envelope to generated noise so that one can use it to resonate. -->
+create_input(self,time,noise_data,amp=0.5,attack=0.1)
+
+<!-- This method is for applying a non-linear jet-flue interaction simulator. -->
+nonlinearity(self,input_array)
+
+<!-- This method is for simulating the acoustics of the organ pipe and writing it to a .wav file (make to sure add '.wav' or ".wav" to the end of your filename). -->
+simulate(self,filename,duration,init_amp=0.5,fs=44100,seed=41)

organ_pipe_waveguide-0.0.1/src/organ_pipe_waveguide.egg-info/SOURCES.txt

@@ -0,0 +1,10 @@
+LICENSE
+README.md
+pyproject.toml
+src/organ_pipe_waveguide/FluePipe.py
+src/organ_pipe_waveguide/__init__.py
+src/organ_pipe_waveguide.egg-info/PKG-INFO
+src/organ_pipe_waveguide.egg-info/SOURCES.txt
+src/organ_pipe_waveguide.egg-info/dependency_links.txt
+src/organ_pipe_waveguide.egg-info/requires.txt
+src/organ_pipe_waveguide.egg-info/top_level.txt

organ_pipe_waveguide-0.0.1/src/organ_pipe_waveguide.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

organ_pipe_waveguide-0.0.1/src/organ_pipe_waveguide.egg-info/requires.txt

@@ -0,0 +1,2 @@
+numpy>=2.1.3
+scipy>=1.17.1

organ_pipe_waveguide-0.0.1/src/organ_pipe_waveguide.egg-info/top_level.txt

@@ -0,0 +1 @@
+organ_pipe_waveguide