bounce.c 46 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889909192939495969798991001011021031041051061071081091101111121131141151161171181191201211221231241251261271281291301311321331341351361371381391401411421431441451461471481491501511521531541551561571581591601611621631641651661671681691701711721731741751761771781791801811821831841851861871881891901911921931941951961971981992002012022032042052062072082092102112122132142152162172182192202212222232242252262272282292302312322332342352362372382392402412422432442452462472482492502512522532542552562572582592602612622632642652662672682692702712722732742752762772782792802812822832842852862872882892902912922932942952962972982993003013023033043053063073083093103113123133143153163173183193203213223233243253263273283293303313323333343353363373383393403413423433443453463473483493503513523533543553563573583593603613623633643653663673683693703713723733743753763773783793803813823833843853863873883893903913923933943953963973983994004014024034044054064074084094104114124134144154164174184194204214224234244254264274284294304314324334344354364374384394404414424434444454464474484494504514524534544554564574584594604614624634644654664674684694704714724734744754764774784794804814824834844854864874884894904914924934944954964974984995005015025035045055065075085095105115125135145155165175185195205215225235245255265275285295305315325335345355365375385395405415425435445455465475485495505515525535545555565575585595605615625635645655665675685695705715725735745755765775785795805815825835845855865875885895905915925935945955965975985996006016026036046056066076086096106116126136146156166176186196206216226236246256266276286296306316326336346356366376386396406416426436446456466476486496506516526536546556566576586596606616626636646656666676686696706716726736746756766776786796806816826836846856866876886896906916926936946956966976986997007017027037047057067077087097107117127137147157167177187197207217227237247257267277287297307317327337347357367377387397407417427437447457467477487497507517527537547557567577587597607617627637647657667677687697707717727737747757767777787797807817827837847857867877887897907917927937947957967977987998008018028038048058068078088098108118128138148158168178188198208218228238248258268278288298308318328338348358368378388398408418428438448458468478488498508518528538548558568578588598608618628638648658668678688698708718728738748758768778788798808818828838848858868878888898908918928938948958968978988999009019029039049059069079089099109119129139149159169179189199209219229239249259269279289299309319329339349359369379389399409419429439449459469479489499509519529539549559569579589599609619629639649659669679689699709719729739749759769779789799809819829839849859869879889899909919929939949959969979989991000100110021003100410051006100710081009101010111012101310141015101610171018101910201021102210231024102510261027102810291030103110321033103410351036103710381039104010411042104310441045104610471048104910501051105210531054105510561057105810591060106110621063106410651066106710681069107010711072107310741075107610771078107910801081108210831084108510861087108810891090109110921093109410951096109710981099110011011102110311041105110611071108110911101111111211131114111511161117111811191120112111221123112411251126112711281129113011311132113311341135113611371138113911401141114211431144114511461147114811491150115111521153115411551156115711581159116011611162116311641165116611671168116911701171117211731174117511761177117811791180118111821183118411851186118711881189119011911192119311941195119611971198119912001201120212031204120512061207120812091210121112121213121412151216121712181219122012211222122312241225122612271228122912301231123212331234123512361237123812391240124112421243124412451246124712481249125012511252125312541255125612571258125912601261126212631264126512661267126812691270127112721273127412751276127712781279128012811282128312841285128612871288128912901291129212931294129512961297129812991300130113021303130413051306130713081309131013111312131313141315131613171318131913201321132213231324132513261327132813291330133113321333133413351336133713381339134013411342134313441345134613471348134913501351135213531354135513561357135813591360136113621363136413651366136713681369137013711372137313741375137613771378137913801381138213831384138513861387138813891390139113921393139413951396139713981399140014011402140314041405140614071408140914101411141214131414141514161417
  1. #include <stdio.h>
  2. #include <stdlib.h>
  3. #include <structures.h>
  4. #include <tkglobals.h>
  5. #include <pnames.h>
  6. #include <filetype.h>
  7. #include <processno.h>
  8. #include <modeno.h>
  9. #include <logic.h>
  10. #include <globcon.h>
  11. #include <cdpmain.h>
  12. #include <math.h>
  13. #include <mixxcon.h>
  14. #include <osbind.h>
  15. #include <standalone.h>
  16. #include <science.h>
  17. #include <ctype.h>
  18. #include <sfsys.h>
  19. #include <string.h>
  20. #include <srates.h>
  21. #define BOUNCESPLICE 3 // 3mS splices for cuts of src in bounce
  22. #define BOUNCEMINDUR 0.151 // Minimum length of bounce unit , to allow splice at start OR end
  23. #define CUTBNC 0
  24. #define KEEPEND 1
  25. #ifdef unix
  26. #define round(x) lround((x))
  27. #endif
  28. char errstr[2400];
  29. int anal_infiles = 1;
  30. int sloom = 0;
  31. int sloombatch = 0;
  32. const char* cdp_version = "6.1.0";
  33. //CDP LIB REPLACEMENTS
  34. static int check_bounce_param_validity_and_consistency(dataptr dz);
  35. static int setup_bounce_application(dataptr dz);
  36. static int parse_sloom_data(int argc,char *argv[],char ***cmdline,int *cmdlinecnt,dataptr dz);
  37. static int parse_infile_and_check_type(char **cmdline,dataptr dz);
  38. static int setup_bounce_param_ranges_and_defaults(dataptr dz);
  39. static int handle_the_outfile(int *cmdlinecnt,char ***cmdline,dataptr dz);
  40. static int setup_and_init_input_param_activity(dataptr dz,int tipc);
  41. static int setup_input_param_defaultval_stores(int tipc,aplptr ap);
  42. static int establish_application(dataptr dz);
  43. static int initialise_vflags(dataptr dz);
  44. static int setup_parameter_storage_and_constants(int storage_cnt,dataptr dz);
  45. static int initialise_is_int_and_no_brk_constants(int storage_cnt,dataptr dz);
  46. static int mark_parameter_types(dataptr dz,aplptr ap);
  47. static int assign_file_data_storage(int infilecnt,dataptr dz);
  48. static int get_tk_cmdline_word(int *cmdlinecnt,char ***cmdline,char *q);
  49. static int get_the_process_no(char *prog_identifier_from_cmdline,dataptr dz);
  50. static int create_bounce_sndbufs(dataptr dz);
  51. static int setup_and_init_input_brktable_constants(dataptr dz,int brkcnt);
  52. static int bounce_param_preprocess(dataptr dz);
  53. static int handle_buffer_output(int passno,double *maxsamp,int *maxpos,int fadelen,int *fadecnt,int *start_of_buf,int bnc1len,double normdrop,int samps_to_write,int *opos,dataptr dz);
  54. static int normalisation_fade(int fadelen, int *fadecnt, int start_of_buf, int absfadestt, double normdrop, dataptr dz);
  55. static int bounce(dataptr dz);
  56. /**************************************** MAIN *********************************************/
  57. int main(int argc,char *argv[])
  58. {
  59. int exit_status;
  60. dataptr dz = NULL;
  61. char **cmdline;
  62. int cmdlinecnt;
  63. int n;
  64. aplptr ap;
  65. int is_launched = FALSE;
  66. if(argc==2 && (strcmp(argv[1],"--version") == 0)) {
  67. fprintf(stdout,"%s\n",cdp_version);
  68. fflush(stdout);
  69. return 0;
  70. }
  71. /* CHECK FOR SOUNDLOOM */
  72. if((sloom = sound_loom_in_use(&argc,&argv)) > 1) {
  73. sloom = 0;
  74. sloombatch = 1;
  75. }
  76. if(sflinit("cdp")){
  77. sfperror("cdp: initialisation\n");
  78. return(FAILED);
  79. }
  80. /* SET UP THE PRINCIPLE DATASTRUCTURE */
  81. if((exit_status = establish_datastructure(&dz))<0) { // CDP LIB
  82. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  83. return(FAILED);
  84. }
  85. if(!sloom) {
  86. if(argc == 1) {
  87. usage1();
  88. return(FAILED);
  89. } else if(argc == 2) {
  90. usage2(argv[1]);
  91. return(FAILED);
  92. }
  93. }
  94. if(!sloom) {
  95. if((exit_status = make_initial_cmdline_check(&argc,&argv))<0) { // CDP LIB
  96. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  97. return(FAILED);
  98. }
  99. cmdline = argv;
  100. cmdlinecnt = argc;
  101. if((get_the_process_no(argv[0],dz))<0)
  102. return(FAILED);
  103. cmdline++;
  104. cmdlinecnt--;
  105. dz->maxmode = 0;
  106. // setup_particular_application =
  107. if((exit_status = setup_bounce_application(dz))<0) {
  108. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  109. return(FAILED);
  110. }
  111. if((exit_status = count_and_allocate_for_infiles(cmdlinecnt,cmdline,dz))<0) { // CDP LIB
  112. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  113. return(FAILED);
  114. }
  115. } else {
  116. //parse_TK_data() =
  117. if((exit_status = parse_sloom_data(argc,argv,&cmdline,&cmdlinecnt,dz))<0) {
  118. exit_status = print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  119. return(exit_status);
  120. }
  121. }
  122. ap = dz->application;
  123. // parse_infile_and_hone_type() =
  124. if((exit_status = parse_infile_and_check_type(cmdline,dz))<0) {
  125. exit_status = print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  126. return(FAILED);
  127. }
  128. // setup_param_ranges_and_defaults() =
  129. if((exit_status = setup_bounce_param_ranges_and_defaults(dz))<0) {
  130. exit_status = print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  131. return(FAILED);
  132. }
  133. // open_first_infile CDP LIB
  134. if((exit_status = open_first_infile(cmdline[0],dz))<0) {
  135. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  136. return(FAILED);
  137. }
  138. cmdlinecnt--;
  139. cmdline++;
  140. // handle_extra_infiles() : redundant
  141. // handle_outfile() =
  142. if((exit_status = handle_the_outfile(&cmdlinecnt,&cmdline,dz))<0) {
  143. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  144. return(FAILED);
  145. }
  146. // handle_formants() redundant
  147. // handle_formant_quiksearch() redundant
  148. // handle_special_data() redundant
  149. if((exit_status = read_parameters_and_flags(&cmdline,&cmdlinecnt,dz))<0) { // CDP LIB
  150. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  151. return(FAILED);
  152. }
  153. // check_param_validity_and_consistency....
  154. if((exit_status = check_bounce_param_validity_and_consistency(dz))<0) {
  155. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  156. return(FAILED);
  157. }
  158. is_launched = TRUE;
  159. dz->bufcnt = 3;
  160. if((dz->sampbuf = (float **)malloc(sizeof(float *) * (dz->bufcnt+1)))==NULL) {
  161. sprintf(errstr,"INSUFFICIENT MEMORY establishing sample buffers.\n");
  162. return(MEMORY_ERROR);
  163. }
  164. if((dz->sbufptr = (float **)malloc(sizeof(float *) * dz->bufcnt))==NULL) {
  165. sprintf(errstr,"INSUFFICIENT MEMORY establishing sample buffer pointers.\n");
  166. return(MEMORY_ERROR);
  167. }
  168. for(n = 0;n <dz->bufcnt; n++)
  169. dz->sampbuf[n] = dz->sbufptr[n] = (float *)0;
  170. dz->sampbuf[n] = (float *)0;
  171. if((exit_status = create_bounce_sndbufs(dz))<0) { // CDP LIB
  172. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  173. return(FAILED);
  174. }
  175. //param_preprocess() redundant
  176. if((exit_status = bounce_param_preprocess(dz))<0) {
  177. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  178. return(FAILED);
  179. }
  180. //spec_process_file =
  181. if((exit_status = bounce(dz))<0) {
  182. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  183. return(FAILED);
  184. }
  185. if((exit_status = complete_output(dz))<0) { // CDP LIB
  186. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  187. return(FAILED);
  188. }
  189. exit_status = print_messages_and_close_sndfiles(FINISHED,is_launched,dz); // CDP LIB
  190. free(dz);
  191. return(SUCCEEDED);
  192. }
  193. /**********************************************
  194. REPLACED CDP LIB FUNCTIONS
  195. **********************************************/
  196. /****************************** SET_PARAM_DATA *********************************/
  197. int set_param_data(aplptr ap, int special_data,int maxparamcnt,int paramcnt,char *paramlist)
  198. {
  199. ap->special_data = (char)special_data;
  200. ap->param_cnt = (char)paramcnt;
  201. ap->max_param_cnt = (char)maxparamcnt;
  202. if(ap->max_param_cnt>0) {
  203. if((ap->param_list = (char *)malloc((size_t)(ap->max_param_cnt+1)))==NULL) {
  204. sprintf(errstr,"INSUFFICIENT MEMORY: for param_list\n");
  205. return(MEMORY_ERROR);
  206. }
  207. strcpy(ap->param_list,paramlist);
  208. }
  209. return(FINISHED);
  210. }
  211. /****************************** SET_VFLGS *********************************/
  212. int set_vflgs
  213. (aplptr ap,char *optflags,int optcnt,char *optlist,char *varflags,int vflagcnt, int vparamcnt,char *varlist)
  214. {
  215. ap->option_cnt = (char) optcnt; /*RWD added cast */
  216. if(optcnt) {
  217. if((ap->option_list = (char *)malloc((size_t)(optcnt+1)))==NULL) {
  218. sprintf(errstr,"INSUFFICIENT MEMORY: for option_list\n");
  219. return(MEMORY_ERROR);
  220. }
  221. strcpy(ap->option_list,optlist);
  222. if((ap->option_flags = (char *)malloc((size_t)(optcnt+1)))==NULL) {
  223. sprintf(errstr,"INSUFFICIENT MEMORY: for option_flags\n");
  224. return(MEMORY_ERROR);
  225. }
  226. strcpy(ap->option_flags,optflags);
  227. }
  228. ap->vflag_cnt = (char) vflagcnt;
  229. ap->variant_param_cnt = (char) vparamcnt;
  230. if(vflagcnt) {
  231. if((ap->variant_list = (char *)malloc((size_t)(vflagcnt+1)))==NULL) {
  232. sprintf(errstr,"INSUFFICIENT MEMORY: for variant_list\n");
  233. return(MEMORY_ERROR);
  234. }
  235. strcpy(ap->variant_list,varlist);
  236. if((ap->variant_flags = (char *)malloc((size_t)(vflagcnt+1)))==NULL) {
  237. sprintf(errstr,"INSUFFICIENT MEMORY: for variant_flags\n");
  238. return(MEMORY_ERROR);
  239. }
  240. strcpy(ap->variant_flags,varflags);
  241. }
  242. return(FINISHED);
  243. }
  244. /***************************** APPLICATION_INIT **************************/
  245. int application_init(dataptr dz)
  246. {
  247. int exit_status;
  248. int storage_cnt;
  249. int tipc, brkcnt;
  250. aplptr ap = dz->application;
  251. if(ap->vflag_cnt>0)
  252. initialise_vflags(dz);
  253. tipc = ap->max_param_cnt + ap->option_cnt + ap->variant_param_cnt;
  254. ap->total_input_param_cnt = (char)tipc;
  255. if(tipc>0) {
  256. if((exit_status = setup_input_param_range_stores(tipc,ap))<0)
  257. return(exit_status);
  258. if((exit_status = setup_input_param_defaultval_stores(tipc,ap))<0)
  259. return(exit_status);
  260. if((exit_status = setup_and_init_input_param_activity(dz,tipc))<0)
  261. return(exit_status);
  262. }
  263. brkcnt = tipc;
  264. //THERE ARE NO INPUTFILE brktables USED IN THIS PROCESS
  265. if(brkcnt>0) {
  266. if((exit_status = setup_and_init_input_brktable_constants(dz,brkcnt))<0)
  267. return(exit_status);
  268. }
  269. if((storage_cnt = tipc + ap->internal_param_cnt)>0) {
  270. if((exit_status = setup_parameter_storage_and_constants(storage_cnt,dz))<0)
  271. return(exit_status);
  272. if((exit_status = initialise_is_int_and_no_brk_constants(storage_cnt,dz))<0)
  273. return(exit_status);
  274. }
  275. if((exit_status = mark_parameter_types(dz,ap))<0)
  276. return(exit_status);
  277. // establish_infile_constants() replaced by
  278. dz->infilecnt = 1;
  279. //establish_bufptrs_and_extra_buffers():
  280. return(FINISHED);
  281. }
  282. /********************** SETUP_PARAMETER_STORAGE_AND_CONSTANTS ********************/
  283. /* RWD mallo changed to calloc; helps debug verison run as release! */
  284. int setup_parameter_storage_and_constants(int storage_cnt,dataptr dz)
  285. {
  286. if((dz->param = (double *)calloc(storage_cnt, sizeof(double)))==NULL) {
  287. sprintf(errstr,"setup_parameter_storage_and_constants(): 1\n");
  288. return(MEMORY_ERROR);
  289. }
  290. if((dz->iparam = (int *)calloc(storage_cnt, sizeof(int) ))==NULL) {
  291. sprintf(errstr,"setup_parameter_storage_and_constants(): 2\n");
  292. return(MEMORY_ERROR);
  293. }
  294. if((dz->is_int = (char *)calloc(storage_cnt, sizeof(char)))==NULL) {
  295. sprintf(errstr,"setup_parameter_storage_and_constants(): 3\n");
  296. return(MEMORY_ERROR);
  297. }
  298. if((dz->no_brk = (char *)calloc(storage_cnt, sizeof(char)))==NULL) {
  299. sprintf(errstr,"setup_parameter_storage_and_constants(): 5\n");
  300. return(MEMORY_ERROR);
  301. }
  302. return(FINISHED);
  303. }
  304. /************** INITIALISE_IS_INT_AND_NO_BRK_CONSTANTS *****************/
  305. int initialise_is_int_and_no_brk_constants(int storage_cnt,dataptr dz)
  306. {
  307. int n;
  308. for(n=0;n<storage_cnt;n++) {
  309. dz->is_int[n] = (char)0;
  310. dz->no_brk[n] = (char)0;
  311. }
  312. return(FINISHED);
  313. }
  314. /***************************** MARK_PARAMETER_TYPES **************************/
  315. int mark_parameter_types(dataptr dz,aplptr ap)
  316. {
  317. int n, m; /* PARAMS */
  318. for(n=0;n<ap->max_param_cnt;n++) {
  319. switch(ap->param_list[n]) {
  320. case('0'): break; /* dz->is_active[n] = 0 is default */
  321. case('i'): dz->is_active[n] = (char)1; dz->is_int[n] = (char)1;dz->no_brk[n] = (char)1; break;
  322. case('I'): dz->is_active[n] = (char)1; dz->is_int[n] = (char)1; break;
  323. case('d'): dz->is_active[n] = (char)1; dz->no_brk[n] = (char)1; break;
  324. case('D'): dz->is_active[n] = (char)1; /* normal case: double val or brkpnt file */ break;
  325. default:
  326. sprintf(errstr,"Programming error: invalid parameter type in mark_parameter_types()\n");
  327. return(PROGRAM_ERROR);
  328. }
  329. } /* OPTIONS */
  330. for(n=0,m=ap->max_param_cnt;n<ap->option_cnt;n++,m++) {
  331. switch(ap->option_list[n]) {
  332. case('i'): dz->is_active[m] = (char)1; dz->is_int[m] = (char)1; dz->no_brk[m] = (char)1; break;
  333. case('I'): dz->is_active[m] = (char)1; dz->is_int[m] = (char)1; break;
  334. case('d'): dz->is_active[m] = (char)1; dz->no_brk[m] = (char)1; break;
  335. case('D'): dz->is_active[m] = (char)1; /* normal case: double val or brkpnt file */ break;
  336. default:
  337. sprintf(errstr,"Programming error: invalid option type in mark_parameter_types()\n");
  338. return(PROGRAM_ERROR);
  339. }
  340. } /* VARIANTS */
  341. for(n=0,m=ap->max_param_cnt + ap->option_cnt;n < ap->variant_param_cnt; n++, m++) {
  342. switch(ap->variant_list[n]) {
  343. case('0'): break;
  344. case('i'): dz->is_active[m] = (char)1; dz->is_int[m] = (char)1; dz->no_brk[m] = (char)1; break;
  345. case('I'): dz->is_active[m] = (char)1; dz->is_int[m] = (char)1; break;
  346. case('d'): dz->is_active[m] = (char)1; dz->no_brk[m] = (char)1; break;
  347. case('D'): dz->is_active[m] = (char)1; /* normal case: double val or brkpnt file */ break;
  348. default:
  349. sprintf(errstr,"Programming error: invalid variant type in mark_parameter_types()\n");
  350. return(PROGRAM_ERROR);
  351. }
  352. } /* INTERNAL */
  353. for(n=0,
  354. m=ap->max_param_cnt + ap->option_cnt + ap->variant_param_cnt; n<ap->internal_param_cnt; n++,m++) {
  355. switch(ap->internal_param_list[n]) {
  356. case('0'): break; /* dummy variables: variables not used: but important for internal paream numbering!! */
  357. case('i'): dz->is_int[m] = (char)1; dz->no_brk[m] = (char)1; break;
  358. case('d'): dz->no_brk[m] = (char)1; break;
  359. default:
  360. sprintf(errstr,"Programming error: invalid internal param type in mark_parameter_types()\n");
  361. return(PROGRAM_ERROR);
  362. }
  363. }
  364. return(FINISHED);
  365. }
  366. /************************ HANDLE_THE_OUTFILE *********************/
  367. int handle_the_outfile(int *cmdlinecnt,char ***cmdline,dataptr dz)
  368. {
  369. int exit_status;
  370. char *filename = (*cmdline)[0];
  371. if(filename[0]=='-' && filename[1]=='f') {
  372. dz->floatsam_output = 1;
  373. dz->true_outfile_stype = SAMP_FLOAT;
  374. filename+= 2;
  375. }
  376. if(!sloom) {
  377. if(file_has_invalid_startchar(filename) || value_is_numeric(filename)) {
  378. sprintf(errstr,"Outfile name %s has invalid start character(s) or looks too much like a number.\n",filename);
  379. return(DATA_ERROR);
  380. }
  381. }
  382. strcpy(dz->outfilename,filename);
  383. if((exit_status = create_sized_outfile(filename,dz))<0)
  384. return(exit_status);
  385. (*cmdline)++;
  386. (*cmdlinecnt)--;
  387. return(FINISHED);
  388. }
  389. /***************************** ESTABLISH_APPLICATION **************************/
  390. int establish_application(dataptr dz)
  391. {
  392. aplptr ap;
  393. if((dz->application = (aplptr)malloc(sizeof (struct applic)))==NULL) {
  394. sprintf(errstr,"establish_application()\n");
  395. return(MEMORY_ERROR);
  396. }
  397. ap = dz->application;
  398. memset((char *)ap,0,sizeof(struct applic));
  399. return(FINISHED);
  400. }
  401. /************************* INITIALISE_VFLAGS *************************/
  402. int initialise_vflags(dataptr dz)
  403. {
  404. int n;
  405. if((dz->vflag = (char *)malloc(dz->application->vflag_cnt * sizeof(char)))==NULL) {
  406. sprintf(errstr,"INSUFFICIENT MEMORY: vflag store,\n");
  407. return(MEMORY_ERROR);
  408. }
  409. for(n=0;n<dz->application->vflag_cnt;n++)
  410. dz->vflag[n] = FALSE;
  411. return FINISHED;
  412. }
  413. /************************* SETUP_INPUT_PARAM_DEFAULTVALS *************************/
  414. int setup_input_param_defaultval_stores(int tipc,aplptr ap)
  415. {
  416. int n;
  417. if((ap->default_val = (double *)malloc(tipc * sizeof(double)))==NULL) {
  418. sprintf(errstr,"INSUFFICIENT MEMORY for application default values store\n");
  419. return(MEMORY_ERROR);
  420. }
  421. for(n=0;n<tipc;n++)
  422. ap->default_val[n] = 0.0;
  423. return(FINISHED);
  424. }
  425. /***************************** SETUP_AND_INIT_INPUT_PARAM_ACTIVITY **************************/
  426. int setup_and_init_input_param_activity(dataptr dz,int tipc)
  427. {
  428. int n;
  429. if((dz->is_active = (char *)malloc((size_t)tipc))==NULL) {
  430. sprintf(errstr,"setup_and_init_input_param_activity()\n");
  431. return(MEMORY_ERROR);
  432. }
  433. for(n=0;n<tipc;n++)
  434. dz->is_active[n] = (char)0;
  435. return(FINISHED);
  436. }
  437. /************************* SETUP_BOUNCE_APPLICATION *******************/
  438. int setup_bounce_application(dataptr dz)
  439. {
  440. int exit_status;
  441. aplptr ap;
  442. if((exit_status = establish_application(dz))<0) // GLOBAL
  443. return(FAILED);
  444. ap = dz->application;
  445. // SEE parstruct FOR EXPLANATION of next 2 functions
  446. if((exit_status = set_param_data(ap,0 ,5,5,"idddd"))<0)
  447. return(FAILED);
  448. if((exit_status = set_vflgs(ap,"s",1,"d","ce",2,0,"00"))<0)
  449. return(FAILED);
  450. // set_legal_infile_structure -->
  451. dz->has_otherfile = FALSE;
  452. // assign_process_logic -->
  453. dz->input_data_type = SNDFILES_ONLY;
  454. dz->process_type = UNEQUAL_SNDFILE;
  455. dz->outfiletype = SNDFILE_OUT;
  456. return application_init(dz); //GLOBAL
  457. }
  458. /************************* PARSE_INFILE_AND_CHECK_TYPE *******************/
  459. int parse_infile_and_check_type(char **cmdline,dataptr dz)
  460. {
  461. int exit_status;
  462. infileptr infile_info;
  463. if(!sloom) {
  464. if((infile_info = (infileptr)malloc(sizeof(struct filedata)))==NULL) {
  465. sprintf(errstr,"INSUFFICIENT MEMORY for infile structure to test file data.");
  466. return(MEMORY_ERROR);
  467. } else if((exit_status = cdparse(cmdline[0],infile_info))<0) {
  468. sprintf(errstr,"Failed to parse input file %s\n",cmdline[0]);
  469. return(PROGRAM_ERROR);
  470. } else if(infile_info->filetype != SNDFILE) {
  471. sprintf(errstr,"File %s is not of correct type\n",cmdline[0]);
  472. return(DATA_ERROR);
  473. } else if((exit_status = copy_parse_info_to_main_structure(infile_info,dz))<0) {
  474. sprintf(errstr,"Failed to copy file parsing information\n");
  475. return(PROGRAM_ERROR);
  476. }
  477. free(infile_info);
  478. }
  479. return(FINISHED);
  480. }
  481. /************************* SETUP_BOUNCE_PARAM_RANGES_AND_DEFAULTS *******************/
  482. int setup_bounce_param_ranges_and_defaults(dataptr dz)
  483. {
  484. int exit_status;
  485. aplptr ap = dz->application;
  486. // set_param_ranges()
  487. ap->total_input_param_cnt = (char)(ap->max_param_cnt + ap->option_cnt + ap->variant_param_cnt);
  488. // NB total_input_param_cnt is > 0 !!!
  489. if((exit_status = setup_input_param_range_stores(ap->total_input_param_cnt,ap))<0)
  490. return(FAILED);
  491. // get_param_ranges()
  492. ap->lo[BNC_NUMBER] = 1;
  493. ap->hi[BNC_NUMBER] = 100;
  494. ap->default_val[BNC_NUMBER] = 8;
  495. ap->lo[BNC_STTSTEP] = 0.04;
  496. ap->hi[BNC_STTSTEP] = 10.0;
  497. ap->default_val[BNC_STTSTEP] = 1.0;
  498. ap->lo[BNC_SHORTEN] = 0.1;
  499. ap->hi[BNC_SHORTEN] = 1.0;
  500. ap->default_val[BNC_SHORTEN] = 0.92;
  501. ap->lo[BNC_ENDLEV] = 0.0;
  502. ap->hi[BNC_ENDLEV] = 1.0;
  503. ap->default_val[BNC_ENDLEV] = 0.05;
  504. ap->lo[BNC_LEVWRP] = .01;
  505. ap->hi[BNC_LEVWRP] = 100.0;
  506. ap->default_val[BNC_LEVWRP] = 1.0;
  507. ap->lo[BNC_MINDUR] = 0.0;
  508. ap->hi[BNC_MINDUR] = 1.0;
  509. ap->default_val[BNC_MINDUR] = 0.02;
  510. dz->maxmode = 0;
  511. if(!sloom)
  512. put_default_vals_in_all_params(dz);
  513. return(FINISHED);
  514. }
  515. /********************************* PARSE_SLOOM_DATA *********************************/
  516. int parse_sloom_data(int argc,char *argv[],char ***cmdline,int *cmdlinecnt,dataptr dz)
  517. {
  518. int exit_status;
  519. int cnt = 1, infilecnt;
  520. int filesize, insams, inbrksize;
  521. double dummy;
  522. int true_cnt = 0;
  523. aplptr ap;
  524. while(cnt<=PRE_CMDLINE_DATACNT) {
  525. if(cnt > argc) {
  526. sprintf(errstr,"Insufficient data sent from TK\n");
  527. return(DATA_ERROR);
  528. }
  529. switch(cnt) {
  530. case(1):
  531. if(sscanf(argv[cnt],"%d",&dz->process)!=1) {
  532. sprintf(errstr,"Cannot read process no. sent from TK\n");
  533. return(DATA_ERROR);
  534. }
  535. break;
  536. case(2):
  537. if(sscanf(argv[cnt],"%d",&dz->mode)!=1) {
  538. sprintf(errstr,"Cannot read mode no. sent from TK\n");
  539. return(DATA_ERROR);
  540. }
  541. if(dz->mode > 0)
  542. dz->mode--;
  543. //setup_particular_application() =
  544. if((exit_status = setup_bounce_application(dz))<0)
  545. return(exit_status);
  546. ap = dz->application;
  547. break;
  548. case(3):
  549. if(sscanf(argv[cnt],"%d",&infilecnt)!=1) {
  550. sprintf(errstr,"Cannot read infilecnt sent from TK\n");
  551. return(DATA_ERROR);
  552. }
  553. if(infilecnt < 1) {
  554. true_cnt = cnt + 1;
  555. cnt = PRE_CMDLINE_DATACNT; /* force exit from loop after assign_file_data_storage */
  556. }
  557. if((exit_status = assign_file_data_storage(infilecnt,dz))<0)
  558. return(exit_status);
  559. break;
  560. case(INPUT_FILETYPE+4):
  561. if(sscanf(argv[cnt],"%d",&dz->infile->filetype)!=1) {
  562. sprintf(errstr,"Cannot read filetype sent from TK (%s)\n",argv[cnt]);
  563. return(DATA_ERROR);
  564. }
  565. break;
  566. case(INPUT_FILESIZE+4):
  567. if(sscanf(argv[cnt],"%d",&filesize)!=1) {
  568. sprintf(errstr,"Cannot read infilesize sent from TK\n");
  569. return(DATA_ERROR);
  570. }
  571. dz->insams[0] = filesize;
  572. break;
  573. case(INPUT_INSAMS+4):
  574. if(sscanf(argv[cnt],"%d",&insams)!=1) {
  575. sprintf(errstr,"Cannot read insams sent from TK\n");
  576. return(DATA_ERROR);
  577. }
  578. dz->insams[0] = insams;
  579. break;
  580. case(INPUT_SRATE+4):
  581. if(sscanf(argv[cnt],"%d",&dz->infile->srate)!=1) {
  582. sprintf(errstr,"Cannot read srate sent from TK\n");
  583. return(DATA_ERROR);
  584. }
  585. break;
  586. case(INPUT_CHANNELS+4):
  587. if(sscanf(argv[cnt],"%d",&dz->infile->channels)!=1) {
  588. sprintf(errstr,"Cannot read channels sent from TK\n");
  589. return(DATA_ERROR);
  590. }
  591. break;
  592. case(INPUT_STYPE+4):
  593. if(sscanf(argv[cnt],"%d",&dz->infile->stype)!=1) {
  594. sprintf(errstr,"Cannot read stype sent from TK\n");
  595. return(DATA_ERROR);
  596. }
  597. break;
  598. case(INPUT_ORIGSTYPE+4):
  599. if(sscanf(argv[cnt],"%d",&dz->infile->origstype)!=1) {
  600. sprintf(errstr,"Cannot read origstype sent from TK\n");
  601. return(DATA_ERROR);
  602. }
  603. break;
  604. case(INPUT_ORIGRATE+4):
  605. if(sscanf(argv[cnt],"%d",&dz->infile->origrate)!=1) {
  606. sprintf(errstr,"Cannot read origrate sent from TK\n");
  607. return(DATA_ERROR);
  608. }
  609. break;
  610. case(INPUT_MLEN+4):
  611. if(sscanf(argv[cnt],"%d",&dz->infile->Mlen)!=1) {
  612. sprintf(errstr,"Cannot read Mlen sent from TK\n");
  613. return(DATA_ERROR);
  614. }
  615. break;
  616. case(INPUT_DFAC+4):
  617. if(sscanf(argv[cnt],"%d",&dz->infile->Dfac)!=1) {
  618. sprintf(errstr,"Cannot read Dfac sent from TK\n");
  619. return(DATA_ERROR);
  620. }
  621. break;
  622. case(INPUT_ORIGCHANS+4):
  623. if(sscanf(argv[cnt],"%d",&dz->infile->origchans)!=1) {
  624. sprintf(errstr,"Cannot read origchans sent from TK\n");
  625. return(DATA_ERROR);
  626. }
  627. break;
  628. case(INPUT_SPECENVCNT+4):
  629. if(sscanf(argv[cnt],"%d",&dz->infile->specenvcnt)!=1) {
  630. sprintf(errstr,"Cannot read specenvcnt sent from TK\n");
  631. return(DATA_ERROR);
  632. }
  633. dz->specenvcnt = dz->infile->specenvcnt;
  634. break;
  635. case(INPUT_WANTED+4):
  636. if(sscanf(argv[cnt],"%d",&dz->wanted)!=1) {
  637. sprintf(errstr,"Cannot read wanted sent from TK\n");
  638. return(DATA_ERROR);
  639. }
  640. break;
  641. case(INPUT_WLENGTH+4):
  642. if(sscanf(argv[cnt],"%d",&dz->wlength)!=1) {
  643. sprintf(errstr,"Cannot read wlength sent from TK\n");
  644. return(DATA_ERROR);
  645. }
  646. break;
  647. case(INPUT_OUT_CHANS+4):
  648. if(sscanf(argv[cnt],"%d",&dz->out_chans)!=1) {
  649. sprintf(errstr,"Cannot read out_chans sent from TK\n");
  650. return(DATA_ERROR);
  651. }
  652. break;
  653. /* RWD these chanegs to samps - tk will have to deal with that! */
  654. case(INPUT_DESCRIPTOR_BYTES+4):
  655. if(sscanf(argv[cnt],"%d",&dz->descriptor_samps)!=1) {
  656. sprintf(errstr,"Cannot read descriptor_samps sent from TK\n");
  657. return(DATA_ERROR);
  658. }
  659. break;
  660. case(INPUT_IS_TRANSPOS+4):
  661. if(sscanf(argv[cnt],"%d",&dz->is_transpos)!=1) {
  662. sprintf(errstr,"Cannot read is_transpos sent from TK\n");
  663. return(DATA_ERROR);
  664. }
  665. break;
  666. case(INPUT_COULD_BE_TRANSPOS+4):
  667. if(sscanf(argv[cnt],"%d",&dz->could_be_transpos)!=1) {
  668. sprintf(errstr,"Cannot read could_be_transpos sent from TK\n");
  669. return(DATA_ERROR);
  670. }
  671. break;
  672. case(INPUT_COULD_BE_PITCH+4):
  673. if(sscanf(argv[cnt],"%d",&dz->could_be_pitch)!=1) {
  674. sprintf(errstr,"Cannot read could_be_pitch sent from TK\n");
  675. return(DATA_ERROR);
  676. }
  677. break;
  678. case(INPUT_DIFFERENT_SRATES+4):
  679. if(sscanf(argv[cnt],"%d",&dz->different_srates)!=1) {
  680. sprintf(errstr,"Cannot read different_srates sent from TK\n");
  681. return(DATA_ERROR);
  682. }
  683. break;
  684. case(INPUT_DUPLICATE_SNDS+4):
  685. if(sscanf(argv[cnt],"%d",&dz->duplicate_snds)!=1) {
  686. sprintf(errstr,"Cannot read duplicate_snds sent from TK\n");
  687. return(DATA_ERROR);
  688. }
  689. break;
  690. case(INPUT_BRKSIZE+4):
  691. if(sscanf(argv[cnt],"%d",&inbrksize)!=1) {
  692. sprintf(errstr,"Cannot read brksize sent from TK\n");
  693. return(DATA_ERROR);
  694. }
  695. if(inbrksize > 0) {
  696. switch(dz->input_data_type) {
  697. case(WORDLIST_ONLY):
  698. break;
  699. case(PITCH_AND_PITCH):
  700. case(PITCH_AND_TRANSPOS):
  701. case(TRANSPOS_AND_TRANSPOS):
  702. dz->tempsize = inbrksize;
  703. break;
  704. case(BRKFILES_ONLY):
  705. case(UNRANGED_BRKFILE_ONLY):
  706. case(DB_BRKFILES_ONLY):
  707. case(ALL_FILES):
  708. case(ANY_NUMBER_OF_ANY_FILES):
  709. if(dz->extrabrkno < 0) {
  710. sprintf(errstr,"Storage location number for brktable not established by CDP.\n");
  711. return(DATA_ERROR);
  712. }
  713. if(dz->brksize == NULL) {
  714. sprintf(errstr,"CDP has not established storage space for input brktable.\n");
  715. return(PROGRAM_ERROR);
  716. }
  717. dz->brksize[dz->extrabrkno] = inbrksize;
  718. break;
  719. default:
  720. sprintf(errstr,"TK sent brktablesize > 0 for input_data_type [%d] not using brktables.\n",
  721. dz->input_data_type);
  722. return(PROGRAM_ERROR);
  723. }
  724. break;
  725. }
  726. break;
  727. case(INPUT_NUMSIZE+4):
  728. if(sscanf(argv[cnt],"%d",&dz->numsize)!=1) {
  729. sprintf(errstr,"Cannot read numsize sent from TK\n");
  730. return(DATA_ERROR);
  731. }
  732. break;
  733. case(INPUT_LINECNT+4):
  734. if(sscanf(argv[cnt],"%d",&dz->linecnt)!=1) {
  735. sprintf(errstr,"Cannot read linecnt sent from TK\n");
  736. return(DATA_ERROR);
  737. }
  738. break;
  739. case(INPUT_ALL_WORDS+4):
  740. if(sscanf(argv[cnt],"%d",&dz->all_words)!=1) {
  741. sprintf(errstr,"Cannot read all_words sent from TK\n");
  742. return(DATA_ERROR);
  743. }
  744. break;
  745. case(INPUT_ARATE+4):
  746. if(sscanf(argv[cnt],"%f",&dz->infile->arate)!=1) {
  747. sprintf(errstr,"Cannot read arate sent from TK\n");
  748. return(DATA_ERROR);
  749. }
  750. break;
  751. case(INPUT_FRAMETIME+4):
  752. if(sscanf(argv[cnt],"%lf",&dummy)!=1) {
  753. sprintf(errstr,"Cannot read frametime sent from TK\n");
  754. return(DATA_ERROR);
  755. }
  756. dz->frametime = (float)dummy;
  757. break;
  758. case(INPUT_WINDOW_SIZE+4):
  759. if(sscanf(argv[cnt],"%f",&dz->infile->window_size)!=1) {
  760. sprintf(errstr,"Cannot read window_size sent from TK\n");
  761. return(DATA_ERROR);
  762. }
  763. break;
  764. case(INPUT_NYQUIST+4):
  765. if(sscanf(argv[cnt],"%lf",&dz->nyquist)!=1) {
  766. sprintf(errstr,"Cannot read nyquist sent from TK\n");
  767. return(DATA_ERROR);
  768. }
  769. break;
  770. case(INPUT_DURATION+4):
  771. if(sscanf(argv[cnt],"%lf",&dz->duration)!=1) {
  772. sprintf(errstr,"Cannot read duration sent from TK\n");
  773. return(DATA_ERROR);
  774. }
  775. break;
  776. case(INPUT_MINBRK+4):
  777. if(sscanf(argv[cnt],"%lf",&dz->minbrk)!=1) {
  778. sprintf(errstr,"Cannot read minbrk sent from TK\n");
  779. return(DATA_ERROR);
  780. }
  781. break;
  782. case(INPUT_MAXBRK+4):
  783. if(sscanf(argv[cnt],"%lf",&dz->maxbrk)!=1) {
  784. sprintf(errstr,"Cannot read maxbrk sent from TK\n");
  785. return(DATA_ERROR);
  786. }
  787. break;
  788. case(INPUT_MINNUM+4):
  789. if(sscanf(argv[cnt],"%lf",&dz->minnum)!=1) {
  790. sprintf(errstr,"Cannot read minnum sent from TK\n");
  791. return(DATA_ERROR);
  792. }
  793. break;
  794. case(INPUT_MAXNUM+4):
  795. if(sscanf(argv[cnt],"%lf",&dz->maxnum)!=1) {
  796. sprintf(errstr,"Cannot read maxnum sent from TK\n");
  797. return(DATA_ERROR);
  798. }
  799. break;
  800. default:
  801. sprintf(errstr,"case switch item missing: parse_sloom_data()\n");
  802. return(PROGRAM_ERROR);
  803. }
  804. cnt++;
  805. }
  806. if(cnt!=PRE_CMDLINE_DATACNT+1) {
  807. sprintf(errstr,"Insufficient pre-cmdline params sent from TK\n");
  808. return(DATA_ERROR);
  809. }
  810. if(true_cnt)
  811. cnt = true_cnt;
  812. *cmdlinecnt = 0;
  813. while(cnt < argc) {
  814. if((exit_status = get_tk_cmdline_word(cmdlinecnt,cmdline,argv[cnt]))<0)
  815. return(exit_status);
  816. cnt++;
  817. }
  818. return(FINISHED);
  819. }
  820. /********************************* GET_TK_CMDLINE_WORD *********************************/
  821. int get_tk_cmdline_word(int *cmdlinecnt,char ***cmdline,char *q)
  822. {
  823. if(*cmdlinecnt==0) {
  824. if((*cmdline = (char **)malloc(sizeof(char *)))==NULL) {
  825. sprintf(errstr,"INSUFFICIENT MEMORY for TK cmdline array.\n");
  826. return(MEMORY_ERROR);
  827. }
  828. } else {
  829. if((*cmdline = (char **)realloc(*cmdline,((*cmdlinecnt)+1) * sizeof(char *)))==NULL) {
  830. sprintf(errstr,"INSUFFICIENT MEMORY for TK cmdline array.\n");
  831. return(MEMORY_ERROR);
  832. }
  833. }
  834. if(((*cmdline)[*cmdlinecnt] = (char *)malloc((strlen(q) + 1) * sizeof(char)))==NULL) {
  835. sprintf(errstr,"INSUFFICIENT MEMORY for TK cmdline item %d.\n",(*cmdlinecnt)+1);
  836. return(MEMORY_ERROR);
  837. }
  838. strcpy((*cmdline)[*cmdlinecnt],q);
  839. (*cmdlinecnt)++;
  840. return(FINISHED);
  841. }
  842. /****************************** ASSIGN_FILE_DATA_STORAGE *********************************/
  843. int assign_file_data_storage(int infilecnt,dataptr dz)
  844. {
  845. int exit_status;
  846. int no_sndfile_system_files = FALSE;
  847. dz->infilecnt = infilecnt;
  848. if((exit_status = allocate_filespace(dz))<0)
  849. return(exit_status);
  850. if(no_sndfile_system_files)
  851. dz->infilecnt = 0;
  852. return(FINISHED);
  853. }
  854. /************************* redundant functions: to ensure libs compile OK *******************/
  855. int assign_process_logic(dataptr dz)
  856. {
  857. return(FINISHED);
  858. }
  859. void set_legal_infile_structure(dataptr dz)
  860. {}
  861. int set_legal_internalparam_structure(int process,int mode,aplptr ap)
  862. {
  863. return(FINISHED);
  864. }
  865. int setup_internal_arrays_and_array_pointers(dataptr dz)
  866. {
  867. return(FINISHED);
  868. }
  869. int establish_bufptrs_and_extra_buffers(dataptr dz)
  870. {
  871. return(FINISHED);
  872. }
  873. int read_special_data(char *str,dataptr dz)
  874. {
  875. return(FINISHED);
  876. }
  877. int inner_loop
  878. (int *peakscore,int *descnt,int *in_start_portion,int *least,int *pitchcnt,int windows_in_buf,dataptr dz)
  879. {
  880. return(FINISHED);
  881. }
  882. int get_process_no(char *prog_identifier_from_cmdline,dataptr dz)
  883. {
  884. return(FINISHED);
  885. }
  886. /******************************** USAGE1 ********************************/
  887. int usage1(void)
  888. {
  889. usage2("bounce");
  890. return(USAGE_ONLY);
  891. }
  892. /**************************** CHECK_BOUNCE_PARAM_VALIDITY_AND_CONSISTENCY *****************************/
  893. int check_bounce_param_validity_and_consistency(dataptr dz)
  894. {
  895. return FINISHED;
  896. }
  897. /******************************** DBTOLEVEL ***********************/
  898. double dbtolevel(double val)
  899. {
  900. int isneg = 0;
  901. if(flteq(val,0.0))
  902. return(1.0);
  903. if(val < 0.0) {
  904. val = -val;
  905. isneg = 1;
  906. }
  907. val /= 20.0;
  908. val = pow(10.0,val);
  909. if(isneg)
  910. val = 1.0/val;
  911. return(val);
  912. }
  913. /********************************************************************************************/
  914. int get_the_process_no(char *prog_identifier_from_cmdline,dataptr dz)
  915. {
  916. if(!strcmp(prog_identifier_from_cmdline,"bounce")) dz->process = BOUNCE;
  917. else {
  918. sprintf(errstr,"Unknown program identification string '%s'\n",prog_identifier_from_cmdline);
  919. return(USAGE_ONLY);
  920. }
  921. return(FINISHED);
  922. }
  923. /******************************** SETUP_AND_INIT_INPUT_BRKTABLE_CONSTANTS ********************************/
  924. int setup_and_init_input_brktable_constants(dataptr dz,int brkcnt)
  925. {
  926. int n;
  927. if((dz->brk = (double **)malloc(brkcnt * sizeof(double *)))==NULL) {
  928. sprintf(errstr,"setup_and_init_input_brktable_constants(): 1\n");
  929. return(MEMORY_ERROR);
  930. }
  931. if((dz->brkptr = (double **)malloc(brkcnt * sizeof(double *)))==NULL) {
  932. sprintf(errstr,"setup_and_init_input_brktable_constants(): 6\n");
  933. return(MEMORY_ERROR);
  934. }
  935. if((dz->brksize = (int *)malloc(brkcnt * sizeof(int)))==NULL) {
  936. sprintf(errstr,"setup_and_init_input_brktable_constants(): 2\n");
  937. return(MEMORY_ERROR);
  938. }
  939. if((dz->firstval = (double *)malloc(brkcnt * sizeof(double)))==NULL) {
  940. sprintf(errstr,"setup_and_init_input_brktable_constants(): 3\n");
  941. return(MEMORY_ERROR);
  942. }
  943. if((dz->lastind = (double *)malloc(brkcnt * sizeof(double)))==NULL) {
  944. sprintf(errstr,"setup_and_init_input_brktable_constants(): 4\n");
  945. return(MEMORY_ERROR);
  946. }
  947. if((dz->lastval = (double *)malloc(brkcnt * sizeof(double)))==NULL) {
  948. sprintf(errstr,"setup_and_init_input_brktable_constants(): 5\n");
  949. return(MEMORY_ERROR);
  950. }
  951. if((dz->brkinit = (int *)malloc(brkcnt * sizeof(int)))==NULL) {
  952. sprintf(errstr,"setup_and_init_input_brktable_constants(): 7\n");
  953. return(MEMORY_ERROR);
  954. }
  955. for(n=0;n<brkcnt;n++) {
  956. dz->brk[n] = NULL;
  957. dz->brkptr[n] = NULL;
  958. dz->brkinit[n] = 0;
  959. dz->brksize[n] = 0;
  960. }
  961. return(FINISHED);
  962. }
  963. /******************************** USAGE2 ********************************/
  964. int usage2(char *str)
  965. {
  966. if(!strcmp(str,"bounce")) {
  967. fprintf(stderr,
  968. "USAGE:\n"
  969. "bounce bounce inf outf count startgap shorten endlevel ewarp [-smin] [-e] [-c]\n"
  970. "\n"
  971. "\"Bounce\" a sound (Accelerating repeats, decaying in level).\n"
  972. "\n"
  973. "COUNT Number of bounces.\n"
  974. "STARTGAP Gap between src start and 1st bounce. (0.04 to 10)\n"
  975. "SHORTEN Bounce gap reduction (multiplier) from 1 bounce to next.\n"
  976. "ENDLEVEL Loudness of last bounce as fraction of src level. (0 to 1)\n"
  977. "EWARP Decay warp. Decresc at start : >1.0 greater : < 1.0 less.\n"
  978. "-sMIN Shrink bounced elements in same proportion as accel. Min = min dur.\n"
  979. " (Value of zero turns off shrinkage).\n"
  980. "-c If repeating elements overlap, cut to avoid this overlap.\n"
  981. "-e Shrink elements by trimming start. (Default, by trimming end).\n"
  982. "\n");
  983. } else
  984. fprintf(stdout,"Unknown option '%s'\n",str);
  985. return(USAGE_ONLY);
  986. }
  987. int usage3(char *str1,char *str2)
  988. {
  989. fprintf(stderr,"Insufficient parameters on command line.\n");
  990. return(USAGE_ONLY);
  991. }
  992. /******************************** GATE ********************************/
  993. int bounce(dataptr dz)
  994. {
  995. int exit_status, chans = dz->infile->channels, passno, bouncecnt, numberlessone, endsampstepset, levcnt;
  996. double levchange, levpos, thislevel, normaliser, maxsamp, splcincr, normdrop = 0.0, lastratio, thisratio, srate = (double)dz->infile->srate;
  997. int gpsampdur, gpsplen, absplen, opos, start_of_buf, maxpos = 0, startsplice, endsplice, sampstep = 0, endsampstep = 0, n, m, k, minsampdur;
  998. int fadelen = 0, fadecnt = 0, enddur, absbouncepos;
  999. float *ibuf = dz->sampbuf[0], *obuf = dz->sampbuf[1];
  1000. int *steps;
  1001. double *levels;
  1002. int minstep;
  1003. if((steps = (int *)malloc((dz->iparam[BNC_NUMBER] + 2) * sizeof(int)))==NULL) { // Step from orig to qist bounce, for BNC_NUMBER bounces
  1004. sprintf(errstr,"INSUFFICIENT MEMORY TO STORE LENGTHS OF BOUNCE-STEPS.\n"); // + a pseudo-step (= end of last bounce if it is to be cut)
  1005. return(MEMORY_ERROR);
  1006. }
  1007. if((levels = (double *)malloc((dz->iparam[BNC_NUMBER] + 1) * sizeof(double)))==NULL) {
  1008. sprintf(errstr,"INSUFFICIENT MEMORY TO STORE LENGTHS OF BOUNCE-STEPS.\n");
  1009. return(MEMORY_ERROR);
  1010. }
  1011. minsampdur = (int)round(dz->param[BNC_MINDUR] * dz->infile->srate) * chans;
  1012. numberlessone = dz->iparam[BNC_NUMBER] - 1;
  1013. gpsplen = (int)round(BOUNCESPLICE * MS_TO_SECS * (double)dz->infile->srate);
  1014. minstep = (gpsplen/2) * chans;
  1015. absplen = gpsplen * chans;
  1016. gpsampdur = dz->insams[0]/chans;
  1017. normaliser = 1.0;
  1018. bouncecnt = 0;
  1019. levchange = 1.0 - dz->param[BNC_ENDLEV];
  1020. levels[0] = 1.0;
  1021. steps[0] = (int)round(dz->param[BNC_STTSTEP] * srate) * chans;
  1022. if(sloom) // Pre-calc total output len (for progress-bar display)
  1023. dz->tempsize = steps[0];
  1024. lastratio = dz->param[BNC_STTSTEP];
  1025. bouncecnt++;
  1026. levcnt = dz->iparam[BNC_NUMBER] - 1;
  1027. while(bouncecnt <= dz->iparam[BNC_NUMBER]) {
  1028. levpos = (double)levcnt/dz->iparam[BNC_NUMBER]; // levpos = Where we are in change of level, linearly interpd e.g. 2/3 1/3 0 for 3 bounces
  1029. levpos = pow(levpos,dz->param[BNC_LEVWRP]); // Warped at (for pow-decay)
  1030. levels[bouncecnt] = levchange * levpos; // i.e. endlevel+2/3change endlevel+1/3change endlevel
  1031. levels[bouncecnt] += dz->param[BNC_ENDLEV];
  1032. thisratio = lastratio * dz->param[BNC_SHORTEN];
  1033. steps[bouncecnt] = (int)round(thisratio * srate) * chans; // Distance to next bounce
  1034. if(steps[bouncecnt] < minstep) {
  1035. dz->iparam[BNC_NUMBER] = bouncecnt;
  1036. fprintf(stdout,"INFO: Bounce Steps get too short for splices at bounce %d : curtailing.\n",bouncecnt);
  1037. fflush(stdout);
  1038. break;
  1039. }
  1040. lastratio = thisratio;
  1041. if(sloom)
  1042. dz->tempsize += steps[bouncecnt];
  1043. bouncecnt++;
  1044. levcnt--;
  1045. }
  1046. thisratio = lastratio * dz->param[BNC_SHORTEN];
  1047. steps[bouncecnt] = (int)round(gpsampdur * thisratio) * chans; // Even last event has a distance-to-next-bounce
  1048. if(sloom && !dz->vflag[CUTBNC]) { // If src not cut
  1049. dz->tempsize -= steps[bouncecnt-1]; // Go back to start of last bounced-element
  1050. dz->tempsize += dz->insams[0]; // and add duration of src (duration of element)
  1051. }
  1052. for(passno = 0;passno < 2; passno++) {
  1053. start_of_buf = 0;
  1054. absbouncepos = 0;
  1055. if(sloom)
  1056. display_virtual_time(0,dz);
  1057. memset((char *)obuf,0,dz->buflen * 2 * sizeof(float)); // Clear outbuf & ovflwbuf
  1058. sndseekEx(dz->ifd[0],0,0);
  1059. dz->total_samps_read = 0;
  1060. dz->samps_left = dz->insams[0];
  1061. if((exit_status = read_samps(ibuf,dz))<0)
  1062. return(exit_status);
  1063. dz->total_samps_written = 0;
  1064. maxsamp = 0.0;
  1065. opos = 0; // Write src, before bouncing
  1066. if(dz->vflag[CUTBNC] && dz->duration > dz->param[BNC_STTSTEP]) {
  1067. endsplice = steps[0];
  1068. startsplice = endsplice - absplen;
  1069. memcpy((char *)obuf,(char *)ibuf,startsplice * sizeof(float));
  1070. n = startsplice;
  1071. k = gpsplen - 1;
  1072. while(k >= 0) {
  1073. splcincr = (double)k/(double)gpsplen;
  1074. for(m = 0; m < chans; m++) {
  1075. obuf[n] = (float)(ibuf[n] * splcincr);
  1076. n++;
  1077. }
  1078. k--;
  1079. }
  1080. } else
  1081. memcpy((char *)obuf,(char *)ibuf,dz->insams[0] * sizeof(float));
  1082. bouncecnt = 1;
  1083. while(bouncecnt <= dz->iparam[BNC_NUMBER]) {
  1084. thislevel = levels[bouncecnt];
  1085. sampstep = steps[bouncecnt-1];
  1086. absbouncepos += sampstep;
  1087. opos = absbouncepos - start_of_buf;
  1088. while(opos >= dz->buflen) { // If outbuffer filled
  1089. if ((exit_status = handle_buffer_output(passno,&maxsamp,&maxpos,fadelen,&fadecnt,&start_of_buf,steps[0],normdrop,dz->buflen,&opos,dz)) < 0)
  1090. return exit_status;
  1091. }
  1092. if(dz->vflag[CUTBNC] || dz->param[BNC_MINDUR] > 0.0) { // If sampstep to next bounce is less than src len, need cut bounced copy to fit
  1093. endsampstepset = 0;
  1094. if(dz->vflag[CUTBNC]) {
  1095. endsampstep = steps[bouncecnt]; // Actual sampstep in samples = max length of current bounced-item
  1096. endsampstepset = 1;
  1097. }
  1098. if(dz->param[BNC_MINDUR] > 0.0) { // If bounced items are shrunk, shirnk proprtionately to shrinkage in bounce-steps
  1099. enddur = (int)round(((double)steps[bouncecnt]/(double)steps[0]) * (double)gpsampdur) * chans;
  1100. enddur = max(enddur,minsampdur); // but not to a duration less than the MIn allowed
  1101. if(endsampstepset)
  1102. endsampstep = min(enddur,endsampstep); // length of bounced item is shorter of 2 possible curtails
  1103. else
  1104. endsampstep = enddur;
  1105. }
  1106. if(endsampstep < dz->insams[0]) { // Overlap
  1107. if(dz->vflag[KEEPEND]) {
  1108. startsplice = dz->insams[0] - endsampstep;
  1109. endsplice = startsplice + absplen; // Splice off start
  1110. k = 0;
  1111. n = startsplice;
  1112. /* TEST */
  1113. if(dz->total_samps_written + opos > 132496) {
  1114. fprintf(stdout,"AAA at %lf\n",(dz->total_samps_written + opos)/(double)dz->infile->srate);
  1115. fflush(stdout);
  1116. }
  1117. /* TEST */
  1118. while(k < gpsplen) {
  1119. splcincr = (double)k/(double)gpsplen;
  1120. for(m = 0; m < chans; m++) {
  1121. obuf[opos] = (float)(obuf[opos] + (ibuf[n++] * thislevel * splcincr));
  1122. opos++;
  1123. }
  1124. k++;
  1125. /* TEST */
  1126. if(dz->total_samps_written + opos > 132496) {
  1127. fprintf(stdout,"BBB at %lf\n",(dz->total_samps_written + opos)/(double)dz->infile->srate);
  1128. fflush(stdout);
  1129. }
  1130. /* TEST */
  1131. }
  1132. while(n < dz->insams[0])
  1133. /* TEST */
  1134. if(dz->total_samps_written + opos > 132496) {
  1135. fprintf(stdout,"CCC at %lf\n",(dz->total_samps_written + opos)/(double)dz->infile->srate);
  1136. fflush(stdout);
  1137. }
  1138. /* TEST */
  1139. obuf[opos++] = (float)(ibuf[n++] * thislevel);
  1140. } else {
  1141. /* TEST */
  1142. if(dz->total_samps_written + opos > 132496) {
  1143. fprintf(stdout,"DDD at %lf\n",(dz->total_samps_written + opos)/(double)dz->infile->srate);
  1144. fflush(stdout);
  1145. }
  1146. /* TEST */
  1147. startsplice = endsampstep - absplen; // Splice off end
  1148. endsplice = endsampstep;
  1149. for(n=0;n<startsplice;n++) {
  1150. obuf[opos] = (float)(obuf[opos] + (ibuf[n] * thislevel));
  1151. opos++;
  1152. }
  1153. /* TEST */
  1154. if(dz->total_samps_written + opos > 132496) {
  1155. fprintf(stdout,"EEE bouncecnt = %d at %lf\n",bouncecnt,(dz->total_samps_written + opos)/(double)dz->infile->srate);
  1156. fflush(stdout);
  1157. }
  1158. /* TEST */
  1159. k = gpsplen - 1;
  1160. while(k >= 0) {
  1161. splcincr = (double)k/(double)gpsplen;
  1162. for(m = 0; m < chans; m++) {
  1163. obuf[opos] = (float)(obuf[opos] + (ibuf[n++] * thislevel * splcincr));
  1164. opos++;
  1165. }
  1166. k--;
  1167. }
  1168. /* TEST */
  1169. if(dz->total_samps_written + opos > 132496) {
  1170. fprintf(stdout,"FFF bouncecnt = %d at %lf\n",bouncecnt,(dz->total_samps_written + opos)/(double)dz->infile->srate);
  1171. fflush(stdout);
  1172. }
  1173. /* TEST */
  1174. }
  1175. } else { // No overlap, straight copy
  1176. /* TEST */
  1177. if(dz->total_samps_written + opos > 132496) {
  1178. fprintf(stdout,"GGG\n");
  1179. fflush(stdout);
  1180. }
  1181. /* TEST */
  1182. for(n=0;n<dz->insams[0];n++) {
  1183. obuf[opos] = (float)(obuf[opos] + (ibuf[n] * thislevel));
  1184. opos++;
  1185. }
  1186. }
  1187. } else { // As overlap possible, ADD new copy to output
  1188. /* TEST */
  1189. if(dz->total_samps_written + opos > 132496) {
  1190. fprintf(stdout,"HHH\n");
  1191. fflush(stdout);
  1192. }
  1193. /* TEST */
  1194. for(n=0;n<dz->insams[0];n++) {
  1195. for(n=0;n<dz->insams[0];n++) {
  1196. obuf[opos] = (float)(obuf[opos] + (ibuf[n] * thislevel));
  1197. opos++;
  1198. }
  1199. }
  1200. }
  1201. bouncecnt++;
  1202. }
  1203. if(opos > 0) {
  1204. /* TEST */
  1205. if(dz->total_samps_written + opos > 132496) {
  1206. fprintf(stdout,"III bouncecnt = %d at %lf\n",bouncecnt,(dz->total_samps_written + opos)/(double)dz->infile->srate);
  1207. fflush(stdout);
  1208. }
  1209. /* TEST */
  1210. while ((exit_status = handle_buffer_output(passno,&maxsamp,&maxpos,fadelen,&fadecnt,&start_of_buf,steps[0],normdrop,opos,&opos,dz)) < 0)
  1211. return exit_status;
  1212. }
  1213. if(passno == 0 && maxsamp > 0.95) {
  1214. normaliser = 0.95/maxsamp;
  1215. normdrop = 1.0 - normaliser;
  1216. fadelen = maxpos - steps[0];
  1217. }
  1218. }
  1219. return FINISHED;
  1220. }
  1221. /*************************** CREATE_BOUNCE_SNDBUFS **************************/
  1222. int create_bounce_sndbufs(dataptr dz)
  1223. {
  1224. int n;
  1225. int bigbufsize;
  1226. int framesize, secsize;
  1227. framesize = F_SECSIZE * dz->infile->channels * 2;
  1228. if(dz->sbufptr == 0 || dz->sampbuf==0) {
  1229. sprintf(errstr,"buffer pointers not allocated: create_sndbufs()\n");
  1230. return(PROGRAM_ERROR);
  1231. }
  1232. dz->buflen = dz->insams[0];
  1233. secsize = dz->buflen/framesize;
  1234. secsize++;
  1235. dz->buflen = secsize * framesize;
  1236. bigbufsize = dz->buflen * sizeof(float);
  1237. if(dz->buflen < 0 || bigbufsize < 0) {
  1238. sprintf(errstr,"INPUT FILE too large for this process.\n");
  1239. return(PROGRAM_ERROR);
  1240. }
  1241. if((dz->bigbuf = (float *)malloc(bigbufsize * dz->bufcnt)) == NULL) {
  1242. sprintf(errstr,"INSUFFICIENT MEMORY to create sound buffers.\n");
  1243. return(PROGRAM_ERROR);
  1244. }
  1245. for(n=0;n<dz->bufcnt;n++)
  1246. dz->sbufptr[n] = dz->sampbuf[n] = dz->bigbuf + (dz->buflen * n);
  1247. dz->sampbuf[n] = dz->bigbuf + (dz->buflen * n);
  1248. return(FINISHED);
  1249. }
  1250. /*************************** BOUNCE_PARAM_PREPROCESS **************************/
  1251. int bounce_param_preprocess(dataptr dz)
  1252. {
  1253. if(dz->param[BNC_MINDUR] > 0.0 && dz->param[BNC_MINDUR] < dz->application->default_val[BNC_MINDUR]) {
  1254. sprintf(errstr,"MIN DURATION OF SHRUNK BOUNCED ELEMENT CANNOT BE LESS THAN %.3lf UNLESS SET TO ZERO (SWITCHES OFF SHRINKAGE)\n",BOUNCEMINDUR);
  1255. return DATA_ERROR;
  1256. }
  1257. return FINISHED;
  1258. }
  1259. /*************************** NORMALISATION_FADE **************************/
  1260. int normalisation_fade(int fadelen, int *fadecnt, int start_of_buf, int absfadestt, double normdrop, dataptr dz)
  1261. {
  1262. int fadestt, fadeend, n;
  1263. double thisdrop, thislevl;
  1264. float *obuf = dz->sampbuf[1];
  1265. int fcnt = *fadecnt;
  1266. // if buffer ends after firstbounce-starts (where fade starts)
  1267. if((dz->total_samps_read > absfadestt) && (fcnt < fadelen)) {
  1268. fcnt = *fadecnt; // && not et at end of fade (fcnt < fadelen)
  1269. fadestt = absfadestt - start_of_buf; // Start of fade relative to start of buffer
  1270. fadestt = max(0,fadestt); // but not before start of buffer.
  1271. fadeend = fadelen - fcnt; // Remaining length of fade
  1272. fadeend += fadestt; // End of fade relative to current fade start
  1273. fadeend = min(fadeend,dz->buflen); // but not beyond end of buffer
  1274. for(n=fadestt;n<fadeend;n++) { // increment through the total drop in level required
  1275. thisdrop = normdrop * (double)fcnt/(double)fadelen;
  1276. thislevl = 1.0 - thisdrop; // Actual level is 1.0 minus current level-drop
  1277. obuf[n] = (float)(obuf[n] * thislevl); // set level required
  1278. fcnt++; // and count where we are in fade
  1279. }
  1280. *fadecnt = fcnt;
  1281. }
  1282. return FINISHED;
  1283. }
  1284. /**************************************** HANDLE_BUFFER_OUTPUT ********************************/
  1285. int handle_buffer_output(int passno,double *maxsamp,int *maxpos,int fadelen,int *fadecnt,int *start_of_buf,int bnc1len,double normdrop,int samps_to_write,int *opos,dataptr dz)
  1286. {
  1287. int exit_status;
  1288. double absval;
  1289. int n;
  1290. float *obuf = dz->sampbuf[1], *ovflwbuf = dz->sampbuf[2];
  1291. if(passno == 0) {
  1292. for(n=0;n<dz->buflen;n++) {
  1293. absval = fabs(obuf[n]);
  1294. if(absval > *maxsamp) {
  1295. *maxsamp = absval;
  1296. *maxpos = *start_of_buf + n;
  1297. }
  1298. }
  1299. if(sloom) {
  1300. dz->total_samps_written += dz->buflen;
  1301. dz->process = GREV;
  1302. display_virtual_time(dz->total_samps_written,dz);
  1303. dz->process = BOUNCE;
  1304. }
  1305. } else {
  1306. if(fadelen) // If normalisation required, do normalisation fade
  1307. normalisation_fade(fadelen,fadecnt,*start_of_buf,bnc1len,normdrop,dz);
  1308. dz->process = GREV; // then write to output
  1309. if((exit_status = write_samps(obuf,samps_to_write,dz))<0)
  1310. return(exit_status);
  1311. dz->process = BOUNCE;
  1312. }
  1313. memcpy((char *)obuf,(char *)ovflwbuf,dz->buflen * sizeof(float));
  1314. memset((char *)ovflwbuf,0,dz->buflen * sizeof(float));
  1315. *opos -= dz->buflen;
  1316. *start_of_buf += dz->buflen;
  1317. return FINISHED;
  1318. }