/dports/science/R-cran-eco/eco/src/ |
H A D | gibbsBase.c | 122 for (i=0; i< n_samp; i++) { in cBaseeco() 144 W[(n_samp+i)][0]=0.0001; in cBaseeco() 147 W[(n_samp+i)][0]=0.9999; in cBaseeco() 149 Wstar[(n_samp+i)][0]=log(W[(n_samp+i)][0])-log(1-W[(n_samp+i)][0]); in cBaseeco() 162 Wstar[(n_samp+x1_samp+i)][1]=log(W[(n_samp+x1_samp+i)][1])-log(1-W[(n_samp+x1_samp+i)][1]); in cBaseeco() 213 for (i=0;i<n_samp;i++){ in cBaseeco() 234 W[n_samp+i][1]=exp(Wstar[n_samp+i][1])/(1+exp(Wstar[n_samp+i][1])); in cBaseeco() 244 W[n_samp+x1_samp+i][0]=exp(Wstar[n_samp+x1_samp+i][0])/(1+exp(Wstar[n_samp+x1_samp+i][0])); in cBaseeco() 288 FreeMatrix(X, n_samp); in cBaseeco() 294 FreeMatrix(W1g, n_samp); in cBaseeco() [all …]
|
H A D | gibbsXBase.c | 144 if (W[(n_samp+i)][0]==0) W[(n_samp+i)][0]=0.0001; in cBaseecoX() 145 if (W[(n_samp+i)][0]==1) W[(n_samp+i)][0]=0.9999; in cBaseecoX() 146 Wstar[(n_samp+i)][0]=log(W[(n_samp+i)][0])-log(1-W[(n_samp+i)][0]); in cBaseecoX() 153 if (W[(n_samp+x1_samp+i)][1]==0) W[(n_samp+x1_samp+i)][1]=0.0001; in cBaseecoX() 154 if (W[(n_samp+x1_samp+i)][1]==1) W[(n_samp+x1_samp+i)][1]=0.9999; in cBaseecoX() 155 Wstar[(n_samp+x1_samp+i)][1]=log(W[(n_samp+x1_samp+i)][1])-log(1-W[(n_samp+x1_samp+i)][1]); in cBaseecoX() 228 W[n_samp+i][1]=exp(Wstar[n_samp+i][1])/(1+exp(Wstar[n_samp+i][1])); in cBaseecoX() 238 W[n_samp+x1_samp+i][0]=exp(Wstar[n_samp+x1_samp+i][0])/(1+exp(Wstar[n_samp+x1_samp+i][0])); in cBaseecoX() 283 FreeMatrix(X, n_samp); in cBaseecoX() 288 FreeMatrix(W1g, n_samp); in cBaseecoX() [all …]
|
H A D | gibbsDP.c | 175 if (W[(n_samp+i)][0]==0) in cDPeco() 178 if (W[(n_samp+i)][0]==1) in cDPeco() 181 Wstar[(n_samp+i)][0]=log(W[(n_samp+i)][0])-log(1-W[(n_samp+i)][0]); in cDPeco() 195 Wstar[(n_samp+x1_samp+i)][1]=log(W[(n_samp+x1_samp+i)][1])-log(1-W[(n_samp+x1_samp+i)][1]); in cDPeco() 282 …dtemp=mu[n_samp+i][1]+Sigma[n_samp+i][0][1]/Sigma[n_samp+i][0][0]*(Wstar[n_samp+i][0]-mu[n_samp+i]… in cDPeco() 283 …dtemp1=Sigma[n_samp+i][1][1]*(1-Sigma[n_samp+i][0][1]*Sigma[n_samp+i][0][1]/(Sigma[n_samp+i][0][0]… in cDPeco() 286 W[n_samp+i][1]=exp(Wstar[n_samp+i][1])/(1+exp(Wstar[n_samp+i][1])); in cDPeco() 292 …dtemp=mu[n_samp+x1_samp+i][0]+Sigma[n_samp+x1_samp+i][0][1]/Sigma[n_samp+x1_samp+i][1][1]*(Wstar[n… in cDPeco() 293 …emp1=Sigma[n_samp+x1_samp+i][0][0]*(1-Sigma[n_samp+x1_samp+i][0][1]*Sigma[n_samp+x1_samp+i][0][1]/… in cDPeco() 296 W[n_samp+x1_samp+i][0]=exp(Wstar[n_samp+x1_samp+i][0])/(1+exp(Wstar[n_samp+x1_samp+i][0])); in cDPeco() [all …]
|
H A D | gibbsZBase.c | 198 if (W[(n_samp+i)][0]==0) W[(n_samp+i)][0]=0.0001; in cBaseecoZ() 199 if (W[(n_samp+i)][0]==1) W[(n_samp+i)][0]=0.9999; in cBaseecoZ() 200 Wstar[(n_samp+i)][0]=log(W[(n_samp+i)][0])-log(1-W[(n_samp+i)][0]); in cBaseecoZ() 208 Wstar[(n_samp+x1_samp+i)][1]=log(W[(n_samp+x1_samp+i)][1])-log(1-W[(n_samp+x1_samp+i)][1]); in cBaseecoZ() 269 dtemp=mu[n_samp+i][1]+Sigma[0][1]/Sigma[0][0]*(Wstar[n_samp+i][0]-mu[n_samp+i][0]); in cBaseecoZ() 274 W[n_samp+i][1]=exp(Wstar[n_samp+i][1])/(1+exp(Wstar[n_samp+i][1])); in cBaseecoZ() 275 Z[(i+n_samp)*n_dim][n_cov]=Wstar[(i+n_samp)][0]; in cBaseecoZ() 276 Z[(i+n_samp)*n_dim+1][n_cov]=Wstar[(i+n_samp)][1]; in cBaseecoZ() 282 …dtemp=mu[n_samp+x1_samp+i][0]+Sigma[0][1]/Sigma[1][1]*(Wstar[n_samp+x1_samp+i][1]-mu[n_samp+x1_sam… in cBaseecoZ() 286 W[n_samp+x1_samp+i][0]=exp(Wstar[n_samp+x1_samp+i][0])/(1+exp(Wstar[n_samp+x1_samp+i][0])); in cBaseecoZ() [all …]
|
H A D | gibbsXDP.c | 176 if (W[(n_samp+i)][0]==0) W[(n_samp+i)][0]=0.0001; in cDPecoX() 177 if (W[(n_samp+i)][0]==1) W[(n_samp+i)][0]=0.9999; in cDPecoX() 178 Wstar[(n_samp+i)][0]=log(W[(n_samp+i)][0])-log(1-W[(n_samp+i)][0]); in cDPecoX() 184 if (W[(n_samp+x1_samp+i)][1]==0) W[(n_samp+x1_samp+i)][1]=0.0001; in cDPecoX() 185 if (W[(n_samp+x1_samp+i)][1]==1) W[(n_samp+x1_samp+i)][1]=0.9999; in cDPecoX() 186 Wstar[(n_samp+x1_samp+i)][1]=log(W[(n_samp+x1_samp+i)][1])-log(1-W[(n_samp+x1_samp+i)][1]); in cDPecoX() 257 if (i<n_samp) in cDPecoX() 276 if (*x1==1 && i>=n_samp && i<(n_samp+x1_samp)) { in cDPecoX() 285 if (*x0==1 && i>=(n_samp+x1_samp) && i<(n_samp+x1_samp+x0_samp)) { in cDPecoX() 444 FreeMatrix(W1g, n_samp); in cDPecoX() [all …]
|
H A D | gibbsBase2C.c | 89 for (i = 0; i < n_samp; i++) in cBase2C() 102 for (i = 0; i < n_samp; i++) in cBase2C() 106 for (i = 0; i < n_samp; i++) in cBase2C() 112 for (i = 0; i < n_samp; i++) { in cBase2C() 151 for (i = 0; i < n_samp; i++){ in cBase2C() 171 for(i = 0; i < n_samp; i++) in cBase2C() 194 FreeMatrix(X, n_samp); in cBase2C() 195 FreeMatrix(W, n_samp); in cBase2C() 196 FreeMatrix(Wstar, n_samp); in cBase2C() 197 FreeMatrix(minU, n_samp); in cBase2C() [all …]
|
H A D | bayes.c | 31 int n_samp, /* sample size */ in NIWupdate() argument 43 for (i=0; i<n_samp; i++) in NIWupdate() 45 Ybar[j] /= n_samp; in NIWupdate() 54 mun[j] = (tau0*mu0[j]+n_samp*Ybar[j])/(tau0+n_samp); in NIWupdate() 57 Sn[j][k] += (tau0*n_samp)*(Ybar[j]-mu0[j])*(Ybar[k]-mu0[k])/(tau0+n_samp); in NIWupdate() 58 for (i=0; i<n_samp; i++) in NIWupdate() 64 rWish(InvSigma, mtemp, nu0+n_samp, n_dim); in NIWupdate() 69 mtemp[j][k] = Sigma[j][k]/(tau0+n_samp); in NIWupdate()
|
H A D | gibbsBaseRC.c | 97 for (i = 0; i < n_samp; i++) in cBaseRC() 103 for (i = 0; i < n_samp; i++) in cBaseRC() 110 for (i = 0; i < n_samp; i++) in cBaseRC() 129 for (i = 0; i < n_samp; i++) { in cBaseRC() 181 for (i = 0; i < n_samp; i++) { in cBaseRC() 236 nu0, S0, n_samp, n_dim); in cBaseRC() 249 for(i = 0; i < n_samp; i++) in cBaseRC() 273 FreeMatrix(X, n_samp); in cBaseRC() 274 FreeMatrix(Y, n_samp); in cBaseRC() 275 Free3DMatrix(W, n_samp, n_dim); in cBaseRC() [all …]
|
H A D | gibbsEM.c | 251 if(i<n_samp) { in cEMeco() 366 n_samp=setP->n_samp; in ecoEStep() 431 for (i=n_samp; i<n_samp+s_samp; i++) { in ecoEStep() 443 for (i=n_samp+s_samp; i<n_samp+s_samp+x1_samp; i++) { in ecoEStep() 454 for (i=n_samp+s_samp+x1_samp; i<n_samp+s_samp+x1_samp+x0_samp; i++) { in ecoEStep() 1165 for (i=n_samp; i<n_samp+s_samp; i++) { in readData() 1180 for (i=n_samp+s_samp; i<n_samp+s_samp+x1_samp; i++) { in readData() 1186 for (i=n_samp+s_samp+x1_samp; i<n_samp+s_samp+x1_samp+x0_samp; i++) { in readData() 1203 int s_max=fmin2(n_samp+x1_samp+x0_samp+s_samp,n_samp+x1_samp+x0_samp+5); in readData() 1420 for (i=0;i<n_samp;i++){ in gridEStep() [all …]
|
H A D | preDP.c | 25 int n_samp = *pin_samp; /* sample size */ in preDP() local 44 for(i=0; i<n_samp; i++) { in preDP()
|
/dports/audio/codec2/codec2-1.0.3/src/ |
H A D | sine.c | 579 int n_samp = c2const->n_samp; in make_synthesis_window() local 588 for(i=n_samp/2-tw; i<n_samp/2+tw; win+=1.0/(2*tw), i++ ) in make_synthesis_window() 590 for(i=n_samp/2+tw; i<3*n_samp/2-tw; i++) in make_synthesis_window() 593 for(i=3*n_samp/2-tw; i<3*n_samp/2+tw; win-=1.0/(2*tw), i++) in make_synthesis_window() 595 for(i=3*n_samp/2+tw; i<2*n_samp; i++) in make_synthesis_window() 612 int n_samp, in synthesise() argument 626 for(i=0; i<n_samp-1; i++) { in synthesise() 627 Sn_[i] = Sn_[i+n_samp]; in synthesise() 629 Sn_[n_samp-1] = 0.0; in synthesise() 665 for(i=n_samp-1,j=0; i<2*n_samp; i++,j++) in synthesise() [all …]
|
H A D | codec2.c | 146 int n_samp = c2->n_samp = c2->c2const.n_samp; in codec2_create() local 175 for(i=0; i<2*n_samp; i++) in codec2_create() 214 for(i=0; i<BPF_N+4*c2->n_samp; i++) in codec2_create() 2047 for(i=0; i<c2->n_samp; i++) { in synthesise_one_frame() 2051 ear_protection(c2->Sn_, c2->n_samp); in synthesise_one_frame() 2053 for(i=0; i<c2->n_samp; i++) { in synthesise_one_frame() 2081 int n_samp = c2->n_samp; in analyse_one_frame() local 2086 for(i=0; i<m_pitch-n_samp; i++) in analyse_one_frame() 2087 c2->Sn[i] = c2->Sn[i+n_samp]; in analyse_one_frame() 2088 for(i=0; i<n_samp; i++) in analyse_one_frame() [all …]
|
/dports/audio/codec2/codec2-1.0.3/unittest/ |
H A D | tnewamp1.c | 45 int n_samp = c2const.n_samp; in main() local 47 short buf[n_samp]; /* input/output buffer */ in main() 124 assert(fread(buf,sizeof(short),n_samp,fin) == n_samp); in main() 128 for(i=0; i<m_pitch-n_samp; i++) { in main() 129 Sn[i] = Sn[i+n_samp]; in main() 131 for(i=0; i<n_samp; i++) { in main() 132 Sn[i+m_pitch-n_samp] = buf[i]; in main() 137 nlp(nlp_states, Sn, n_samp, &pitch, Sw, W, &prev_f0); in main()
|
/dports/math/combblas/CombBLAS_beta_16_2/psort-1.0/include/psort/ |
H A D | psort_samplesort.h | 205 long n_samp = s*k; in parallel_samplesort() local 206 _ValueType *s_splitters = new _ValueType[n_samp]; in parallel_samplesort() 207 sample_splitters(first, n_loc, s_splitters, n_samp); in parallel_samplesort() 213 MPI_Gather(s_splitters, n_samp, MPI_valueType, in parallel_samplesort() 216 _ValueType *gath_splitters = new _ValueType[nproc * n_samp]; in parallel_samplesort() 217 MPI_Gather(s_splitters, n_samp, MPI_valueType, in parallel_samplesort() 218 gath_splitters, n_samp, MPI_valueType, in parallel_samplesort() 222 sort(gath_splitters, gath_splitters + (nproc * n_samp), comp); in parallel_samplesort() 225 double stride = (double) (nproc * n_samp) / (nproc * k); in parallel_samplesort()
|
/dports/graphics/gmt/gmt-6.3.0/src/segy/ |
H A D | segy2grd.c | 290 uint64_t ij, ij0, n_samp = 0, isamp; in GMT_segy2grd() local 447 if ((n_samp = segy_samp_rd (header)) != 0) n_samp = Ctrl->L.value; in GMT_segy2grd() 450 if ((n_samp - ij0) > (uint64_t)Grid->header->n_rows) n_samp = Grid->header->n_rows + ij0; in GMT_segy2grd() 455 for (isamp = 0; isamp < n_samp; ++isamp) { in GMT_segy2grd() 462 …for (ij = ij0; ij < n_samp ; ij++) { /* n*idy is index of first sample to be included in the grid… in GMT_segy2grd() 518 if ((n_samp = segy_samp_rd (header)) != 0) n_samp = Ctrl->L.value; in GMT_segy2grd() 523 for (isamp = 0; isamp < n_samp; ++isamp) { in GMT_segy2grd() 534 for (isamp = 0; isamp< n_samp; ++isamp) { in GMT_segy2grd()
|
H A D | pssegy.c | 333 GMT_LOCAL float pssegy_rms (float *data, uint32_t n_samp) { in pssegy_rms() argument 337 for (ix = 0; ix < n_samp; ix++) in pssegy_rms() 339 sumsq /= ((double) n_samp); in pssegy_rms() 461 …plot_trace (struct GMT_CTRL *GMT, float *data, double dy, double x0, int n_samp, int do_fill, int … in pssegy_plot_trace() argument 468 for (iy = 1; iy < n_samp; iy++) { /* loop over samples on trace - refer to pairs iy-1, iy */ in pssegy_plot_trace() 485 uint32_t n_samp = 0, n_tracelist = 0; in GMT_pssegy() local 700 if ((n_samp = segy_samp_rd (header)) != 0) n_samp = Ctrl->L.value; in GMT_pssegy() 705 for (iy = 0; iy < n_samp; ++iy) { in GMT_pssegy() 713 scale = pssegy_rms (data, n_samp); in GMT_pssegy() 716 for (iy = 0; iy < n_samp; iy++) { /* scale bias and clip each sample in the trace */ in GMT_pssegy() [all …]
|
/dports/biology/py-biom-format/biom-format-2.1.10/biom/cli/ |
H A D | table_head.py | 30 def head(input_fp, output_fp, n_obs, n_samp): argument 40 table = load_table(input_fp).head(n=n_obs, m=n_samp)
|
/dports/devel/aarch64-none-elf-gcc/gcc-8.4.0/gcc/testsuite/g++.dg/graphite/ |
H A D | pr41305.C | 321 foo (int n_samp) 323 vector <double> xi (n_samp, 0); 324 for (int n = 0; n <n_samp; ++n) 326 vector <double> cos_w_n (n_samp); 329 vector <double> cos_wd (n_samp);
|
/dports/lang/gcc6-aux/gcc-6-20180516/gcc/testsuite/g++.dg/graphite/ |
H A D | pr41305.C | 321 foo (int n_samp) 323 vector <double> xi (n_samp, 0); 324 for (int n = 0; n <n_samp; ++n) 326 vector <double> cos_w_n (n_samp); 329 vector <double> cos_wd (n_samp);
|
/dports/lang/gcc9/gcc-9.4.0/gcc/testsuite/g++.dg/graphite/ |
H A D | pr41305.C | 321 foo (int n_samp) 323 vector <double> xi (n_samp, 0); 324 for (int n = 0; n <n_samp; ++n) 326 vector <double> cos_w_n (n_samp); 329 vector <double> cos_wd (n_samp);
|
/dports/devel/arm-none-eabi-gcc492/gcc-4.9.2/gcc/testsuite/g++.dg/graphite/ |
H A D | pr41305.C | 321 foo (int n_samp) 323 vector <double> xi (n_samp, 0); 324 for (int n = 0; n <n_samp; ++n) 326 vector <double> cos_w_n (n_samp); 329 vector <double> cos_wd (n_samp);
|
/dports/devel/avr-gcc/gcc-10.2.0/gcc/testsuite/g++.dg/graphite/ |
H A D | pr41305.C | 321 foo (int n_samp) 323 vector <double> xi (n_samp, 0); 324 for (int n = 0; n <n_samp; ++n) 326 vector <double> cos_w_n (n_samp); 329 vector <double> cos_wd (n_samp);
|
/dports/devel/riscv64-gcc/gcc-8.3.0/gcc/testsuite/g++.dg/graphite/ |
H A D | pr41305.C | 321 foo (int n_samp) 323 vector <double> xi (n_samp, 0); 324 for (int n = 0; n <n_samp; ++n) 326 vector <double> cos_w_n (n_samp); 329 vector <double> cos_wd (n_samp);
|
/dports/lang/gnat_util/gcc-6-20180516/gcc/testsuite/g++.dg/graphite/ |
H A D | pr41305.C | 321 foo (int n_samp) 323 vector <double> xi (n_samp, 0); 324 for (int n = 0; n <n_samp; ++n) 326 vector <double> cos_w_n (n_samp); 329 vector <double> cos_wd (n_samp);
|
/dports/lang/gcc48/gcc-4.8.5/gcc/testsuite/g++.dg/graphite/ |
H A D | pr41305.C | 321 foo (int n_samp) 323 vector <double> xi (n_samp, 0); 324 for (int n = 0; n <n_samp; ++n) 326 vector <double> cos_w_n (n_samp); 329 vector <double> cos_wd (n_samp);
|