Lines Matching refs:ki
166 int ki,kj,kk,kh; /* Counters. */ local
202 for (ki=0; ki<4; ki++) lder[ki] = 2;
204 for (ki=0; ki<=kder; ki++) knmb += ki + 1;
209 for (ki=0; ki<5; ki++)
211 stpos[ki] = (double)0.0;
212 stpos[5+ki] = (double)1.0;
214 for (ki=0; ki<4; ki++)
216 stcrt2[ki] = (double)0.0;
217 stcrt2[4+ki] = (double)1.0;
219 for (ki=0; ki<7; ki++)
221 stcrt1[ki] = (double)0.0;
222 stcrt1[7+ki] = (double)1.0;
234 for (ki=0; ki<icurv; ki++)
238 if (vboundc[ki]->idim != kdim) goto err102;
243 s1221(vboundc[2*ki+1],1,tpar,&kleft,sdum,&kstat);
246 memcopy(stwist+ki*kdim,sdum+kdim,kdim,double);
252 for (ki=0; ki<kdim; ki++)
256 stwist2[(kj*4+2)*kdim+ki] = stwist[kk*kdim+ki]/(double)4.0;
272 for (ki=1; ki<icurv; ki++)
279 memcopy(sder+ki*knmb*kdim,sder,kdim,DOUBLE);
283 memcopy(sder+(ki*knmb+1)*kdim,stang+(ki-1)*kdim,kdim,DOUBLE);
284 memcopy(sder+(ki*knmb+2)*kdim,stang+ki*kdim,kdim,DOUBLE);
290 memcopy(sder+(ki*knmb+3)*kdim,sder+((ki-1)*knmb+5)*kdim,kdim,DOUBLE);
294 sh1461_s9coef(stang+(icurv-1)*kdim,stang,stang+(ki-1)*kdim,
295 stang+ki*kdim,kdim,&tro00,&tro01,&tro10,&tro11,&kstat);
301 sder+(ki*knmb+4)*kdim,&kstat);
306 sder+(ki*knmb+5)*kdim,&kstat);
311 for (ki=0; ki<icurv; ki++)
317 memcopy(spos,sder+ki*knmb*kdim,kdim,DOUBLE);
318 memcopy(spos+kdim,sder+(ki*knmb+2)*kdim,kdim,DOUBLE);
319 memcopy(spos+2*kdim,sder+(ki*knmb+5)*kdim,kdim,DOUBLE);
320 memcopy(scrt2,sder+(ki*knmb+1)*kdim,kdim,DOUBLE);
321 memcopy(scrt2+kdim,sder+(ki*knmb+4)*kdim,kdim,DOUBLE);
322 memcopy(stwist2+4*ki*kdim,sder+(ki*knmb+4)*kdim,kdim,DOUBLE);
327 kj = (ki < icurv-1) ? ki+1 : 0;
346 stwist2[(ki*4+3)*kdim+kh] = scrt2[2*kdim+kh] = -sdum2[kdim+kh]/(double)4.0;
347 stwist2[(kj*4+1)*kdim+kh] = -stwist2[(ki*4+3)*kdim+kh];
360 sh1461_s9chcoor(sblend,4,spos,5,scrt2,4,sder+(knmb*ki+1)*kdim,sder+(knmb*ki+2)*kdim,
370 if ((qbound[8*ki+6] = newCurve(5,5,stpos,spos,1,kdim,1)) == SISL_NULL)
372 if ((qbound[8*ki+7] = newCurve(4,4,stcrt2,scrt2,1,kdim,1)) == SISL_NULL)
377 s1710(qc,tpar,qbound+8*ki+4,qbound+8*kj+2,&kstat);
380 s1710(vboundc[2*kj+1],tpar,qbound+8*ki+5,qbound+8*kj+3,&kstat);
385 s1706(qbound[8*ki+4]);
386 s1706(qbound[8*ki+5]);
390 for (kk=0; kk<kdim*(qbound[8*ki+5]->in); kk++)
391 *(qbound[8*ki+5]->ecoef+kk) *= (double)0.5;
400 kk = (kh > 3) ? 8*ki : 8*kj;
411 for (ki=0; ki<icurv; ki++)
415 s1401(qbound+8*ki,stwist2+4*ki*kdim,vsurf+ki,&kstat);
453 for (ki=0; ki<8*icurv; ki++)
454 if (qbound[ki]) freeCurve(qbound[ki]);
496 int ki; /* Index. */ local
506 for (ki=0; ki<idim; ki++)
508 econd[idim+ki] = ONE_THIRD*econd[idim+ki] + econd[ki];
509 econd[2*idim+ki] = ONE_THIRD*econd[2*idim+ki] + econd[3*idim+ki];
517 for (ki=0; ki<idim; ki++)
519 econd[idim+ki] = ONE_FOURTH*econd[idim+ki] + econd[ki];
520 econd[2*idim+ki] = econd[2*idim+ki]/(double)12.0
521 + (double)2.0*econd[idim+ki] - econd[ki];
522 econd[3*idim+ki] = ONE_FOURTH*econd[3*idim+ki] + econd[4*idim+ki];
739 int ki; /* Counters. */ local
775 for (ki=0; ki<4; ki++) sblend2[ki] = tmy1*eblend[ki] - (double)1.0;
779 for (ki=0; ki<4*idim; ki++)
780 scoef[ki] = (double)4.0*(epos[idim+ki] - epos[ki]);
792 for (ki=0; ki<7*idim; ki++)
793 ecrt2[ki] = (tlambda*sc1[ki] + sc2[ki]);
849 int ki; /* Index. */ local
855 for (ki=0; ki<idim; ki++)
859 ecoefnew[ki] = eblend[0]*ecoef[ki];
860 ecoefnew[idim+ki] = (eblend[0]*ecoef[idim+ki]
861 + eblend[1]*ecoef[ki])/(double)2.0;
862 ecoefnew[2*idim+ki] = (eblend[0]*ecoef[2*idim+ki]
863 + (double)3.0*eblend[1]*ecoef[idim+ki]
864 + eblend[2]*ecoef[ki])/(double)5.0;
865 ecoefnew[3*idim+ki] = (eblend[0]*ecoef[3*idim+ki] + eblend[3]*ecoef[ki]
866 + (double)9.0*(eblend[1]*ecoef[2*idim+ki] +
867 eblend[2]*ecoef[idim+ki]))/(double)20.0;
868 ecoefnew[4*idim+ki] = (eblend[1]*ecoef[3*idim+ki]
869 + (double)3.0*eblend[2]*ecoef[2*idim+ki]
870 + eblend[3]*ecoef[idim+ki])/(double)5.0;
871 ecoefnew[5*idim+ki] = (eblend[2]*ecoef[3*idim+ki]
872 + eblend[3]*ecoef[2*idim+ki])/(double)2.0;
873 ecoefnew[6*idim+ki] = eblend[3]*ecoef[3*idim+ki];