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