quirk.c 40 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889909192939495969798991001011021031041051061071081091101111121131141151161171181191201211221231241251261271281291301311321331341351361371381391401411421431441451461471481491501511521531541551561571581591601611621631641651661671681691701711721731741751761771781791801811821831841851861871881891901911921931941951961971981992002012022032042052062072082092102112122132142152162172182192202212222232242252262272282292302312322332342352362372382392402412422432442452462472482492502512522532542552562572582592602612622632642652662672682692702712722732742752762772782792802812822832842852862872882892902912922932942952962972982993003013023033043053063073083093103113123133143153163173183193203213223233243253263273283293303313323333343353363373383393403413423433443453463473483493503513523533543553563573583593603613623633643653663673683693703713723733743753763773783793803813823833843853863873883893903913923933943953963973983994004014024034044054064074084094104114124134144154164174184194204214224234244254264274284294304314324334344354364374384394404414424434444454464474484494504514524534544554564574584594604614624634644654664674684694704714724734744754764774784794804814824834844854864874884894904914924934944954964974984995005015025035045055065075085095105115125135145155165175185195205215225235245255265275285295305315325335345355365375385395405415425435445455465475485495505515525535545555565575585595605615625635645655665675685695705715725735745755765775785795805815825835845855865875885895905915925935945955965975985996006016026036046056066076086096106116126136146156166176186196206216226236246256266276286296306316326336346356366376386396406416426436446456466476486496506516526536546556566576586596606616626636646656666676686696706716726736746756766776786796806816826836846856866876886896906916926936946956966976986997007017027037047057067077087097107117127137147157167177187197207217227237247257267277287297307317327337347357367377387397407417427437447457467477487497507517527537547557567577587597607617627637647657667677687697707717727737747757767777787797807817827837847857867877887897907917927937947957967977987998008018028038048058068078088098108118128138148158168178188198208218228238248258268278288298308318328338348358368378388398408418428438448458468478488498508518528538548558568578588598608618628638648658668678688698708718728738748758768778788798808818828838848858868878888898908918928938948958968978988999009019029039049059069079089099109119129139149159169179189199209219229239249259269279289299309319329339349359369379389399409419429439449459469479489499509519529539549559569579589599609619629639649659669679689699709719729739749759769779789799809819829839849859869879889899909919929939949959969979989991000100110021003100410051006100710081009101010111012101310141015101610171018101910201021102210231024102510261027102810291030103110321033103410351036103710381039104010411042104310441045104610471048104910501051105210531054105510561057105810591060106110621063106410651066106710681069107010711072107310741075107610771078107910801081108210831084108510861087108810891090109110921093109410951096109710981099110011011102110311041105110611071108110911101111111211131114111511161117111811191120112111221123112411251126112711281129113011311132113311341135113611371138113911401141114211431144114511461147114811491150115111521153115411551156115711581159116011611162116311641165116611671168116911701171117211731174117511761177117811791180118111821183118411851186118711881189119011911192119311941195119611971198119912001201120212031204120512061207120812091210121112121213121412151216121712181219122012211222122312241225122612271228122912301231123212331234123512361237123812391240124112421243124412451246124712481249125012511252125312541255125612571258125912601261126212631264126512661267126812691270127112721273127412751276127712781279128012811282128312841285128612871288128912901291129212931294129512961297129812991300130113021303130413051306130713081309131013111312131313141315131613171318131913201321132213231324
  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 <ctype.h>
  17. #include <sfsys.h>
  18. #include <string.h>
  19. #include <srates.h>
  20. #ifdef unix
  21. #define round(x) lround((x))
  22. #endif
  23. char errstr[2400];
  24. int anal_infiles = 1;
  25. int sloom = 0;
  26. int sloombatch = 0;
  27. const char* cdp_version = "7.0.0";
  28. //CDP LIB REPLACEMENTS
  29. static int quirk_param_preprocess(dataptr dz);
  30. static int setup_quirk_application(dataptr dz);
  31. static int parse_sloom_data(int argc,char *argv[],char ***cmdline,int *cmdlinecnt,dataptr dz);
  32. static int parse_infile_and_check_type(char **cmdline,dataptr dz);
  33. static int setup_quirk_param_ranges_and_defaults(dataptr dz);
  34. static int handle_the_outfile(int *cmdlinecnt,char ***cmdline,dataptr dz);
  35. static int setup_and_init_input_param_activity(dataptr dz,int tipc);
  36. static int setup_input_param_defaultval_stores(int tipc,aplptr ap);
  37. static int establish_application(dataptr dz);
  38. static int initialise_vflags(dataptr dz);
  39. static int setup_parameter_storage_and_constants(int storage_cnt,dataptr dz);
  40. static int initialise_is_int_and_no_brk_constants(int storage_cnt,dataptr dz);
  41. static int mark_parameter_types(dataptr dz,aplptr ap);
  42. static int assign_file_data_storage(int infilecnt,dataptr dz);
  43. static int get_tk_cmdline_word(int *cmdlinecnt,char ***cmdline,char *q);
  44. static int get_the_process_no(char *prog_identifier_from_cmdline,dataptr dz);
  45. static int get_the_mode_from_cmdline(char *str,dataptr dz);
  46. static int setup_and_init_input_brktable_constants(dataptr dz,int brkcnt);
  47. static int quirk(dataptr dz);
  48. /**************************************** MAIN *********************************************/
  49. int main(int argc,char *argv[])
  50. {
  51. int exit_status;
  52. dataptr dz = NULL;
  53. char **cmdline;
  54. int cmdlinecnt;
  55. int n;
  56. aplptr ap;
  57. int is_launched = FALSE;
  58. if(argc==2 && (strcmp(argv[1],"--version") == 0)) {
  59. fprintf(stdout,"%s\n",cdp_version);
  60. fflush(stdout);
  61. return 0;
  62. }
  63. /* CHECK FOR SOUNDLOOM */
  64. if((sloom = sound_loom_in_use(&argc,&argv)) > 1) {
  65. sloom = 0;
  66. sloombatch = 1;
  67. }
  68. if(sflinit("cdp")){
  69. sfperror("cdp: initialisation\n");
  70. return(FAILED);
  71. }
  72. /* SET UP THE PRINCIPLE DATASTRUCTURE */
  73. if((exit_status = establish_datastructure(&dz))<0) { // CDP LIB
  74. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  75. return(FAILED);
  76. }
  77. if(!sloom) {
  78. if(argc == 1) {
  79. usage1();
  80. return(FAILED);
  81. } else if(argc == 2) {
  82. usage2(argv[1]);
  83. return(FAILED);
  84. }
  85. }
  86. if(!sloom) {
  87. if((exit_status = make_initial_cmdline_check(&argc,&argv))<0) { // CDP LIB
  88. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  89. return(FAILED);
  90. }
  91. cmdline = argv;
  92. cmdlinecnt = argc;
  93. if((get_the_process_no(argv[0],dz))<0)
  94. return(FAILED);
  95. cmdline++;
  96. cmdlinecnt--;
  97. dz->maxmode = 2;
  98. if((exit_status = get_the_mode_from_cmdline(cmdline[0],dz))<0) {
  99. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  100. return(exit_status);
  101. }
  102. cmdline++;
  103. cmdlinecnt--;
  104. // setup_particular_application =
  105. if((exit_status = setup_quirk_application(dz))<0) {
  106. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  107. return(FAILED);
  108. }
  109. if((exit_status = count_and_allocate_for_infiles(cmdlinecnt,cmdline,dz))<0) { // CDP LIB
  110. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  111. return(FAILED);
  112. }
  113. } else {
  114. //parse_TK_data() =
  115. if((exit_status = parse_sloom_data(argc,argv,&cmdline,&cmdlinecnt,dz))<0) {
  116. exit_status = print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  117. return(exit_status);
  118. }
  119. }
  120. ap = dz->application;
  121. // parse_infile_and_hone_type() =
  122. if((exit_status = parse_infile_and_check_type(cmdline,dz))<0) {
  123. exit_status = print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  124. return(FAILED);
  125. }
  126. // setup_param_ranges_and_defaults() =
  127. if((exit_status = setup_quirk_param_ranges_and_defaults(dz))<0) {
  128. exit_status = print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  129. return(FAILED);
  130. }
  131. // open_first_infile CDP LIB
  132. if((exit_status = open_first_infile(cmdline[0],dz))<0) {
  133. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  134. return(FAILED);
  135. }
  136. cmdlinecnt--;
  137. cmdline++;
  138. // handle_extra_infiles() : redundant
  139. // handle_outfile() =
  140. if((exit_status = handle_the_outfile(&cmdlinecnt,&cmdline,dz))<0) {
  141. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  142. return(FAILED);
  143. }
  144. // handle_formants() redundant
  145. // handle_formant_quiksearch() redundant
  146. // handle_special_data() redundant
  147. if((exit_status = read_parameters_and_flags(&cmdline,&cmdlinecnt,dz))<0) { // CDP LIB
  148. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  149. return(FAILED);
  150. }
  151. is_launched = TRUE;
  152. dz->bufcnt = 1;
  153. if((dz->sampbuf = (float **)malloc(sizeof(float *) * (dz->bufcnt+1)))==NULL) {
  154. sprintf(errstr,"INSUFFICIENT MEMORY establishing sample buffers.\n");
  155. return(MEMORY_ERROR);
  156. }
  157. if((dz->sbufptr = (float **)malloc(sizeof(float *) * dz->bufcnt))==NULL) {
  158. sprintf(errstr,"INSUFFICIENT MEMORY establishing sample buffer pointers.\n");
  159. return(MEMORY_ERROR);
  160. }
  161. for(n = 0;n <dz->bufcnt; n++)
  162. dz->sampbuf[n] = dz->sbufptr[n] = (float *)0;
  163. dz->sampbuf[n] = (float *)0;
  164. if((exit_status = create_sndbufs(dz))<0) { // CDP LIB
  165. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  166. return(FAILED);
  167. }
  168. // check_param_validity_and_consistency....
  169. if((exit_status = quirk_param_preprocess(dz))<0) {
  170. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  171. return(FAILED);
  172. }
  173. //param_preprocess() redundant
  174. //spec_process_file =
  175. if((exit_status = quirk(dz))<0) {
  176. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  177. return(FAILED);
  178. }
  179. if((exit_status = complete_output(dz))<0) { // CDP LIB
  180. print_messages_and_close_sndfiles(exit_status,is_launched,dz);
  181. return(FAILED);
  182. }
  183. exit_status = print_messages_and_close_sndfiles(FINISHED,is_launched,dz); // CDP LIB
  184. free(dz);
  185. return(SUCCEEDED);
  186. }
  187. /**********************************************
  188. REPLACED CDP LIB FUNCTIONS
  189. **********************************************/
  190. /****************************** SET_PARAM_DATA *********************************/
  191. int set_param_data(aplptr ap, int special_data,int maxparamcnt,int paramcnt,char *paramlist)
  192. {
  193. ap->special_data = (char)special_data;
  194. ap->param_cnt = (char)paramcnt;
  195. ap->max_param_cnt = (char)maxparamcnt;
  196. if(ap->max_param_cnt>0) {
  197. if((ap->param_list = (char *)malloc((size_t)(ap->max_param_cnt+1)))==NULL) {
  198. sprintf(errstr,"INSUFFICIENT MEMORY: for param_list\n");
  199. return(MEMORY_ERROR);
  200. }
  201. strcpy(ap->param_list,paramlist);
  202. }
  203. return(FINISHED);
  204. }
  205. /****************************** SET_VFLGS *********************************/
  206. int set_vflgs
  207. (aplptr ap,char *optflags,int optcnt,char *optlist,char *varflags,int vflagcnt, int vparamcnt,char *varlist)
  208. {
  209. ap->option_cnt = (char) optcnt; /*RWD added cast */
  210. if(optcnt) {
  211. if((ap->option_list = (char *)malloc((size_t)(optcnt+1)))==NULL) {
  212. sprintf(errstr,"INSUFFICIENT MEMORY: for option_list\n");
  213. return(MEMORY_ERROR);
  214. }
  215. strcpy(ap->option_list,optlist);
  216. if((ap->option_flags = (char *)malloc((size_t)(optcnt+1)))==NULL) {
  217. sprintf(errstr,"INSUFFICIENT MEMORY: for option_flags\n");
  218. return(MEMORY_ERROR);
  219. }
  220. strcpy(ap->option_flags,optflags);
  221. }
  222. ap->vflag_cnt = (char) vflagcnt;
  223. ap->variant_param_cnt = (char) vparamcnt;
  224. if(vflagcnt) {
  225. if((ap->variant_list = (char *)malloc((size_t)(vflagcnt+1)))==NULL) {
  226. sprintf(errstr,"INSUFFICIENT MEMORY: for variant_list\n");
  227. return(MEMORY_ERROR);
  228. }
  229. strcpy(ap->variant_list,varlist);
  230. if((ap->variant_flags = (char *)malloc((size_t)(vflagcnt+1)))==NULL) {
  231. sprintf(errstr,"INSUFFICIENT MEMORY: for variant_flags\n");
  232. return(MEMORY_ERROR);
  233. }
  234. strcpy(ap->variant_flags,varflags);
  235. }
  236. return(FINISHED);
  237. }
  238. /***************************** APPLICATION_INIT **************************/
  239. int application_init(dataptr dz)
  240. {
  241. int exit_status;
  242. int storage_cnt;
  243. int tipc, brkcnt;
  244. aplptr ap = dz->application;
  245. if(ap->vflag_cnt>0)
  246. initialise_vflags(dz);
  247. tipc = ap->max_param_cnt + ap->option_cnt + ap->variant_param_cnt;
  248. ap->total_input_param_cnt = (char)tipc;
  249. if(tipc>0) {
  250. if((exit_status = setup_input_param_range_stores(tipc,ap))<0)
  251. return(exit_status);
  252. if((exit_status = setup_input_param_defaultval_stores(tipc,ap))<0)
  253. return(exit_status);
  254. if((exit_status = setup_and_init_input_param_activity(dz,tipc))<0)
  255. return(exit_status);
  256. }
  257. brkcnt = tipc;
  258. //THERE ARE NO INPUTFILE brktables USED IN THIS PROCESS
  259. if(brkcnt>0) {
  260. if((exit_status = setup_and_init_input_brktable_constants(dz,brkcnt))<0)
  261. return(exit_status);
  262. }
  263. if((storage_cnt = tipc + ap->internal_param_cnt)>0) {
  264. if((exit_status = setup_parameter_storage_and_constants(storage_cnt,dz))<0)
  265. return(exit_status);
  266. if((exit_status = initialise_is_int_and_no_brk_constants(storage_cnt,dz))<0)
  267. return(exit_status);
  268. }
  269. if((exit_status = mark_parameter_types(dz,ap))<0)
  270. return(exit_status);
  271. // establish_infile_constants() replaced by
  272. dz->infilecnt = 1;
  273. //establish_bufptrs_and_extra_buffers():
  274. return(FINISHED);
  275. }
  276. /********************** SETUP_PARAMETER_STORAGE_AND_CONSTANTS ********************/
  277. /* RWD mallo changed to calloc; helps debug verison run as release! */
  278. int setup_parameter_storage_and_constants(int storage_cnt,dataptr dz)
  279. {
  280. if((dz->param = (double *)calloc(storage_cnt, sizeof(double)))==NULL) {
  281. sprintf(errstr,"setup_parameter_storage_and_constants(): 1\n");
  282. return(MEMORY_ERROR);
  283. }
  284. if((dz->iparam = (int *)calloc(storage_cnt, sizeof(int) ))==NULL) {
  285. sprintf(errstr,"setup_parameter_storage_and_constants(): 2\n");
  286. return(MEMORY_ERROR);
  287. }
  288. if((dz->is_int = (char *)calloc(storage_cnt, sizeof(char)))==NULL) {
  289. sprintf(errstr,"setup_parameter_storage_and_constants(): 3\n");
  290. return(MEMORY_ERROR);
  291. }
  292. if((dz->no_brk = (char *)calloc(storage_cnt, sizeof(char)))==NULL) {
  293. sprintf(errstr,"setup_parameter_storage_and_constants(): 5\n");
  294. return(MEMORY_ERROR);
  295. }
  296. return(FINISHED);
  297. }
  298. /************** INITIALISE_IS_INT_AND_NO_BRK_CONSTANTS *****************/
  299. int initialise_is_int_and_no_brk_constants(int storage_cnt,dataptr dz)
  300. {
  301. int n;
  302. for(n=0;n<storage_cnt;n++) {
  303. dz->is_int[n] = (char)0;
  304. dz->no_brk[n] = (char)0;
  305. }
  306. return(FINISHED);
  307. }
  308. /***************************** MARK_PARAMETER_TYPES **************************/
  309. int mark_parameter_types(dataptr dz,aplptr ap)
  310. {
  311. int n, m; /* PARAMS */
  312. for(n=0;n<ap->max_param_cnt;n++) {
  313. switch(ap->param_list[n]) {
  314. case('0'): break; /* dz->is_active[n] = 0 is default */
  315. case('i'): dz->is_active[n] = (char)1; dz->is_int[n] = (char)1;dz->no_brk[n] = (char)1; break;
  316. case('I'): dz->is_active[n] = (char)1; dz->is_int[n] = (char)1; break;
  317. case('d'): dz->is_active[n] = (char)1; dz->no_brk[n] = (char)1; break;
  318. case('D'): dz->is_active[n] = (char)1; /* normal case: double val or brkpnt file */ break;
  319. default:
  320. sprintf(errstr,"Programming error: invalid parameter type in mark_parameter_types()\n");
  321. return(PROGRAM_ERROR);
  322. }
  323. } /* OPTIONS */
  324. for(n=0,m=ap->max_param_cnt;n<ap->option_cnt;n++,m++) {
  325. switch(ap->option_list[n]) {
  326. case('i'): dz->is_active[m] = (char)1; dz->is_int[m] = (char)1; dz->no_brk[m] = (char)1; break;
  327. case('I'): dz->is_active[m] = (char)1; dz->is_int[m] = (char)1; break;
  328. case('d'): dz->is_active[m] = (char)1; dz->no_brk[m] = (char)1; break;
  329. case('D'): dz->is_active[m] = (char)1; /* normal case: double val or brkpnt file */ break;
  330. default:
  331. sprintf(errstr,"Programming error: invalid option type in mark_parameter_types()\n");
  332. return(PROGRAM_ERROR);
  333. }
  334. } /* VARIANTS */
  335. for(n=0,m=ap->max_param_cnt + ap->option_cnt;n < ap->variant_param_cnt; n++, m++) {
  336. switch(ap->variant_list[n]) {
  337. case('0'): break;
  338. case('i'): dz->is_active[m] = (char)1; dz->is_int[m] = (char)1; dz->no_brk[m] = (char)1; break;
  339. case('I'): dz->is_active[m] = (char)1; dz->is_int[m] = (char)1; break;
  340. case('d'): dz->is_active[m] = (char)1; dz->no_brk[m] = (char)1; break;
  341. case('D'): dz->is_active[m] = (char)1; /* normal case: double val or brkpnt file */ break;
  342. default:
  343. sprintf(errstr,"Programming error: invalid variant type in mark_parameter_types()\n");
  344. return(PROGRAM_ERROR);
  345. }
  346. } /* INTERNAL */
  347. for(n=0,
  348. m=ap->max_param_cnt + ap->option_cnt + ap->variant_param_cnt; n<ap->internal_param_cnt; n++,m++) {
  349. switch(ap->internal_param_list[n]) {
  350. case('0'): break; /* dummy variables: variables not used: but important for internal paream numbering!! */
  351. case('i'): dz->is_int[m] = (char)1; dz->no_brk[m] = (char)1; break;
  352. case('d'): dz->no_brk[m] = (char)1; break;
  353. default:
  354. sprintf(errstr,"Programming error: invalid internal param type in mark_parameter_types()\n");
  355. return(PROGRAM_ERROR);
  356. }
  357. }
  358. return(FINISHED);
  359. }
  360. /************************ HANDLE_THE_OUTFILE *********************/
  361. int handle_the_outfile(int *cmdlinecnt,char ***cmdline,dataptr dz)
  362. {
  363. int exit_status;
  364. char *filename = (*cmdline)[0];
  365. if(filename[0]=='-' && filename[1]=='f') {
  366. dz->floatsam_output = 1;
  367. dz->true_outfile_stype = SAMP_FLOAT;
  368. filename+= 2;
  369. }
  370. if(!sloom) {
  371. if(file_has_invalid_startchar(filename) || value_is_numeric(filename)) {
  372. sprintf(errstr,"Outfile name %s has invalid start character(s) or looks too much like a number.\n",filename);
  373. return(DATA_ERROR);
  374. }
  375. }
  376. strcpy(dz->outfilename,filename);
  377. if((exit_status = create_sized_outfile(filename,dz))<0)
  378. return(exit_status);
  379. (*cmdline)++;
  380. (*cmdlinecnt)--;
  381. return(FINISHED);
  382. }
  383. /***************************** ESTABLISH_APPLICATION **************************/
  384. int establish_application(dataptr dz)
  385. {
  386. aplptr ap;
  387. if((dz->application = (aplptr)malloc(sizeof (struct applic)))==NULL) {
  388. sprintf(errstr,"establish_application()\n");
  389. return(MEMORY_ERROR);
  390. }
  391. ap = dz->application;
  392. memset((char *)ap,0,sizeof(struct applic));
  393. return(FINISHED);
  394. }
  395. /************************* INITIALISE_VFLAGS *************************/
  396. int initialise_vflags(dataptr dz)
  397. {
  398. int n;
  399. if((dz->vflag = (char *)malloc(dz->application->vflag_cnt * sizeof(char)))==NULL) {
  400. sprintf(errstr,"INSUFFICIENT MEMORY: vflag store,\n");
  401. return(MEMORY_ERROR);
  402. }
  403. for(n=0;n<dz->application->vflag_cnt;n++)
  404. dz->vflag[n] = FALSE;
  405. return FINISHED;
  406. }
  407. /************************* SETUP_INPUT_PARAM_DEFAULTVALS *************************/
  408. int setup_input_param_defaultval_stores(int tipc,aplptr ap)
  409. {
  410. int n;
  411. if((ap->default_val = (double *)malloc(tipc * sizeof(double)))==NULL) {
  412. sprintf(errstr,"INSUFFICIENT MEMORY for application default values store\n");
  413. return(MEMORY_ERROR);
  414. }
  415. for(n=0;n<tipc;n++)
  416. ap->default_val[n] = 0.0;
  417. return(FINISHED);
  418. }
  419. /***************************** SETUP_AND_INIT_INPUT_PARAM_ACTIVITY **************************/
  420. int setup_and_init_input_param_activity(dataptr dz,int tipc)
  421. {
  422. int n;
  423. if((dz->is_active = (char *)malloc((size_t)tipc))==NULL) {
  424. sprintf(errstr,"setup_and_init_input_param_activity()\n");
  425. return(MEMORY_ERROR);
  426. }
  427. for(n=0;n<tipc;n++)
  428. dz->is_active[n] = (char)0;
  429. return(FINISHED);
  430. }
  431. /************************* SETUP_QUIRK_APPLICATION *******************/
  432. int setup_quirk_application(dataptr dz)
  433. {
  434. int exit_status;
  435. aplptr ap;
  436. if((exit_status = establish_application(dz))<0) // GLOBAL
  437. return(FAILED);
  438. ap = dz->application;
  439. // SEE parstruct FOR EXPLANATION of next 2 functions
  440. if((exit_status = set_param_data(ap,0 ,1,1,"d"))<0)
  441. return(FAILED);
  442. if((exit_status = set_vflgs(ap,"",0,"","",0,0,""))<0)
  443. return(FAILED);
  444. // set_legal_infile_structure -->
  445. dz->has_otherfile = FALSE;
  446. // assign_process_logic -->
  447. dz->input_data_type = SNDFILES_ONLY;
  448. dz->process_type = EQUAL_SNDFILE;
  449. dz->outfiletype = SNDFILE_OUT;
  450. return application_init(dz); //GLOBAL
  451. }
  452. /************************* PARSE_INFILE_AND_CHECK_TYPE *******************/
  453. int parse_infile_and_check_type(char **cmdline,dataptr dz)
  454. {
  455. int exit_status;
  456. infileptr infile_info;
  457. if(!sloom) {
  458. if((infile_info = (infileptr)malloc(sizeof(struct filedata)))==NULL) {
  459. sprintf(errstr,"INSUFFICIENT MEMORY for infile structure to test file data.");
  460. return(MEMORY_ERROR);
  461. } else if((exit_status = cdparse(cmdline[0],infile_info))<0) {
  462. sprintf(errstr,"Failed to parse input file %s\n",cmdline[0]);
  463. return(PROGRAM_ERROR);
  464. } else if(infile_info->filetype != SNDFILE) {
  465. sprintf(errstr,"File %s is not of correct type\n",cmdline[0]);
  466. return(DATA_ERROR);
  467. } else if(infile_info->channels != 1) {
  468. sprintf(errstr,"File %s is not of correct type (must be mono)\n",cmdline[0]);
  469. return(DATA_ERROR);
  470. } else if((exit_status = copy_parse_info_to_main_structure(infile_info,dz))<0) {
  471. sprintf(errstr,"Failed to copy file parsing information\n");
  472. return(PROGRAM_ERROR);
  473. }
  474. free(infile_info);
  475. }
  476. return(FINISHED);
  477. }
  478. /************************* SETUP_QUIRK_PARAM_RANGES_AND_DEFAULTS *******************/
  479. int setup_quirk_param_ranges_and_defaults(dataptr dz)
  480. {
  481. int exit_status;
  482. aplptr ap = dz->application;
  483. // set_param_ranges()
  484. ap->total_input_param_cnt = (char)(ap->max_param_cnt + ap->option_cnt + ap->variant_param_cnt);
  485. // NB total_input_param_cnt is > 0 !!!
  486. if((exit_status = setup_input_param_range_stores(ap->total_input_param_cnt,ap))<0)
  487. return(FAILED);
  488. // get_param_ranges()
  489. ap->lo[0] = 0.01;
  490. ap->hi[0] = 100;
  491. ap->default_val[0] = 1;
  492. dz->maxmode = 2;
  493. if(!sloom)
  494. put_default_vals_in_all_params(dz);
  495. return(FINISHED);
  496. }
  497. /********************************* PARSE_SLOOM_DATA *********************************/
  498. int parse_sloom_data(int argc,char *argv[],char ***cmdline,int *cmdlinecnt,dataptr dz)
  499. {
  500. int exit_status;
  501. int cnt = 1, infilecnt;
  502. int filesize, insams, inbrksize;
  503. double dummy;
  504. int true_cnt = 0;
  505. aplptr ap;
  506. while(cnt<=PRE_CMDLINE_DATACNT) {
  507. if(cnt > argc) {
  508. sprintf(errstr,"Insufficient data sent from TK\n");
  509. return(DATA_ERROR);
  510. }
  511. switch(cnt) {
  512. case(1):
  513. if(sscanf(argv[cnt],"%d",&dz->process)!=1) {
  514. sprintf(errstr,"Cannot read process no. sent from TK\n");
  515. return(DATA_ERROR);
  516. }
  517. break;
  518. case(2):
  519. if(sscanf(argv[cnt],"%d",&dz->mode)!=1) {
  520. sprintf(errstr,"Cannot read mode no. sent from TK\n");
  521. return(DATA_ERROR);
  522. }
  523. if(dz->mode > 0)
  524. dz->mode--;
  525. //setup_particular_application() =
  526. if((exit_status = setup_quirk_application(dz))<0)
  527. return(exit_status);
  528. ap = dz->application;
  529. break;
  530. case(3):
  531. if(sscanf(argv[cnt],"%d",&infilecnt)!=1) {
  532. sprintf(errstr,"Cannot read infilecnt sent from TK\n");
  533. return(DATA_ERROR);
  534. }
  535. if(infilecnt < 1) {
  536. true_cnt = cnt + 1;
  537. cnt = PRE_CMDLINE_DATACNT; /* force exit from loop after assign_file_data_storage */
  538. }
  539. if((exit_status = assign_file_data_storage(infilecnt,dz))<0)
  540. return(exit_status);
  541. break;
  542. case(INPUT_FILETYPE+4):
  543. if(sscanf(argv[cnt],"%d",&dz->infile->filetype)!=1) {
  544. sprintf(errstr,"Cannot read filetype sent from TK (%s)\n",argv[cnt]);
  545. return(DATA_ERROR);
  546. }
  547. break;
  548. case(INPUT_FILESIZE+4):
  549. if(sscanf(argv[cnt],"%d",&filesize)!=1) {
  550. sprintf(errstr,"Cannot read infilesize sent from TK\n");
  551. return(DATA_ERROR);
  552. }
  553. dz->insams[0] = filesize;
  554. break;
  555. case(INPUT_INSAMS+4):
  556. if(sscanf(argv[cnt],"%d",&insams)!=1) {
  557. sprintf(errstr,"Cannot read insams sent from TK\n");
  558. return(DATA_ERROR);
  559. }
  560. dz->insams[0] = insams;
  561. break;
  562. case(INPUT_SRATE+4):
  563. if(sscanf(argv[cnt],"%d",&dz->infile->srate)!=1) {
  564. sprintf(errstr,"Cannot read srate sent from TK\n");
  565. return(DATA_ERROR);
  566. }
  567. break;
  568. case(INPUT_CHANNELS+4):
  569. if(sscanf(argv[cnt],"%d",&dz->infile->channels)!=1) {
  570. sprintf(errstr,"Cannot read channels sent from TK\n");
  571. return(DATA_ERROR);
  572. }
  573. break;
  574. case(INPUT_STYPE+4):
  575. if(sscanf(argv[cnt],"%d",&dz->infile->stype)!=1) {
  576. sprintf(errstr,"Cannot read stype sent from TK\n");
  577. return(DATA_ERROR);
  578. }
  579. break;
  580. case(INPUT_ORIGSTYPE+4):
  581. if(sscanf(argv[cnt],"%d",&dz->infile->origstype)!=1) {
  582. sprintf(errstr,"Cannot read origstype sent from TK\n");
  583. return(DATA_ERROR);
  584. }
  585. break;
  586. case(INPUT_ORIGRATE+4):
  587. if(sscanf(argv[cnt],"%d",&dz->infile->origrate)!=1) {
  588. sprintf(errstr,"Cannot read origrate sent from TK\n");
  589. return(DATA_ERROR);
  590. }
  591. break;
  592. case(INPUT_MLEN+4):
  593. if(sscanf(argv[cnt],"%d",&dz->infile->Mlen)!=1) {
  594. sprintf(errstr,"Cannot read Mlen sent from TK\n");
  595. return(DATA_ERROR);
  596. }
  597. break;
  598. case(INPUT_DFAC+4):
  599. if(sscanf(argv[cnt],"%d",&dz->infile->Dfac)!=1) {
  600. sprintf(errstr,"Cannot read Dfac sent from TK\n");
  601. return(DATA_ERROR);
  602. }
  603. break;
  604. case(INPUT_ORIGCHANS+4):
  605. if(sscanf(argv[cnt],"%d",&dz->infile->origchans)!=1) {
  606. sprintf(errstr,"Cannot read origchans sent from TK\n");
  607. return(DATA_ERROR);
  608. }
  609. break;
  610. case(INPUT_SPECENVCNT+4):
  611. if(sscanf(argv[cnt],"%d",&dz->infile->specenvcnt)!=1) {
  612. sprintf(errstr,"Cannot read specenvcnt sent from TK\n");
  613. return(DATA_ERROR);
  614. }
  615. dz->specenvcnt = dz->infile->specenvcnt;
  616. break;
  617. case(INPUT_WANTED+4):
  618. if(sscanf(argv[cnt],"%d",&dz->wanted)!=1) {
  619. sprintf(errstr,"Cannot read wanted sent from TK\n");
  620. return(DATA_ERROR);
  621. }
  622. break;
  623. case(INPUT_WLENGTH+4):
  624. if(sscanf(argv[cnt],"%d",&dz->wlength)!=1) {
  625. sprintf(errstr,"Cannot read wlength sent from TK\n");
  626. return(DATA_ERROR);
  627. }
  628. break;
  629. case(INPUT_OUT_CHANS+4):
  630. if(sscanf(argv[cnt],"%d",&dz->out_chans)!=1) {
  631. sprintf(errstr,"Cannot read out_chans sent from TK\n");
  632. return(DATA_ERROR);
  633. }
  634. break;
  635. /* RWD these chanegs to samps - tk will have to deal with that! */
  636. case(INPUT_DESCRIPTOR_BYTES+4):
  637. if(sscanf(argv[cnt],"%d",&dz->descriptor_samps)!=1) {
  638. sprintf(errstr,"Cannot read descriptor_samps sent from TK\n");
  639. return(DATA_ERROR);
  640. }
  641. break;
  642. case(INPUT_IS_TRANSPOS+4):
  643. if(sscanf(argv[cnt],"%d",&dz->is_transpos)!=1) {
  644. sprintf(errstr,"Cannot read is_transpos sent from TK\n");
  645. return(DATA_ERROR);
  646. }
  647. break;
  648. case(INPUT_COULD_BE_TRANSPOS+4):
  649. if(sscanf(argv[cnt],"%d",&dz->could_be_transpos)!=1) {
  650. sprintf(errstr,"Cannot read could_be_transpos sent from TK\n");
  651. return(DATA_ERROR);
  652. }
  653. break;
  654. case(INPUT_COULD_BE_PITCH+4):
  655. if(sscanf(argv[cnt],"%d",&dz->could_be_pitch)!=1) {
  656. sprintf(errstr,"Cannot read could_be_pitch sent from TK\n");
  657. return(DATA_ERROR);
  658. }
  659. break;
  660. case(INPUT_DIFFERENT_SRATES+4):
  661. if(sscanf(argv[cnt],"%d",&dz->different_srates)!=1) {
  662. sprintf(errstr,"Cannot read different_srates sent from TK\n");
  663. return(DATA_ERROR);
  664. }
  665. break;
  666. case(INPUT_DUPLICATE_SNDS+4):
  667. if(sscanf(argv[cnt],"%d",&dz->duplicate_snds)!=1) {
  668. sprintf(errstr,"Cannot read duplicate_snds sent from TK\n");
  669. return(DATA_ERROR);
  670. }
  671. break;
  672. case(INPUT_BRKSIZE+4):
  673. if(sscanf(argv[cnt],"%d",&inbrksize)!=1) {
  674. sprintf(errstr,"Cannot read brksize sent from TK\n");
  675. return(DATA_ERROR);
  676. }
  677. if(inbrksize > 0) {
  678. switch(dz->input_data_type) {
  679. case(WORDLIST_ONLY):
  680. break;
  681. case(PITCH_AND_PITCH):
  682. case(PITCH_AND_TRANSPOS):
  683. case(TRANSPOS_AND_TRANSPOS):
  684. dz->tempsize = inbrksize;
  685. break;
  686. case(BRKFILES_ONLY):
  687. case(UNRANGED_BRKFILE_ONLY):
  688. case(DB_BRKFILES_ONLY):
  689. case(ALL_FILES):
  690. case(ANY_NUMBER_OF_ANY_FILES):
  691. if(dz->extrabrkno < 0) {
  692. sprintf(errstr,"Storage location number for brktable not established by CDP.\n");
  693. return(DATA_ERROR);
  694. }
  695. if(dz->brksize == NULL) {
  696. sprintf(errstr,"CDP has not established storage space for input brktable.\n");
  697. return(PROGRAM_ERROR);
  698. }
  699. dz->brksize[dz->extrabrkno] = inbrksize;
  700. break;
  701. default:
  702. sprintf(errstr,"TK sent brktablesize > 0 for input_data_type [%d] not using brktables.\n",
  703. dz->input_data_type);
  704. return(PROGRAM_ERROR);
  705. }
  706. break;
  707. }
  708. break;
  709. case(INPUT_NUMSIZE+4):
  710. if(sscanf(argv[cnt],"%d",&dz->numsize)!=1) {
  711. sprintf(errstr,"Cannot read numsize sent from TK\n");
  712. return(DATA_ERROR);
  713. }
  714. break;
  715. case(INPUT_LINECNT+4):
  716. if(sscanf(argv[cnt],"%d",&dz->linecnt)!=1) {
  717. sprintf(errstr,"Cannot read linecnt sent from TK\n");
  718. return(DATA_ERROR);
  719. }
  720. break;
  721. case(INPUT_ALL_WORDS+4):
  722. if(sscanf(argv[cnt],"%d",&dz->all_words)!=1) {
  723. sprintf(errstr,"Cannot read all_words sent from TK\n");
  724. return(DATA_ERROR);
  725. }
  726. break;
  727. case(INPUT_ARATE+4):
  728. if(sscanf(argv[cnt],"%f",&dz->infile->arate)!=1) {
  729. sprintf(errstr,"Cannot read arate sent from TK\n");
  730. return(DATA_ERROR);
  731. }
  732. break;
  733. case(INPUT_FRAMETIME+4):
  734. if(sscanf(argv[cnt],"%lf",&dummy)!=1) {
  735. sprintf(errstr,"Cannot read frametime sent from TK\n");
  736. return(DATA_ERROR);
  737. }
  738. dz->frametime = (float)dummy;
  739. break;
  740. case(INPUT_WINDOW_SIZE+4):
  741. if(sscanf(argv[cnt],"%f",&dz->infile->window_size)!=1) {
  742. sprintf(errstr,"Cannot read window_size sent from TK\n");
  743. return(DATA_ERROR);
  744. }
  745. break;
  746. case(INPUT_NYQUIST+4):
  747. if(sscanf(argv[cnt],"%lf",&dz->nyquist)!=1) {
  748. sprintf(errstr,"Cannot read nyquist sent from TK\n");
  749. return(DATA_ERROR);
  750. }
  751. break;
  752. case(INPUT_DURATION+4):
  753. if(sscanf(argv[cnt],"%lf",&dz->duration)!=1) {
  754. sprintf(errstr,"Cannot read duration sent from TK\n");
  755. return(DATA_ERROR);
  756. }
  757. break;
  758. case(INPUT_MINBRK+4):
  759. if(sscanf(argv[cnt],"%lf",&dz->minbrk)!=1) {
  760. sprintf(errstr,"Cannot read minbrk sent from TK\n");
  761. return(DATA_ERROR);
  762. }
  763. break;
  764. case(INPUT_MAXBRK+4):
  765. if(sscanf(argv[cnt],"%lf",&dz->maxbrk)!=1) {
  766. sprintf(errstr,"Cannot read maxbrk sent from TK\n");
  767. return(DATA_ERROR);
  768. }
  769. break;
  770. case(INPUT_MINNUM+4):
  771. if(sscanf(argv[cnt],"%lf",&dz->minnum)!=1) {
  772. sprintf(errstr,"Cannot read minnum sent from TK\n");
  773. return(DATA_ERROR);
  774. }
  775. break;
  776. case(INPUT_MAXNUM+4):
  777. if(sscanf(argv[cnt],"%lf",&dz->maxnum)!=1) {
  778. sprintf(errstr,"Cannot read maxnum sent from TK\n");
  779. return(DATA_ERROR);
  780. }
  781. break;
  782. default:
  783. sprintf(errstr,"case switch item missing: parse_sloom_data()\n");
  784. return(PROGRAM_ERROR);
  785. }
  786. cnt++;
  787. }
  788. if(cnt!=PRE_CMDLINE_DATACNT+1) {
  789. sprintf(errstr,"Insufficient pre-cmdline params sent from TK\n");
  790. return(DATA_ERROR);
  791. }
  792. if(true_cnt)
  793. cnt = true_cnt;
  794. *cmdlinecnt = 0;
  795. while(cnt < argc) {
  796. if((exit_status = get_tk_cmdline_word(cmdlinecnt,cmdline,argv[cnt]))<0)
  797. return(exit_status);
  798. cnt++;
  799. }
  800. return(FINISHED);
  801. }
  802. /********************************* GET_TK_CMDLINE_WORD *********************************/
  803. int get_tk_cmdline_word(int *cmdlinecnt,char ***cmdline,char *q)
  804. {
  805. if(*cmdlinecnt==0) {
  806. if((*cmdline = (char **)malloc(sizeof(char *)))==NULL) {
  807. sprintf(errstr,"INSUFFICIENT MEMORY for TK cmdline array.\n");
  808. return(MEMORY_ERROR);
  809. }
  810. } else {
  811. if((*cmdline = (char **)realloc(*cmdline,((*cmdlinecnt)+1) * sizeof(char *)))==NULL) {
  812. sprintf(errstr,"INSUFFICIENT MEMORY for TK cmdline array.\n");
  813. return(MEMORY_ERROR);
  814. }
  815. }
  816. if(((*cmdline)[*cmdlinecnt] = (char *)malloc((strlen(q) + 1) * sizeof(char)))==NULL) {
  817. sprintf(errstr,"INSUFFICIENT MEMORY for TK cmdline item %d.\n",(*cmdlinecnt)+1);
  818. return(MEMORY_ERROR);
  819. }
  820. strcpy((*cmdline)[*cmdlinecnt],q);
  821. (*cmdlinecnt)++;
  822. return(FINISHED);
  823. }
  824. /****************************** ASSIGN_FILE_DATA_STORAGE *********************************/
  825. int assign_file_data_storage(int infilecnt,dataptr dz)
  826. {
  827. int exit_status;
  828. int no_sndfile_system_files = FALSE;
  829. dz->infilecnt = infilecnt;
  830. if((exit_status = allocate_filespace(dz))<0)
  831. return(exit_status);
  832. if(no_sndfile_system_files)
  833. dz->infilecnt = 0;
  834. return(FINISHED);
  835. }
  836. /************************* redundant functions: to ensure libs compile OK *******************/
  837. int assign_process_logic(dataptr dz)
  838. {
  839. return(FINISHED);
  840. }
  841. void set_legal_infile_structure(dataptr dz)
  842. {}
  843. int set_legal_internalparam_structure(int process,int mode,aplptr ap)
  844. {
  845. return(FINISHED);
  846. }
  847. int setup_internal_arrays_and_array_pointers(dataptr dz)
  848. {
  849. return(FINISHED);
  850. }
  851. int establish_bufptrs_and_extra_buffers(dataptr dz)
  852. {
  853. return(FINISHED);
  854. }
  855. int read_special_data(char *str,dataptr dz)
  856. {
  857. return(FINISHED);
  858. }
  859. int inner_loop
  860. (int *peakscore,int *descnt,int *in_start_portion,int *least,int *pitchcnt,int windows_in_buf,dataptr dz)
  861. {
  862. return(FINISHED);
  863. }
  864. int get_process_no(char *prog_identifier_from_cmdline,dataptr dz)
  865. {
  866. return(FINISHED);
  867. }
  868. /******************************** USAGE1 ********************************/
  869. int usage1(void)
  870. {
  871. usage2("quirk");
  872. return(USAGE_ONLY);
  873. }
  874. /********************************************************************************************/
  875. int get_the_process_no(char *prog_identifier_from_cmdline,dataptr dz)
  876. {
  877. if(!strcmp(prog_identifier_from_cmdline,"quirk")) dz->process = QUIRK;
  878. else {
  879. sprintf(errstr,"Unknown program identification string '%s'\n",prog_identifier_from_cmdline);
  880. return(USAGE_ONLY);
  881. }
  882. return(FINISHED);
  883. }
  884. /******************************** SETUP_AND_INIT_INPUT_BRKTABLE_CONSTANTS ********************************/
  885. int setup_and_init_input_brktable_constants(dataptr dz,int brkcnt)
  886. {
  887. int n;
  888. if((dz->brk = (double **)malloc(brkcnt * sizeof(double *)))==NULL) {
  889. sprintf(errstr,"setup_and_init_input_brktable_constants(): 1\n");
  890. return(MEMORY_ERROR);
  891. }
  892. if((dz->brkptr = (double **)malloc(brkcnt * sizeof(double *)))==NULL) {
  893. sprintf(errstr,"setup_and_init_input_brktable_constants(): 6\n");
  894. return(MEMORY_ERROR);
  895. }
  896. if((dz->brksize = (int *)malloc(brkcnt * sizeof(int)))==NULL) {
  897. sprintf(errstr,"setup_and_init_input_brktable_constants(): 2\n");
  898. return(MEMORY_ERROR);
  899. }
  900. if((dz->firstval = (double *)malloc(brkcnt * sizeof(double)))==NULL) {
  901. sprintf(errstr,"setup_and_init_input_brktable_constants(): 3\n");
  902. return(MEMORY_ERROR);
  903. }
  904. if((dz->lastind = (double *)malloc(brkcnt * sizeof(double)))==NULL) {
  905. sprintf(errstr,"setup_and_init_input_brktable_constants(): 4\n");
  906. return(MEMORY_ERROR);
  907. }
  908. if((dz->lastval = (double *)malloc(brkcnt * sizeof(double)))==NULL) {
  909. sprintf(errstr,"setup_and_init_input_brktable_constants(): 5\n");
  910. return(MEMORY_ERROR);
  911. }
  912. if((dz->brkinit = (int *)malloc(brkcnt * sizeof(int)))==NULL) {
  913. sprintf(errstr,"setup_and_init_input_brktable_constants(): 7\n");
  914. return(MEMORY_ERROR);
  915. }
  916. for(n=0;n<brkcnt;n++) {
  917. dz->brk[n] = NULL;
  918. dz->brkptr[n] = NULL;
  919. dz->brkinit[n] = 0;
  920. dz->brksize[n] = 0;
  921. }
  922. return(FINISHED);
  923. }
  924. /******************************** USAGE2 ********************************/
  925. int usage2(char *str)
  926. {
  927. if(!strcmp(str,"quirk")) {
  928. fprintf(stderr,
  929. "USAGE:\n"
  930. "quirk quirk 1-2 infile outfile powfac\n"
  931. "\n"
  932. "Distort signal by raising sample values to a power.\n"
  933. "\n"
  934. "POWFAC Range 0.01 to 100.\n"
  935. " < 1: exaggerates signal contour.\n"
  936. " > 1: smooths signal contour.\n"
  937. "\n"
  938. "MODE 1 Applies power factor over amplitude range of individual half-wavesets.\n"
  939. "MODE 2 Applies power factor over amplitude range of entire signal.\n"
  940. "\n");
  941. } else
  942. fprintf(stdout,"Unknown option '%s'\n",str);
  943. return(USAGE_ONLY);
  944. }
  945. int usage3(char *str1,char *str2)
  946. {
  947. fprintf(stderr,"Insufficient parameters on command line.\n");
  948. return(USAGE_ONLY);
  949. }
  950. /****************************** QUIRK_PARAM_PREPROCESS *********************************/
  951. int quirk_param_preprocess(dataptr dz)
  952. {
  953. int exit_status;
  954. int n, cnt = 0;
  955. float *ibuf = dz->sampbuf[0], maxsamp;
  956. int halfwaveset_cnt = 0;
  957. int phase, initial_phase = 0;
  958. int *wavesetsampcnt;
  959. float *wavesetmaxsamp;
  960. dz->samps_left = dz->insams[0];
  961. // DETERMINE INITIAL PHASE
  962. while(dz->samps_left > 0) {
  963. if((exit_status = read_samps(ibuf,dz))<0)
  964. return(exit_status);
  965. n = 0;
  966. while(n < dz->ssampsread) {
  967. if(ibuf[n] > 0.0) {
  968. initial_phase = 1;
  969. break;
  970. } else if(ibuf[n] < 0.0) {
  971. initial_phase = -1;
  972. break;
  973. }
  974. n++;
  975. }
  976. if(initial_phase != 0)
  977. break;
  978. }
  979. if(initial_phase == 0) {
  980. fprintf(stderr,"No signal found in soundfile: Cannot proceed\n");
  981. return DATA_ERROR;
  982. }
  983. phase = initial_phase;
  984. if((sndseekEx(dz->ifd[0],0,0) < 0)){
  985. sprintf(errstr,"sndseek failed\n");
  986. return SYSTEM_ERROR;
  987. }
  988. dz->samps_left = dz->insams[0];
  989. dz->total_samps_read = 0;
  990. while(dz->samps_left > 0) {
  991. if((exit_status = read_samps(ibuf,dz))<0)
  992. return(exit_status);
  993. n = 0;
  994. while(n < dz->ssampsread) {
  995. switch(initial_phase) {
  996. case(1):
  997. if(phase == initial_phase) {
  998. if(ibuf[n] < 0.0) { // i.e. phase has changed from 1 (initial_phase) to -1
  999. halfwaveset_cnt++; // Count half-waveset
  1000. phase = -phase; // Switch phase
  1001. }
  1002. } else {
  1003. if(ibuf[n] >= 0.0) { // i.e. phase has changed from -1 to 1
  1004. halfwaveset_cnt++; // Count half-waveset
  1005. phase = -phase; // Switch phase
  1006. }
  1007. }
  1008. break;
  1009. case(-1):
  1010. if(phase == initial_phase) {
  1011. if(ibuf[n] >= 0.0) { // i.e. phase has changed from -1 (initial_phase) to 1
  1012. halfwaveset_cnt++; // Count half-waveset
  1013. phase = -phase; // Switch phase
  1014. }
  1015. } else {
  1016. if(ibuf[n] < 0.0) { // i.e. phase has changed from 1 to -1
  1017. halfwaveset_cnt++; // Count half-waveset
  1018. phase = -phase; // Switch phase
  1019. }
  1020. }
  1021. break;
  1022. }
  1023. n++;
  1024. }
  1025. }
  1026. if(cnt > 0)
  1027. halfwaveset_cnt++; // Include any partial waveset at file end (i.e. signal doesn't go to zero)
  1028. dz->itemcnt = halfwaveset_cnt;
  1029. if((sndseekEx(dz->ifd[0],0,0) < 0)){
  1030. sprintf(errstr,"sndseek failed\n");
  1031. return SYSTEM_ERROR;
  1032. }
  1033. dz->samps_left = dz->insams[0];
  1034. dz->total_samps_read = 0;
  1035. dz->total_samps_written = 0;
  1036. // establish array to store waveset sizes and maxsamples
  1037. if((dz->iparray = (int **)malloc(sizeof(int *)))==NULL) {
  1038. sprintf(errstr,"INSUFFICIENT MEMORY for waveset samplecnts storage.\n");
  1039. return(MEMORY_ERROR);
  1040. }
  1041. if((dz->iparray[0] = (int *)malloc(dz->itemcnt * sizeof(int)))==NULL) {
  1042. sprintf(errstr,"INSUFFICIENT MEMORY for waveset samplecnts storage.\n");
  1043. return(MEMORY_ERROR);
  1044. }
  1045. if((dz->fptr = (float **)malloc(sizeof(float *)))==NULL) {
  1046. sprintf(errstr,"INSUFFICIENT MEMORY for waveset maxamps store.\n");
  1047. return(MEMORY_ERROR);
  1048. }
  1049. if((dz->fptr[0] = (float *)malloc(dz->itemcnt * sizeof(float)))==NULL) {
  1050. sprintf(errstr,"INSUFFICIENT MEMORY for waveset maxamps store.\n");
  1051. return(MEMORY_ERROR);
  1052. }
  1053. wavesetsampcnt = dz->iparray[0];
  1054. wavesetmaxsamp = dz->fptr[0];
  1055. halfwaveset_cnt = 0;
  1056. maxsamp = 0.0;
  1057. phase = initial_phase;
  1058. while(dz->samps_left > 0) {
  1059. if((exit_status = read_samps(ibuf,dz))<0)
  1060. return(exit_status);
  1061. n = 0;
  1062. while(n < dz->ssampsread) {
  1063. switch(initial_phase) {
  1064. case(1):
  1065. if(phase == initial_phase) {
  1066. if(ibuf[n] >= 0.0) {
  1067. if(ibuf[n] > maxsamp)
  1068. maxsamp = ibuf[n];
  1069. cnt++;
  1070. }
  1071. else { // i.e. phase has changed from 1 (initial_phase) to -1
  1072. wavesetsampcnt[halfwaveset_cnt] = cnt;
  1073. wavesetmaxsamp[halfwaveset_cnt] = maxsamp;
  1074. cnt = 1;
  1075. maxsamp = ibuf[n];
  1076. halfwaveset_cnt++;
  1077. phase = -phase; // Switch phase
  1078. }
  1079. } else {
  1080. if(ibuf[n] < 0.0) {
  1081. if(ibuf[n] < maxsamp)
  1082. maxsamp = ibuf[n];
  1083. cnt++;
  1084. } else { // i.e. phase has changed from -1 to 1
  1085. wavesetsampcnt[halfwaveset_cnt] = cnt;
  1086. wavesetmaxsamp[halfwaveset_cnt] = maxsamp;
  1087. cnt = 1;
  1088. maxsamp = ibuf[n];
  1089. halfwaveset_cnt++;
  1090. phase = -phase; // Switch phase
  1091. }
  1092. }
  1093. break;
  1094. case(-1):
  1095. if(phase == initial_phase) {
  1096. if(ibuf[n] < 0.0) {
  1097. if(ibuf[n] < maxsamp)
  1098. maxsamp = ibuf[n];
  1099. cnt++;
  1100. } else { // i.e. phase has changed from -1 (initial_phase) to 1
  1101. wavesetsampcnt[halfwaveset_cnt] = cnt;
  1102. wavesetmaxsamp[halfwaveset_cnt] = maxsamp;
  1103. cnt = 1;
  1104. maxsamp = ibuf[n];
  1105. halfwaveset_cnt++;
  1106. phase = -phase; // Switch phase
  1107. }
  1108. } else {
  1109. if(ibuf[n] >= 0.0) {
  1110. if(ibuf[n] > maxsamp)
  1111. maxsamp = ibuf[n];
  1112. cnt++;
  1113. } else { // i.e. phase has changed from 1 to -1
  1114. wavesetsampcnt[halfwaveset_cnt] = cnt;
  1115. wavesetmaxsamp[halfwaveset_cnt] = maxsamp;
  1116. cnt = 1;
  1117. maxsamp = ibuf[n];
  1118. halfwaveset_cnt++;
  1119. phase = -phase; // Switch phase
  1120. }
  1121. }
  1122. break;
  1123. }
  1124. n++;
  1125. }
  1126. }
  1127. if(cnt > 0) {
  1128. wavesetsampcnt[halfwaveset_cnt] = cnt;
  1129. wavesetmaxsamp[halfwaveset_cnt] = maxsamp;
  1130. }
  1131. if((sndseekEx(dz->ifd[0],0,0) < 0)){
  1132. sprintf(errstr,"sndseek failed\n");
  1133. return SYSTEM_ERROR;
  1134. }
  1135. dz->samps_left = dz->insams[0];
  1136. dz->total_samps_read = 0;
  1137. dz->total_samps_written = 0;
  1138. return FINISHED;
  1139. }
  1140. /****************************** QUIRK *********************************/
  1141. int quirk(dataptr dz)
  1142. {
  1143. int exit_status;
  1144. float *buf = dz->sampbuf[0];
  1145. int *wavesetsampcnt = dz->iparray[0];
  1146. float *wavesetmaxsamp = dz->fptr[0], localmaxsamp;
  1147. double powfac = dz->param[0], samp, maxsamp;
  1148. int bufpos = 0, thiscnt, k, n;
  1149. switch(dz->mode) {
  1150. case(0): // Power factor applied over amplitude range of halfcycle
  1151. if((exit_status = read_samps(buf,dz))<0)
  1152. return(exit_status);
  1153. for(k = 0; k < dz->itemcnt;k++) {
  1154. thiscnt = wavesetsampcnt[k];
  1155. localmaxsamp = wavesetmaxsamp[k];
  1156. for(n = 0; n < thiscnt; n++) {
  1157. if(localmaxsamp != 0.0) {
  1158. samp = buf[bufpos]/localmaxsamp; // +ve/+ve = +ve OR -ve/-ve = +ve : Range 0 to 1
  1159. samp = pow(samp,powfac); // Range 0 to 1
  1160. samp *= localmaxsamp; // Range 0 to local-maxsamp(+ve or -ve)
  1161. buf[bufpos] = (float)samp;
  1162. }
  1163. if(++bufpos >= dz->buflen) {
  1164. if((exit_status = write_samps(buf,dz->buflen,dz))<0)
  1165. return(exit_status);
  1166. if((exit_status = read_samps(buf,dz))<0)
  1167. return(exit_status);
  1168. bufpos = 0;
  1169. }
  1170. }
  1171. }
  1172. break;
  1173. case(1): // Power factor applied over maximum range of entire sound
  1174. maxsamp = 0.0;
  1175. for(k = 0; k < dz->itemcnt;k++) {
  1176. if(fabs(wavesetmaxsamp[k]) > maxsamp)
  1177. maxsamp = fabs(wavesetmaxsamp[k]); // Find maxsamp in file
  1178. }
  1179. if((exit_status = read_samps(buf,dz))<0)
  1180. return(exit_status);
  1181. for(k = 0; k < dz->itemcnt;k++) {
  1182. thiscnt = wavesetsampcnt[k];
  1183. localmaxsamp = wavesetmaxsamp[k];
  1184. for(n = 0; n < thiscnt; n++) {
  1185. samp = fabs(buf[bufpos])/maxsamp; // range minsamp to maxsamp --> range 0 to 1
  1186. samp = pow(samp,powfac); // range 0 to 1
  1187. samp *= maxsamp; // range 0 to maxsamp
  1188. if(localmaxsamp < 0.0)
  1189. samp = -samp;
  1190. buf[bufpos] = (float)samp;
  1191. if(++bufpos >= dz->buflen) {
  1192. if((exit_status = write_samps(buf,dz->buflen,dz))<0)
  1193. return(exit_status);
  1194. if((exit_status = read_samps(buf,dz))<0)
  1195. return(exit_status);
  1196. bufpos = 0;
  1197. }
  1198. }
  1199. }
  1200. break;
  1201. }
  1202. if(bufpos > 0) {
  1203. if((exit_status = write_samps(buf,bufpos,dz))<0)
  1204. return(exit_status);
  1205. }
  1206. return FINISHED;
  1207. }
  1208. /****************************** GET_MODE *********************************/
  1209. int get_the_mode_from_cmdline(char *str,dataptr dz)
  1210. {
  1211. char temp[200], *p;
  1212. if(sscanf(str,"%s",temp)!=1) {
  1213. sprintf(errstr,"Cannot read mode of program.\n");
  1214. return(USAGE_ONLY);
  1215. }
  1216. p = temp + strlen(temp) - 1;
  1217. while(p >= temp) {
  1218. if(!isdigit(*p)) {
  1219. fprintf(stderr,"Invalid mode of program entered.\n");
  1220. return(USAGE_ONLY);
  1221. }
  1222. p--;
  1223. }
  1224. if(sscanf(str,"%d",&dz->mode)!=1) {
  1225. fprintf(stderr,"Cannot read mode of program.\n");
  1226. return(USAGE_ONLY);
  1227. }
  1228. if(dz->mode <= 0 || dz->mode > dz->maxmode) {
  1229. fprintf(stderr,"Program mode value [%d] is out of range [1 - %d].\n",dz->mode,dz->maxmode);
  1230. return(USAGE_ONLY);
  1231. }
  1232. dz->mode--; /* CHANGE TO INTERNAL REPRESENTATION OF MODE NO */
  1233. return(FINISHED);
  1234. }