k                 242 src/CalculiX.h              ITG *inewton,ITG *nset,ITG *ifreebody,ITG *k));
k                2148 src/CalculiX.h              double *rwork,ITG *lrw,ITG *k,ITG *jout,double *rpar,
k                  82 src/arpack.c     ITG *inum=NULL,k,ido,ldz,iparam[11],ipntr[14],lworkl,ngraph=1,im,
k                 161 src/arpack.c         for(k=1;k<=*nbody;k++){
k                 163 src/arpack.c   			     iendset,ialset,&inewton,nset,&ifreebody,&k));
k                 211 src/arpack.c   		  for(k=mi[2]**ne;k<mi[2]*(*ne+*nslavs);k++) ielorien[k]=0;
k                 215 src/arpack.c   	      for(k=mi[2]**ne;k<mi[2]*(*ne+*nslavs);k++) ielmat[k]=1;
k                 219 src/arpack.c   		  for(k=*nstate_*mi[0]**ne;k<*nstate_*mi[0]*(*ne+*nslavs);k++){
k                 220 src/arpack.c   		      xstate[k]=0.;
k                 268 src/arpack.c   		  for(k=mi[2]**ne;k<mi[2]*(*ne+*nintpoint);k++) ielorien[k]=0;
k                 271 src/arpack.c   	      for(k=mi[2]**ne;k<mi[2]*(*ne+*nintpoint);k++) ielmat[k]=1;
k                 295 src/arpack.c   		      for(k=*nstate_*mi[0]*ne0;
k                 296 src/arpack.c   			  k<*nstate_*mi[0]*(ne0+maxprevcontel);++k){
k                 297 src/arpack.c   			  xstateini[k]=xstate[k];
k                 302 src/arpack.c   		  for(k=*nstate_*mi[0]*ne0;k<*nstate_*mi[0]*(ne0+*nintpoint);k++){
k                 303 src/arpack.c   		      xstate[k]=0.;
k                 325 src/arpack.c   		  for(k=0;k<*nstate_*mi[0]*(ne0+*nintpoint);++k){
k                 326 src/arpack.c   		      xstateini[k]=xstate[k];
k                 347 src/arpack.c         for(k=0;k<*nstate_*mi[0]*(ne0+*nslavs);++k){
k                 348 src/arpack.c         xstateini[k]=xstate[k];
k                 494 src/arpack.c       NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                 839 src/arpack.c       for(k=0;k<6*mi[0]*ne0;k++){eme[k]=0.;}
k                 842 src/arpack.c       for(k=0;k<neq[1];++k)
k                 843 src/arpack.c         temp_array[k]=0.;
k                 849 src/arpack.c       for(k=0;k<neq[1];++k)
k                 850 src/arpack.c         sum+=z[lint+k]*temp_array[k];
k                 851 src/arpack.c       for(k=0;k<neq[1];++k)
k                 852 src/arpack.c         z[lint+k]=z[lint+k]/sqrt(sum);
k                 947 src/arpack.c   	for(k=0;k<*nstate_*mi[0]*(ne0+maxprevcontel);++k){
k                 948 src/arpack.c   	    xstate[k]=xstateini[k];
k                1003 src/arpack.c     for(k=0;k<6*mi[0]*ne0;k++){eme[k]=0.;}
k                  72 src/arpackbu.c   ITG *inum=NULL,k,ido,dz,iparam[11],ipntr[11],lworkl,im,nasym=0,
k                 121 src/arpackbu.c       for(k=1;k<=*nbody;k++){
k                 123 src/arpackbu.c 			     iendset,ialset,&inewton,nset,&ifreebody,&k));
k                 231 src/arpackbu.c   for(k=0;k<neq[0];++k){
k                 232 src/arpackbu.c       b[k]=fext[k]-f[k];
k                 241 src/arpackbu.c     NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                 357 src/arpackbu.c   for(k=0;k<mt**nk;++k){
k                 358 src/arpackbu.c     vold[k]=v[k];
k                 644 src/arpackbu.c     for(k=0;k<6*mi[0]**ne;k++){eme[k]=0.;}
k                  83 src/arpackcs.c   ITG *inum=NULL,k,ido,ldz,iparam[11],ipntr[14],lworkl,idir,nherm=1,
k                 164 src/arpackcs.c       for(k=1;k<=*nbody;k++){
k                 166 src/arpackcs.c 			     iendset,ialset,&inewton,nset,&ifreebody,&k));
k                 214 src/arpackcs.c 		  for(k=mi[2]**ne;k<mi[2]*(*ne+*nslavs);k++) ielorien[k]=0;
k                 218 src/arpackcs.c 	      for(k=mi[2]**ne;k<mi[2]*(*ne+*nslavs);k++) ielmat[k]=1;
k                 222 src/arpackcs.c 		  for(k=*nstate_*mi[0]**ne;k<*nstate_*mi[0]*(*ne+*nslavs);k++){
k                 223 src/arpackcs.c 		      xstate[k]=0.;
k                 271 src/arpackcs.c 		  for(k=mi[2]**ne;k<mi[2]*(*ne+*nintpoint);k++) ielorien[k]=0;
k                 274 src/arpackcs.c 	      for(k=mi[2]**ne;k<mi[2]*(*ne+*nintpoint);k++) ielmat[k]=1;
k                 298 src/arpackcs.c 		      for(k=*nstate_*mi[0]*ne0;
k                 299 src/arpackcs.c 			  k<*nstate_*mi[0]*(ne0+maxprevcontel);++k){
k                 300 src/arpackcs.c 			  xstateini[k]=xstate[k];
k                 305 src/arpackcs.c 		  for(k=*nstate_*mi[0]*ne0;k<*nstate_*mi[0]*(ne0+*nintpoint);k++){
k                 306 src/arpackcs.c 		      xstate[k]=0.;
k                 328 src/arpackcs.c 		  for(k=0;k<*nstate_*mi[0]*(ne0+*nintpoint);++k){
k                 329 src/arpackcs.c 		      xstateini[k]=xstate[k];
k                 350 src/arpackcs.c       for(k=0;k<*nstate_*mi[0]*(ne0+*nslavs);++k){
k                 351 src/arpackcs.c 	  xstateini[k]=xstate[k];
k                 558 src/arpackcs.c 	  NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                 635 src/arpackcs.c       for(k=0;k<11;k++)iparam[k]=0;
k                 921 src/arpackcs.c 	  for(k=0;k<neq2;k++){
k                 922 src/arpackcs.c 	      amp=z[2*j*neq[1]+k]*z[2*j*neq[1]+k]+
k                 923 src/arpackcs.c 		  z[2*j*neq[1]+neq2+k]*z[2*j*neq[1]+neq2+k];
k                 924 src/arpackcs.c 	      if(amp>ampmax){ampmax=amp;kmax1=k;}
k                 928 src/arpackcs.c 	  for(k=0;k<neq2;k++){
k                 929 src/arpackcs.c 	      amp=z[(2*j+1)*neq[1]+k]*z[(2*j+1)*neq[1]+k]+
k                 930 src/arpackcs.c 		  z[(2*j+1)*neq[1]+neq2+k]*z[(2*j+1)*neq[1]+neq2+k];
k                 931 src/arpackcs.c 	      if(amp>ampmax){ampmax=amp;kmax2=k;}
k                1244 src/arpackcs.c 	  for(k=0;k<neq[1];k+=neq[1]/2){
k                1248 src/arpackcs.c 	      if(k==0) {kk=0;kkv=0;kk6=0;kkx=0;if(*nprint>0)FORTRAN(writere,());}
k                1285 src/arpackcs.c 			  if(k==0) {
k                1304 src/arpackcs.c 		    f,&fn[kkv],nactdof,&iout,qa,vold,&z[lint+k],
k                1323 src/arpackcs.c 		    f,&fn[kkv],nactdof,&iout,qa,vold,&z[lint+k],
k                2235 src/arpackcs.c 	for(k=0;k<*nstate_*mi[0]*(ne0+maxprevcontel);++k){
k                2236 src/arpackcs.c 	    xstate[k]=xstateini[k];
k                  42 src/calcresidual.c     ITG j,k,mt=mi[1]+1;
k                  48 src/calcresidual.c 	for(k=0;k<neq[1];++k){
k                  49 src/calcresidual.c 	    b[k]=fext[k]-f[k];
k                  56 src/calcresidual.c 	for(k=0;k<*nk;++k){
k                  57 src/calcresidual.c 	    if(nactdof[mt*k]!=0){
k                  58 src/calcresidual.c 		aux2[nactdof[mt*k]-1]=(vold[mt*k]-vini[mt*k])/(*dtime);}
k                  60 src/calcresidual.c 		if(nactdof[mt*k+j]!=0){aux2[nactdof[mt*k+j]-1]=accold[mt*k+j];}
k                  69 src/calcresidual.c 	for(k=0;k<neq[0];++k){
k                  70 src/calcresidual.c 	    b[k]=scal1*(fext[k]-f[k])-*alpha*(fextini[k]-fini[k])-b[k];
k                  72 src/calcresidual.c 	for(k=neq[0];k<neq[1];++k){
k                  73 src/calcresidual.c 	    b[k]=fext[k]-f[k]-b[k];
k                  79 src/calcresidual.c 	    for(k=0;k<*nk;++k){
k                  80 src/calcresidual.c 		if(nactdof[mt*k]!=0){aux2[nactdof[mt*k]-1]=0.;}
k                  82 src/calcresidual.c 		    if(nactdof[mt*k+j]!=0){
k                  83 src/calcresidual.c 			aux2[nactdof[mt*k+j]-1]=veold[mt*k+j];}
k                  91 src/calcresidual.c 	    for(k=0;k<neq[0];++k){
k                  92 src/calcresidual.c 		b[k]-=scal1*cv[k]-*alpha*cvini[k];
k                 100 src/calcresidual.c 	for(k=0;k<*nk;++k){
k                 101 src/calcresidual.c 	    if(nactdof[mt*k]!=0){
k                 102 src/calcresidual.c 		aux2[nactdof[mt*k]-1]=(vold[mt*k]-vini[mt*k])/(*dtime);}
k                 104 src/calcresidual.c 		if(nactdof[mt*k+j]!=0){aux2[nactdof[mt*k+j]-1]=accold[mt*k+j];}
k                 108 src/calcresidual.c 	for(k=0;k<neq[0];++k){
k                 109 src/calcresidual.c 	    b[k]=scal1*(fext[k]-f[k])-*alpha*(fextini[k]-fini[k])
k                 110 src/calcresidual.c 		-adb[k]*aux2[k];
k                 112 src/calcresidual.c 	for(k=neq[0];k<neq[1];++k){
k                 113 src/calcresidual.c 	    b[k]=fext[k]-f[k]-adb[k]*aux2[k];
k                  41 src/calcresidual_em.c     ITG j,k,mt=mi[1]+1,jstart;
k                  46 src/calcresidual_em.c 	for(k=0;k<neq[1];++k){
k                  47 src/calcresidual_em.c 	    b[k]=fext[k]-f[k];
k                  63 src/calcresidual_em.c 	for(k=0;k<*nk;++k){
k                  65 src/calcresidual_em.c 		if(nactdof[mt*k+j]!=0){aux2[nactdof[mt*k+j]-1]=(vold[mt*k+j]-vini[mt*k+j])/(*dtime);}
k                  77 src/calcresidual_em.c 	for(k=0;k<neq[1];++k){
k                  78 src/calcresidual_em.c 	    b[k]=fext[k]-f[k]-b[k];
k                  79 src/ccx_2.8.c    k,nzs[3],nmpc_=0,nload_=0,nforc_=0,istep,istat,nboun_=0,nintpoint=0,
k                 618 src/ccx_2.8.c  		for(k=0;k<6;k++){
k                 619 src/ccx_2.8.c  		    sti[6*mi[0]*i+6*j+k]=prestr[6*mi[0]*i+6*j+k];
k                 628 src/ccx_2.8.c  		for(k=0;k<6;k++){
k                 629 src/ccx_2.8.c  		    eme[6*mi[0]*i+6*j+k]=prestr[6*mi[0]*i+6*j+k];
k                1302 src/ccx_2.8.c  	  for(k=0;k<nk;++k){
k                1303 src/ccx_2.8.c  	      t1[k]=t0[k];
k                  56 src/checkconvergence.c 	ngraph=1,k,*ipneigh=NULL,*neigh=NULL,*inum=NULL,id,istart,iend,inew,
k                 167 src/checkconvergence.c 		for(k=0;k<*nk;++k){t1[k]=vold[mt*k];}
k                 206 src/checkconvergence.c 		    NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                 328 src/checkconvergence.c 	    NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                 460 src/checkconvergence.c 		    NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                 511 src/checkconvergence.c 			NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                 538 src/checkconvergence.c 			NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                 607 src/checkconvergence.c 			NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                 631 src/checkconvergence.c 			NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                  77 src/complexfreq.c   ITG nev,i,j,k,idof,*inum=NULL,*ipobody=NULL,inewton=0,id,
k                 351 src/complexfreq.c 	      for(k=0;k<neq[1];k++){
k                 352 src/complexfreq.c 		  z[j*neq[1]+k]=z[i*neq[1]+k];
k                 489 src/complexfreq.c       for(k=1;k<=*nbody;k++){
k                 491 src/complexfreq.c 			     iendset,ialset,&inewton,nset,&ifreebody,&k));
k                 531 src/complexfreq.c 	      for(k=0;k<neq[1];k++){
k                 532 src/complexfreq.c 		  cc[i*nev+j]+=z[j*neq[1]+k]*zc[i*neq[1]+k];
k                 638 src/complexfreq.c       for(k=0;k<nev;k++){
k                 639 src/complexfreq.c 	  eigxr[k]=xx[2*l*nev+2*k];
k                 640 src/complexfreq.c 	  eigxi[k]=xx[2*l*nev+2*k+1];
k                 656 src/complexfreq.c 	      for(k=0;k<nev;k++){
k                 661 src/complexfreq.c 		      (xx[2*l*nev+2*k]*z[(long long)2*k*neqact+i]
k                 662 src/complexfreq.c 		       -xx[2*l*nev+2*k+1]*z[(long long)(2*k+1)*neqact+i]);
k                 667 src/complexfreq.c 		      (xx[2*l*nev+2*k]*z[(long long)(2*k+1)*neqact+i]
k                 668 src/complexfreq.c 		       +xx[2*l*nev+2*k+1]*z[(long long)2*k*neqact+i]);
k                 694 src/complexfreq.c 	  for(k=0;k<neq[1];k++){
k                 695 src/complexfreq.c 	      sum+=zz[l*neq[1]+k]*z[k];
k                 699 src/complexfreq.c 	  for(k=0;k<neq[1];k++){
k                 700 src/complexfreq.c 	      zz[l*neq[1]+k]/=sum;
k                 721 src/complexfreq.c 		  for(k=0;k<neq[1];k++){
k                 722 src/complexfreq.c 		      xmr[l*nev+m]+=zz[m*neq[1]+k]*z[k];
k                 727 src/complexfreq.c 	      for(k=0;k<neqact;k++){zi[neqact+k]=-zz[2*l*neqact+k];}
k                 731 src/complexfreq.c 		  for(k=0;k<neq[1];k++){
k                 732 src/complexfreq.c 		      xmi[l*nev+m]+=zz[m*neq[1]+k]*z[k];
k                 772 src/complexfreq.c 		  for(k=0;k<nev;k++){
k                 773 src/complexfreq.c 		      zz[(2*l+j)*neq[1]+i]+=xx[2*l*nev+2*k+j]*
k                 774 src/complexfreq.c                             z[(long long)k*neq[1]+i];
k                 804 src/complexfreq.c 	  for(k=0;k<neq[1];k++){
k                 805 src/complexfreq.c 	      sum+=zz[2*l*neq[1]+k]*z[k];
k                 812 src/complexfreq.c 	  for(k=0;k<neq[1];k++){
k                 813 src/complexfreq.c 	      sum+=zz[(2*l+1)*neq[1]+k]*z[k];
k                 817 src/complexfreq.c 	  for(k=0;k<2*neq[1];k++){
k                 818 src/complexfreq.c 	      zz[2*l*neq[1]+k]/=sum;
k                 845 src/complexfreq.c 		  for(k=0;k<neq[1];k++){
k                 846 src/complexfreq.c 		      xmr[l*nev+m]+=zz[2*m*neq[1]+k]*z[k];
k                 848 src/complexfreq.c 		  for(k=0;k<neq[1];k++){
k                 849 src/complexfreq.c 		      xmi[l*nev+m]-=zz[(2*m+1)*neq[1]+k]*z[k];
k                 861 src/complexfreq.c 		  for(k=0;k<neq[1];k++){
k                 862 src/complexfreq.c 		      xmr[l*nev+m]+=zz[(2*m+1)*neq[1]+k]*z[k];
k                 864 src/complexfreq.c 		  for(k=0;k<neq[1];k++){
k                 865 src/complexfreq.c 		      xmi[l*nev+m]+=zz[2*m*neq[1]+k]*z[k];
k                1272 src/complexfreq.c     for(k=0;k<*mcs;k++){
k                1273 src/complexfreq.c 	theta=nm[j]*2.*pi/cs[17*k];
k                1274 src/complexfreq.c 	cs[17*k+14]=cos(theta);
k                1275 src/complexfreq.c 	cs[17*k+15]=sin(theta);
k                1287 src/complexfreq.c     for(k=0;k<neq[1];k+=neq[1]/2){
k                1291 src/complexfreq.c       if(k==0) {kk=0;kkv=0;kk6=0;kkx=0;if(*nprint>0)FORTRAN(writere,());}
k                1328 src/complexfreq.c             if(k==0) {
k                1347 src/complexfreq.c 	    f,&fn[kkv],nactdof,&iout,qa,vold,&zz[lint+k],
k                1366 src/complexfreq.c             f,&fn[kkv],nactdof,&iout,qa,vold,&zz[lint+k],
k                  35 src/dfdbj.c      ITG j,j1,jdof,kdof,k,k1,l,id,index,ist,id1,ist1,index1,id2,ist2,index2,
k                  62 src/dfdbj.c    		  k=*nactcont-1;
k                  63 src/dfdbj.c    		  l=k-1;
k                  64 src/dfdbj.c    		  while(k>id){
k                  65 src/dfdbj.c    		      ikactcont[k]=ikactcont[l];
k                  66 src/dfdbj.c    		      ilactcont[k--]=ilactcont[l--];
k                  78 src/dfdbj.c    	      for(k=0; k<*nope; k++){
k                  80 src/dfdbj.c    		      sl=s[(3*k+k1)*100+i3];
k                  81 src/dfdbj.c    		      kdof=nactdof[mt*(konl[k]-1)+k1+1];
k                 100 src/dfdbj.c    			  kdof=8*(konl[k]-1)+k1+1;
k                 171 src/dfdbj.c    				  k=*nactcont-1;
k                 172 src/dfdbj.c    				  l=k-1;
k                 174 src/dfdbj.c    				      ikactcont[k]=ikactcont[l];
k                 175 src/dfdbj.c    				      ilactcont[k--]=ilactcont[l--];
k                 176 src/dfdbj.c    				  }while(k>id);
k                 187 src/dfdbj.c    			      for(k=0; k<*nope; k++){
k                 189 src/dfdbj.c    				      sl=s[(3*k+k1)*100+i3];
k                 190 src/dfdbj.c    				      kdof=nactdof[mt*(konl[k]-1)+k1+1];
k                 210 src/dfdbj.c    					  kdof=8*(konl[k]-1)+k1+1;
k                  77 src/dyna.c       ITG nev,i,j,k,idof,*inum=NULL,*ipobody=NULL,inewton=0,id,
k                 679 src/dyna.c     	      for(k=0;k<neq[1];k++){
k                 680 src/dyna.c     		  cc[i*nev+j]+=z[j*neq[1]+k]*zc[i*neq[1]+k];
k                 730 src/dyna.c     	  for(k=mi[2]**ne;k<mi[2]*(*ne+*nslavs);k++) ielorien[k]=0;
k                 733 src/dyna.c           for(k=mi[2]**ne;k<mi[2]*(*ne+*nslavs);k++) ielmat[k]=1;
k                 741 src/dyna.c     	  for(k=*nstate_*mi[0]**ne;k<*nstate_*mi[0]*(*ne+*nslavs);k++){
k                 742 src/dyna.c     	      xstate[k]=0.;
k                 747 src/dyna.c     	  for(k=0;k<*nstate_*mi[0]*(*ne+*nslavs);++k){
k                 748 src/dyna.c     	      xstateini[k]=xstate[k];
k                 827 src/dyna.c           for(k=0;k<nev;k++){
k                 828 src/dyna.c     	cd[k]+=z[k*neq[1]+i]*temp_array2[i];
k                 847 src/dyna.c           for(k=0;k<nev;k++){
k                 848 src/dyna.c     	cv[k]+=z[k*neq[1]+i]*temp_array2[i];
k                 864 src/dyna.c       for(k=0;k<7**nbody;k++){xbodyact[k]=xbody[k];}
k                 865 src/dyna.c       for(k=0;k<7**nbody;k++){xbodydiff[k]=xbody[k];}
k                 876 src/dyna.c           for(k=1;k<=*nbody;k++){
k                 878 src/dyna.c     			     iendset,ialset,&inewton,nset,&ifreebody,&k));
k                1218 src/dyna.c     	    for(k=0;k<*nstate_*mi[0]*(ne0+*nslavs);++k){
k                1219 src/dyna.c     		xstateini[k]=xstate[k];
k                1763 src/dyna.c     		for(k=0;k<ne0;k++){
k                1764 src/dyna.c     		    if(ipkon[k]<-1){
k                1767 src/dyna.c     		    }else if(ipkon[k]!=-1){
k                1768 src/dyna.c     			ipkon[k]=-ipkon[k]-2;
k                1771 src/dyna.c     		for(k=0;k<nmdelem;k++){
k                1772 src/dyna.c     		    ielem=imdelem[k]-1;
k                1802 src/dyna.c     		for(k=0;k<ne0;k++){
k                1803 src/dyna.c     		    if(ipkon[k]<-1){ipkon[k]=-2-ipkon[k];}
k                1809 src/dyna.c     	  for(k=0;k<6*mi[0]*ne0;++k){
k                1810 src/dyna.c     	    sti[k]=stx[k];
k                1915 src/dyna.c       for(k=0;k<*nboun;++k){xboun[k]=xbounact[k];}
k                1916 src/dyna.c       for(k=0;k<*nforc;++k){xforc[k]=xforcact[k];}
k                1917 src/dyna.c       for(k=0;k<2**nload;++k){xload[k]=xloadact[k];}
k                1918 src/dyna.c       for(k=0;k<7**nbody;k=k+7){xbody[k]=xbodyact[k];}
k                1920 src/dyna.c         for(k=0;k<*nk;++k){t1[k]=t1act[k];}
k                 100 src/dynacont.c   ITG i,j,k,l,init,*itg=NULL,ntg=0,maxlenmpc,icascade=0,loop,
k                 138 src/dynacont.c   for(k=0;k<7**nbody;k++){xbodydiff[k]=xbody[k];}
k                 217 src/dynacont.c 	  k=*nactcont-1;
k                 218 src/dynacont.c 	  l=k-1;
k                 219 src/dynacont.c 	  while(k>id){
k                 220 src/dynacont.c 	      ikactcont[k--]=ikactcont[l--];
k                 268 src/dynacont.c 	for(k=0;k<*nstate_*mi[0]*(*ne0+*nslavs);++k){
k                 269 src/dynacont.c 	  xstate[k]=xstateini[k];
k                 769 src/dynacont.c 	  for(k=0; k<*nactcont; k++){
k                 770 src/dynacont.c 	      i1=ikactcont[k];
k                 771 src/dynacont.c 	      i2=(ilactcont[k]-1)**nev;
k                 781 src/dynacont.c 	  for(k=0; k<*nactcont; k++){
k                 782 src/dynacont.c 	      i1=ikactcont[k];
k                 783 src/dynacont.c 	      i2=(ilactcont[k]-1)**nev;
k                  84 src/electromagnetics.c   ITG *inum=NULL,k,iout=0,icntrl,iinc=0,jprint=0,iit=-1,jnz=0,
k                 149 src/electromagnetics.c   for(k=0;k<3;k++){
k                 150 src/electromagnetics.c       strcpy1(&jobnamef[k*132],&jobnamec[k*132],132);
k                 204 src/electromagnetics.c   for(k=0;k<*nboun;++k){xbounact[k]=xbounold[k];}
k                 209 src/electromagnetics.c   for(k=0;k<7**nbody;k++){xbodyact[k]=xbody[k];}
k                 216 src/electromagnetics.c       for(k=1;k<=*nbody;k++){
k                 218 src/electromagnetics.c 			     iendset,ialset,&inewton,nset,&ifreebody,&k));
k                 479 src/electromagnetics.c   for(k=0;k<neq[1];++k){
k                 480 src/electromagnetics.c       b[k]=fext[k]-f[k];
k                 718 src/electromagnetics.c       for(k=index;k<memmpc_;k++){
k                 719 src/electromagnetics.c 	  nodempc[3*k-1]=k+1;
k                 775 src/electromagnetics.c       for(k=0;k<*nk;++k){t1act[k]=t1old[k];}
k                 793 src/electromagnetics.c 	  for(k=0;k<*nboun;++k){xbounini[k]=xbounact[k];}
k                 795 src/electromagnetics.c 	      for(k=0;k<*nk;++k){t1ini[k]=t1act[k];}
k                 797 src/electromagnetics.c 	  for(k=0;k<neq[1];++k){
k                 798 src/electromagnetics.c 	      fini[k]=f[k];
k                 801 src/electromagnetics.c 	      for(k=0;k<mt**nk;++k){
k                 802 src/electromagnetics.c 		  veini[k]=veold[k];
k                 804 src/electromagnetics.c 	      for(k=0;k<neq[1];++k){
k                 805 src/electromagnetics.c 		  fextini[k]=fext[k];
k                 999 src/electromagnetics.c 	      for(k=0;k<*nk;++k){t1act[k]=vold[mt*k];}
k                1044 src/electromagnetics.c 	      NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                1064 src/electromagnetics.c 		  for(k=0;k<neq[0];++k){
k                1065 src/electromagnetics.c 		      ad[k]=adb[k]/dtime+ad[k];
k                1067 src/electromagnetics.c 		  for(k=0;k<nzs[0];++k){
k                1068 src/electromagnetics.c 		      au[k]=aub[k]/dtime+au[k];
k                1074 src/electromagnetics.c 		      for(k=nzs[2];k<nzs[2]+nzs[0];++k){
k                1075 src/electromagnetics.c 			  au[k]=aub[k]/dtime+au[k];
k                1083 src/electromagnetics.c 		  for(k=neq[0];k<neq[1];++k){
k                1084 src/electromagnetics.c 		      ad[k]=adb[k]/dtime+ad[k];
k                1086 src/electromagnetics.c 		  for(k=nzs[0];k<nzs[1];++k){
k                1087 src/electromagnetics.c 		      au[k]=aub[k]/dtime+au[k];
k                1093 src/electromagnetics.c 		      for(k=nzs[2]+nzs[0];k<nzs[2]+nzs[1];++k){
k                1094 src/electromagnetics.c 			  au[k]=aub[k]/dtime+au[k];
k                1222 src/electromagnetics.c 	  for(k=0;k<2;++k){
k                1223 src/electromagnetics.c 	      ram2[k]=ram1[k];
k                1224 src/electromagnetics.c 	      ram1[k]=ram[k];
k                1225 src/electromagnetics.c 	      ram[k]=0.;
k                1231 src/electromagnetics.c 	      for(k=neq[0];k<neq[1];++k){
k                1232 src/electromagnetics.c 		  err=fabs(b[k]);
k                1233 src/electromagnetics.c 		  if(err>ram[1]){ram[1]=err;ram[3]=k+0.5;}
k                1295 src/electromagnetics.c 	  for(k=0;k<*nboun;++k){xbounact[k]=xbounini[k];}
k                1297 src/electromagnetics.c 	      for(k=0;k<*nk;++k){t1act[k]=t1ini[k];}
k                1299 src/electromagnetics.c 	  for(k=0;k<neq[1];++k){
k                1300 src/electromagnetics.c 	      f[k]=fini[k];
k                1303 src/electromagnetics.c 	      for(k=0;k<mt**nk;++k){
k                1304 src/electromagnetics.c 		  veold[k]=veini[k];
k                1306 src/electromagnetics.c 	      for(k=0;k<neq[1];++k){
k                1307 src/electromagnetics.c 		  fext[k]=fextini[k];
k                1481 src/electromagnetics.c       for(k=0;k<mt**nk;++k){veold[k]=0.;}
k                1488 src/electromagnetics.c   for(k=0;k<*nboun;++k){
k                1493 src/electromagnetics.c       if(ndirboun[k]==0){
k                1500 src/electromagnetics.c       }else if((ndirboun[k]>0)&&(ndirboun[k]<4)){
k                1501 src/electromagnetics.c 	  node=nodeboun[k];
k                1509 src/electromagnetics.c       xbounold[k]=xbounact[k];
k                1511 src/electromagnetics.c   for(k=0;k<*nforc;++k){xforcold[k]=xforcact[k];}
k                1512 src/electromagnetics.c   for(k=0;k<2**nload;++k){xloadold[k]=xloadact[k];}
k                1513 src/electromagnetics.c   for(k=0;k<7**nbody;k=k+7){xbodyold[k]=xbodyact[k];}
k                1515 src/electromagnetics.c       for(k=0;k<*nk;++k){t1old[k]=t1act[k];}
k                1516 src/electromagnetics.c       for(k=0;k<*nk;++k){vold[mt*k]=t1act[k];}
k                1519 src/electromagnetics.c       for(k=0;k<*nk;++k){t1[k]=vold[mt*k];}
k                1521 src/electromagnetics.c 	  for(k=0;k<*nk;++k){t1old[k]=t1act[k];}
k                  63 src/exo.c        ITG i,j,k,l,m,n,o,indexe,nemax,nlayer,noutloc,iset,iselect,ncomp,nope,
k                 182 src/exo.c      	  for(k=0;k<mi[2];k++){
k                 183 src/exo.c      	    if(ielmat[i*mi[2]+k]==0) break;
k                 187 src/exo.c      	  for(k=0;k<nlayer;k++){
k                 273 src/exo.c      	  for(k=0;k<mi[2];k++){
k                 274 src/exo.c      	    if(ielmat[i*mi[2]+k]==0) break;
k                 277 src/exo.c      	  for(k=0;k<nlayer;k++){
k                 397 src/exo.c      	    for(k=0;k<mi[2];k++){
k                 398 src/exo.c      	      if(ielmat[i*mi[2]+k]==0) break;
k                 413 src/exo.c            k=0; o=0;
k                 421 src/exo.c      	    for(m=0;m<12;m++){connect[k++] = node_map_inv[kon[indexe+m]-1];}
k                 422 src/exo.c      	    for(m=16;m<20;m++){connect[k++] = node_map_inv[kon[indexe+m]-1];}
k                 423 src/exo.c      	    for(m=12;m<16;m++){connect[k++] = node_map_inv[kon[indexe+m]-1];}
k                 425 src/exo.c      	    for(m=0;m<9;m++){connect[k++] = node_map_inv[kon[indexe+m]-1];}
k                 426 src/exo.c      	    for(m=12;m<15;m++){connect[k++] = node_map_inv[kon[indexe+m]-1];}
k                 427 src/exo.c      	    for(m=9;m<12;m++){connect[k++] = node_map_inv[kon[indexe+m]-1];}
k                 435 src/exo.c      	      for(m=0;m<12;m++){connect[k++] = node_map_inv[kon[indexe+28+20*n+m]-1];}
k                 436 src/exo.c      	      for(m=16;m<20;m++){connect[k++] = node_map_inv[kon[indexe+28+20*n+m]-1];}
k                 437 src/exo.c      	      for(m=12;m<16;m++){connect[k++] = node_map_inv[kon[indexe+28+20*n+m]-1];}
k                 441 src/exo.c      	      connect[k++] = node_map_inv[kon[indexe+20+j]-1];
k                 445 src/exo.c      	      connect[k++] = node_map_inv[kon[indexe+j]-1];
k                  47 src/exoselect.c   ITG i,j,k,l,m,o,nksegment;
k                  90 src/exoselect.c       for(k=istartset[*iset-1]-1;k<iendset[*iset-1];k++){
k                  91 src/exoselect.c 	if(ialset[k]>0){
k                  95 src/exoselect.c 	    i=ialset[k]+l*nksegment-1;
k                 117 src/exoselect.c 	  l=ialset[k-2];
k                 119 src/exoselect.c 	    l-=ialset[k];
k                 120 src/exoselect.c 	    if(l>=ialset[k-1]) break;
k                 152 src/exoselect.c       k=j;
k                 154 src/exoselect.c       k=5;
k                 156 src/exoselect.c       k=4;
k                 157 src/exoselect.c     }else{k=j;}
k                 159 src/exoselect.c     int errr = ex_put_nodal_var (exoid, time_step, k+1+countvar, nout, nodal_var_vals);
k                  32 src/exosetfind.c   ITG i,j,k,l,n,s,e,gen;
k                  72 src/exosetfind.c 	    for (k=ialset[j]; k<=ialset[j+1]; k-=ialset[j+2]){
k                  73 src/exosetfind.c 	      set_nums[n++]=exoset_check(k-1, node_map_inv, nk, &dropped, &unidentified);
k                  35 src/exovector.c   ITG i,j,k,l,m,n,ii,jj,kk;
k                  84 src/exovector.c       for(k=istartset[*iset-1]-1;k<iendset[*iset-1];k++){
k                  85 src/exovector.c 	if(ialset[k]>0){
k                  88 src/exovector.c 	    i=ialset[k]+l*nksegment-1;
k                 101 src/exovector.c 	  l=ialset[k-2];
k                 103 src/exovector.c 	    l-=ialset[k];
k                 104 src/exovector.c 	    if(l>=ialset[k-1]) break;
k                  71 src/expand.c       ITG *inum=NULL,k,idir,lfin,j,iout=0,index,inode,id,i,idof,im,
k                 467 src/expand.c   	for(k=0;k<2*neqh;k+=neqh){
k                 471 src/expand.c   	    if(k==0) {kk=0;kkv=0;kk6=0;}
k                 503 src/expand.c   			if(k==0) {
k                 523 src/expand.c   	      f,&fn[kkv],nactdof,&iout,qa,vold,&z[lint+k],
k                 611 src/expand.c   		for(k=0;k<*nsectors;k++){
k                 613 src/expand.c   		    idof=nactdof[mt*(k**nk+i)+j1]-1;
k                 619 src/expand.c   				    zdof[(long long)j**nzdof+id-1]=vt[k*mt**nk+mt*i+j1];
k                 621 src/expand.c   				    zdof[(long long)2*j**nzdof+id-1]=vt[k*mt**nk+mt*i+j1];
k                 622 src/expand.c   				    zdof[(long long)(2*j+1)**nzdof+id-1]=vti[k*mt**nk+mt*i+j1];
k                 728 src/expand.c   		    for(k=0;k<numnodes;k++){
k                 729 src/expand.c   			sum+=coefright[k];
k                 731 src/expand.c   		    for(k=0;k<numnodes;k++){
k                 732 src/expand.c   			coefright[k]/=sum;
k                 739 src/expand.c   		for(k=0;k<numnodes;k++){
k                 741 src/expand.c   		    nodempc[3**mpcend]=noderight[k]+ileft**nk;
k                 744 src/expand.c   		    coefmpc[*mpcend]=-coefright[k];
k                 794 src/expand.c   				for(k=0;k<3;k++){
k                 795 src/expand.c   				    if(fabs(a[3*k+idir-1]-coef)<1.e-10){
k                 798 src/expand.c   					coef=a[3*k+idir-1];
k                 799 src/expand.c   					locdir=k;
k                 802 src/expand.c   				    if(fabs(a[3*k+idir-1]+coef)<1.e-10){
k                 805 src/expand.c   					coef=-a[3*k+idir-1];
k                 806 src/expand.c   					locdir=k;
k                 888 src/expand.c   	      for(k=0; k<iendset[j]-istartset[j]+1; k++){
k                 889 src/expand.c   		ialset[tint-1]=ialset[istartset[j]-1+k];
k                 894 src/expand.c   		for(k=0; k<iendset[j]-istartset[j]+1; k++){
k                 895 src/expand.c   		  ialset[tint-1]=(ialset[istartset[j]-1+k]!=-1)?ialset[istartset[j]-1+k]+*nk*l:-1;
k                 907 src/expand.c   	      for(k=0; k<iendset[j]-istartset[j]+1; k++){
k                 908 src/expand.c   		ialset[tint-1]=ialset[istartset[j]-1+k];
k                 913 src/expand.c   		for(k=0; k<iendset[j]-istartset[j]+1; k++){
k                 914 src/expand.c   		  tint2=((ITG)(ialset[istartset[j]-1+k]))/10;
k                 915 src/expand.c   		  ialset[tint-1]=(ialset[istartset[j]-1+k]!=-1)?(tint2+*ne*l)*10+(ialset[istartset[j]-1+k]-(tint2*10)):-1;
k                  64 src/frd.c        ITG null,one,i,j,k,indexe,nemax,nlayer,noutloc,iset,iselect,ncomp,nope,
k                 339 src/frd.c      	  for(k=0;k<mi[2];k++){
k                 340 src/frd.c      	    if(ielmat[i*mi[2]+k]==0) break;
k                 343 src/frd.c      	  for(k=0;k<nlayer;k++){
k                 348 src/frd.c      	      for(j=0;j<10;j++)fprintf(f1,"%10" ITGFORMAT "",kon[indexe+28+20*k+j]);
k                 350 src/frd.c      	      for(j=10;j<12;j++)fprintf(f1,"%10" ITGFORMAT "",kon[indexe+28+20*k+j]);
k                 351 src/frd.c      	      for(j=16;j<19;j++)fprintf(f1,"%10" ITGFORMAT "",kon[indexe+28+20*k+j]);
k                 352 src/frd.c      	      for(j=19;j<20;j++)fprintf(f1,"%10" ITGFORMAT "",kon[indexe+28+20*k+j]);
k                 353 src/frd.c      	      for(j=12;j<16;j++)fprintf(f1,"%10" ITGFORMAT "",kon[indexe+28+20*k+j]);
k                 360 src/frd.c      	      for(j=0;j<10;j++){iw=(int)kon[indexe+28+20*k+j];
k                 362 src/frd.c      	      for(j=10;j<12;j++){iw=(int)kon[indexe+28+20*k+j];
k                 364 src/frd.c      	      for(j=16;j<19;j++){iw=(int)kon[indexe+28+20*k+j];
k                 366 src/frd.c      	      for(j=19;j<20;j++){iw=(int)kon[indexe+28+20*k+j];
k                 368 src/frd.c      	      for(j=12;j<16;j++){iw=(int)kon[indexe+28+20*k+j];
k                  43 src/frdselect.c   ITG i,j,k,l,m,n,nksegment;
k                 116 src/frdselect.c     for(k=istartset[*iset-1]-1;k<iendset[*iset-1];k++){
k                 117 src/frdselect.c       if(ialset[k]>0){
k                 119 src/frdselect.c 	  i=ialset[k]+l*nksegment-1;
k                 184 src/frdselect.c 	l=ialset[k-2];
k                 186 src/frdselect.c 	  l-=ialset[k];
k                 187 src/frdselect.c 	  if(l>=ialset[k-1]) break;
k                  28 src/frdset.c     ITG j,k;
k                  73 src/frdset.c   	k=ialset[j-2];
k                  75 src/frdset.c   	  k=k-ialset[j];
k                  76 src/frdset.c   	  if(k>=ialset[j-1]) break;
k                  78 src/frdset.c   	    if(inum[k-1]<0) (*noutloc)++;
k                  80 src/frdset.c   	    if(inum[k-1]>0) (*noutloc)++;
k                  82 src/frdset.c   	    if(inum[k-1]!=0) (*noutloc)++;
k                  29 src/frdvector.c   ITG i,k,l,m,nksegment;
k                  87 src/frdvector.c     for(k=istartset[*iset-1]-1;k<iendset[*iset-1];k++){
k                  88 src/frdvector.c       if(ialset[k]>0){
k                  90 src/frdvector.c 	  i=ialset[k]+l*nksegment-1;
k                 121 src/frdvector.c 	l=ialset[k-2];
k                 123 src/frdvector.c 	  l-=ialset[k];
k                 124 src/frdvector.c 	  if(l>=ialset[k-1]) break;
k                  45 src/inicont.c        *ilmpc=NULL,*nodempc=NULL,nmpc_,i,j,k,ncone;
k                  72 src/linstatic.c   ITG *inum=NULL,k,*icol=NULL,*irow=NULL,ielas,icmd=0,iinc=1,nasym=0,i,j,ic,ir,
k                 121 src/linstatic.c   for(k=0;k<*nboun;++k){xbounact[k]=xbounold[k];}
k                 126 src/linstatic.c   for(k=0;k<7**nbody;k++){xbodyact[k]=xbody[k];}
k                 129 src/linstatic.c     for(k=0;k<*nk;++k){t1act[k]=t1old[k];}
k                 137 src/linstatic.c       for(k=1;k<=*nbody;k++){
k                 139 src/linstatic.c 			     iendset,ialset,&inewton,nset,&ifreebody,&k));
k                 227 src/linstatic.c   for(k=0;k<*neq;++k){
k                 228 src/linstatic.c       b[k]=fext[k]-f[k];
k                 515 src/linstatic.c     NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                 537 src/linstatic.c   for(k=0;k<*nboun;++k){xbounold[k]=xbounact[k];}
k                 538 src/linstatic.c   for(k=0;k<*nforc;++k){xforcold[k]=xforcact[k];}
k                 539 src/linstatic.c   for(k=0;k<2**nload;++k){xloadold[k]=xloadact[k];}
k                 540 src/linstatic.c   for(k=0;k<7**nbody;k=k+7){xbodyold[k]=xbodyact[k];}
k                 542 src/linstatic.c     for(k=0;k<*nk;++k){t1old[k]=t1act[k];}
k                 543 src/linstatic.c     for(k=0;k<*nk;++k){vold[mt*k]=t1act[k];}
k                  41 src/mastruct.c   ITG i,j,k,l,jj,ll,id,index,jdof1,jdof2,idof1,idof2,mpc1,mpc2,id1,id2,
k                  91 src/mastruct.c 	      for(k=1;k<4;++k){
k                  92 src/mastruct.c 		  nactdof[mt*node+k]=1;
k                 249 src/mastruct.c 	k=jj-mi[1]*j;
k                 252 src/mastruct.c 	jdof1=nactdof[mt*(node1-1)+k+1];
k                 274 src/mastruct.c 	      idof2=8*node1+k-7;}
k                 313 src/mastruct.c 	    idof1=8*node1+k-7;
k                 560 src/mastruct.c 		for(k=0;k<4;++k){
k                 561 src/mastruct.c 		    if(nactdof[mt*j+k]==i+1){
k                 563 src/mastruct.c 			idof1=k;
k                 631 src/mastruct.c     k=0;
k                 637 src/mastruct.c 	if(k!=mast1[i]){
k                 638 src/mastruct.c 	  for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
k                 639 src/mastruct.c 	  k=mast1[i];
k                 641 src/mastruct.c 	++icol[k-1];
k                 645 src/mastruct.c     for(l=k;l<neq[1]+1;++l){jq[l]=nmast+1;}
k                 674 src/mastruct.c 	    if(k!=irow[i]){
k                 675 src/mastruct.c 		for(l=k;l<irow[i];++l){jq[l]=i+1-isubtract;}
k                 676 src/mastruct.c 		k=irow[i];
k                 678 src/mastruct.c 	    ++icol[k-1];
k                 681 src/mastruct.c 	for(l=k;l<neq[2]+1;++l){jq[l]=nmastboun+1;}
k                  41 src/mastructcs.c   ITG i,j,k,l,jj,ll,id,index,jdof1,jdof2,idof1,idof2,mpc1,mpc2,id1,id2,
k                  93 src/mastructcs.c       for(k=1;k<4;++k){
k                  94 src/mastructcs.c 	nactdof[mt*node+k]=1;
k                 178 src/mastructcs.c       k=jj-3*j;
k                 181 src/mastructcs.c       jdof1=nactdof[mt*(node1-1)+k+1];
k                 205 src/mastructcs.c 	    idof2=8*node1+k-7;}
k                 260 src/mastructcs.c 	  idof1=8*node1+k-7;
k                 503 src/mastructcs.c   k=0;
k                 509 src/mastructcs.c       if(k!=mast1[i]){
k                 510 src/mastructcs.c 	for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
k                 511 src/mastructcs.c 	k=mast1[i];
k                 513 src/mastructcs.c       ++icol[k-1];
k                 517 src/mastructcs.c   for(l=k;l<neq[0]+1;++l){jq[l]=nmast+1;}
k                  40 src/mastructem.c   ITG i,j,k,l,jj,ll,id,index,jdof1,jdof2,idof1,idof2,mpc1,mpc2,id1,id2,
k                  93 src/mastructem.c 	      for(k=jmin;k<jmax;++k){
k                  94 src/mastructem.c 		  nactdof[mt*node+k]=1;
k                 234 src/mastructem.c 	k=jj-mi[1]*j;
k                 237 src/mastructem.c 	jdof1=nactdof[mt*(node1-1)+k+1];
k                 259 src/mastructem.c 	      idof2=8*node1+k-7;}
k                 289 src/mastructem.c 	    idof1=8*node1+k-7;
k                 526 src/mastructem.c 		for(k=0;k<4;++k){
k                 527 src/mastructem.c 		    if(nactdof[mt*j+k]==i+1){
k                 529 src/mastructem.c 			idof1=k;
k                 585 src/mastructem.c     k=0;
k                 591 src/mastructem.c 	if(k!=mast1[i]){
k                 592 src/mastructem.c 	  for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
k                 593 src/mastructem.c 	  k=mast1[i];
k                 595 src/mastructem.c 	++icol[k-1];
k                 599 src/mastructem.c     for(l=k;l<neq[1]+1;++l){jq[l]=nmast+1;}
k                  32 src/mastructf.c   ITG i,j,k,l,index,idof1,idof2,node1,isubtract,nmast,ifree=0,istart,istartold,
k                 124 src/mastructf.c   k=0;
k                 130 src/mastructf.c 	  if(k!=mast1[i]){
k                 131 src/mastructf.c 	      for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
k                 132 src/mastructf.c 	      k=mast1[i];
k                 134 src/mastructf.c 	  ++icol[k-1];
k                 138 src/mastructf.c   for(l=k;l<*neq+1;++l){jq[l]=nmast+1;}
k                  35 src/mastructrad.c   ITG three=3,i,j,k,l,ii,jj,icav,jcav,*mast1rad=NULL,*irowrad=NULL,
k                 116 src/mastructrad.c   k=0;
k                 122 src/mastructrad.c 	  if(k!=mast1rad[i]){
k                 123 src/mastructrad.c 	      for(l=k;l<mast1rad[i];++l){jqrad[l]=i+1-isubtract;}
k                 124 src/mastructrad.c 	      k=mast1rad[i];
k                 126 src/mastructrad.c 	  ++icolrad[k-1];
k                 130 src/mastructrad.c   for(l=k;l<*ntr+1;++l){jqrad[l]=nmast+1;}
k                  38 src/matrixstorage.c   long long *ipoint=NULL,k;
k                  91 src/matrixstorage.c     k=0;
k                  93 src/matrixstorage.c       ai[k]=i+1;
k                  94 src/matrixstorage.c       aj[k]=i+1;
k                  95 src/matrixstorage.c       aa[k]=ad[i];
k                  96 src/matrixstorage.c       k++;
k                 101 src/matrixstorage.c         ai[k]=i+1;
k                 102 src/matrixstorage.c         aj[k]=irow[l];
k                 103 src/matrixstorage.c         aa[k]=au[l];
k                 104 src/matrixstorage.c         k++;l++;
k                 132 src/matrixstorage.c     k=0;
k                 134 src/matrixstorage.c       ai[k]=i+1;
k                 135 src/matrixstorage.c       aj[k]=i+1;
k                 136 src/matrixstorage.c       aa[k]=ad[i];
k                 137 src/matrixstorage.c       k++;
k                 142 src/matrixstorage.c         ai[k]=i+1;
k                 143 src/matrixstorage.c         aj[k]=irow[l];
k                 144 src/matrixstorage.c         aa[k]=au[l];
k                 145 src/matrixstorage.c         k++;
k                 146 src/matrixstorage.c         ai[k]=irow[l];
k                 147 src/matrixstorage.c         aj[k]=i+1;
k                 148 src/matrixstorage.c         aa[k]=au[l];
k                 149 src/matrixstorage.c         k++;
k                 156 src/matrixstorage.c     k=0;
k                 158 src/matrixstorage.c       jref=aj[k];
k                 159 src/matrixstorage.c       kstart=k;
k                 161 src/matrixstorage.c         k++;
k                 162 src/matrixstorage.c 	if(k==ndim2) break;
k                 163 src/matrixstorage.c         if(aj[k]!=jref) break;
k                 165 src/matrixstorage.c       klen=k-kstart;
k                 192 src/matrixstorage.c       k=(long long)i3c*neq3+i3l;
k                 193 src/matrixstorage.c       FORTRAN(nidentll,(ipoint,&k,&npoint,&id));
k                 196 src/matrixstorage.c         ipoint[npoint-1]=k;
k                 199 src/matrixstorage.c       else if(ipoint[id-1]!=k){
k                 207 src/matrixstorage.c         ipoint[npoint-1]=k;
k                 260 src/matrixstorage.c     k=-1;
k                 270 src/matrixstorage.c         k++;
k                 271 src/matrixstorage.c         ai[k]=il;
k                 272 src/matrixstorage.c         aj[k]=ic;
k                 273 src/matrixstorage.c         aa[k]=aa3[9*i+j];
k                 281 src/matrixstorage.c   k=0;
k                 283 src/matrixstorage.c     jref=aj[k];
k                 284 src/matrixstorage.c     kstart=k;
k                 286 src/matrixstorage.c       k++;
k                 287 src/matrixstorage.c       if(aj[k]!=jref) break;
k                 289 src/matrixstorage.c     klen=k-kstart;
k                 321 src/matrixstorage.c     k=0;
k                 323 src/matrixstorage.c       ai[k]=i+1;
k                 324 src/matrixstorage.c       aj[k]=i+1;
k                 325 src/matrixstorage.c       aa[k]=adb[i];
k                 326 src/matrixstorage.c       k++;
k                 331 src/matrixstorage.c         ai[k]=i+1;
k                 332 src/matrixstorage.c         aj[k]=irow[l];
k                 333 src/matrixstorage.c         aa[k]=aub[l];
k                 334 src/matrixstorage.c         k++;l++;
k                 356 src/matrixstorage.c     k=0;
k                 358 src/matrixstorage.c       ai[k]=i+1;
k                 359 src/matrixstorage.c       aj[k]=i+1;
k                 360 src/matrixstorage.c       aa[k]=adb[i];
k                 361 src/matrixstorage.c       k++;
k                 366 src/matrixstorage.c         ai[k]=i+1;
k                 367 src/matrixstorage.c         aj[k]=irow[l];
k                 368 src/matrixstorage.c         aa[k]=aub[l];
k                 369 src/matrixstorage.c         k++;
k                 370 src/matrixstorage.c         ai[k]=irow[l];
k                 371 src/matrixstorage.c         aj[k]=i+1;
k                 372 src/matrixstorage.c         aa[k]=aub[l];
k                 373 src/matrixstorage.c         k++;
k                 380 src/matrixstorage.c     k=0;
k                 382 src/matrixstorage.c       jref=aj[k];
k                 383 src/matrixstorage.c       kstart=k;
k                 385 src/matrixstorage.c         k++;
k                 386 src/matrixstorage.c 	if(k==ndim2) break;
k                 387 src/matrixstorage.c         if(aj[k]!=jref) break;
k                 389 src/matrixstorage.c       klen=k-kstart;
k                 416 src/matrixstorage.c       k=(long long)i3c*neq3+i3l;
k                 417 src/matrixstorage.c       FORTRAN(nidentll,(ipoint,&k,&npoint,&id));
k                 420 src/matrixstorage.c         ipoint[npoint-1]=k;
k                 423 src/matrixstorage.c       else if(ipoint[id-1]!=k){
k                 431 src/matrixstorage.c         ipoint[npoint-1]=k;
k                 484 src/matrixstorage.c     k=-1;
k                 494 src/matrixstorage.c         k++;
k                 495 src/matrixstorage.c         ai[k]=il;
k                 496 src/matrixstorage.c         aj[k]=ic;
k                 497 src/matrixstorage.c         aa[k]=aa3[9*i+j];
k                 505 src/matrixstorage.c   k=0;
k                 507 src/matrixstorage.c     jref=aj[k];
k                 508 src/matrixstorage.c     kstart=k;
k                 510 src/matrixstorage.c       k++;
k                 511 src/matrixstorage.c       if(aj[k]!=jref) break;
k                 513 src/matrixstorage.c     klen=k-kstart;
k                  85 src/nonlingeo.c   ITG *inum=NULL,k,iout=0,icntrl,iinc=0,jprint=0,iit=-1,jnz=0,
k                 195 src/nonlingeo.c   for(k=0;k<3;k++){
k                 196 src/nonlingeo.c       strcpy1(&jobnamef[k*132],&jobnamec[k*132],132);
k                 281 src/nonlingeo.c   for(k=0;k<*nboun;++k){xbounact[k]=xbounold[k];}
k                 286 src/nonlingeo.c   for(k=0;k<7**nbody;k++){xbodyact[k]=xbody[k];}
k                 293 src/nonlingeo.c       for(k=1;k<=*nbody;k++){
k                 295 src/nonlingeo.c 			     iendset,ialset,&inewton,nset,&ifreebody,&k));
k                 451 src/nonlingeo.c 	      for(k=mi[2]**ne;k<mi[2]*(*ne+*nslavs);k++) ielorien[k]=0;
k                 455 src/nonlingeo.c 	  for(k=mi[2]**ne;k<mi[2]*(*ne+*nslavs);k++) ielmat[k]=1;
k                 459 src/nonlingeo.c 	      for(k=*nstate_*mi[0]**ne;k<*nstate_*mi[0]*(*ne+*nslavs);k++){
k                 460 src/nonlingeo.c 		  xstate[k]=0.;
k                 483 src/nonlingeo.c       for(k=0;k<3*memmpc_;k++){nodempcref[k]=nodempc[k];}
k                 485 src/nonlingeo.c       for(k=0;k<memmpc_;k++){coefmpcref[k]=coefmpc[k];}
k                 491 src/nonlingeo.c       for(k=0;k<*nk;++k){t1act[k]=t1old[k];}
k                 521 src/nonlingeo.c       for(k=0;k<*nstate_*mi[0]*(*ne+*nslavs);++k){
k                 522 src/nonlingeo.c 	  xstateini[k]=xstate[k];
k                 612 src/nonlingeo.c 	  for(k=0;k<3*memmpcref_;k++){nodempc[k]=nodempcref[k];}
k                 614 src/nonlingeo.c 	  for(k=0;k<memmpcref_;k++){coefmpc[k]=coefmpcref[k];}
k                 626 src/nonlingeo.c 	  for(k=0;k<3*memmpc_;k++){nodempcref[k]=nodempc[k];}
k                 628 src/nonlingeo.c 	  for(k=0;k<memmpc_;k++){coefmpcref[k]=coefmpc[k];}
k                 742 src/nonlingeo.c       for(k=0;k<neq[0];++k){
k                 743 src/nonlingeo.c 	  b[k]=fext[k]-f[k];
k                 754 src/nonlingeo.c 	  for(k=0;k<neq[0];++k){
k                 755 src/nonlingeo.c 	      ad[k]=adb[k]+scal1*ad[k];
k                 757 src/nonlingeo.c 	  for(k=0;k<nzs[0];++k){
k                 758 src/nonlingeo.c 	      au[k]=aub[k]+scal1*au[k];
k                 801 src/nonlingeo.c 	  for(k=0;k<neq[0];++k){
k                 802 src/nonlingeo.c 	      b[k]=(fext[k]-f[k])/adb[k];
k                 808 src/nonlingeo.c       for(k=neq[0];k<neq[1];++k){
k                 809 src/nonlingeo.c 	  b[k]=0.;
k                 844 src/nonlingeo.c 	  for(k=0;k<6*mi[0]*ne0;++k){
k                 845 src/nonlingeo.c 	      sti[k]=stx[k];
k                 883 src/nonlingeo.c 	  for(k=0;k<*nboun;++k){xbounini[k]=xbounact[k];}
k                 885 src/nonlingeo.c 	      for(k=0;k<*nk;++k){t1ini[k]=t1act[k];}
k                 887 src/nonlingeo.c 	  for(k=0;k<neq[1];++k){
k                 888 src/nonlingeo.c 	      fini[k]=f[k];
k                 891 src/nonlingeo.c 	      for(k=0;k<mt**nk;++k){
k                 892 src/nonlingeo.c 		  veini[k]=veold[k];
k                 893 src/nonlingeo.c 		  accini[k]=accold[k];
k                 895 src/nonlingeo.c 	      for(k=0;k<neq[1];++k){
k                 897 src/nonlingeo.c 		  fextini[k]=fext[k];
k                 898 src/nonlingeo.c 		  cvini[k]=cv[k];
k                 902 src/nonlingeo.c 	      for(k=0;k<6*mi[0]*ne0;++k){
k                 903 src/nonlingeo.c 		  stiini[k]=sti[k];
k                 904 src/nonlingeo.c 		  emeini[k]=eme[k];
k                 908 src/nonlingeo.c 	      for(k=0;k<mi[0]*ne0;++k){enerini[k]=ener[k];}
k                 912 src/nonlingeo.c 		  for(k=0;k<*nstate_*mi[0]*(ne0+*nslavs);++k){
k                 913 src/nonlingeo.c 		      xstateini[k]=xstate[k];
k                 997 src/nonlingeo.c 	  for(k=0;k<3*memmpcref_;k++){nodempc[k]=nodempcref[k];}
k                 999 src/nonlingeo.c 	  for(k=0;k<memmpcref_;k++){coefmpc[k]=coefmpcref[k];}
k                1052 src/nonlingeo.c 		      for(k=mi[2]**ne;k<mi[2]*(*ne+*nintpoint);k++) ielorien[k]=0;
k                1055 src/nonlingeo.c 		  for(k=mi[2]**ne;k<mi[2]*(*ne+*nintpoint);k++) ielmat[k]=1;
k                1103 src/nonlingeo.c 			  for(k=*nstate_*mi[0]*ne0;
k                1104 src/nonlingeo.c                                 k<*nstate_*mi[0]*(ne0+maxprevcontel);++k){
k                1105 src/nonlingeo.c 			      xstateini[k]=xstate[k];
k                1110 src/nonlingeo.c 		      for(k=*nstate_*mi[0]*ne0;k<*nstate_*mi[0]*(ne0+*nintpoint);k++){
k                1111 src/nonlingeo.c 			  xstate[k]=0.;
k                1133 src/nonlingeo.c 		      for(k=0;k<*nstate_*mi[0]*(ne0+*nintpoint);++k){
k                1134 src/nonlingeo.c 			  xstateini[k]=xstate[k];
k                1155 src/nonlingeo.c 	  for(k=0;k<3*memmpc_;k++){nodempcref[k]=nodempc[k];}
k                1157 src/nonlingeo.c 	  for(k=0;k<memmpc_;k++){coefmpcref[k]=coefmpc[k];}
k                1191 src/nonlingeo.c 	  for(k=0;k<*nboun;++k){
k                1192 src/nonlingeo.c 	      err=fabs(xbounact[k]-xbounini[k]);
k                1249 src/nonlingeo.c 	  for(k=0;k<neq[1];++k){b[k]=f[k];}
k                1273 src/nonlingeo.c 	  for(k=0;k<neq[1];++k){f[k]=f[k]+b[k];}
k                1301 src/nonlingeo.c 	      for(k=0;k<6*mi[0]*ne0;++k){
k                1302 src/nonlingeo.c 		  sti[k]=stx[k];
k                1323 src/nonlingeo.c 	for(k=0;k<nzs[1]-nzs[0];k++) iruc[k]=irow[k+nzs[0]]-neq[0];
k                1365 src/nonlingeo.c 	      for(k=0;k<3*memmpcref_;k++){nodempc[k]=nodempcref[k];}
k                1367 src/nonlingeo.c 	      for(k=0;k<memmpcref_;k++){coefmpc[k]=coefmpcref[k];}
k                1420 src/nonlingeo.c 	      for(k=0;k<*nk;++k){t1act[k]=vold[mt*k];}
k                1433 src/nonlingeo.c 	      for(k=0;k<3*memmpc_;k++){nodempcref[k]=nodempc[k];}
k                1435 src/nonlingeo.c 	      for(k=0;k<memmpc_;k++){coefmpcref[k]=coefmpc[k];}
k                1584 src/nonlingeo.c 	  for(k=0;k<neq[0];k++){adc[k]=alpham*adb[k]+betam*ad[k];}
k                1587 src/nonlingeo.c 	      for(k=0;k<nzs[0];k++){auc[k]=alpham*aub[k]+betam*au[k];}
k                1590 src/nonlingeo.c 	      for(k=0;k<2*nzs[0];k++){auc[k]=alpham*aub[k]+betam*au[k];}
k                1626 src/nonlingeo.c 	  NNEW(inum,ITG,*nk);for(k=0;k<*nk;k++) inum[k]=1;
k                1658 src/nonlingeo.c 		  for(k=0;k<neq[0];++k){
k                1659 src/nonlingeo.c 		      ad[k]=adb[k]+scal1*ad[k];
k                1661 src/nonlingeo.c 		  for(k=0;k<nzs[0];++k){
k                1662 src/nonlingeo.c 		      au[k]=aub[k]+scal1*au[k];
k                1668 src/nonlingeo.c 		      for(k=nzs[2];k<nzs[2]+nzs[0];++k){
k                1669 src/nonlingeo.c 			  au[k]=aub[k]+scal1*au[k];
k                1677 src/nonlingeo.c 		      for(k=0;k<neq[0];++k){
k                1678 src/nonlingeo.c 			  ad[k]+=scal1*adc[k];
k                1680 src/nonlingeo.c 		      for(k=0;k<nzs[0];++k){
k                1681 src/nonlingeo.c 			  au[k]+=scal1*auc[k];
k                1687 src/nonlingeo.c 			  for(k=nzs[2];k<nzs[2]+nzs[0];++k){
k                1688 src/nonlingeo.c 			      au[k]+=scal1*auc[k];
k                1698 src/nonlingeo.c 		  for(k=neq[0];k<neq[1];++k){
k                1699 src/nonlingeo.c 		      ad[k]=adb[k]/dtime+ad[k];
k                1701 src/nonlingeo.c 		  for(k=nzs[0];k<nzs[1];++k){
k                1702 src/nonlingeo.c 		      au[k]=aub[k]/dtime+au[k];
k                1708 src/nonlingeo.c 		      for(k=nzs[2]+nzs[0];k<nzs[2]+nzs[1];++k){
k                1709 src/nonlingeo.c 			  au[k]=aub[k]/dtime+au[k];
k                1805 src/nonlingeo.c 	      for(k=0;k<neq[0];++k){
k                1806 src/nonlingeo.c 		  b[k]=b[k]/adb[k];
k                1810 src/nonlingeo.c 	      for(k=neq[0];k<neq[1];++k){
k                1811 src/nonlingeo.c 		  b[k]=b[k]*dtime/adb[k];
k                1861 src/nonlingeo.c 	  for(k=0;k<6*mi[0]*ne0;++k){
k                1862 src/nonlingeo.c 	      sti[k]=stx[k];
k                1887 src/nonlingeo.c       for(k=0;k<2;++k){
k                1888 src/nonlingeo.c 	  ram2[k]=ram1[k];
k                1889 src/nonlingeo.c 	  ram1[k]=ram[k];
k                1890 src/nonlingeo.c 	  ram[k]=0.;
k                1893 src/nonlingeo.c 	  for(k=0;k<neq[0];++k){
k                1894 src/nonlingeo.c 	      err=fabs(b[k]);
k                1895 src/nonlingeo.c 	      if(err>ram[0]){ram[0]=err;ram[2]=k+0.5;}
k                1899 src/nonlingeo.c 	  for(k=neq[0];k<neq[1];++k){
k                1900 src/nonlingeo.c 	      err=fabs(b[k]);
k                1901 src/nonlingeo.c 	      if(err>ram[1]){ram[1]=err;ram[3]=k+0.5;}
k                1908 src/nonlingeo.c 	  for(k=4;k<8;++k){
k                1909 src/nonlingeo.c 	      ram2[k]=ram1[k];
k                1910 src/nonlingeo.c 	      ram1[k]=ram[k];
k                2021 src/nonlingeo.c       for(k=0;k<*nboun;++k){xbounact[k]=xbounini[k];}
k                2023 src/nonlingeo.c 	for(k=0;k<*nk;++k){t1act[k]=t1ini[k];}
k                2025 src/nonlingeo.c       for(k=0;k<neq[1];++k){
k                2026 src/nonlingeo.c 	  f[k]=fini[k];
k                2029 src/nonlingeo.c 	for(k=0;k<mt**nk;++k){
k                2030 src/nonlingeo.c 	  veold[k]=veini[k];
k                2031 src/nonlingeo.c 	  accold[k]=accini[k];
k                2033 src/nonlingeo.c 	for(k=0;k<neq[1];++k){
k                2035 src/nonlingeo.c 	  fext[k]=fextini[k];
k                2036 src/nonlingeo.c 	  cv[k]=cvini[k];
k                2040 src/nonlingeo.c 	  for(k=0;k<6*mi[0]*ne0;++k){
k                2041 src/nonlingeo.c 	      sti[k]=stiini[k];
k                2042 src/nonlingeo.c 	      eme[k]=emeini[k];
k                2046 src/nonlingeo.c 	  for(k=0;k<mi[0]*ne0;++k){ener[k]=enerini[k];}
k                2048 src/nonlingeo.c       for(k=0;k<*nstate_*mi[0]*(ne0+maxprevcontel);++k){
k                2049 src/nonlingeo.c 	  xstate[k]=xstateini[k];
k                2305 src/nonlingeo.c     for(k=0;k<mt**nk;++k){veold[k]=0.;}
k                2312 src/nonlingeo.c   for(k=0;k<*nboun;++k){
k                2317 src/nonlingeo.c       if(ndirboun[k]==0){
k                2324 src/nonlingeo.c       }else if((ndirboun[k]>0)&&(ndirboun[k]<4)){
k                2325 src/nonlingeo.c 	  node=nodeboun[k];
k                2333 src/nonlingeo.c       xbounold[k]=xbounact[k];
k                2335 src/nonlingeo.c   for(k=0;k<*nforc;++k){xforcold[k]=xforcact[k];}
k                2336 src/nonlingeo.c   for(k=0;k<2**nload;++k){xloadold[k]=xloadact[k];}
k                2337 src/nonlingeo.c   for(k=0;k<7**nbody;k=k+7){xbodyold[k]=xbodyact[k];}
k                2339 src/nonlingeo.c     for(k=0;k<*nk;++k){t1old[k]=t1act[k];}
k                2340 src/nonlingeo.c     for(k=0;k<*nk;++k){vold[mt*k]=t1act[k];}
k                2343 src/nonlingeo.c     for(k=0;k<*nk;++k){t1[k]=vold[mt*k];}
k                2345 src/nonlingeo.c 	for(k=0;k<*nk;++k){t1old[k]=t1act[k];}
k                2389 src/nonlingeo.c       for(k=0;k<3*memmpcref_;k++){nodempc[k]=nodempcref[k];}
k                2391 src/nonlingeo.c       for(k=0;k<memmpcref_;k++){coefmpc[k]=coefmpcref[k];}
k                  40 src/pardiso.c    ITG i,j,k,l,maxfct=1,mnum=1,phase=12,nrhs=1,*perm=NULL,mtype,
k                  93 src/pardiso.c        k=ndim;
k                 100 src/pardiso.c  		  icolpardiso[--k]=irow[--l];
k                 101 src/pardiso.c  		  aupardiso[k]=au[l];
k                 103 src/pardiso.c  	      pointers[i]=k--;
k                 104 src/pardiso.c  	      icolpardiso[k]=i+1;
k                 105 src/pardiso.c  	      aupardiso[k]=ad[i];
k                 112 src/pardiso.c  		  icolpardiso[--k]=irow[--l];
k                 113 src/pardiso.c  		  aupardiso[k]=au[l]-*sigma*aub[l];
k                 115 src/pardiso.c  	      pointers[i]=k--;
k                 116 src/pardiso.c  	      icolpardiso[k]=i+1;
k                 117 src/pardiso.c  	      aupardiso[k]=ad[i]-*sigma*adb[i];
k                 135 src/pardiso.c  	  k=0;
k                 139 src/pardiso.c  		  if(au[k]>1.e-12||au[k]<-1.e-12){
k                 141 src/pardiso.c  		      irowpardiso[k2]=irow[k];
k                 142 src/pardiso.c  		      aupardiso[k2]=au[k];
k                 145 src/pardiso.c  		  k++;	      
k                 166 src/pardiso.c  	  k=0;
k                 170 src/pardiso.c  	      kstart=k;
k                 172 src/pardiso.c  		  if(irowpardiso[k]!=j ){
k                 173 src/pardiso.c  		      n=k-kstart;		  
k                 176 src/pardiso.c  		      pointers[i+1]=k+1;
k                 179 src/pardiso.c  		      if(k+1==ndim){
k                 180 src/pardiso.c  			  n=k-kstart+1;	  
k                 185 src/pardiso.c  			  k++;	       
k                 207 src/pardiso.c  	  k=0;
k                 210 src/pardiso.c  		  icolpardiso[k]=i+1;
k                 211 src/pardiso.c  		  irowpardiso[k]=irow[k];
k                 212 src/pardiso.c  		  aupardiso[k]=au[k];
k                 213 src/pardiso.c  		  k++;
k                 226 src/pardiso.c  	  k=0;
k                 230 src/pardiso.c  	      kstart=k;
k                 235 src/pardiso.c  		  if(irowpardiso[k]!=j){
k                 236 src/pardiso.c  		      n=k-kstart;
k                 239 src/pardiso.c  		      pointers[i+1]=k+1;
k                 245 src/pardiso.c  		      if(k+1==ndim){
k                 246 src/pardiso.c  			  n=k-kstart+1;
k                 254 src/pardiso.c  			  k++;
k                 267 src/pardiso.c  	  k=ndim;
k                 269 src/pardiso.c  	      l=k+1;
k                 271 src/pardiso.c  		  icolpardiso[--k]=irow[j-1];
k                 272 src/pardiso.c  		  aupardiso[k]=au[j+*nzs3-1];
k                 274 src/pardiso.c  	      icolpardiso[--k]=i+1;
k                 275 src/pardiso.c  	      aupardiso[k]=ad[i];
k                 277 src/pardiso.c  		  icolpardiso[--k]=icolpardiso[j-1];
k                 278 src/pardiso.c  		  aupardiso[k]=aupardiso[j-1];
k                  39 src/pardiso_as.c   ITG i,j,k,maxfct=1,mnum=1,mtype=11,phase=12,nrhs=1,*perm=NULL,
k                  80 src/pardiso_as.c   k=0;
k                  91 src/pardiso_as.c       pointersas[i]=k+1;
k                 102 src/pardiso_as.c 	    irowpardisoas[k]=j+1;
k                 103 src/pardiso_as.c 	    aupardisoas[k]=au[index];
k                 104 src/pardiso_as.c 	    k++;
k                 111 src/pardiso_as.c       irowpardisoas[k]=i+1;
k                 112 src/pardiso_as.c       aupardisoas[k]=ad[i];
k                 113 src/pardiso_as.c       k++;
k                 118 src/pardiso_as.c 	irowpardisoas[k]=irow[j-1];
k                 119 src/pardiso_as.c 	aupardisoas[k]=au[j+*nzs-1];
k                 120 src/pardiso_as.c 	k++;
k                 123 src/pardiso_as.c     pointersas[*neq]=k+1;
k                 126 src/pardiso_as.c       pointersas[i]=k+1;
k                 137 src/pardiso_as.c 	    irowpardisoas[k]=j+1;
k                 138 src/pardiso_as.c 	    aupardisoas[k]=au[index]-*sigma*aub[index];
k                 139 src/pardiso_as.c 	    k++;
k                 146 src/pardiso_as.c       irowpardisoas[k]=i+1;
k                 147 src/pardiso_as.c       aupardisoas[k]=ad[i]-*sigma*adb[i];
k                 148 src/pardiso_as.c       k++;
k                 153 src/pardiso_as.c 	irowpardisoas[k]=irow[j-1];
k                 154 src/pardiso_as.c 	aupardisoas[k]=au[j+*nzs-1]-*sigma*aub[j+*nzs-1];
k                 155 src/pardiso_as.c 	k++;
k                 158 src/pardiso_as.c     pointersas[*neq]=k+1;
k                 151 src/pcgsolver.c   ITG			i=0, k=0, ncg=0,iam,ier=0;
k                 190 src/pcgsolver.c   for (k=1; k<=*niter; k++, ncg++)
k                 208 src/pcgsolver.c       if (k!=1 && (ram<=c1*qam))	break;
k                 209 src/pcgsolver.c       if (k!=1)
k                 240 src/pcgsolver.c   if(k==*niter){
k                 326 src/pcgsolver.c   ITG				i=0, j=0, jlo=0, jup=0, k=0;
k                 344 src/pcgsolver.c 	  k = ia[j];
k                 345 src/pcgsolver.c 	  z[i] += A[j]*p[k];
k                 346 src/pcgsolver.c 	  z[k] += A[j]*p[i];
k                 377 src/pcgsolver.c 	ITG				i=0, j=0, jlo=0, jup=0, k=0, klo=0, kup=0, l=0, llo=0, lup=0;
k                 402 src/pcgsolver.c 			for (k=klo; k<=kup; k++)
k                 404 src/pcgsolver.c 				m = ia[k];
k                 411 src/pcgsolver.c 					C[k] -= C[j]*C[l];
k                 507 src/pcgsolver.c 	ITG			i=0, k=0, ncg=0,iam;
k                 530 src/pcgsolver.c 	for (k=1; k<=(*niter); k++, ncg++)
k                 539 src/pcgsolver.c 		if (k!=1 && (ram<=c1*qam)) break;
k                 541 src/pcgsolver.c 		if (k!=1)
k                 565 src/pcgsolver.c 	if(k==*niter){
k                  38 src/prediction.c     ITG j,k,mt=mi[1]+1;
k                  49 src/prediction.c 	    for(k=0;k<*nk;++k){
k                  51 src/prediction.c 		    dextrapol=*dtime*veold[mt*k+j]+scal1*accold[mt*k+j];
k                  52 src/prediction.c 		    if((fabs(dextrapol)>uam[0])&&(nactdof[mt*k+j]>0)) {uam[0]=fabs(dextrapol);}
k                  53 src/prediction.c 		    v[mt*k+j]=vold[mt*k+j]+dextrapol;
k                  54 src/prediction.c 		    veold[mt*k+j]=veold[mt*k+j]+scal2*accold[mt*k+j];
k                  55 src/prediction.c 		    accold[mt*k+j]=0.;
k                  59 src/prediction.c 	    for(k=0;k<*nk;++k){
k                  61 src/prediction.c 		    v[mt*k+j]=vold[mt*k+j];
k                  64 src/prediction.c 	    for(k=0;k<*nk;++k){
k                  65 src/prediction.c 		dextrapol=*dtime*veold[mt*k];
k                  67 src/prediction.c 		if((fabs(dextrapol)>uam[1])&&(nactdof[mt*k]>0)) {uam[1]=fabs(dextrapol);}
k                  68 src/prediction.c 		v[mt*k]+=dextrapol;
k                  71 src/prediction.c 	    for(k=0;k<*nk;++k){
k                  73 src/prediction.c 		    dextrapol=*dtime*veold[mt*k+j]+scal1*accold[mt*k+j];
k                  76 src/prediction.c 			if((fabs(dextrapol)>uam[1])&&(nactdof[mt*k]>0)) {uam[1]=fabs(dextrapol);}
k                  78 src/prediction.c 			if((fabs(dextrapol)>uam[0])&&(nactdof[mt*k+j]>0)) {uam[0]=fabs(dextrapol);}
k                  80 src/prediction.c 		    v[mt*k+j]=vold[mt*k+j]+dextrapol;
k                  81 src/prediction.c 		    veold[mt*k+j]=veold[mt*k+j]+scal2*accold[mt*k+j];
k                  82 src/prediction.c 		    accold[mt*k+j]=0.;
k                  94 src/prediction.c 		for(k=0;k<*nk;++k){
k                  97 src/prediction.c 			    dextrapol=*dtime*veold[mt*k+j];
k                  98 src/prediction.c 			    if((fabs(dextrapol)>uam[0])&&(nactdof[mt*k+j]>0)) {uam[0]=fabs(dextrapol);}	
k                  99 src/prediction.c 			    v[mt*k+j]=vold[mt*k+j]+dextrapol;
k                 101 src/prediction.c 			    v[mt*k+j]=vold[mt*k+j];
k                 106 src/prediction.c 		for(k=0;k<*nk;++k){
k                 108 src/prediction.c 			v[mt*k+j]=vold[mt*k+j];
k                 111 src/prediction.c 		for(k=0;k<*nk;++k){
k                 113 src/prediction.c 			dextrapol=*dtime*veold[mt*k];
k                 115 src/prediction.c 			if((fabs(dextrapol)>uam[1])&&(nactdof[mt*k]>0)) {uam[1]=fabs(dextrapol);}	
k                 116 src/prediction.c 			v[mt*k]+=dextrapol;
k                 120 src/prediction.c 		for(k=0;k<*nk;++k){
k                 123 src/prediction.c 			    dextrapol=*dtime*veold[mt*k+j];
k                 126 src/prediction.c 				if((fabs(dextrapol)>uam[1])&&(nactdof[mt*k+j]>0)) {uam[1]=fabs(dextrapol);}
k                 128 src/prediction.c 				if((fabs(dextrapol)>uam[0])&&(nactdof[mt*k+j]>0)) {uam[0]=fabs(dextrapol);}
k                 130 src/prediction.c 			    v[mt*k+j]=vold[mt*k+j]+dextrapol;
k                 132 src/prediction.c 			    v[mt*k+j]=vold[mt*k+j];
k                 139 src/prediction.c 	    for(k=0;k<*nk;++k){
k                 141 src/prediction.c 		    v[mt*k+j]=vold[mt*k+j];
k                  38 src/prediction_em.c     ITG j,k,mt=mi[1]+1,jstart;
k                  51 src/prediction_em.c 	for(k=0;k<*nk;++k){
k                  53 src/prediction_em.c 		dextrapol=*dtime*veold[mt*k+j];
k                  56 src/prediction_em.c 		    if((fabs(dextrapol)>uam[1])&&(nactdof[mt*k]>0)) {uam[1]=fabs(dextrapol);}
k                  58 src/prediction_em.c 		v[mt*k+j]=vold[mt*k+j]+dextrapol;
k                  68 src/prediction_em.c 	    for(k=0;k<*nk;++k){
k                  71 src/prediction_em.c 			dextrapol=*dtime*veold[mt*k+j];
k                  74 src/prediction_em.c 			    if((fabs(dextrapol)>uam[1])&&(nactdof[mt*k]>0)) {uam[1]=fabs(dextrapol);}
k                  76 src/prediction_em.c 			v[mt*k+j]=vold[mt*k+j]+dextrapol;
k                  78 src/prediction_em.c 			v[mt*k+j]=vold[mt*k+j];
k                  84 src/prediction_em.c 	    for(k=0;k<*nk;++k){
k                  86 src/prediction_em.c 		    v[mt*k+j]=vold[mt*k+j];
k                  26 src/preiter.c    ITG precFlg,niter=5000000,ndim,i,j,k,ier,*icol=NULL,*irow=NULL,
k                  62 src/preiter.c    k=*nzs;
k                  65 src/preiter.c        icol[--k]=i+1;
k                  69 src/preiter.c    k=*nzs;
k                  72 src/preiter.c      au[k]=ad[i];
k                  73 src/preiter.c      irow[k]=++j;
k                  74 src/preiter.c      icol[k]=j;
k                  75 src/preiter.c      ++k;
k                  35 src/radcyc.c        node,i1,i2,nope,iel,indexe,j,k,ielset,node1,node2,node3,l,jj;
k                 120 src/radcyc.c       for(k=1;k<idtie;k++){
k                 121 src/radcyc.c         j=i+k**ntri;
k                 122 src/radcyc.c         kontri[4*j]=node1+k**nk;
k                 123 src/radcyc.c         kontri[4*j+1]=node2+k**nk;
k                 124 src/radcyc.c         kontri[4*j+2]=node3+k**nk;
k                  36 src/readinput.c   ITG i,j,k,n,in=0,nlinemax=100000,irestartread,irestartstep,
k                  90 src/readinput.c       k=0;
k                  95 src/readinput.c 	      if((buff[i]=='\0')||(buff[i]=='\n')||(buff[i]=='\r')||(k==1320)) break;
k                  97 src/readinput.c 	      buff[k]=buff[i];
k                  98 src/readinput.c 	      k++;
k                 103 src/readinput.c 	      if((buff[i]=='\0')||(buff[i]=='\n')||(buff[i]=='\r')||(k==1320)) break;
k                 104 src/readinput.c 	      buff[k]=buff[i];
k                 105 src/readinput.c 	      k++;
k                 111 src/readinput.c       if(k==0) continue;
k                 137 src/readinput.c 		  }while(j<k);
k                 142 src/readinput.c 	  }while(j<k);
k                 154 src/readinput.c 	  lincludefn=k;
k                 179 src/readinput.c       if(ipoinpc[*nline-1]+k>ncharmax){
k                 186 src/readinput.c       for(j=0;j<k;j++){
k                 189 src/readinput.c       ipoinpc[*nline]=ipoinpc[*nline-1]+k;
k                 196 src/readinput.c 	  strcpy1(&buff[k]," ",1);
k                  34 src/sgi.c        ITG i,j,k,l,*pointers=NULL,method;
k                  49 src/sgi.c        k=ndim;
k                  56 src/sgi.c      	irowsgi[--k]=irow[--l]-1;
k                  57 src/sgi.c      	ausgi[k]=au[l];
k                  59 src/sgi.c            pointers[i]=--k;
k                  60 src/sgi.c            irowsgi[k]=i;
k                  61 src/sgi.c            ausgi[k]=ad[i];
k                  68 src/sgi.c      	irowsgi[--k]=irow[--l]-1;
k                  69 src/sgi.c      	ausgi[k]=au[l]-*sigma*aub[l];
k                  71 src/sgi.c            pointers[i]=--k;
k                  72 src/sgi.c            irowsgi[k]=i;
k                  73 src/sgi.c            ausgi[k]=ad[i]-*sigma*adb[i];
k                  75 src/steadystate.c   ITG nev,i,j,k, *inum=NULL,*ipobody=NULL,inewton=0,nsectors,im,
k                 861 src/steadystate.c       for(k=0;k<7**nbody;k++){xbodyact[k]=xbody[k];}
k                 872 src/steadystate.c 	  for(k=1;k<=*nbody;k++){
k                 874 src/steadystate.c 				 iendset,ialset,&inewton,nset,&ifreebody,&k));
k                 995 src/steadystate.c 		      for(k=0;k<neq[1];k++){
k                 996 src/steadystate.c 			  cc[i*nev+j]+=z[(long long)j*neq[1]+k]*zc[i*neq[1]+k];
k                1895 src/steadystate.c       for(k=0;k<*nboun;++k){xboun[k]=xbounact[k];}
k                1896 src/steadystate.c       for(k=0;k<*nforc;++k){xforc[k]=xforcact[k];}
k                1897 src/steadystate.c       for(k=0;k<2**nload;++k){xload[k]=xloadact[k];}
k                1898 src/steadystate.c       for(k=0;k<7**nbody;k=k+7){xbody[k]=xbodyact[k];}
k                1900 src/steadystate.c 	  for(k=0;k<*nk;++k){t1[k]=t1act[k];}
k                2196 src/steadystate.c 	  for(k=0;k<7**nbody;k++){xbodyr[k]=xbody[k];}
k                2209 src/steadystate.c 	      for(k=1;k<=*nbody;k++){
k                2211 src/steadystate.c 			    iendset,ialset,&inewton,nset,&ifreebody,&k));
k                2271 src/steadystate.c 			  for(k=0;k<neq[1];k++){
k                2272 src/steadystate.c 			      cc[i*nev+j]+=z[(long long)j*neq[1]+k]*zc[i*neq[1]+k];
k                  41 src/storecontactdof.c   ITG j,j1,jdof,id,k,l,ist,index,node,ndir,*ikactcont=*ikactcontp;
k                  61 src/storecontactdof.c 		  k=*nactcont-1;
k                  62 src/storecontactdof.c 		  l=k-1;
k                  63 src/storecontactdof.c 		  while(k>id){
k                  64 src/storecontactdof.c 		      ikactcont[k--]=ikactcont[l--];
k                  99 src/storecontactdof.c 				  k=*nactcont-1;
k                 100 src/storecontactdof.c 				  l=k-1;
k                 101 src/storecontactdof.c 				  while(k>id){
k                 102 src/storecontactdof.c 				      ikactcont[k--]=ikactcont[l--];
k                  45 src/tau.c        ITG i,j,k,l,*irow=NULL;
k                  65 src/tau.c        k=ndim;
k                  72 src/tau.c      	irowtau[--k]=irow[--l]-1;
k                  73 src/tau.c      	autau[k]=au[l];
k                  75 src/tau.c            pointtau[i]=--k;
k                  76 src/tau.c            irowtau[k]=i;
k                  77 src/tau.c            autau[k]=ad[i];
k                  84 src/tau.c      	irowtau[--k]=irow[--l]-1;
k                  85 src/tau.c      	autau[k]=au[l]-*sigma*aub[l];
k                  87 src/tau.c            pointtau[i]=--k;
k                  88 src/tau.c            irowtau[k]=i;
k                  89 src/tau.c            autau[k]=ad[i]-*sigma*adb[i];
k                  36 src/tiedcontact.c   ITG *itietri=NULL,*koncont=NULL,nconf,i,k,*nx=NULL,im,
k                 149 src/tiedcontact.c   for(k=index;k<*memmpc_;k++){
k                 150 src/tiedcontact.c       nodempc[3*k-1]=k+1;