static void Rerror (struct d9r_io *io, int_16 tag, const char *string, int_16 code) { struct d9r_tag_metadata *md = d9r_tag_metadata (io, tag); if (md->aux != (void *)0) { struct d9c_tag_status *mds = (struct d9c_tag_status *)(md->aux); struct d9c_status *status = (struct d9c_status *)(io->aux); switch (mds->code) { case d9c_attaching: status->code = d9c_error; break; case d9c_walking_read: case d9c_walking_create: case d9c_walking_write: case d9c_opening_read: case d9c_opening_write: case d9c_ready_read: case d9c_ready_write: case d9c_ready_write_working: io_finish (mds->io); kill_fid (io, mds->fid); default: break; } if (status->error != (void *)0) { status->error (io, string, status->aux); } } }
int main(void) { mat *m; size_t dim[2]; m = mat_create(6,6); io_init(); io_set_fmt(IO_ASCII); mat_load(m,dim,"ex_io.dat:m"); printf("%dx%d matrix loaded from %s:\n",(int)(dim[0]),(int)(dim[1]),FNAME); mat_print(m,"% .15e"); io_set_fmt(IO_XML); mat_save("ex_io.xml:m",'w',m); io_finish(); mat_destroy(m); return EXIT_SUCCESS; }
int main(int argc, char* argv[]) { /* parsing arguments */ /********************************************/ double latspac_nu; qcd_options *opt; strbuf name[3],manf_name,unit,sink[3],source[3]; size_t binsize,nboot,propdim[2],nt; char *cpt1,*cpt2; corr_no lastc; fit_model *fm_pt; opt = qcd_arg_parse(argc,argv,A_PLOT|A_SAVE_RS|A_LOAD_RG|A_PROP_NAME\ |A_PROP_LOAD|A_LATSPAC|A_QCOMP|A_FIT,2,0); cpt1 = strtok(opt->ss,"/"); printf("%s\n",cpt1); if (cpt1) { cpt2 = strchr(cpt1,':'); if (cpt2 == NULL) { fprintf(stderr,"error: sink/source option %s is invalid\n",\ cpt1); exit(EXIT_FAILURE); } strbufcpy(source[PP],cpt2+1); *cpt2 = '\0'; strbufcpy(sink[PP],cpt1); } else { fprintf(stderr,"error: sink/source option %s is invalid\n",\ opt->ss); exit(EXIT_FAILURE); } cpt1 = strtok(NULL,"/"); if (cpt1) { cpt2 = strchr(cpt1,':'); if (cpt2 == NULL) { fprintf(stderr,"error: sink/source option %s is invalid\n",\ cpt1); exit(EXIT_FAILURE); } strbufcpy(source[AP],cpt2+1); *cpt2 = '\0'; strbufcpy(sink[AP],cpt1); } else { strbufcpy(sink[AP],sink[PP]); strbufcpy(source[AP],source[PP]); } cpt1 = strtok(NULL,"/"); if (cpt1) { cpt2 = strchr(cpt1,':'); if (cpt2 == NULL) { fprintf(stderr,"error: sink/source option %s is invalid\n",\ cpt1); exit(EXIT_FAILURE); } strbufcpy(source[AA],cpt2+1); *cpt2 = '\0'; strbufcpy(sink[AA],cpt1); } else { strbufcpy(sink[AA],sink[PP]); strbufcpy(source[AA],source[PP]); } cpt1 = strtok(opt->channel[0],"/"); if (cpt1) { sprintf(name[PP],"%s_%s_%s_%s",cpt1,opt->quark[0],sink[PP],source[PP]); } else { fprintf(stderr,"error: channel option %s is invalid\n",\ opt->channel[0]); exit(EXIT_FAILURE); } cpt1 = strtok(NULL,"/"); if (cpt1) { sprintf(name[AP],"%s_%s_%s_%s",cpt1,opt->quark[0],sink[AP],source[AP]); } else { fprintf(stderr,"error: channel option %s is invalid\n",\ opt->channel[0]); exit(EXIT_FAILURE); } cpt1 = strtok(NULL,"/"); if (cpt1) { sprintf(name[AA],"%s_%s_%s_%s",cpt1,opt->quark[0],sink[AA],source[AA]); lastc = AA; fm_pt = &fm_pseudosc3; } else { lastc = AP; fm_pt = &fm_pseudosc2; } strbufcpy(manf_name,opt->manf_name); binsize = opt->binsize; nboot = opt->nboot; latspac_nu = opt->latspac_nu; if (opt->have_latspac) { strbufcpy(unit," (MeV)"); } else { strbufcpy(unit,""); } latan_set_verb(opt->latan_verb); minimizer_set_alg(opt->minimizer); mat_ar_loadbin(NULL,propdim,manf_name,name[PP],1); nt = propdim[0]; io_set_fmt(opt->latan_fmt); io_init(); /* loading datas */ /********************************************/ size_t ndat,nbdat; mat **prop[3]; corr_no c; ndat = (size_t)get_nfile(manf_name); nbdat = ndat/binsize + ((ndat%binsize == 0) ? 0 : 1); for (c=PP;c<=lastc;c++) { prop[c] = mat_ar_create(nbdat,propdim[0],propdim[1]); qcd_printf(opt,"-- loading %s datas from %s...\n",name[c],manf_name); mat_ar_loadbin(prop[c],NULL,manf_name,name[c],binsize); } /* propagator */ /********************************************/ rs_sample *s_mprop[3]; mat *mprop[3],*sigmprop[3]; for (c=PP;c<=lastc;c++) { s_mprop[c] = rs_sample_create(propdim[0],propdim[1],nboot); sigmprop[c] = mat_create(propdim[0],propdim[1]); qcd_printf(opt,"-- resampling %s mean propagator...\n",name[c]); randgen_set_state(opt->state); resample(s_mprop[c],prop[c],nbdat,&rs_mean,BOOT,NULL); mprop[c] = rs_sample_pt_cent_val(s_mprop[c]); rs_sample_varp(sigmprop[c],s_mprop[c]); mat_eqsqrt(sigmprop[c]); } /* effective mass */ /********************************************/ rs_sample *s_effmass[3]; mat *tem,*em[3],*sigem[3]; size_t emdim[2]; get_effmass_size(emdim,mprop[PP],1,EM_ACOSH); tem = mat_create(emdim[0],1); for (c=PP;c<=lastc;c++) { s_effmass[c] = rs_sample_create(emdim[0],emdim[1],nboot); sigem[c] = mat_create(emdim[0],emdim[1]); qcd_printf(opt,"-- resampling %s effective mass...\n",name[c]); rs_sample_effmass(s_effmass[c],tem,s_mprop[c],1,EM_ACOSH); em[c] = rs_sample_pt_cent_val(s_effmass[c]); rs_sample_varp(sigem[c],s_effmass[c]); mat_eqsqrt(sigem[c]); } /* fit mass */ /********************************************/ fit_data *d; rs_sample *s_fit; mat *fit,*limit,*sigfit,*scanres_t,*scanres_chi2,*scanres_mass,\ *scanres_masserr; size_t npar,nti,tibeg,range[2],ta; size_t i,j; strbuf buf,range_info,latan_path; double pref_i,mass_i; d = fit_data_create(nt,1,(size_t)(lastc+1)); tibeg = (size_t)(opt->range[0][0]); range[0] = 0; range[1] = 0; npar = fit_model_get_npar(fm_pt,&nt); s_fit = rs_sample_create(npar,1,nboot); fit = rs_sample_pt_cent_val(s_fit); limit = mat_create(npar,2); sigfit = mat_create(npar,1); /** print operation **/ if (!opt->do_range_scan) { qcd_printf(opt,"-- fitting and resampling...\n"); } else { qcd_printf(opt,"-- scanning ranges [ti,%u] from ti= %u\n",\ opt->range[0][1],opt->range[0][0]); opt->nmanrange = 1; } /** check ranges **/ strbufcpy(range_info,""); qcd_printf(opt,"%-20s: ","corrected range(s)"); fit_data_fit_all_points(d,false); for (i=0;i<opt->nmanrange;i++) { sprintf(buf,"_%u_%u",opt->range[i][0],opt->range[i][1]); strbufcat(range_info,buf); range[0] = (size_t)(opt->range[i][0]); range[1] = (size_t)(opt->range[i][1]); correct_range(range,mprop[PP],sigmprop[PP],opt); fit_data_fit_range(d,range[0],range[1],true); } qcd_printf(opt,"\n"); nti = MAX(range[1] - 1 - tibeg,1); scanres_t = mat_create(nti,1); scanres_chi2 = mat_create(nti,1); scanres_mass = mat_create(nti,1); scanres_masserr = mat_create(nti,1); /** set model **/ fit_data_set_model(d,fm_pt,&nt); /** set correlation filter **/ if (opt->corr == NO_COR) { for (i=0;i<nt;i++) for (j=0;j<nt;j++) { if (i != j) { fit_data_set_data_cor(d,i,j,false); } } } /** set initial parameter values **/ ta = nt/8-(size_t)(mat_get(tem,0,0)); mass_i = mat_get(em[PP],ta,0); if (latan_isnan(mass_i)) { mass_i = 0.3; } mat_set(fit,0,0,mass_i); pref_i = mat_get(mprop[PP],ta,0)/(exp(-mass_i*ta)+exp(-mass_i*(nt-ta))); if (latan_isnan(pref_i)) { pref_i = 1.0; } mat_set(fit,1,0,sqrt(pref_i)); mat_set(fit,2,0,sqrt(pref_i)); qcd_printf(opt,"%-22smass= %e -- prefactor_0= %e\n","initial parameters: ",\ mat_get(fit,0,0),pref_i); /** set parameter limits **/ mat_cst(limit,latan_nan()); mat_set(limit,0,0,0.0); mat_set(limit,1,0,0.0); mat_set(limit,2,0,0.0); /** positive AP correlator **/ rs_sample_eqabs(s_mprop[AP]); /** set x **/ for (i=0;i<nt;i++) { fit_data_set_x(d,i,0,(double)(i)-opt->tshift); } /** regular correlator fit... **/ if (!opt->do_range_scan) { latan_set_warn(false); rs_data_fit(s_fit,limit,NULL,s_mprop,d,NO_COR,NULL); latan_set_warn(true); rs_data_fit(s_fit,limit,NULL,s_mprop,d,opt->corr,NULL); rs_sample_varp(sigfit,s_fit); mat_eqsqrt(sigfit); if (fit_data_get_chi2pdof(d) > 2.0) { fprintf(stderr,"warning: bad final fit (chi^2/dof= %.2e)\n",\ fit_data_get_chi2pdof(d)); } qcd_printf(opt,"-- results:\n"); qcd_printf(opt,"%-10s= %.8f +/- %.8e %s\n","mass",\ mat_get(fit,0,0)/latspac_nu, \ mat_get(sigfit,0,0)/latspac_nu,unit); qcd_printf(opt,"%-10s= %.8f +/- %.8e %s\n","decay",\ mat_get(fit,1,0)/latspac_nu, \ mat_get(sigfit,1,0)/latspac_nu,unit); qcd_printf(opt,"%-10s= %.8f +/- %.8e %s\n","norm",\ mat_get(fit,2,0)/latspac_nu, \ mat_get(sigfit,2,0)/latspac_nu,unit); qcd_printf(opt,"%-10s= %d\n","dof",fit_data_get_dof(d)); qcd_printf(opt,"%-10s= %e\n","chi^2/dof",fit_data_get_chi2pdof(d)); if (opt->do_save_rs_sample) { sprintf(latan_path,"%s_pseudosc_fit%s_%s.boot:%s_pseudosc_fit%s_%s",\ opt->quark[0],range_info,manf_name,opt->quark[0],\ range_info,manf_name); rs_sample_save_subsamp(latan_path,'w',s_fit,0,0,1,0); } } /** ...or fit range scanning **/ else { qcd_printf(opt,"\n%-5s %-12s a*M_%-8s %-12s","ti/a","chi^2/dof",\ opt->quark[0],"error"); for (i=tibeg;i<range[1]-1;i++) { latan_set_warn(false); rs_data_fit(s_fit,limit,NULL,s_mprop,d,NO_COR,NULL); latan_set_warn(true); rs_data_fit(s_fit,limit,NULL,s_mprop,d,opt->corr,NULL); rs_sample_varp(sigfit,s_fit); mat_eqsqrt(sigfit); mat_set(scanres_t,i-tibeg,0,(double)(i)); mat_set(scanres_chi2,i-tibeg,0,fit_data_get_chi2pdof(d)); mat_set(scanres_mass,i-tibeg,0,mat_get(fit,0,0)); mat_set(scanres_masserr,i-tibeg,0,mat_get(sigfit,0,0)); qcd_printf(opt,"\n% -4d % -.5e % -.5e % -.5e",(int)(i),\ fit_data_get_chi2pdof(d),mat_get(fit,0,0), \ mat_get(sigfit,0,0)); fit_data_fit_point(d,i,false); } qcd_printf(opt,"\n\n"); } /* plot */ /********************************************/ if (opt->do_plot) { mat *mbuf,*em_i,*sigem_i,*par,*ft[3],*comp[3]; plot *p; strbuf key,dirname,color; size_t maxt,t,npoint; double dmaxt,nmass; mbuf = mat_create(1,1); em_i = mat_create(nrow(em[PP]),1); sigem_i = mat_create(nrow(em[PP]),1); par = mat_create(2,1); if (!opt->do_range_scan) { maxt = nt; dmaxt = (double)maxt; npoint = fit_data_fit_point_num(d); /** chi^2 plot **/ p = plot_create(); i = 0; for (c=PP;c<=lastc;c++) { ft[c] = mat_create(npoint,1); comp[c] = mat_create(npoint,1); for (t=0;t<nt;t++) { if (fit_data_is_fit_point(d,t)) { mat_set(ft[c],i%npoint,0,(double)(t)+0.33*(double)(c)); mat_set(comp[c],i%npoint,0,mat_get(d->chi2_comp,i,0)); i++; } } } plot_set_scale_manual(p,-1.0,dmaxt,-5.0,5.0); plot_add_plot(p,"0.0 lt -1 lc rgb 'black' notitle",""); plot_add_plot(p,"1.0 lt -1 lc rgb 'black' notitle",""); plot_add_plot(p,"-1.0 lt -1 lc rgb 'black' notitle",""); plot_add_plot(p,"2.0 lt -1 lc rgb 'dark-gray' notitle",""); plot_add_plot(p,"-2.0 lt -1 lc rgb 'dark-gray' notitle",""); plot_add_plot(p,"3.0 lt -1 lc rgb 'gray' notitle",""); plot_add_plot(p,"-3.0 lt -1 lc rgb 'gray' notitle",""); plot_add_plot(p,"4.0 lt -1 lc rgb 'light-gray' notitle",""); plot_add_plot(p,"-4.0 lt -1 lc rgb 'light-gray' notitle",""); plot_set_ylabel(p,"standard deviations"); for (c=PP;c<=lastc;c++) { sprintf(color,"%d",(int)(c)+1); plot_add_points(p,ft[c],comp[c],c_name[c],color,"impulses"); } plot_disp(p); if (opt->do_save_plot) { sprintf(dirname,"%s_dev",opt->save_plot_dir); plot_save(dirname,p); } plot_destroy(p); for (c=PP;c<=lastc;c++) { /** propagator plot **/ p = plot_create(); fit_data_fit_all_points(d,true); plot_set_scale_ylog(p); plot_set_scale_xmanual(p,0,dmaxt); sprintf(key,"%s %s propagator",opt->quark[0],c_name[c]); mat_eqabs(mprop[c]); plot_add_fit(p,d,c,mbuf,0,fit,0,dmaxt,1000,false,\ PF_FIT|PF_DATA,key,"","rgb 'red'","rgb 'red'"); plot_disp(p); if (opt->do_save_plot) { sprintf(dirname,"%s_prop_%s",opt->save_plot_dir,c_name[c]); plot_save(dirname,p); } plot_destroy(p); /** effective mass plot **/ p = plot_create(); mat_eqmuls(em[c],1.0/latspac_nu); mat_eqmuls(sigem[c],1.0/latspac_nu); nmass = mat_get(fit,0,0)/latspac_nu; plot_add_hlineerr(p,nmass, mat_get(sigfit,0,0)/latspac_nu,\ "rgb 'red'"); sprintf(key,"%s %s effective mass",opt->quark[0],c_name[c]); plot_add_dat(p,tem,em[c],NULL,sigem[c],key,"rgb 'blue'"); plot_disp(p); if (opt->do_save_plot) { sprintf(dirname,"%s_em_%s",opt->save_plot_dir,c_name[c]); plot_save(dirname,p); } plot_destroy(p); } } else { /* chi^2 plot */ p = plot_create(); plot_set_scale_manual(p,0,(double)(nt/2),0,5.0); plot_add_hline(p,1.0,"rgb 'black'"); plot_add_dat(p,scanres_t,scanres_chi2,NULL,NULL,"chi^2/dof",\ "rgb 'blue'"); plot_disp(p); if (opt->do_save_plot) { sprintf(dirname,"%s_chi2",opt->save_plot_dir); plot_save(dirname,p); } plot_destroy(p); /* mass plot */ p = plot_create(); plot_set_scale_xmanual(p,0,(double)(nt/2)); sprintf(key,"a*M_%s",opt->quark[0]); plot_add_dat(p,scanres_t,scanres_mass,NULL,scanres_masserr,key,\ "rgb 'red'"); plot_disp(p); if (opt->do_save_plot) { sprintf(dirname,"%s_mass",opt->save_plot_dir); plot_save(dirname,p); } plot_destroy(p); } mat_destroy(em_i); mat_destroy(sigem_i); mat_destroy(mbuf); mat_destroy(par); for (c=PP;c<=lastc;c++) { mat_destroy(ft[c]); mat_destroy(comp[c]); } } /* desallocation */ /********************************************/ free(opt); io_finish(); mat_ar_destroy(prop[0],nbdat); mat_ar_destroy(prop[1],nbdat); for (c=PP;c<=lastc;c++) { rs_sample_destroy(s_mprop[c]); mat_destroy(sigmprop[c]); rs_sample_destroy(s_effmass[c]); mat_destroy(sigem[c]); } mat_destroy(tem); fit_data_destroy(d); rs_sample_destroy(s_fit); mat_destroy(limit); mat_destroy(sigfit); mat_destroy(scanres_t); mat_destroy(scanres_chi2); mat_destroy(scanres_mass); mat_destroy(scanres_masserr); return EXIT_SUCCESS; }
int main(int argc, char *argv[]) { int i,j; rs_sample *s1,*s2,*res; size_t s1_dim[2],s2_dim[2],s1_nsample,s2_nsample; mat *sig; bool do_save_res, show_usage; strbuf out_elname,out_fname,in_elname[2],in_fname[2],in_path[2],out_path; char *spt[2]; io_fmt_no fmt; /* argument parsing */ i = 1; j = 0; show_usage = false; do_save_res = false; fmt = io_get_fmt(); if (argc <= 1) { show_usage = true; } else { while (i < argc) { if (strcmp(argv[i],"-o") == 0) { if (i == argc - 1) { show_usage = true; break; } else { strbufcpy(out_path,argv[i+1]); do_save_res = true; i += 2; } } else if (strcmp(argv[i],"-f") == 0) { if (i == argc - 1) { show_usage = true; break; } else { if (strcmp(argv[i+1],"xml") == 0) { fmt = IO_XML; } else if (strcmp(argv[i+1],"ascii") == 0) { fmt = IO_ASCII; } else { fprintf(stderr,"error: format %s unknown\n",argv[i+1]); return EXIT_FAILURE; } i += 2; } } else { if (j < 2) { strbufcpy(in_path[j],argv[i]); j++; i++; } else { show_usage = true; break; } } } } if (show_usage) { fprintf(stderr,"usage: %s <in sample 1> <in sample 2> [-o <out sample>] [-f {ascii|xml}]\n",\ argv[0]); return EXIT_FAILURE; } /* I/O init */ io_set_fmt(fmt); io_init(); /* getting sizes */ rs_sample_load(NULL,&s1_nsample,s1_dim,in_path[0]); rs_sample_load(NULL,&s2_nsample,s2_dim,in_path[1]); /* allocation */ s1 = rs_sample_create(s1_dim[0],s1_dim[1],s1_nsample); s2 = rs_sample_create(s2_dim[0],s2_dim[1],s2_nsample); res = rs_sample_create(s1_dim[0],s1_dim[1],s1_nsample); sig = mat_create(s1_dim[0],s1_dim[1]); /* loading samples */ printf("-- loading sample from %s...\n",in_path[0]); rs_sample_load(s1,NULL,NULL,in_path[0]); printf("-- loading sample from %s...\n",in_path[1]); rs_sample_load(s2,NULL,NULL,in_path[1]); /* multiplying samples */ printf("-- executing operation on samples...\n"); rs_sample_binop(res,s1,s2,&BINOP); /* result output */ rs_sample_varp(sig,res); mat_eqsqrt(sig); printf("central value:\n"); mat_print(rs_sample_pt_cent_val(res),"%e"); printf("standard deviation:\n"); mat_print(sig,"%e"); if (do_save_res) { get_elname(out_fname,out_elname,out_path); if (strlen(out_elname) == 0) { get_elname(in_fname[0],in_elname[0],in_path[0]); get_elname(in_fname[1],in_elname[1],in_path[1]); spt[0] = (strlen(in_elname[0]) == 0) ? in_fname[0] : in_elname[0]; spt[1] = (strlen(in_elname[1]) == 0) ? in_fname[1] : in_elname[1]; sprintf(out_path,"%s%c%s_%s_%s",out_fname,LATAN_PATH_SEP,argv[0],\ spt[0],spt[1]); } rs_sample_save(out_path,'w',res); } /* desallocation */ rs_sample_destroy(s1); rs_sample_destroy(s2); rs_sample_destroy(res); mat_destroy(sig); /* I/O finish */ io_finish(); return EXIT_SUCCESS; }
int main(int argc, char *argv[]) { rs_sample *s1,*res,*buf; size_t s1_dim[2],res_dim[2],s1_nsample,a,b,s,k,l; int i,j; mat *sig; bool do_save_res, show_usage, have_s; char *spt; strbuf out_elname,out_fname,in_elname,in_fname,in_path,out_path; io_fmt_no fmt; /* argument parsing */ i = 1; j = 0; a = 0; b = 0; s = 1; show_usage = false; do_save_res = false; have_s = false; fmt = io_get_fmt(); if (argc <= 3) { show_usage = true; } else { while (i < argc) { if (strcmp(argv[i],"-o") == 0) { if (i == argc - 1) { show_usage = true; break; } else { strbufcpy(out_path,argv[i+1]); do_save_res = true; i += 2; } } else if (strcmp(argv[i],"-f") == 0) { if (i == argc - 1) { show_usage = true; break; } else { if (strcmp(argv[i+1],"xml") == 0) { fmt = IO_XML; } else if (strcmp(argv[i+1],"ascii") == 0) { fmt = IO_ASCII; } else { fprintf(stderr,"error: format %s unknown\n",argv[i+1]); return EXIT_FAILURE; } i += 2; } } else if (strcmp(argv[i],"-s") == 0) { if (i == argc - 1) { show_usage = true; break; } else { s = (size_t)atol(argv[i+1]); have_s = true; i += 2; } } else { if (j == 0) { strbufcpy(in_path,argv[i]); j++; i++; } else if (j == 1) { a = (size_t)atol(argv[i]); j++; i++; } else if (j == 2) { b = (size_t)atol(argv[i]); j++; i++; } else { show_usage = true; break; } } } } if (show_usage) { fprintf(stderr,"usage: %s <in sample> <a> <b> [-o <out sample>] [-f {ascii|xml}]\n",\ argv[0]); return EXIT_FAILURE; } /* I/O init */ io_set_fmt(fmt); io_init(); /* getting sizes */ rs_sample_load(NULL,&s1_nsample,s1_dim,in_path); l = b - a + 1; s = (have_s) ? s : l; res_dim[0] = s; res_dim[1] = s1_dim[1]; /* allocation */ s1 = rs_sample_create(s1_dim[0],s1_dim[1],s1_nsample); buf = rs_sample_create(1,s1_dim[1],s1_nsample); res = rs_sample_create(res_dim[0],res_dim[1],s1_nsample); sig = mat_create(res_dim[0],res_dim[1]); /* loading samples */ printf("-- loading sample from %s...\n",in_path); rs_sample_load(s1,NULL,NULL,in_path); /* multiplying samples */ printf("-- taking subsample [%d,%d]...\n",(int)a,(int)b); for (k=0;k<s;k++) { rs_sample_get_subsamp(buf,s1,a+(k%l),0,a+(k%l),s1_dim[1]-1); rs_sample_set_subsamp(res,buf,k,0,k,s1_dim[1]-1); } /* result output */ rs_sample_varp(sig,res); mat_eqsqrt(sig); printf("central value:\n"); mat_print(rs_sample_pt_cent_val(res),"%e"); printf("standard deviation:\n"); mat_print(sig,"%e"); if (do_save_res) { get_elname(out_fname,out_elname,out_path); if (strlen(out_elname) == 0) { get_elname(in_fname,in_elname,in_path); spt = (strlen(in_elname) == 0) ? in_fname : in_elname; sprintf(out_path,"%s%c%s_%s",out_fname,LATAN_PATH_SEP,argv[0],spt); } rs_sample_save(out_path,'w',res); } /* desallocation */ rs_sample_destroy(s1); rs_sample_destroy(res); rs_sample_destroy(buf); mat_destroy(sig); /* I/O finish */ io_finish(); return EXIT_SUCCESS; }