1 #include "stdio.h"
2 #ifndef mips
3 #include "stdlib.h"
4 #endif
5 #include "xlisp.h"
6 #include "sound.h"
7
8 #include "falloc.h"
9 #include "cext.h"
10 #include "instrclarall.h"
11
12 void clarinet_all_free(snd_susp_type a_susp);
13
14
15 typedef struct clarinet_all_susp_struct {
16 snd_susp_node susp;
17 int64_t terminate_cnt;
18 sound_type breath_env;
19 int breath_env_cnt;
20 sample_block_values_type breath_env_ptr;
21 sound_type freq_env;
22 int freq_env_cnt;
23 sample_block_values_type freq_env_ptr;
24 sound_type reed_stiffness;
25 int reed_stiffness_cnt;
26 sample_block_values_type reed_stiffness_ptr;
27 sound_type noise_env;
28 int noise_env_cnt;
29 sample_block_values_type noise_env_ptr;
30
31 struct instr *clar;
32 double frequency;
33 float breath_scale;
34 float reed_scale;
35 float noise_scale;
36 } clarinet_all_susp_node, *clarinet_all_susp_type;
37
38 #include "instr.h"
39 #include "upsample.h"
40
41
clarinet_all_nsnn_fetch(snd_susp_type a_susp,snd_list_type snd_list)42 void clarinet_all_nsnn_fetch(snd_susp_type a_susp, snd_list_type snd_list)
43 {
44 clarinet_all_susp_type susp = (clarinet_all_susp_type) a_susp;
45 int cnt = 0; /* how many samples computed */
46 int togo;
47 int n;
48 sample_block_type out;
49 register sample_block_values_type out_ptr;
50
51 register sample_block_values_type out_ptr_reg;
52
53 register struct instr * clar_reg;
54 register double frequency_reg;
55 register float breath_scale_reg;
56 register float reed_scale_reg;
57 register float noise_scale_reg;
58 register sample_block_values_type noise_env_ptr_reg;
59 register sample_block_values_type reed_stiffness_ptr_reg;
60 register sample_type freq_env_scale_reg = susp->freq_env->scale;
61 register sample_block_values_type freq_env_ptr_reg;
62 register sample_block_values_type breath_env_ptr_reg;
63 falloc_sample_block(out, "clarinet_all_nsnn_fetch");
64 out_ptr = out->samples;
65 snd_list->block = out;
66
67 while (cnt < max_sample_block_len) { /* outer loop */
68 /* first compute how many samples to generate in inner loop: */
69 /* don't overflow the output sample block: */
70 togo = max_sample_block_len - cnt;
71
72 /* don't run past the breath_env input sample block: */
73 susp_check_term_samples(breath_env, breath_env_ptr, breath_env_cnt);
74 togo = min(togo, susp->breath_env_cnt);
75
76 /* don't run past the freq_env input sample block: */
77 susp_check_samples(freq_env, freq_env_ptr, freq_env_cnt);
78 togo = min(togo, susp->freq_env_cnt);
79
80 /* don't run past the reed_stiffness input sample block: */
81 susp_check_samples(reed_stiffness, reed_stiffness_ptr, reed_stiffness_cnt);
82 togo = min(togo, susp->reed_stiffness_cnt);
83
84 /* don't run past the noise_env input sample block: */
85 susp_check_samples(noise_env, noise_env_ptr, noise_env_cnt);
86 togo = min(togo, susp->noise_env_cnt);
87
88 /* don't run past terminate time */
89 if (susp->terminate_cnt != UNKNOWN &&
90 susp->terminate_cnt <= susp->susp.current + cnt + togo) {
91 togo = (int) (susp->terminate_cnt - (susp->susp.current + cnt));
92 if (togo < 0) togo = 0; /* avoids rounding errros */
93 if (togo == 0) break;
94 }
95
96 n = togo;
97 clar_reg = susp->clar;
98 frequency_reg = susp->frequency;
99 breath_scale_reg = susp->breath_scale;
100 reed_scale_reg = susp->reed_scale;
101 noise_scale_reg = susp->noise_scale;
102 noise_env_ptr_reg = susp->noise_env_ptr;
103 reed_stiffness_ptr_reg = susp->reed_stiffness_ptr;
104 freq_env_ptr_reg = susp->freq_env_ptr;
105 breath_env_ptr_reg = susp->breath_env_ptr;
106 out_ptr_reg = out_ptr;
107 if (n) do { /* the inner sample computation loop */
108 controlChange(clar_reg, 128, breath_scale_reg * *breath_env_ptr_reg++);
109 controlChange(clar_reg, 2, reed_scale_reg * *reed_stiffness_ptr_reg++);
110 controlChange(clar_reg, 4, noise_scale_reg * *noise_env_ptr_reg++);
111 setFrequency(clar_reg, frequency_reg + (freq_env_scale_reg * *freq_env_ptr_reg++));
112 *out_ptr_reg++ = (sample_type) tick(clar_reg);
113 } while (--n); /* inner loop */
114
115 susp->clar = clar_reg;
116 /* using noise_env_ptr_reg is a bad idea on RS/6000: */
117 susp->noise_env_ptr += togo;
118 /* using reed_stiffness_ptr_reg is a bad idea on RS/6000: */
119 susp->reed_stiffness_ptr += togo;
120 /* using freq_env_ptr_reg is a bad idea on RS/6000: */
121 susp->freq_env_ptr += togo;
122 /* using breath_env_ptr_reg is a bad idea on RS/6000: */
123 susp->breath_env_ptr += togo;
124 out_ptr += togo;
125 susp_took(breath_env_cnt, togo);
126 susp_took(freq_env_cnt, togo);
127 susp_took(reed_stiffness_cnt, togo);
128 susp_took(noise_env_cnt, togo);
129 cnt += togo;
130 } /* outer loop */
131
132 /* test for termination */
133 if (togo == 0 && cnt == 0) {
134 snd_list_terminate(snd_list);
135 } else {
136 snd_list->block_len = cnt;
137 susp->susp.current += cnt;
138 }
139 } /* clarinet_all_nsnn_fetch */
140
141
clarinet_all_toss_fetch(snd_susp_type a_susp,snd_list_type snd_list)142 void clarinet_all_toss_fetch(snd_susp_type a_susp, snd_list_type snd_list)
143 {
144 clarinet_all_susp_type susp = (clarinet_all_susp_type) a_susp;
145 time_type final_time = susp->susp.t0;
146 int n;
147
148 /* fetch samples from breath_env up to final_time for this block of zeros */
149 while ((ROUNDBIG((final_time - susp->breath_env->t0) * susp->breath_env->sr)) >=
150 susp->breath_env->current)
151 susp_get_samples(breath_env, breath_env_ptr, breath_env_cnt);
152 /* fetch samples from freq_env up to final_time for this block of zeros */
153 while ((ROUNDBIG((final_time - susp->freq_env->t0) * susp->freq_env->sr)) >=
154 susp->freq_env->current)
155 susp_get_samples(freq_env, freq_env_ptr, freq_env_cnt);
156 /* fetch samples from reed_stiffness up to final_time for this block of zeros */
157 while ((ROUNDBIG((final_time - susp->reed_stiffness->t0) * susp->reed_stiffness->sr)) >=
158 susp->reed_stiffness->current)
159 susp_get_samples(reed_stiffness, reed_stiffness_ptr, reed_stiffness_cnt);
160 /* fetch samples from noise_env up to final_time for this block of zeros */
161 while ((ROUNDBIG((final_time - susp->noise_env->t0) * susp->noise_env->sr)) >=
162 susp->noise_env->current)
163 susp_get_samples(noise_env, noise_env_ptr, noise_env_cnt);
164 /* convert to normal processing when we hit final_count */
165 /* we want each signal positioned at final_time */
166 n = (int) ROUNDBIG((final_time - susp->breath_env->t0) * susp->breath_env->sr -
167 (susp->breath_env->current - susp->breath_env_cnt));
168 susp->breath_env_ptr += n;
169 susp_took(breath_env_cnt, n);
170 n = (int) ROUNDBIG((final_time - susp->freq_env->t0) * susp->freq_env->sr -
171 (susp->freq_env->current - susp->freq_env_cnt));
172 susp->freq_env_ptr += n;
173 susp_took(freq_env_cnt, n);
174 n = (int) ROUNDBIG((final_time - susp->reed_stiffness->t0) * susp->reed_stiffness->sr -
175 (susp->reed_stiffness->current - susp->reed_stiffness_cnt));
176 susp->reed_stiffness_ptr += n;
177 susp_took(reed_stiffness_cnt, n);
178 n = (int) ROUNDBIG((final_time - susp->noise_env->t0) * susp->noise_env->sr -
179 (susp->noise_env->current - susp->noise_env_cnt));
180 susp->noise_env_ptr += n;
181 susp_took(noise_env_cnt, n);
182 susp->susp.fetch = susp->susp.keep_fetch;
183 (*(susp->susp.fetch))(a_susp, snd_list);
184 }
185
186
clarinet_all_mark(snd_susp_type a_susp)187 void clarinet_all_mark(snd_susp_type a_susp)
188 {
189 clarinet_all_susp_type susp = (clarinet_all_susp_type) a_susp;
190 sound_xlmark(susp->breath_env);
191 sound_xlmark(susp->freq_env);
192 sound_xlmark(susp->reed_stiffness);
193 sound_xlmark(susp->noise_env);
194 }
195
196
clarinet_all_free(snd_susp_type a_susp)197 void clarinet_all_free(snd_susp_type a_susp)
198 {
199 clarinet_all_susp_type susp = (clarinet_all_susp_type) a_susp;
200 deleteInstrument(susp->clar);
201 sound_unref(susp->breath_env);
202 sound_unref(susp->freq_env);
203 sound_unref(susp->reed_stiffness);
204 sound_unref(susp->noise_env);
205 ffree_generic(susp, sizeof(clarinet_all_susp_node), "clarinet_all_free");
206 }
207
208
clarinet_all_print_tree(snd_susp_type a_susp,int n)209 void clarinet_all_print_tree(snd_susp_type a_susp, int n)
210 {
211 clarinet_all_susp_type susp = (clarinet_all_susp_type) a_susp;
212 indent(n);
213 stdputstr("breath_env:");
214 sound_print_tree_1(susp->breath_env, n);
215
216 indent(n);
217 stdputstr("freq_env:");
218 sound_print_tree_1(susp->freq_env, n);
219
220 indent(n);
221 stdputstr("reed_stiffness:");
222 sound_print_tree_1(susp->reed_stiffness, n);
223
224 indent(n);
225 stdputstr("noise_env:");
226 sound_print_tree_1(susp->noise_env, n);
227 }
228
229
snd_make_clarinet_all(double freq,sound_type breath_env,sound_type freq_env,double vibrato_freq,double vibrato_gain,sound_type reed_stiffness,sound_type noise_env,rate_type sr)230 sound_type snd_make_clarinet_all(double freq, sound_type breath_env, sound_type freq_env, double vibrato_freq, double vibrato_gain, sound_type reed_stiffness, sound_type noise_env, rate_type sr)
231 {
232 register clarinet_all_susp_type susp;
233 /* sr specified as input parameter */
234 time_type t0 = breath_env->t0;
235 sample_type scale_factor = 1.0F;
236 time_type t0_min = t0;
237 falloc_generic(susp, clarinet_all_susp_node, "snd_make_clarinet_all");
238 susp->clar = initInstrument(CLARINET, ROUND32(sr));
239 noteOn(susp->clar, freq, 1.0);
240 controlChange(susp->clar, 11, CLAR_CONTROL_CHANGE_CONST * vibrato_freq);
241 controlChange(susp->clar, 1, CLAR_CONTROL_CHANGE_CONST * vibrato_gain);;
242 susp->frequency = freq;
243 susp->breath_scale = breath_env->scale * CLAR_CONTROL_CHANGE_CONST;
244 susp->reed_scale = reed_stiffness->scale * CLAR_CONTROL_CHANGE_CONST;
245 susp->noise_scale = noise_env->scale * CLAR_CONTROL_CHANGE_CONST;
246
247 /* make sure no sample rate is too high */
248 if (breath_env->sr > sr) {
249 sound_unref(breath_env);
250 snd_badsr();
251 } else if (breath_env->sr < sr) breath_env = snd_make_up(sr, breath_env);
252 if (freq_env->sr > sr) {
253 sound_unref(freq_env);
254 snd_badsr();
255 } else if (freq_env->sr < sr) freq_env = snd_make_up(sr, freq_env);
256 if (reed_stiffness->sr > sr) {
257 sound_unref(reed_stiffness);
258 snd_badsr();
259 } else if (reed_stiffness->sr < sr) reed_stiffness = snd_make_up(sr, reed_stiffness);
260 if (noise_env->sr > sr) {
261 sound_unref(noise_env);
262 snd_badsr();
263 } else if (noise_env->sr < sr) noise_env = snd_make_up(sr, noise_env);
264 susp->susp.fetch = clarinet_all_nsnn_fetch;
265 susp->terminate_cnt = UNKNOWN;
266 /* handle unequal start times, if any */
267 if (t0 < breath_env->t0) sound_prepend_zeros(breath_env, t0);
268 if (t0 < freq_env->t0) sound_prepend_zeros(freq_env, t0);
269 if (t0 < reed_stiffness->t0) sound_prepend_zeros(reed_stiffness, t0);
270 if (t0 < noise_env->t0) sound_prepend_zeros(noise_env, t0);
271 /* minimum start time over all inputs: */
272 t0_min = min(breath_env->t0, min(freq_env->t0, min(reed_stiffness->t0, min(noise_env->t0, t0))));
273 /* how many samples to toss before t0: */
274 susp->susp.toss_cnt = (long) ((t0 - t0_min) * sr + 0.5);
275 if (susp->susp.toss_cnt > 0) {
276 susp->susp.keep_fetch = susp->susp.fetch;
277 susp->susp.fetch = clarinet_all_toss_fetch;
278 }
279
280 /* initialize susp state */
281 susp->susp.free = clarinet_all_free;
282 susp->susp.sr = sr;
283 susp->susp.t0 = t0;
284 susp->susp.mark = clarinet_all_mark;
285 susp->susp.print_tree = clarinet_all_print_tree;
286 susp->susp.name = "clarinet_all";
287 susp->susp.log_stop_cnt = UNKNOWN;
288 susp->susp.current = 0;
289 susp->breath_env = breath_env;
290 susp->breath_env_cnt = 0;
291 susp->freq_env = freq_env;
292 susp->freq_env_cnt = 0;
293 susp->reed_stiffness = reed_stiffness;
294 susp->reed_stiffness_cnt = 0;
295 susp->noise_env = noise_env;
296 susp->noise_env_cnt = 0;
297 return sound_create((snd_susp_type)susp, t0, sr, scale_factor);
298 }
299
300
snd_clarinet_all(double freq,sound_type breath_env,sound_type freq_env,double vibrato_freq,double vibrato_gain,sound_type reed_stiffness,sound_type noise_env,rate_type sr)301 sound_type snd_clarinet_all(double freq, sound_type breath_env, sound_type freq_env, double vibrato_freq, double vibrato_gain, sound_type reed_stiffness, sound_type noise_env, rate_type sr)
302 {
303 sound_type breath_env_copy = sound_copy(breath_env);
304 sound_type freq_env_copy = sound_copy(freq_env);
305 sound_type reed_stiffness_copy = sound_copy(reed_stiffness);
306 sound_type noise_env_copy = sound_copy(noise_env);
307 return snd_make_clarinet_all(freq, breath_env_copy, freq_env_copy, vibrato_freq, vibrato_gain, reed_stiffness_copy, noise_env_copy, sr);
308 }
309