i                  55 src/CalculiX.h void FORTRAN(addimdnodecload,(ITG *nodeforc,ITG *i,ITG *imdnode, 
i                  62 src/CalculiX.h              ITG *kon,char *lakon,ITG *i,ITG *imdnode,ITG *nmdnode,
i                  69 src/CalculiX.h 	     ITG *mi,ITG *izdof,ITG *nzdof,ITG *i,ITG *iznode,ITG *nznode,
i                  74 src/CalculiX.h 	     ITG *mi,ITG *i,ITG *iznode,ITG *nznode,ITG *nk,
i                 238 src/CalculiX.h void *biotsavartmt(ITG *i);
i                 344 src/CalculiX.h void *calcviewmt(ITG *i);
i                 601 src/CalculiX.h void *ddotc1mt(ITG *i);
i                1487 src/CalculiX.h void FORTRAN(mpcrem,(ITG *i,ITG *mpcfree,ITG *nodempc,ITG *nmpc,ITG *ikmpc,
i                1868 src/CalculiX.h void *resultsemmt(ITG *i);
i                1916 src/CalculiX.h void *resultsmechmt(ITG *i);
i                1952 src/CalculiX.h void *resultsthermemmt(ITG *i);
i                1954 src/CalculiX.h void *resultsthermmt(ITG *i);
i                2054 src/CalculiX.h         ITG *i,ITG *l,ITG *ntri));
i                  84 src/arpackcs.c     info,rvec=1,*select=NULL,lfin,j,lint,iout=1,nm,index,inode,id,i,idof,
i                 427 src/arpackcs.c       for(i=0;i<*nk;i++) inocs[i]=-1;
i                 428 src/arpackcs.c       for(i=0;i<*ne;i++) ielcs[i]=-1;
i                 431 src/arpackcs.c   for(i=0;i<*mcs;i++){
i                 432 src/arpackcs.c       is=cs[17*i+4];
i                 434 src/arpackcs.c       ielset=cs[17*i+12];
i                 440 src/arpackcs.c 	      ielcs[iel]=i;
i                 455 src/arpackcs.c 		  inocs[node]=i;
i                 464 src/arpackcs.c 		  ielcs[iel]=i;
i                 474 src/arpackcs.c 		      inocs[node]=i;
i                1104 src/arpackcs.c 	  for(i=1;i<is;i++){
i                1106 src/arpackcs.c 	      theta=i*2.*pi/cs[17*jj];
i                1110 src/arpackcs.c 		      cot[3*l+i*3**nk]=cot[3*l];
i                1111 src/arpackcs.c 		      cot[1+3*l+i*3**nk]=cot[1+3*l]+theta;
i                1112 src/arpackcs.c 		      cot[2+3*l+i*3**nk]=cot[2+3*l];
i                1113 src/arpackcs.c 		      if(*ntrans>0){inotrt[2*l+i*2**nk]=inotrt[2*l];}
i                1116 src/arpackcs.c 	      for(l=0;l<*nkon;l++){kont[l+i**nkon]=kon[l]+i**nk;}
i                1120 src/arpackcs.c 			  ipkont[l+i**ne]=ipkon[l]+i**nkon;
i                1121 src/arpackcs.c 			  ielmatt[mi[2]*(l+i**ne)]=ielmat[mi[2]*l];
i                1124 src/arpackcs.c 			      lakont[l2+i*8**ne]=lakon[l2];
i                1246 src/arpackcs.c 	      for(i=0;i<6*mi[0]*ne0;i++){eme[i]=0.;}
i                1255 src/arpackcs.c 	      for(i=0;i<*nmpc;i++){
i                1256 src/arpackcs.c 		  index=ipompc[i]-1;
i                1267 src/arpackcs.c 		      if(strcmp1(&labmpc[20*i],"CYCLIC")==0){
i                1268 src/arpackcs.c 			  icomplex=atoi(&labmpc[20*i+6]);}
i                1269 src/arpackcs.c 		      else if(strcmp1(&labmpc[20*i],"SUBCYCLIC")==0){
i                1349 src/arpackcs.c 	      for(i=0;i<*nset;i++){
i                1350 src/arpackcs.c 		  if(strcmp1(&set[81*i],"RAYN")==0){
i                1351 src/arpackcs.c 		      iray=ialset[istartset[i]-1];
i                1440 src/arpackcs.c 	      for(i=0;i<*nset;i++){
i                1441 src/arpackcs.c 		  if(strcmp1(&set[81*i],"STRESSDOMAINN")==0){
i                1442 src/arpackcs.c 		      ielset=i+1;
i                1593 src/arpackcs.c       for(i=0;i<*nset;i++){
i                1594 src/arpackcs.c 	if(strcmp1(&set[81*i],"STRAINDOMAINN")==0){
i                1595 src/arpackcs.c 	  ielset=i+1;
i                1774 src/arpackcs.c       for(i=1;i<is;i++){
i                1776 src/arpackcs.c         for(l=0;l<*nk;l++){inumt[l+i**nk]=inum[l];}
i                1778 src/arpackcs.c         theta=i*nm*2.*pi/cs[17*jj];
i                1794 src/arpackcs.c 			  vt[l+mt**nk*i]=v[l];
i                1801 src/arpackcs.c                 vt[l+mt**nk*i]=ctl*v[l]-stl*v[l+mt**nk];
i                1822 src/arpackcs.c 			  vt[l+mt**nk*(i+ngraph)]=v[l+mt**nk];
i                1829 src/arpackcs.c                 vt[l+mt**nk*(i+ngraph)]=stl*v[l]+ctl*v[l+mt**nk];
i                1837 src/arpackcs.c 	      if(inocs[l]==jj) t1t[l+*nk*i]=t1[l];
i                1853 src/arpackcs.c 			  stnt[l+6**nk*i]=stn[l];
i                1860 src/arpackcs.c                 stnt[l+6**nk*i]=ctl*stn[l]-stl*stn[l+6**nk];
i                1881 src/arpackcs.c 			  stnt[l+6**nk*(i+ngraph)]=stn[l+6**nk];
i                1888 src/arpackcs.c                 stnt[l+6**nk*(i+ngraph)]=stl*stn[l]+ctl*stn[l+6**nk];
i                1906 src/arpackcs.c 			  eent[l+6**nk*i]=een[l];
i                1913 src/arpackcs.c                 eent[l+6**nk*i]=ctl*een[l]-stl*een[l+6**nk];
i                1934 src/arpackcs.c 			  eent[l+6**nk*(i+ngraph)]=een[l+6**nk];
i                1941 src/arpackcs.c                 eent[l+6**nk*(i+ngraph)]=stl*een[l]+ctl*een[l+6**nk];
i                1961 src/arpackcs.c 			  emnt[l+6**nk*i]=emn[l];
i                1968 src/arpackcs.c                 emnt[l+6**nk*i]=ctl*emn[l]-stl*emn[l+6**nk];
i                1989 src/arpackcs.c 			  emnt[l+6**nk*(i+ngraph)]=emn[l+6**nk];
i                1996 src/arpackcs.c                 emnt[l+6**nk*(i+ngraph)]=stl*emn[l]+ctl*emn[l+6**nk];
i                2017 src/arpackcs.c 			  cdnt[l+6**nk*i]=cdn[l];
i                2024 src/arpackcs.c                 cdnt[l+6**nk*i]=ctl*cdn[l]-stl*cdn[l+6**nk];
i                2045 src/arpackcs.c 			  cdnt[l+6**nk*(i+ngraph)]=cdn[l+6**nk];
i                2052 src/arpackcs.c                 cdnt[l+6**nk*(i+ngraph)]=stl*cdn[l]+ctl*cdn[l+6**nk];
i                2070 src/arpackcs.c 			  fnt[l+mt**nk*i]=fn[l];
i                2077 src/arpackcs.c                 fnt[l+mt**nk*i]=ctl*fn[l]-stl*fn[l+mt**nk];
i                2098 src/arpackcs.c 		    fnt[l+mt**nk*(i+ngraph)]=fn[l+mt**nk];
i                2105 src/arpackcs.c                 fnt[l+mt**nk*(i+ngraph)]=stl*fn[l]+ctl*fn[l+mt**nk];
i                2113 src/arpackcs.c             if(inocs[l]==jj) enernt[l+*nk*i]=0.;
i                2325 src/arpackcs.c   for(i=0;i<6*mi[0]*ne0;i++){eme[i]=0.;}
i                  34 src/biosav.c       ITG i,j,*ithread=NULL,nkphi,idelta,isum;
i                  96 src/biosav.c       for(i=0;i<*nk;i++){
i                  97 src/biosav.c   	if(inomat[i]==1) nkphi++;
i                 106 src/biosav.c       for(i=0;i<num_cpus;i++){
i                 107 src/biosav.c   	nkapar[i]=isum;
i                 108 src/biosav.c   	if(i!=num_cpus-1){
i                 113 src/biosav.c   	nkepar[i]=isum-1;
i                 118 src/biosav.c       i=0;
i                 126 src/biosav.c   		nkapar[j]=i;
i                 130 src/biosav.c   		    i++;
i                 131 src/biosav.c   		    if(inomat[i]==1){
i                 141 src/biosav.c   		nkepar[j]=i;
i                 146 src/biosav.c   		    i++;
i                 147 src/biosav.c   		    if(inomat[i]==1){
i                 166 src/biosav.c       for(i=0;i<num_cpus;i++){
i                 167 src/biosav.c   	ithread[i]=i;
i                 168 src/biosav.c   	pthread_create(&tid[i],NULL,(void *)biotsavartmt,(void *)&ithread[i]);
i                 170 src/biosav.c       for(i=0;i<num_cpus;i++)pthread_join(tid[i], NULL);
i                 180 src/biosav.c   void *biotsavartmt(ITG *i){
i                 184 src/biosav.c       nka=nkapar[*i]+1;
i                 185 src/biosav.c       nkb=nkepar[*i]+1;
i                  59 src/cascade.c      ITG i,j,index,id,idof,nterm,idepend,*nodempc=NULL,
i                 106 src/cascade.c      for(i=0;i<*nmpc;i++){
i                 108 src/cascade.c  	    FORTRAN(nident,(ikboun,&ikmpc[i],nboun,&id));}
i                 111 src/cascade.c  	    if(ikboun[id-1]==ikmpc[i]){
i                 112 src/cascade.c  		if(strcmp1(&labmpc[20*i],"FLUID")!=0){
i                 114 src/cascade.c  			   (ikmpc[i])/8+1,ikmpc[i]-8*((ikmpc[i])/8));
i                 117 src/cascade.c  			   (-ikmpc[i])/8+1,-ikmpc[i]-8*((-ikmpc[i])/8));
i                 128 src/cascade.c      for(i=0;i<*nmpc;i++){
i                 136 src/cascade.c  	if((strcmp1(&labmpc[20*i],"                    ")==0) ||
i                 137 src/cascade.c  	   (strcmp1(&labmpc[20*i],"CYCLIC")==0) ||
i                 138 src/cascade.c  	   (strcmp1(&labmpc[20*i],"SUBCYCLIC")==0)||
i                 139 src/cascade.c  	   (strcmp1(&labmpc[20*i],"CONTACT")==0)||
i                 140 src/cascade.c  	   (strcmp1(&labmpc[20*i],"FLUID")==0)||
i                 141 src/cascade.c             (*iperturb<2)) jmpc[i]=0;
i                 145 src/cascade.c  	else if((strcmp1(&labmpc[20*i],"RIGID")==0) ||
i                 146 src/cascade.c  	   (strcmp1(&labmpc[20*i],"KNOT")==0) ||
i                 147 src/cascade.c  	   (strcmp1(&labmpc[20*i],"PLANE")==0) ||
i                 148 src/cascade.c  	   (strcmp1(&labmpc[20*i],"STRAIGHT")==0)||
i                 149 src/cascade.c             (strcmp1(&labmpc[20*i],"ISOCHORIC")==0)) jmpc[i]=1;
i                 154 src/cascade.c  	    jmpc[i]=1;
i                 167 src/cascade.c          for(i=0;i<*nmpc;i++){
i                 169 src/cascade.c  	    if(strcmp1(&labmpc[20*i],"FLUID")!=0){
i                 175 src/cascade.c  	    if(jmpc[i]==1) nl=1;
i                 178 src/cascade.c  	    index=nodempc[3*ipompc[i]-1];
i                 229 src/cascade.c  		    index1=ipompc[i];
i                 256 src/cascade.c  		    	    index1=ipompc[i];
i                 266 src/cascade.c  		    if((strcmp1(&labmpc[20*i],"                    ")==0)&&
i                 268 src/cascade.c  			strcpy1(&labmpc[20*i],"SUBCYCLIC",9);
i                 313 src/cascade.c  	    index1=ipompc[i];
i                 341 src/cascade.c  	    index1=ipompc[i];
i                 384 src/cascade.c  	for(i=*nmpc-1;i>-1;i--){
i                 385 src/cascade.c  	    index=ipompc[i];
i                 404 src/cascade.c  		icoef[2*ifree]=i+1;
i                 420 src/cascade.c  	for(i=0;i<*nmpc;i++){
i                 421 src/cascade.c  	    idof=ikmpc[i];
i                 422 src/cascade.c  	    icolumn=ilmpc[i]-1;
i                 578 src/cascade.c  	for(i=nindep;i>0;i--){
i                 579 src/cascade.c  	    idof=indepdof[i-1];
i                 680 src/cascade.c  	for(i=0;i<*nmpc;i++){
i                 681 src/cascade.c  	    j=ilmpc[i]-1;
i                 682 src/cascade.c  	    idof=ikmpc[i];
i                 710 src/cascade.c      for(i=0;i<*nmpc;i++){
i                 711 src/cascade.c  	index=ipompc[i];
i                  77 src/ccx_2.8.c    nmethod,neq[3]={0,0,0},i,mpcfree=1,mei[4],j,nzl,nam,nbounold=0,
i                 106 src/ccx_2.8.c    for(i=1;i<argc;i++){
i                 107 src/ccx_2.8.c      if(strcmp1(argv[i],"-i")==0) {
i                 108 src/ccx_2.8.c      strcpy(jobnamec,argv[i+1]);strcpy1(jobnamef,argv[i+1],132);jin++;break;}
i                 109 src/ccx_2.8.c      if(strcmp1(argv[i],"-v")==0) {
i                 116 src/ccx_2.8.c    for(i=1;i<argc;i++){
i                 117 src/ccx_2.8.c      if(strcmp1(argv[i],"-o")==0) {
i                 118 src/ccx_2.8.c      strcpy(output,argv[i+1]);break;}
i                 222 src/ccx_2.8.c  	for(i=0;i<ne_;i++) ielprop[i]=-1;
i                 230 src/ccx_2.8.c  	for(i=0;i<2*nkon_;i++) iponor[i]=-1;
i                 256 src/ccx_2.8.c      for(i=0;i<3*memmpc_;i+=3){nodempc[i+2]=i/3+2;}
i                 474 src/ccx_2.8.c      for(i=7*nbodyold;i<7*nbody_;i++) xbodyold[i]=0;
i                 616 src/ccx_2.8.c  	for(i=0;i<ne;i++){
i                 619 src/ccx_2.8.c  		    sti[6*mi[0]*i+6*j+k]=prestr[6*mi[0]*i+6*j+k];
i                 626 src/ccx_2.8.c  	for(i=0;i<ne;i++){
i                 629 src/ccx_2.8.c  		    eme[6*mi[0]*i+6*j+k]=prestr[6*mi[0]*i+6*j+k];
i                 647 src/ccx_2.8.c  	for(i=0;i<nboun;i++){
i                 648 src/ccx_2.8.c  	    if(strcmp1(&typeboun[i],"F")==0) continue;
i                 649 src/ccx_2.8.c  	    if(ndirboun[i]==0){
i                 650 src/ccx_2.8.c  		xbounold[i]=vold[mt*(nodeboun[i]-1)];
i                 659 src/ccx_2.8.c        for(i=0;i<nk;i++) t1old[i]=t0[i];
i                 723 src/ccx_2.8.c  	for(i=nxstate;i<nstate_*mi[0]*(ne+nslavs);i++){xstate[i]=0.;}
i                 726 src/ccx_2.8.c  	for(i=nxstate;i<nstate_*mi[0]*(ne+nintpoint);i++){xstate[i]=0.;}
i                 818 src/ccx_2.8.c      for(i=nforcold;i<nforc;i++) xforcold[i]=0;
i                 821 src/ccx_2.8.c      for(i=2*nloadold;i<2*nload;i++) xloadold[i]=0;
i                 937 src/ccx_2.8.c  	  for(i=2*nload1;i<2*nload;i++)iamload[i]=0;
i                 942 src/ccx_2.8.c        for(i=mi[2]*ne0;i<mi[2]*ne;i++)ielmat[i]=1;
i                 953 src/ccx_2.8.c        for(i=0;i<nload;i++){
i                 954 src/ccx_2.8.c  	  strcpy1(&sideloadtemp[20*i],&sideload[20*i],20);
i                 955 src/ccx_2.8.c  	  if((strcmp1(&sideload[20*i]," ")==0)&&
i                 956 src/ccx_2.8.c  	     (strcmp1(&sideload[20*i+1]," ")!=0)){
i                 957 src/ccx_2.8.c  	      strcpy1(&sideloadtemp[20*i],"F",1);
i                1308 src/ccx_2.8.c        for (i=0;i<nboun;i++) {
i                1309 src/ccx_2.8.c  	  nodebounold[i]=nodeboun[i];
i                1310 src/ccx_2.8.c  	  ndirbounold[i]=ndirboun[i];
i                1319 src/ccx_2.8.c  	  for (i=0;i<nboun;i++) {
i                1320 src/ccx_2.8.c  	      if(iamboun[i]>0){
i                1321 src/ccx_2.8.c  		  if(namta[3*iamboun[i]-1]>0){
i                1322 src/ccx_2.8.c  		      iamboun[i]=0;
i                1323 src/ccx_2.8.c  		      xboun[i]=xbounold[i];}
i                1326 src/ccx_2.8.c  	  for (i=0;i<nforc;i++){
i                1327 src/ccx_2.8.c  	      if(iamforc[i]>0){
i                1328 src/ccx_2.8.c  		  if(namta[3*iamforc[i]-1]>0){
i                1329 src/ccx_2.8.c  		      iamforc[i]=0;
i                1330 src/ccx_2.8.c  		      xforc[i]=xforcold[i];}
i                1333 src/ccx_2.8.c  	  for (i=0;i<2*nload;i++){
i                1334 src/ccx_2.8.c  	      if(iamload[i]>0){
i                1335 src/ccx_2.8.c  		  if(namta[3*iamload[i]-1]>0){
i                1336 src/ccx_2.8.c  		      iamload[i]=0;
i                1337 src/ccx_2.8.c  		      xload[i]=xloadold[i];}
i                1340 src/ccx_2.8.c  	  for (i=1;i<3*nbody;i=i+3){
i                1341 src/ccx_2.8.c  	      if(ibody[i]>0){
i                1342 src/ccx_2.8.c  		  if(namta[3*ibody[i]-1]>0){
i                1343 src/ccx_2.8.c  		      ibody[i]=0;
i                1344 src/ccx_2.8.c  		      xbody[7*(i-1)/3]=xbodyold[7*(i-1)/3];}
i                1348 src/ccx_2.8.c  	      if(iamt1[i]>0){
i                1349 src/ccx_2.8.c  		  if(namta[3*iamt1[i]-1]>0){
i                1350 src/ccx_2.8.c  		      iamt1[i]=0;
i                1351 src/ccx_2.8.c  		      t1[i]=t1old[i];}
i                1375 src/ccx_2.8.c        for(i=nload-1;i>=nload0;i--){
i                1376 src/ccx_2.8.c  	  if(strcmp2(&sideload[20*i],"                    ",20)==0){
i                1377 src/ccx_2.8.c  	      iload=nelemload[2*i+1];
i                  57 src/checkconvergence.c         i,j,mt=mi[1]+1,iexceed;
i                 188 src/checkconvergence.c 	    for(i=0;i<*nk;i++){
i                 190 src/checkconvergence.c 		    veold[mt*i+j]=(vold[mt*i+j]-vini[mt*i+j])/(*dtime);
i                  15 src/compare.c      ITG     i;
i                  17 src/compare.c      i = 0;
i                  18 src/compare.c      while ((str1[i]==str2[i]) && (i<length))
i                  19 src/compare.c          i++;
i                  21 src/compare.c      return i;               /* Return how far we got before a difference occurred, or the variable      */
i                  80 src/compfluid.c       iconvergence=0,symmetryflag,inputformat,i,*inum=NULL,iitf,ifreefa,
i                 625 src/compfluid.c       for(i=0;i<neq;i++){rwork[i]=1./adv[i];}
i                 906 src/compfluid.c 	  for(i=0;i<neq;i++){rwork[i]=1./ad[i];}
i                 966 src/compfluid.c 	  for(i=0;i<*nprint;i++){
i                 967 src/compfluid.c 	      if(strcmp1(&prlab[6*i],"SVF")==0) isti=1;
i                 968 src/compfluid.c 	      if(strcmp1(&prlab[6*i],"HFLF")==0)iqfx=1;
i                  77 src/complexfreq.c   ITG nev,i,j,k,idof,*inum=NULL,*ipobody=NULL,inewton=0,id,
i                 156 src/complexfreq.c   for(i=neq[1];i>0;i--){
i                 157 src/complexfreq.c       if(icol[i-1]>0){
i                 158 src/complexfreq.c 	  *nzl=i;
i                 268 src/complexfreq.c       for(i=0;i<nev;i++){nm[i]=-1;}
i                 309 src/complexfreq.c 	  for(i=nev;i<nev+nevd;i++){nm[i]=nmd;}
i                 342 src/complexfreq.c       for(i=0;i<nev;i++){
i                 343 src/complexfreq.c 	  if((i/2)*2==i){
i                 346 src/complexfreq.c 		  d[j]=d[i];
i                 348 src/complexfreq.c 		  d[2*j]=d[2*i];d[2*j+1]=d[2*i+1];
i                 350 src/complexfreq.c 	      nm[j]=nm[i];
i                 352 src/complexfreq.c 		  z[j*neq[1]+k]=z[i*neq[1]+k];
i                 363 src/complexfreq.c       for(i=0;i<*mcs;i++){
i                 364 src/complexfreq.c 	  if(cs[17*i]>nsectors) nsectors=(ITG)(cs[17*i]+0.5);
i                 386 src/complexfreq.c       for(i=0;i<*nk;i++) inocs[i]=-1;
i                 387 src/complexfreq.c       for(i=0;i<*ne;i++) ielcs[i]=-1;
i                 390 src/complexfreq.c     for(i=0;i<*mcs;i++){
i                 391 src/complexfreq.c       is=cs[17*i+4];
i                 393 src/complexfreq.c       ielset=cs[17*i+12];
i                 399 src/complexfreq.c 	  ielcs[iel]=i;
i                 409 src/complexfreq.c 	    inocs[node]=i;
i                 418 src/complexfreq.c 	    ielcs[iel]=i;
i                 428 src/complexfreq.c 	      inocs[node]=i;
i                 443 src/complexfreq.c       for(i=nev-2;i>-1;i--){
i                 444 src/complexfreq.c 	  if(fabs(d[i])<0.0001*fabs(d[i+1])) setnull=0.;
i                 445 src/complexfreq.c 	  d[i]*=setnull;
i                 449 src/complexfreq.c       for(i=nev-2;i>-1;i--){
i                 450 src/complexfreq.c 	  if(sqrt(d[2*i]*d[2*i]+d[2*i+1]*d[2*i+1])<
i                 451 src/complexfreq.c                  0.0001*sqrt(d[2*i+2]*d[2*i+2]+d[2*i+3]*d[2*i+3])) setnull=0.;
i                 452 src/complexfreq.c 	  d[2*i]*=setnull;
i                 453 src/complexfreq.c 	  d[2*i+1]*=setnull;
i                 466 src/complexfreq.c       for(i=1;i<nev;i++){
i                 467 src/complexfreq.c 	  if(nm[i]==inmd) continue;
i                 468 src/complexfreq.c 	  iendnmd[nmd]=i;
i                 470 src/complexfreq.c 	  istartnmd[nmd]=i+1;
i                 471 src/complexfreq.c 	  inmd=nm[i];
i                 519 src/complexfreq.c       for(i=0;i<nev;i++){
i                 520 src/complexfreq.c 	  FORTRAN(op_corio,(&neq[1],&z[i*neq[1]],&zc[i*neq[1]],adc,auc,
i                 529 src/complexfreq.c       for(i=0;i<nev;i++){
i                 530 src/complexfreq.c 	  for(j=0;j<=i;j++){
i                 532 src/complexfreq.c 		  cc[i*nev+j]+=z[j*neq[1]+k]*zc[i*neq[1]+k];
i                 539 src/complexfreq.c       for(i=0;i<nev;i++){
i                 540 src/complexfreq.c 	  for(j=i+1;j<nev;j++){
i                 541 src/complexfreq.c 	      cc[i*nev+j]=-cc[j*nev+i];
i                 655 src/complexfreq.c 	  for(i=0;i<neqact;i++){
i                 660 src/complexfreq.c 		  zz[(long long)2*l*neqact+i]+=
i                 661 src/complexfreq.c 		      (xx[2*l*nev+2*k]*z[(long long)2*k*neqact+i]
i                 662 src/complexfreq.c 		       -xx[2*l*nev+2*k+1]*z[(long long)(2*k+1)*neqact+i]);
i                 666 src/complexfreq.c 		  zz[(long long)(2*l+1)*neqact+i]+=
i                 667 src/complexfreq.c 		      (xx[2*l*nev+2*k]*z[(long long)(2*k+1)*neqact+i]
i                 668 src/complexfreq.c 		       +xx[2*l*nev+2*k+1]*z[(long long)2*k*neqact+i]);
i                 771 src/complexfreq.c 	      for(i=0;i<neq[1];i++){
i                 773 src/complexfreq.c 		      zz[(2*l+j)*neq[1]+i]+=xx[2*l*nev+2*k+j]*
i                 774 src/complexfreq.c                             z[(long long)k*neq[1]+i];
i                 936 src/complexfreq.c 	  for(i=0;i<nev;i++){
i                 937 src/complexfreq.c 	      eigx[2*i]=eigxx[2*i]*eigxx[2*i]-eigxx[2*i+1]*eigxx[2*i+1];
i                 938 src/complexfreq.c 	      eigx[2*i+1]=2.*eigxx[2*i]*eigxx[2*i+1];
i                1031 src/complexfreq.c 	  for(i=0;i<nev;i++){
i                1032 src/complexfreq.c 	      eigx[2*i]=eigxx[2*i]*eigxx[2*i]-eigxx[2*i+1]*eigxx[2*i+1];
i                1033 src/complexfreq.c 	      eigx[2*i+1]=2.*eigxx[2*i]*eigxx[2*i+1];
i                1164 src/complexfreq.c       for(i=1;i<is;i++){
i                1166 src/complexfreq.c 	theta=i*2.*pi/cs[17*jj];
i                1170 src/complexfreq.c 	    cot[3*l+i*3**nk]=cot[3*l];
i                1171 src/complexfreq.c 	    cot[1+3*l+i*3**nk]=cot[1+3*l]+theta;
i                1172 src/complexfreq.c 	    cot[2+3*l+i*3**nk]=cot[2+3*l];
i                1173 src/complexfreq.c 	    if(*ntrans>0){inotrt[2*l+i*2**nk]=inotrt[2*l];}
i                1176 src/complexfreq.c 	for(l=0;l<*nkon;l++){kont[l+i**nkon]=kon[l]+i**nk;}
i                1180 src/complexfreq.c 	      ipkont[l+i**ne]=ipkon[l]+i**nkon;
i                1181 src/complexfreq.c 	      ielmatt[mi[2]*(l+i**ne)]=ielmat[mi[2]*l];
i                1184 src/complexfreq.c 		lakont[l2+i*8**ne]=lakon[l2];
i                1289 src/complexfreq.c       for(i=0;i<6*mi[0]**ne;i++){eme[i]=0.;}
i                1298 src/complexfreq.c       for(i=0;i<*nmpc;i++){
i                1299 src/complexfreq.c 	index=ipompc[i]-1;
i                1310 src/complexfreq.c 	  if(strcmp1(&labmpc[20*i],"CYCLIC")==0){
i                1311 src/complexfreq.c             icomplex=atoi(&labmpc[20*i+6]);}
i                1312 src/complexfreq.c 	  else if(strcmp1(&labmpc[20*i],"SUBCYCLIC")==0){
i                1430 src/complexfreq.c 	for(i=1;i<is;i++){
i                1432 src/complexfreq.c 	  for(l=0;l<*nk;l++){inumt[l+i**nk]=inum[l];}
i                1434 src/complexfreq.c 	  theta=i*nm[j]*2.*pi/cs[17*jj];
i                1450 src/complexfreq.c 		      vt[l+mt**nk*i]=v[l];
i                1457 src/complexfreq.c 		  vt[l+mt**nk*i]=ctl*v[l]-stl*v[l+mt**nk];
i                1479 src/complexfreq.c 		      vt[l+mt**nk*(i+ngraph)]=v[l+mt**nk];
i                1486 src/complexfreq.c 		  vt[l+mt**nk*(i+ngraph)]=stl*v[l]+ctl*v[l+mt**nk];
i                1494 src/complexfreq.c 	      if(inocs[l]==jj) t1t[l+*nk*i]=t1[l];
i                1510 src/complexfreq.c 		      stnt[l+6**nk*i]=stn[l];
i                1517 src/complexfreq.c 		  stnt[l+6**nk*i]=ctl*stn[l]-stl*stn[l+6**nk];
i                1538 src/complexfreq.c 		      stnt[l+6**nk*(i+ngraph)]=stn[l+6**nk];
i                1545 src/complexfreq.c 		  stnt[l+6**nk*(i+ngraph)]=stl*stn[l]+ctl*stn[l+6**nk];
i                1563 src/complexfreq.c 		      eent[l+6**nk*i]=een[l];
i                1570 src/complexfreq.c 		  eent[l+6**nk*i]=ctl*een[l]-stl*een[l+6**nk];
i                1591 src/complexfreq.c 				  eent[l+6**nk*(i+ngraph)]=een[l+6**nk];
i                1598 src/complexfreq.c 			  eent[l+6**nk*(i+ngraph)]=stl*een[l]+ctl*een[l+6**nk];
i                1616 src/complexfreq.c 		      emnt[l+6**nk*i]=emn[l];
i                1623 src/complexfreq.c 		  emnt[l+6**nk*i]=ctl*emn[l]-stl*emn[l+6**nk];
i                1644 src/complexfreq.c 				  emnt[l+6**nk*(i+ngraph)]=emn[l+6**nk];
i                1651 src/complexfreq.c 			  emnt[l+6**nk*(i+ngraph)]=stl*emn[l]+ctl*emn[l+6**nk];
i                1669 src/complexfreq.c 		      fnt[l+mt**nk*i]=fn[l];
i                1676 src/complexfreq.c 		  fnt[l+mt**nk*i]=ctl*fn[l]-stl*fn[l+mt**nk];
i                1697 src/complexfreq.c 				  fnt[l+mt**nk*(i+ngraph)]=fn[l+mt**nk];
i                1704 src/complexfreq.c 			  fnt[l+mt**nk*(i+ngraph)]=stl*fn[l]+ctl*fn[l+mt**nk];
i                1712 src/complexfreq.c 	      if(inocs[l]==jj) enernt[l+*nk*i]=0.;
i                  44 src/contact.c      ITG i,ntrimax,*nx=NULL,*ny=NULL,*nz=NULL,*ipompc=NULL,*ikmpc=NULL,
i                  71 src/contact.c      for(i=0;i<*ntie;i++){
i                  72 src/contact.c  	if(itietri[2*i+1]-itietri[2*i]+1>ntrimax)
i                  73 src/contact.c  	    ntrimax=itietri[2*i+1]-itietri[2*i]+1;
i                  77 src/dyna.c       ITG nev,i,j,k,idof,*inum=NULL,*ipobody=NULL,inewton=0,id,
i                 186 src/dyna.c       for(i=neq[1];i>0;i--){
i                 187 src/dyna.c           if(icol[i-1]>0){
i                 188 src/dyna.c     	  *nzl=i;
i                 250 src/dyna.c     	  for(i=0;i<*nload;i++){
i                 251 src/dyna.c     	      iload=i+1;
i                 257 src/dyna.c     	  for(i=0;i<*nforc;i++){
i                 258 src/dyna.c     	      iforc=i+1;
i                 298 src/dyna.c       for(i=0;i<nmddof;i++){imddof[i]-=1;}
i                 325 src/dyna.c           for(i=0;i<nev;i++){
i                 326 src/dyna.c     	  if(d[i]>0){d[i]=sqrt(d[i]);}else{d[i]=0.;}
i                 384 src/dyna.c     	  for(i=nev;i<nev+nevd;i++){
i                 385 src/dyna.c     	      if(d[i]>0){d[i]=sqrt(d[i]);}else{d[i]=0.;}
i                 388 src/dyna.c     	  for(i=nev;i<nev+nevd;i++){nm[i]=nmd;}
i                 420 src/dyna.c           for(i=0;i<*mcs;i++){
i                 422 src/dyna.c     	  if(cs[17*i]>nsectors) nsectors=(ITG)(cs[17*i]+0.5);
i                 440 src/dyna.c           for(i=3**nk;i<3**nk*nsectors;i++){co[i]=0.;}
i                 451 src/dyna.c           for(i=*ne;i<*ne*nsectors;i++){ipkon[i]=-1;}
i                 472 src/dyna.c           for(i=0;i<*nmpc;i++){ipompcold[i]=ipompc[i];}
i                 473 src/dyna.c           for(i=0;i<3**mpcend;i++){nodempcold[i]=nodempc[i];}
i                 474 src/dyna.c           for(i=0;i<*mpcend;i++){coefmpcold[i]=coefmpc[i];}
i                 475 src/dyna.c           for(i=0;i<20**nmpc;i++){labmpcold[i]=labmpc[i];}
i                 476 src/dyna.c           for(i=0;i<*nmpc;i++){ikmpcold[i]=ikmpc[i];}
i                 477 src/dyna.c           for(i=0;i<*nmpc;i++){ilmpcold[i]=ilmpc[i];}
i                 495 src/dyna.c           for(i=0; i<*ntie; i++){
i                 496 src/dyna.c     	if(tieset[i*(81*3)+80]=='C'){
i                 500 src/dyna.c     	  memcpy(tchar2,&tieset[i*(81*3)+81],81);
i                 502 src/dyna.c     	  memcpy(tchar3,&tieset[i*(81*3)+81+81],81);
i                 595 src/dyna.c       for(i=0;i<*nmpc;i++){
i                 596 src/dyna.c           if((strcmp1(&labmpc[20*i]," ")!=0)&&
i                 597 src/dyna.c              (strcmp1(&labmpc[20*i],"CONTACT")!=0)&&
i                 598 src/dyna.c              (strcmp1(&labmpc[20*i],"CYCLIC")!=0)&&
i                 599 src/dyna.c              (strcmp1(&labmpc[20*i],"SUBCYCLIC")!=0)){
i                 624 src/dyna.c       for(i=nev-2;i>-1;i--){
i                 625 src/dyna.c           if(fabs(d[i])<0.0001*fabs(d[i+1])) setnull=0.;
i                 626 src/dyna.c           d[i]*=setnull;
i                 632 src/dyna.c       for(i=0;i<*ne;i++){
i                 633 src/dyna.c           if(ipkon[i]<0) continue;
i                 634 src/dyna.c           if(strcmp1(&lakon[i*8],"ED")==0){
i                 668 src/dyna.c           for(i=0;i<nev;i++){
i                 669 src/dyna.c     	  FORTRAN(op,(&neq[1],&z[i*neq[1]],&zc[i*neq[1]],adc,auc,
i                 677 src/dyna.c           for(i=0;i<nev;i++){
i                 678 src/dyna.c     	  for(j=0;j<=i;j++){
i                 680 src/dyna.c     		  cc[i*nev+j]+=z[j*neq[1]+k]*zc[i*neq[1]+k];
i                 687 src/dyna.c           for(i=0;i<nev;i++){
i                 688 src/dyna.c     	  for(j=i;j<nev;j++){
i                 689 src/dyna.c     	      cc[i*nev+j]=cc[j*nev+i];
i                 772 src/dyna.c           for(i=0;i<nev;i++){
i                 773 src/dyna.c     	if(fabs(d[i])>(1.e-10)){
i                 774 src/dyna.c     	  zeta[i]=(alpham+betam*d[i]*d[i])/(2.*d[i]);
i                 779 src/dyna.c     	  zeta[i]=0.;
i                 795 src/dyna.c         for(i=0; i<imax; i++){
i                 796 src/dyna.c           zeta[i]=xmodal[11+i];     
i                 811 src/dyna.c         for(i=0;i<neq[1];i++){temp_array1[i]=0;temp_array2[i]=0;}
i                 815 src/dyna.c         for(i=0;i<*nk;i++){
i                 817 src/dyna.c     	if(nactdof[mt*i+j]!=0){
i                 818 src/dyna.c     	  idof=nactdof[mt*i+j]-1;
i                 819 src/dyna.c     	  temp_array1[idof]=vold[mt*i+j];
i                 826 src/dyna.c         for(i=0;i<neq[1];i++){
i                 828 src/dyna.c     	cd[k]+=z[k*neq[1]+i]*temp_array2[i];
i                 834 src/dyna.c         for(i=0;i<neq[1];i++){temp_array1[i]=0;temp_array2[i]=0;}
i                 835 src/dyna.c         for(i=0;i<*nk;i++){
i                 837 src/dyna.c     	if(nactdof[mt*i+j]!=0){
i                 838 src/dyna.c     	  idof=nactdof[mt*i+j]-1;
i                 839 src/dyna.c     	  temp_array1[idof]=veold[mt*i+j];
i                 846 src/dyna.c         for(i=0;i<neq[1];i++){
i                 848 src/dyna.c     	cv[k]+=z[k*neq[1]+i]*temp_array2[i];
i                 924 src/dyna.c       for(i=0;i<*nboun;i++){
i                 925 src/dyna.c           if(fabs(xboun[i])>1.e-10){
i                1059 src/dyna.c           for(i=ne0;i<*ne;i++){
i                1060 src/dyna.c     	  indexe=ipkon[i];
i                1061 src/dyna.c     	  imat=ielmat[mi[2]*i];
i                1063 src/dyna.c     	  for(j=0;j<8;j++){lakonl[j]=lakon[8*i+j];}
i                1097 src/dyna.c           for(i=0;i<nev;i++){
i                1098 src/dyna.c     	  i2=(long long)i*neq[1];
i                1099 src/dyna.c     	  aamech[i]=0.;
i                1102 src/dyna.c     		  aamech[i]+=z[i2+ikactmech[j]]*b[ikactmech[j]];
i                1106 src/dyna.c     		  aamech[i]+=z[i2+j]*b[j];
i                1109 src/dyna.c     	  aanew[i]=aamech[i];
i                1112 src/dyna.c     		  aanew[i]+=z[i2+ikactcont[j]]*bcont[ikactcont[j]];
i                1117 src/dyna.c           for(i=0;i<nev;i++){aamech[i]=0.;}
i                1122 src/dyna.c     		  for(i=0;i<nev;i++){
i                1123 src/dyna.c     		      aamech[i]+=z[(long long)i*nzdof+id-1]*b[ikactmech[j]];
i                1134 src/dyna.c     		      for(i=0;i<nev;i++){
i                1135 src/dyna.c     			  aanew[i]+=z[(long long)i*nzdof+id-1]*bcont[ikactcont[j]];
i                1163 src/dyna.c           for(i=0;i<*nprint;i++){
i                1164 src/dyna.c     	  if((strcmp1(&prlab[6*i],"S")==0)||
i                1165 src/dyna.c     	     (strcmp1(&prlab[6*i],"E")==0)||
i                1166 src/dyna.c     	     (strcmp1(&prlab[6*i],"PEEQ")==0)||
i                1167 src/dyna.c     	     (strcmp1(&prlab[6*i],"ENER")==0)||
i                1168 src/dyna.c     	     (strcmp1(&prlab[6*i],"SDV")==0)||
i                1169 src/dyna.c     	     (strcmp1(&prlab[6*i],"CDIS")==0)||
i                1170 src/dyna.c     	     (strcmp1(&prlab[6*i],"CSTR")==0)||
i                1171 src/dyna.c     	     (strcmp1(&prlab[6*i],"CELS")==0)||
i                1172 src/dyna.c     	     (strcmp1(&prlab[6*i],"RF")==0)) {intpointvar=1;break;}
i                1180 src/dyna.c           for(i=0;i<*nprint;i++){
i                1181 src/dyna.c     	  if((strcmp1(&prlab[6*i],"HFL")==0)||
i                1182 src/dyna.c     	     (strcmp1(&prlab[6*i],"RFL")==0)) {intpointvar=1;break;}
i                1208 src/dyna.c     	    for(i=0;i<nmdnode;i++){
i                1209 src/dyna.c     		i1=mt*(imdnode[i]-1);
i                1321 src/dyna.c     		  for(i=0;i<nactmech;i++){
i                1322 src/dyna.c     		      if(fabs(b[ikactmech[i]])>bbmax) bbmax=fabs(b[ikactmech[i]]);
i                1325 src/dyna.c     		  for(i=0;i<neq[1];i++){
i                1326 src/dyna.c     		      if(fabs(b[i])>bbmax) bbmax=fabs(b[i]);
i                1334 src/dyna.c     		  for(i=0;i<nactmech;i++){
i                1335 src/dyna.c     		      if(fabs(b[ikactmech[i]]-bold[ikactmech[i]])>bbmax) 
i                1336 src/dyna.c     			  bbmax=fabs(b[ikactmech[i]]-bold[ikactmech[i]]);
i                1339 src/dyna.c     		  for(i=0;i<neq[1];i++){
i                1340 src/dyna.c     		      if(fabs(b[i])>bbmax) bbmax=fabs(b[i]-bold[i]);
i                1347 src/dyna.c     		  for(i=0;i<nactmech;i++){
i                1348 src/dyna.c     		      bold[ikactmech[i]]=b[ikactmech[i]];
i                1436 src/dyna.c     	      for(i=0;i<nactcont;i++){
i                1437 src/dyna.c     		  jdof=ikactcont[i];
i                1451 src/dyna.c     	  for(i=0;i<nev;i++){
i                1452 src/dyna.c     	      i2=(long long)i*neq[1];
i                1453 src/dyna.c     	      aa[i]=aanew[i];
i                1454 src/dyna.c     	      if(iabsload==2){aamech[i]=0.;}
i                1457 src/dyna.c     		      aamech[i]+=z[i2+ikactmech[j]]*b[ikactmech[j]];
i                1461 src/dyna.c     		      aamech[i]+=z[i2+j]*b[j];
i                1465 src/dyna.c     	      aanew[i]=aamech[i];
i                1468 src/dyna.c     		      aanew[i]+=z[i2+ikactcont[j]]*bcont[ikactcont[j]];
i                1472 src/dyna.c     	      bb[i]=(aanew[i]-aa[i])/dtime;
i                1473 src/dyna.c     	      aa[i]=aanew[i]-bb[i]*time;
i                1476 src/dyna.c     	  for(i=0;i<nev;i++){
i                1478 src/dyna.c     	      if(iabsload==2){aamech[i]=0.;}
i                1484 src/dyna.c     		      for(i=0;i<nev;i++){
i                1485 src/dyna.c     			  aamech[i]+=z[(long long)i*nzdof+id-1]*b[ikactmech[j]];
i                1496 src/dyna.c     			  for(i=0;i<nev;i++){
i                1497 src/dyna.c     			      aanew[i]+=z[(long long)i*nzdof+id-1]*bcont[ikactcont[j]];
i                1503 src/dyna.c     	  for(i=0;i<nev;i++){
i                1504 src/dyna.c           	  	bb[i]=(aanew[i]-aa[i])/(dtime);
i                1505 src/dyna.c     	  	aa[i]=aanew[i]-bb[i]*time;
i                1519 src/dyna.c     	      for(i=0;i<liw;i++){iwork[i]=0;}
i                1522 src/dyna.c     	      for(i=0;i<lrw;i++){rwork[i]=0.;}
i                1633 src/dyna.c     	      for(i=0;i<nmddof;i++){
i                1634 src/dyna.c     		  b[imddof[i]]=bmin[imddof[i]];
i                1635 src/dyna.c     		  bp[imddof[i]]=bv[imddof[i]];
i                1644 src/dyna.c     	      for(i=0;i<nmddof;i++){
i                1645 src/dyna.c     		  b[imddof[i]]=0.;
i                1646 src/dyna.c     		  bp[imddof[i]]=0.;
i                1653 src/dyna.c     	      for(i=0;i<neq[1];i++){
i                1655 src/dyna.c     		      b[i]+=bj[j]*z[(long long)j*neq[1]+i];
i                1656 src/dyna.c     		      bp[i]+=bjp[j]*z[(long long)j*neq[1]+i];
i                1660 src/dyna.c     	      for(i=0;i<nmddof;i++){
i                1662 src/dyna.c     		      b[imddof[i]]+=bj[j]*z[(long long)j*neq[1]+imddof[i]];
i                1663 src/dyna.c     		      bp[imddof[i]]+=bjp[j]*z[(long long)j*neq[1]+imddof[i]];
i                1668 src/dyna.c     	  for(i=0;i<nmddof;i++){
i                1669 src/dyna.c     	      FORTRAN(nident,(izdof,&imddof[i],&nzdof,&id));
i                1671 src/dyna.c     		  if(izdof[id-1]==imddof[i]){
i                1673 src/dyna.c     			  b[imddof[i]]+=bj[j]*z[(long long)j*nzdof+id-1];
i                1674 src/dyna.c     			  bp[imddof[i]]+=bjp[j]*z[(long long)j*nzdof+id-1];
i                1839 src/dyna.c     	      for(i=1;i<mt**nk;i=i+mt){
i                1840 src/dyna.c     		  if(fabs(v[i])>resultmax) resultmax=fabs(v[i]);}
i                1841 src/dyna.c     	      for(i=2;i<mt**nk;i=i+mt){
i                1842 src/dyna.c     		  if(fabs(v[i])>resultmax) resultmax=fabs(v[i]);}
i                1843 src/dyna.c     	      for(i=3;i<mt**nk;i=i+mt){
i                1844 src/dyna.c     		  if(fabs(v[i])>resultmax) resultmax=fabs(v[i]);}
i                1846 src/dyna.c     	      for(i=0;i<mt**nk;i=i+mt){
i                1847 src/dyna.c     		  if(fabs(v[i])>resultmax) resultmax=fabs(v[i]);}
i                1872 src/dyna.c           for(i=0;i<neq[1];i++){
i                1874 src/dyna.c     	      b[i]+=bj[j]*z[(long long)j*neq[1]+i];
i                1875 src/dyna.c     	      bp[i]+=bjp[j]*z[(long long)j*neq[1]+i];
i                2037 src/dyna.c           for(i=0;i<*nmpc;i++){ipompc[i]=ipompcold[i];}
i                2038 src/dyna.c           for(i=0;i<3**mpcend;i++){nodempc[i]=nodempcold[i];}
i                2039 src/dyna.c           for(i=0;i<*mpcend;i++){coefmpc[i]=coefmpcold[i];}
i                2040 src/dyna.c           for(i=0;i<20**nmpc;i++){labmpc[i]=labmpcold[i];}
i                2041 src/dyna.c           for(i=0;i<*nmpc;i++){ikmpc[i]=ikmpcold[i];}
i                2042 src/dyna.c           for(i=0;i<*nmpc;i++){ilmpc[i]=ilmpcold[i];}
i                2054 src/dyna.c           for(i=0;i<*nload;i++){
i                2055 src/dyna.c     	  if(nelemload[2*i]<nsectors){
i                2056 src/dyna.c     	      nelemload[2*i]-=*ne*nelemload[2*i+1];
i                2058 src/dyna.c     	      nelemload[2*i]-=*ne*(nelemload[2*i+1]-nsectors);
i                2074 src/dyna.c           for(i=0;i<*nforc;i++){
i                2075 src/dyna.c     	  if(nodeforc[2*i+1]<nsectors){
i                2076 src/dyna.c     	      nodeforc[2*i]-=*nk*nodeforc[2*i+1];
i                2078 src/dyna.c     	      nodeforc[2*i]-=*nk*(nodeforc[2*i+1]-nsectors);
i                 100 src/dynacont.c   ITG i,j,k,l,init,*itg=NULL,ntg=0,maxlenmpc,icascade=0,loop,
i                 170 src/dynacont.c   for(i=*ne0;i<*ne;i++){
i                 171 src/dynacont.c       indexe=ipkon[i];
i                 172 src/dynacont.c       imat=ielmat[mi[2]*i];
i                 174 src/dynacont.c       for(j=0;j<8;j++){lakonl[j]=lakon[8*i+j];}
i                 203 src/dynacont.c   for(i=0;i<nactcont1;i++){
i                 204 src/dynacont.c       jdof=ikactcont1[i];
i                 231 src/dynacont.c       for(i=0;i<*nactcont;i++){
i                 232 src/dynacont.c 	  jdof=ikactcont[i];
i                 257 src/dynacont.c 	  for(i=0;i<*nmdnode;i++){
i                 258 src/dynacont.c 	      i1=mt*(imdnode[i]-1);
i                 275 src/dynacont.c     for(i=0;i<*nev;i++){
i                 276 src/dynacont.c       aa[i]+=bb[i]*(*time-*dtime);
i                 367 src/dynacont.c 	  for(i=0;i<*nev;i++){
i                 368 src/dynacont.c 	      i2=(long long)i*neq[1];
i                 370 src/dynacont.c 	      if(iabsload==2){aamech[i]=0.;}
i                 373 src/dynacont.c 		      aamech[i]+=z[i2+ikactmech[j]]*b[ikactmech[j]];
i                 377 src/dynacont.c 		      aamech[i]+=z[i2+j]*b[j];
i                 382 src/dynacont.c 	  for(i=0;i<*nev;i++){
i                 383 src/dynacont.c 	      if(iabsload==2){aamech[i]=0.;}
i                 389 src/dynacont.c 		      for(i=0;i<*nev;i++){
i                 390 src/dynacont.c 			  aamech[i]+=z[i**nzdof+id-1]*b[ikactmech[j]];
i                 406 src/dynacont.c     for(i=0;i<*nev;i++){
i                 408 src/dynacont.c       aanew[i]=aamech[i];
i                 410 src/dynacont.c       bb[i]=(aanew[i]-aa[i])/(*dtime);
i                 411 src/dynacont.c       aa[i]=aanew[i]-bb[i]**time;
i                 626 src/dynacont.c 	      for(i=0;i<*nmddof;i++){
i                 627 src/dynacont.c 		  b[imddof[i]]=bmin[imddof[i]];
i                 628 src/dynacont.c 		  bp[imddof[i]]=bv[imddof[i]];
i                 637 src/dynacont.c 	      for(i=0;i<*nmddof;i++){
i                 638 src/dynacont.c 		  b[imddof[i]]=0.;
i                 639 src/dynacont.c 		  bp[imddof[i]]=0.;
i                 646 src/dynacont.c 	      for(i=0;i<neq[1];i++){
i                 648 src/dynacont.c 		      b[i]+=bj[j]*z[(long long)j*neq[1]+i];
i                 649 src/dynacont.c 		      bp[i]+=bjp[j]*z[(long long)j*neq[1]+i];
i                 653 src/dynacont.c 	      for(i=0;i<*nmddof;i++){
i                 655 src/dynacont.c 		      b[imddof[i]]+=bj[j]*z[(long long)j*neq[1]+imddof[i]];
i                 656 src/dynacont.c 		      bp[imddof[i]]+=bjp[j]*z[(long long)j*neq[1]+imddof[i]];
i                 661 src/dynacont.c 	  for(i=0;i<*nmddof;i++){
i                 662 src/dynacont.c 	      FORTRAN(nident,(izdof,&imddof[i],nzdof,&id));
i                 664 src/dynacont.c 		  if(izdof[id-1]==imddof[i]){
i                 666 src/dynacont.c 			  b[imddof[i]]+=bj[j]*z[j**nzdof+id-1];
i                 667 src/dynacont.c 			  bp[imddof[i]]+=bjp[j]*z[j**nzdof+id-1];
i                 721 src/dynacont.c       for(i=*ne0;i<*ne;i++){
i                 722 src/dynacont.c 	indexe=ipkon[i];
i                 723 src/dynacont.c 	imat=ielmat[mi[2]*i];
i                 725 src/dynacont.c 	for(j=0;j<8;j++){lakonl[j]=lakon[8*i+j];}
i                 818 src/dynacont.c       for(i=0;i<*nev;i++){
i                 819 src/dynacont.c 	if(fabs(bj[i])>bjmax) bjmax=fabs(bj[i]);
i                 820 src/dynacont.c 	if(fabs(dbj[i])>dbjmax) dbjmax=fabs(dbj[i]);
i                 882 src/dynacont.c   for(i=0;i<*nev;i++){
i                 883 src/dynacont.c     aa[i]+=bb[i]*(*time-*dtime);
i                 891 src/dynacont.c       for(i=0;i<*nev;i++){
i                 892 src/dynacont.c 	  i2=(long long)i*neq[1];
i                 894 src/dynacont.c 	  aanew[i]=aamech[i];
i                 896 src/dynacont.c 	      aanew[i]+=z[i2+ikactcont[j]]*bcont[ikactcont[j]];
i                 899 src/dynacont.c 	  bb[i]=(aanew[i]-aa[i])/(*dtime);
i                 900 src/dynacont.c 	  aa[i]=aanew[i]-bb[i]**time;
i                 908 src/dynacont.c 		  for(i=0;i<*nev;i++){
i                 909 src/dynacont.c 		      aanew[i]+=z[i**nzdof+id-1]*bcont[ikactcont[j]];
i                 914 src/dynacont.c       for(i=0;i<*nev;i++){
i                 915 src/dynacont.c       	  bb[i]=(aanew[i]-aa[i])/(*dtime);
i                 916 src/dynacont.c 	  aa[i]=aanew[i]-bb[i]**time;
i                  48 src/dynboun.c      ITG idiff[3],i,j,ic,ir,im,symmetryflag=0;
i                  76 src/dynboun.c        for(i=0;i<*nboun;i++){
i                  77 src/dynboun.c  	  xbounmin[i]=xbounold[i];
i                  78 src/dynboun.c  	  xbounact[i]=xbounold[i];
i                 106 src/dynboun.c        for(i=0;i<*nboun;i++){
i                 107 src/dynboun.c  	  ic=neq[1]+i;
i                 110 src/dynboun.c  	      bmin[ir]=bmin[ir]-au[j]*xbounmin[i];
i                 139 src/dynboun.c    for(i=0;i<*nboun;i++){
i                 140 src/dynboun.c        if(fabs(xbounact[i]-xbounmin[i])>1.e-10){
i                 146 src/dynboun.c        for(i=0;i<*nboun;i++){
i                 147 src/dynboun.c  	  ic=neq[1]+i;
i                 150 src/dynboun.c  	      bact[ir]=bact[ir]-au[j]*xbounact[i];
i                 179 src/dynboun.c    for(i=0;i<*nboun;i++){
i                 180 src/dynboun.c        if(fabs(xbounplus[i]-xbounact[i])>1.e-10){
i                 186 src/dynboun.c        for(i=0;i<*nboun;i++){
i                 187 src/dynboun.c  	  ic=neq[1]+i;
i                 190 src/dynboun.c  	      bplus[ir]=bplus[ir]-au[j]*xbounplus[i];
i                 220 src/dynboun.c  	  for(i=0;i<neq[1];i++){bplus[i]=bact[i];}
i                 222 src/dynboun.c        for(i=0;i<neq[1];i++){
i                 226 src/dynboun.c  	  bv[i]=(bplus[i]-bmin[i])/deltatime2;
i                 230 src/dynboun.c  	  ba[i]=(bmin[i]-2.*bact[i]+bplus[i])/deltatimesq;
i                 232 src/dynboun.c  	  b1[i]=ba[i]+*alpham*bv[i];
i                 233 src/dynboun.c  	  b2[i]=*betam*bv[i];
i                 235 src/dynboun.c  	  bmin[i]=bact[i];
i                 236 src/dynboun.c  	  bact[i]=bplus[i];
i                 240 src/dynboun.c        for(i=0;i<neq[1];i++){bnew[i]=-bplus[i];}
i                 243 src/dynboun.c  	  for(i=0;i<neq[1];i++){
i                 244 src/dynboun.c  	      bnew[i]-=bplus[i];
i                 245 src/dynboun.c  	      b[i]+=bnew[i];
i                 248 src/dynboun.c  	  for(i=0;i<neq[1];i++){
i                 249 src/dynboun.c  	      bnew[i]-=bplus[i];
i                 250 src/dynboun.c  	      bdiff[i]=bnew[i]-bprev[i];
i                 251 src/dynboun.c  	      b[i]+=bdiff[i];
i                 256 src/dynboun.c  	  for(i=0;i<neq[1];i++){
i                 257 src/dynboun.c  	      bnew[i]-=bplus[i];
i                 258 src/dynboun.c  	      bdiff[i]+=bnew[i]-bprev[i];
i                 259 src/dynboun.c  	      b[i]+=bdiff[i];
i                 271 src/dynboun.c  	  for(i=0;i<neq[1];i++){
i                 272 src/dynboun.c  	      bdiff[i]=-bprev[i];
i                 273 src/dynboun.c  	      b[i]+=bdiff[i];
i                 279 src/dynboun.c  	  for(i=0;i<neq[1];i++){
i                 280 src/dynboun.c  	      bdiff[i]+=-bprev[i];
i                 281 src/dynboun.c  	      b[i]+=bdiff[i];
i                 102 src/electromagnetics.c       coriolis=0,*ipneigh=NULL,*neigh=NULL,i,icfd=0,id,node,networknode,
i                 390 src/electromagnetics.c   for(i=0;i<*ne;i++){
i                 391 src/electromagnetics.c       if(strcmp1(&lakon[8*i+6],"L")!=0){
i                 392 src/electromagnetics.c 	  ipkon[i]=-ipkon[i]-2;
i                 553 src/electromagnetics.c   for(i=0;i<*ne;i++){
i                 554 src/electromagnetics.c       if(strcmp1(&lakon[8*i+6],"L")!=0){
i                 555 src/electromagnetics.c 	  ipkon[i]=-ipkon[i]-2;
i                 556 src/electromagnetics.c 	  indexe=ipkon[i];
i                 557 src/electromagnetics.c 	  if(strcmp1(&lakon[8*i+3],"6")==0){nope=6;}
i                 558 src/electromagnetics.c 	  else if(strcmp1(&lakon[8*i+3],"8")==0){nope=8;}
i                 559 src/electromagnetics.c 	  else if(strcmp1(&lakon[8*i+3],"1")==0){nope=15;}
i                 626 src/electromagnetics.c   for(i=0;i<*ne;i++){
i                 627 src/electromagnetics.c       if(strcmp1(&lakon[8*i+6],"L")==0){
i                 628 src/electromagnetics.c 	  ipkon[i]=-ipkon[i]-2;
i                 687 src/electromagnetics.c   for(i=1;i<3;i++){
i                 688 src/electromagnetics.c       imast=iactive[i];
i                 851 src/electromagnetics.c       for(i=0;i<3**nk;i++){h0[i]/=h0scale;}
i                 862 src/electromagnetics.c       for(i=0;i<3**nk;i++){h0[i]*=h0scale;}
i                 864 src/electromagnetics.c       for(i=0;i<3;i++){cam[i]=0.;}for(i=3;i<5;i++){cam[i]=0.5;}
i                 958 src/electromagnetics.c 	      for(i=0;i<3**nk;i++){h0[i]/=h0scale;}
i                 970 src/electromagnetics.c 	      for(i=0;i<3**nk;i++){h0[i]*=h0scale;}
i                 972 src/electromagnetics.c 	      for(i=0;i<3;i++){cam[i]=0.;}for(i=3;i<5;i++){cam[i]=0.5;}
i                  19 src/elemChecker.c   ITG i, j;
i                  25 src/elemChecker.c   for (i=0; i<sum_e; i++)
i                  27 src/elemChecker.c     if( (elem[i].type == 1)||(elem[i].type == 4) )   /* HEXA */
i                  29 src/elemChecker.c       v_result( &node[elem[i].nod[0]].nx, &node[elem[i].nod[1]].nx, v12);
i                  30 src/elemChecker.c       v_result( &node[elem[i].nod[0]].nx, &node[elem[i].nod[3]].nx, v13);
i                  31 src/elemChecker.c       v_result( &node[elem[i].nod[0]].nx, &node[elem[i].nod[4]].nx, v15);
i                  45 src/elemChecker.c           epuf[j] = elem[i].nod[j];
i                  47 src/elemChecker.c         elem[i].nod[0] = epuf[4];
i                  48 src/elemChecker.c         elem[i].nod[1] = epuf[5];
i                  49 src/elemChecker.c         elem[i].nod[2] = epuf[6];
i                  50 src/elemChecker.c         elem[i].nod[3] = epuf[7];
i                  51 src/elemChecker.c         elem[i].nod[4] = epuf[0];
i                  52 src/elemChecker.c         elem[i].nod[5] = epuf[1];
i                  53 src/elemChecker.c         elem[i].nod[6] = epuf[2];
i                  54 src/elemChecker.c         elem[i].nod[7] = epuf[3];
i                  56 src/elemChecker.c         if (elem[i].type == 4)
i                  60 src/elemChecker.c             epuf[j] = elem[i].nod[j+8];
i                  61 src/elemChecker.c             epuf[j+4] = elem[i].nod[j+16];
i                  63 src/elemChecker.c           elem[i].nod[8] = epuf[4];
i                  64 src/elemChecker.c           elem[i].nod[9] = epuf[5];
i                  65 src/elemChecker.c           elem[i].nod[10] = epuf[6];
i                  66 src/elemChecker.c           elem[i].nod[11] = epuf[7];
i                  67 src/elemChecker.c           elem[i].nod[16] = epuf[0];
i                  68 src/elemChecker.c           elem[i].nod[17] = epuf[1];
i                  69 src/elemChecker.c           elem[i].nod[18] = epuf[2];
i                  70 src/elemChecker.c           elem[i].nod[19] = epuf[3];
i                  76 src/elemChecker.c     else if( (elem[i].type == 2)||(elem[i].type == 5) )   /* PENTA */
i                  78 src/elemChecker.c       v_result( &node[elem[i].nod[0]].nx, &node[elem[i].nod[1]].nx, v12);
i                  79 src/elemChecker.c       v_result( &node[elem[i].nod[0]].nx, &node[elem[i].nod[2]].nx, v13);
i                  80 src/elemChecker.c       v_result( &node[elem[i].nod[0]].nx, &node[elem[i].nod[3]].nx, v15);
i                  94 src/elemChecker.c           epuf[j] = elem[i].nod[j];
i                  96 src/elemChecker.c         elem[i].nod[0] = epuf[3];
i                  97 src/elemChecker.c         elem[i].nod[1] = epuf[4];
i                  98 src/elemChecker.c         elem[i].nod[2] = epuf[5];
i                  99 src/elemChecker.c         elem[i].nod[3] = epuf[0];
i                 100 src/elemChecker.c         elem[i].nod[4] = epuf[1];
i                 101 src/elemChecker.c         elem[i].nod[5] = epuf[2];
i                 103 src/elemChecker.c         if (elem[i].type == 5)
i                 107 src/elemChecker.c             epuf[j] = elem[i].nod[j+6];
i                 108 src/elemChecker.c             epuf[j+3] = elem[i].nod[j+12];
i                 110 src/elemChecker.c           elem[i].nod[6] = epuf[3];
i                 111 src/elemChecker.c           elem[i].nod[7] = epuf[4];
i                 112 src/elemChecker.c           elem[i].nod[8] = epuf[5];
i                 113 src/elemChecker.c           elem[i].nod[12] = epuf[0];
i                 114 src/elemChecker.c           elem[i].nod[13] = epuf[1];
i                 115 src/elemChecker.c           elem[i].nod[14] = epuf[2];
i                 121 src/elemChecker.c     else if( (elem[i].type == 3)||(elem[i].type == 6) )   /* TET */
i                 123 src/elemChecker.c       v_result( &node[elem[i].nod[0]].nx, &node[elem[i].nod[1]].nx, v12);
i                 124 src/elemChecker.c       v_result( &node[elem[i].nod[0]].nx, &node[elem[i].nod[2]].nx, v13);
i                 125 src/elemChecker.c       v_result( &node[elem[i].nod[0]].nx, &node[elem[i].nod[3]].nx, v15);
i                 137 src/elemChecker.c         for (j=0; j<4; j++) epuf[j] = elem[i].nod[j];
i                 138 src/elemChecker.c         elem[i].nod[0] = epuf[1];
i                 139 src/elemChecker.c         elem[i].nod[1] = epuf[2];
i                 140 src/elemChecker.c         elem[i].nod[2] = epuf[3];
i                 141 src/elemChecker.c         elem[i].nod[3] = epuf[0];
i                 143 src/elemChecker.c         if (elem[i].type == 6)
i                 145 src/elemChecker.c           for (j=4; j<10; j++) epuf[j] = elem[i].nod[j];
i                 146 src/elemChecker.c           elem[i].nod[4] = epuf[5];
i                 147 src/elemChecker.c           elem[i].nod[5] = epuf[9];
i                 148 src/elemChecker.c           elem[i].nod[6] = epuf[8];
i                 149 src/elemChecker.c           elem[i].nod[7] = epuf[4];
i                 150 src/elemChecker.c           elem[i].nod[8] = epuf[6];
i                 151 src/elemChecker.c           elem[i].nod[9] = epuf[7];
i                  63 src/exo.c        ITG i,j,k,l,m,n,o,indexe,nemax,nlayer,noutloc,iset,iselect,ncomp,nope,
i                  97 src/exo.c          for(i=0;i<*nk;i++){
i                  98 src/exo.c            if(inum[i]==0) continue;
i                  99 src/exo.c            if(inum[i]>0) noutplus++;
i                 100 src/exo.c            if(inum[i]<0) noutmin++;
i                 120 src/exo.c          for(i=0;i<*nk;i++){
i                 121 src/exo.c            if(inum[i]==0){continue;}
i                 124 src/exo.c            node_map[j] = i+1;
i                 125 src/exo.c            node_map_inv[i] = j+1;
i                 126 src/exo.c            x[j]   = co[3*i];
i                 127 src/exo.c            y[j]   = co[3*i+1];
i                 128 src/exo.c            z[j++] = co[3*i+2];
i                 131 src/exo.c          for(i=0;i<*nk;i++){
i                 132 src/exo.c            node_map[j] = i+1;
i                 133 src/exo.c            node_map_inv[i] = j+1;
i                 134 src/exo.c            x[j]   = co[3*i];
i                 135 src/exo.c            y[j]   = co[3*i+1];
i                 136 src/exo.c            z[j++] = co[3*i+2];
i                 170 src/exo.c            for(i=0;i<*ne0;i++){
i                 171 src/exo.c      	if(ipkon[i]<-1){
i                 173 src/exo.c      	}else if(strcmp1(&lakon[8*i],"ESPRNGC")==0){
i                 175 src/exo.c      	}else if(strcmp1(&lakon[8*i],"ESPRNGF")==0){
i                 177 src/exo.c      	}else if(strcmp1(&lakon[8*i],"DCOUP3D")==0){
i                 179 src/exo.c      	}else if(strcmp2(&lakon[8*i+6],"LC",2)==0){
i                 183 src/exo.c      	    if(ielmat[i*mi[2]+k]==0) break;
i                 242 src/exo.c          for(i=0;i<*ne0;i++){ // For each element.  Composite elements are
i                 245 src/exo.c            if(ipkon[i]<=-1){
i                 247 src/exo.c            }else if(strcmp1(&lakon[8*i],"F")==0){
i                 249 src/exo.c            }else if(strcmp1(&lakon[8*i],"ESPRNGC")==0){
i                 251 src/exo.c            }else if(strcmp1(&lakon[8*i],"ESPRNGF")==0){
i                 253 src/exo.c            }else if(strcmp1(&lakon[8*i],"DCOUP3D")==0){
i                 256 src/exo.c      	indexe=ipkon[i];
i                 259 src/exo.c            elem_map[l] = i+1;
i                 261 src/exo.c            strcpy1(curblk,&lakon[8*i],5);
i                 262 src/exo.c            strcpy1(material,&matname[80*(ielmat[i*mi[2]]-1)],5);
i                 263 src/exo.c            if(strcmp1(&lakon[8*i+3],"2")==0){
i                 265 src/exo.c      	if(((strcmp1(&lakon[8*i+6]," ")==0)||
i                 267 src/exo.c      	    (strcmp1(&lakon[8*i+6],"I")==0))&&
i                 268 src/exo.c      	   (strcmp2(&lakon[8*i+6],"LC",2)!=0)){
i                 270 src/exo.c      	}else if(strcmp2(&lakon[8*i+6],"LC",2)==0){
i                 274 src/exo.c      	    if(ielmat[i*mi[2]+k]==0) break;
i                 279 src/exo.c      	    elem_map[l] = i+1;
i                 282 src/exo.c      	}else if(strcmp1(&lakon[8*i+6],"B")==0){
i                 289 src/exo.c            }else if(strcmp1(&lakon[8*i+3],"8")==0){
i                 290 src/exo.c      	if((strcmp1(&lakon[8*i+6]," ")==0)||
i                 294 src/exo.c      	}else if(strcmp1(&lakon[8*i+6],"B")==0){
i                 296 src/exo.c      	  if(strcmp1(&lakon[8*i+4],"R")==0){
i                 298 src/exo.c      	  }else if(strcmp1(&lakon[8*i+4],"I")==0){
i                 303 src/exo.c      	  if(strcmp1(&lakon[8*i+4],"R")==0){
i                 305 src/exo.c      	  }else if(strcmp1(&lakon[8*i+4],"I")==0){
i                 309 src/exo.c            }else if((strcmp1(&lakon[8*i+3],"10")==0)||
i                 310 src/exo.c      	       (strcmp1(&lakon[8*i+3],"14")==0)){
i                 313 src/exo.c            }else if(strcmp1(&lakon[8*i+3],"4")==0){
i                 316 src/exo.c            }else if(strcmp1(&lakon[8*i+3],"15")==0){
i                 317 src/exo.c      	if((strcmp1(&lakon[8*i+6]," ")==0)||
i                 325 src/exo.c            }else if(strcmp1(&lakon[8*i+3],"6")==0){
i                 326 src/exo.c      	if((strcmp1(&lakon[8*i+6]," ")==0)||
i                 336 src/exo.c            }else if(strcmp1(&lakon[8*i],"D")==0){
i                 347 src/exo.c            }else if((strcmp1(&lakon[8*i],"E")==0)&&
i                 348 src/exo.c      	       (strcmp1(&lakon[8*i+6],"A")==0)){
i                 393 src/exo.c            for(i=0;i<*ne0;i++){
i                 394 src/exo.c      	if(ipkon[i]<0) continue;
i                 398 src/exo.c      	      if(ielmat[i*mi[2]+k]==0) break;
i                 415 src/exo.c            for(i=0;i<*ne0;i++){
i                 416 src/exo.c      	if(ipkon[i]<0) continue;
i                 417 src/exo.c      	indexe=ipkon[i];
i                 431 src/exo.c      	      if(ielmat[i*mi[2]+l]==0) break;
i                 925 src/exo.c      	for(i=*ne-1;i>=0;i--){
i                 926 src/exo.c      	  if((strcmp1(&lakon[8*i+1],"S")!=0)||(strcmp1(&lakon[8*i+6],"C")!=0))
i                 929 src/exo.c      	noutloc=*ne-i-1;
i                 932 src/exo.c      	  for(i=*ne-1;i>=0;i--){
i                 933 src/exo.c      	    if((strcmp1(&lakon[8*i+1],"S")!=0)||(strcmp1(&lakon[8*i+6],"C")!=0))
i                 935 src/exo.c      	    strcpy1(text,&lakon[8*i+7],1);
i                 937 src/exo.c      	    nodes=node_map_inv[kon[ipkon[i]+nope-1]-1]-1;
i                 938 src/exo.c      	    nodal_var_vals[nodes]=stx[6*mi[0]*i+j]; 
i                 956 src/exo.c      	for(i=*ne-1;i>=0;i--){
i                 957 src/exo.c      	  if((strcmp1(&lakon[8*i+1],"S")!=0)||(strcmp1(&lakon[8*i+6],"C")!=0))
i                 960 src/exo.c      	noutloc=*ne-i-1;
i                 963 src/exo.c      	for(i=*ne-1;i>=0;i--){
i                 964 src/exo.c      	  if((strcmp1(&lakon[8*i+1],"S")!=0)||(strcmp1(&lakon[8*i+6],"C")!=0))
i                 966 src/exo.c      	  nope=atoi(&lakon[8*i+7])+1;
i                 967 src/exo.c      	  nodes=node_map_inv[kon[ipkon[i]+nope-1]-1]-1;
i                 968 src/exo.c      	  nodal_var_vals[nodes]=ener[i*mi[0]];
i                 994 src/exo.c      	for(i=0;i<*nstate_;i++){
i                 995 src/exo.c      	  ifieldstate[i]=1;icompstate[i]=i;
i                  47 src/exoselect.c   ITG i,j,k,l,m,o,nksegment;
i                  68 src/exoselect.c       for(i=0;i<*nkcoords;i++){
i                  73 src/exoselect.c 	  if(inum[i]<=0) continue;
i                  75 src/exoselect.c 	  if(inum[i]>=0) continue;
i                  77 src/exoselect.c 	  if(inum[i]==0) continue;
i                  82 src/exoselect.c 	  nodal_var_vals[m]=field1[i*nfield[0]+icomp[j]];
i                  84 src/exoselect.c 	  nodal_var_vals[m]=field2[i*nfield[1]+icomp[j]];
i                  95 src/exoselect.c 	    i=ialset[k]+l*nksegment-1;
i                 101 src/exoselect.c 	      if(inum[i]<=0) continue;
i                 103 src/exoselect.c 	      if(inum[i]>=0) continue;
i                 105 src/exoselect.c 	      if(inum[i]==0) continue;
i                 110 src/exoselect.c 	      nodal_var_vals[m]=field1[i*nfield[0]+icomp[j]];
i                 112 src/exoselect.c 	      nodal_var_vals[m]=field2[i*nfield[1]+icomp[j]];
i                 123 src/exoselect.c 	      i=l+m*nksegment-1;
i                 129 src/exoselect.c 		if(inum[i]<=0) continue;
i                 131 src/exoselect.c 		if(inum[i]>=0) continue;
i                 133 src/exoselect.c 		if(inum[i]==0) continue;
i                 138 src/exoselect.c 		nodal_var_vals[o]=field1[i*nfield[0]+icomp[j]];
i                 140 src/exoselect.c 		nodal_var_vals[o]=field2[i*nfield[1]+icomp[j]];
i                  32 src/exosetfind.c   ITG i,j,k,l,n,s,e,gen;
i                  45 src/exosetfind.c   for (i=0; i<*nset; i++){
i                  50 src/exosetfind.c       s=istartset[i]-1;
i                  51 src/exosetfind.c       e=iendset[i]-1;
i                  98 src/exosetfind.c     strncpy(tmpstr,set+i*81,81);
i                 104 src/exosetfind.c 	errr = ex_put_node_set_param (exoid, i, n, 0); // CURRENTLY NO DISTRIBUTIONS ADDED
i                 106 src/exosetfind.c 	errr = ex_put_node_set       (exoid, i, set_nums);
i                 153 src/exosetfind.c     for (i=0; i<*nset; i++){
i                 154 src/exosetfind.c       strncpy(tmpstr,set+i*81,81);
i                  35 src/exovector.c   ITG i,j,k,l,m,n,ii,jj,kk;
i                  63 src/exovector.c 	for(i=0;i<*nkcoords;i++){
i                  64 src/exovector.c 	  if(inum[i]<=0) continue;
i                  65 src/exovector.c 	  nodal_var_vals[m++]=v[(mi[1]+1)*i+j];
i                  69 src/exovector.c 	for(i=0;i<*nkcoords;i++){
i                  70 src/exovector.c 	  if(inum[i]<=0) continue;
i                  71 src/exovector.c 	  if(inotr[2*i]==0){
i                  72 src/exovector.c 	    nodal_var_vals[m++]=v[(mi[1]+1)*i+j];
i                  74 src/exovector.c 	    ii=(mi[1]+1)*i+1;
i                  75 src/exovector.c 	    jj=(mi[1]+1)*i+2;
i                  76 src/exovector.c 	    kk=(mi[1]+1)*i+3;
i                  77 src/exovector.c 	    FORTRAN(transformatrix,(&trab[7*(inotr[2*i]-1)],&co[3*i],a));
i                  88 src/exovector.c 	    i=ialset[k]+l*nksegment-1;
i                  89 src/exovector.c 	    if(inum[i]<=0) continue;
i                  90 src/exovector.c 	    if((*ntrans==0)||(strcmp1(&filabl[5],"G")==0)||(inotr[2*i]==0)){
i                  91 src/exovector.c 	      nodal_var_vals[m++]=v[(mi[1]+1)*i+j];
i                  93 src/exovector.c 	      FORTRAN(transformatrix,(&trab[7*(inotr[2*i]-1)],&co[3*i],a));
i                  94 src/exovector.c 	      ii=(mi[1]+1)*i+1;
i                  95 src/exovector.c 	      jj=(mi[1]+1)*i+2;
i                  96 src/exovector.c 	      kk=(mi[1]+1)*i+3;
i                 107 src/exovector.c 	      i=l+m*nksegment-1;      
i                 108 src/exovector.c 	      if(inum[i]<=0) continue;
i                 109 src/exovector.c 	      if((*ntrans==0)||(strcmp1(&filabl[5],"G")==0)||(inotr[2*i]==0)){
i                 110 src/exovector.c 		nodal_var_vals[n++]=v[(mi[1]+1)*i+j];
i                 112 src/exovector.c 		FORTRAN(transformatrix,(&trab[7*(inotr[2*i]-1)],&co[3*i],a));
i                 113 src/exovector.c 		ii=(mi[1]+1)*i+1;
i                 114 src/exovector.c 		jj=(mi[1]+1)*i+2;
i                 115 src/exovector.c 		kk=(mi[1]+1)*i+3;
i                  71 src/expand.c       ITG *inum=NULL,k,idir,lfin,j,iout=0,index,inode,id,i,idof,im,
i                 116 src/expand.c       for(i=1;i<87*nlabel;i++) filabt[i]=' ';
i                 130 src/expand.c   	for(i=0;i<*nk;i++) inocs[i]=-1;
i                 131 src/expand.c   	for(i=0;i<*ne;i++) ielcs[i]=-1;
i                 136 src/expand.c       for(i=0;i<*mcs;i++){
i                 137 src/expand.c   	is=cs[17*i];
i                 139 src/expand.c   	ielset=cs[17*i+12];
i                 145 src/expand.c   		ielcs[iel]=i;
i                 162 src/expand.c   			ithermal,&csmass[i]));
i                 166 src/expand.c   		    inocs[node]=i;
i                 175 src/expand.c   		    ielcs[iel]=i;
i                 188 src/expand.c   			ithermal,&csmass[i]));
i                 192 src/expand.c   			inocs[node]=i;
i                 210 src/expand.c       for(i=1;i<*nsectors;i++){
i                 212 src/expand.c   	    imddof[i**nmddof+j]=imddof[j]+i*neqh;
i                 215 src/expand.c   	    imdnode[i**nmdnode+j]=imdnode[j]+i**nk;
i                 218 src/expand.c   	    imdboun[i**nmdboun+j]=imdboun[j]+i**nboun;
i                 221 src/expand.c   	    imdmpc[i**nmdmpc+j]=imdmpc[j]+i**nmpc;
i                 247 src/expand.c   	for(i=1;i<is;i++){
i                 249 src/expand.c   	    theta=i*2.*pi/cs[17*jj];
i                 253 src/expand.c   		    co[3*l+i*3**nk]=co[3*l];
i                 254 src/expand.c   		    co[1+3*l+i*3**nk]=co[1+3*l]+theta;
i                 255 src/expand.c   		    co[2+3*l+i*3**nk]=co[2+3*l];
i                 256 src/expand.c   		    if(*ntrans>0) inotr[2*l+i*2**nk]=inotr[2*l];
i                 259 src/expand.c   	    for(l=0;l<*nkon;l++){kon[l+i**nkon]=kon[l]+i**nk;}
i                 263 src/expand.c   			ipkon[l+i**ne]=ipkon[l]+i**nkon;
i                 264 src/expand.c   			ielmat[mi[2]*(l+i**ne)]=ielmat[mi[2]*l];
i                 265 src/expand.c   			if(*norien>0) ielorien[l+i**ne]=ielorien[l];
i                 268 src/expand.c   			    lakon[l2+i*8**ne]=lakon[l2];
i                 271 src/expand.c   			ipkon[l+i**ne]=ipkon[l];
i                 284 src/expand.c       for(i=1;i<*nsectors;i++){
i                 285 src/expand.c   	lint=i*mt**nk;
i                 288 src/expand.c   		nactdof[lint+j]=nactdof[j]+i*neqh;
i                 298 src/expand.c       for(i=1;i<*nsectors;i++){
i                 300 src/expand.c   	    nodeboun[i**nboun+j]=nodeboun[j]+i**nk;
i                 301 src/expand.c   	    ndirboun[i**nboun+j]=ndirboun[j];
i                 302 src/expand.c   	    xboun[i**nboun+j]=xboun[j];
i                 303 src/expand.c   	    xbounold[i**nboun+j]=xbounold[j];
i                 304 src/expand.c   	    if(*nam>0) iamboun[i**nboun+j]=iamboun[j];
i                 305 src/expand.c   	    ikboun[i**nboun+j]=ikboun[j]+8*i**nk;
i                 306 src/expand.c   	    ilboun[i**nboun+j]=ilboun[j]+i**nboun;
i                 312 src/expand.c       for(i=0;i<*nload;i++){
i                 313 src/expand.c   	if(nelemload[2*i+1]<*nsectors){
i                 314 src/expand.c   	    nelemload[2*i]+=*ne*nelemload[2*i+1];
i                 316 src/expand.c   	    nelemload[2*i]+=*ne*(nelemload[2*i+1]-(*nsectors));
i                 318 src/expand.c   	iload=i+1;
i                 354 src/expand.c       i=0;
i                 355 src/expand.c       while(i<*nforc){
i                 356 src/expand.c   	node=nodeforc[2*i];
i                 367 src/expand.c   	    if(nodeforc[2*i+1]<*nsectors){
i                 368 src/expand.c   		nodeforc[2*i]+=*nk*nodeforc[2*i+1];
i                 370 src/expand.c   		nodeforc[2*i]+=*nk*(nodeforc[2*i+1]-(*nsectors));
i                 372 src/expand.c   	    i++;iforc=i;
i                 383 src/expand.c   	    if((ndirforc[i]!=1)||
i                 384 src/expand.c   	       (nodeforc[2*i+2]!=node)||(ndirforc[i+1]!=2)||
i                 385 src/expand.c   	       (nodeforc[2*i+4]!=node)||(ndirforc[i+2]!=3)||
i                 397 src/expand.c   	    fin[0]=xforc[i];
i                 398 src/expand.c   	    fin[1]=xforc[i+1];
i                 399 src/expand.c   	    fin[2]=xforc[i+2];
i                 405 src/expand.c   	    if(nodeforc[2*i+1]<*nsectors){
i                 406 src/expand.c   		nodeforc[2*i]+=*nk*nodeforc[2*i+1];
i                 408 src/expand.c   		nodeforc[2*i]+=*nk*(nodeforc[2*i+1]-(*nsectors));
i                 410 src/expand.c   	    nodeforc[2*i+2]=nodeforc[2*i];
i                 411 src/expand.c   	    nodeforc[2*i+4]=nodeforc[2*i];
i                 416 src/expand.c   	    node=nodeforc[2*i];
i                 422 src/expand.c   	    xforc[i]=fout[0];
i                 423 src/expand.c   	    xforc[i+1]=fout[1];
i                 424 src/expand.c   	    xforc[i+2]=fout[2];
i                 429 src/expand.c   		i++;iforc=i;
i                 469 src/expand.c   	    for(i=0;i<6*mi[0]**ne;i++){eme[i]=0.;}
i                 473 src/expand.c   	    for(i=0;i<*nmpc;i++){
i                 474 src/expand.c   		index=ipompc[i]-1;
i                 485 src/expand.c   		    if(strcmp1(&labmpc[20*i],"CYCLIC")==0){
i                 486 src/expand.c   			icomplex=atoi(&labmpc[20*i+6]);}
i                 487 src/expand.c   		    else if(strcmp1(&labmpc[20*i],"SUBCYCLIC")==0){
i                 565 src/expand.c   	    for(i=1;i<*nsectors;i++){
i                 567 src/expand.c   		theta=i*nm[j]*2.*pi/cs[17*jj];
i                 583 src/expand.c   				    vt[l+mt**nk*i]=v[l];
i                 584 src/expand.c   				    if(*nherm!=1)vti[l+mt**nk*i]=v[l+mt**nk];
i                 591 src/expand.c   			    vt[l+mt**nk*i]=ctl*v[l]-stl*v[l+mt**nk];
i                 592 src/expand.c   			    if(*nherm!=1)vti[l+mt**nk*i]=stl*v[l]+ctl*v[l+mt**nk];
i                 607 src/expand.c   	    i=iznode[ll]-1;
i                 613 src/expand.c   		    idof=nactdof[mt*(k**nk+i)+j1]-1;
i                 619 src/expand.c   				    zdof[(long long)j**nzdof+id-1]=vt[k*mt**nk+mt*i+j1];
i                 621 src/expand.c   				    zdof[(long long)2*j**nzdof+id-1]=vt[k*mt**nk+mt*i+j1];
i                 622 src/expand.c   				    zdof[(long long)(2*j+1)**nzdof+id-1]=vti[k*mt**nk+mt*i+j1];
i                 641 src/expand.c   	for(i=0;i<*mcs;i++){
i                 642 src/expand.c   	  if (nm[j]==0||(nm[j]==(ITG)((cs[17*i]/2))&&(fmod(cs[17*i],2.)==0.))){
i                 643 src/expand.c   	    sum+=cs[17*i]*csmass[i];
i                 645 src/expand.c   	    sum+=cs[17*i]*csmass[i]/2.;
i                 647 src/expand.c   	  summass+=csmass[i];
i                 659 src/expand.c   	    for(i=0;i<*nzdof;i++){zdof[(long long)j**nzdof+i]/=sum;}
i                 661 src/expand.c   	    for(i=0;i<*nzdof;i++){zdof[(long long)(2*j)**nzdof+i]/=sum;}
i                 662 src/expand.c   	    for(i=0;i<*nzdof;i++){zdof[(long long)(2*j+1)**nzdof+i]/=sum;}
i                 675 src/expand.c   	for(i=0;i<*nev;i++){
i                 677 src/expand.c   		xmr[i**nev+j]/=(sumi[i]*sumi[j]);
i                 678 src/expand.c   		xmi[i**nev+j]/=(sumi[i]*sumi[j]);
i                 689 src/expand.c       for(i=0;i<*nsectors;i++){
i                 690 src/expand.c   	if(i==0){
i                 693 src/expand.c   	    ileft=i-1;
i                 702 src/expand.c   	    ikmpc[*nmpc]=ikmpc[j]+8*i**nk;
i                 703 src/expand.c   	    ilmpc[*nmpc]=ilmpc[j]+i**nmpcold;
i                 735 src/expand.c   		nodempc[3**mpcend]=nodeleft+i**nk;
i                 762 src/expand.c   		    nodenew=node+i**nk;
i                 841 src/expand.c   	for(i=1;i<*nsectors;i++){
i                 842 src/expand.c   	    lint=i**nk;
i                 850 src/expand.c   	    for(i=1;i<*nsectors;i++){
i                 851 src/expand.c   		lint=i**nk;
i                 875 src/expand.c         for(i=0; i<*ntie; i++){
i                 876 src/expand.c   	if(tieset[i*(81*3)+80]=='C'){
i                 877 src/expand.c   	  memcpy(tchar2,&tieset[i*(81*3)+81],81);
i                 879 src/expand.c   	  memcpy(tchar3,&tieset[i*(81*3)+81+81],81);
i                  64 src/frd.c        ITG null,one,i,j,k,indexe,nemax,nlayer,noutloc,iset,iselect,ncomp,nope,
i                 127 src/frd.c            for(i=0;i<*nk;i++){
i                 128 src/frd.c      	  if(inum[i]==0) continue;
i                 130 src/frd.c      	  if(inum[i]>0) noutplus++;
i                 131 src/frd.c      	  if(inum[i]<0) noutmin++;
i                 148 src/frd.c          for(i=0;i<20;i++) newdate[i]=' ';
i                 206 src/frd.c          for(i=0;i<*nmat;i++){
i                 207 src/frd.c      	strcpy1(material,&matname[80*i],58);
i                 208 src/frd.c      	fprintf(f1,"%5sUMAT%5" ITGFORMAT "%58s\n",p1,i+1,material);
i                 222 src/frd.c            for(i=0;i<*nk;i++){
i                 223 src/frd.c      	if(inum[i]==0) continue;
i                 225 src/frd.c      	    fprintf(f1,"%3s%10" ITGFORMAT "%12.5E%12.5E%12.5E\n",m1,i+1,(float)co[3*i],
i                 226 src/frd.c      		  (float)co[3*i+1],(float)co[3*i+2]);
i                 228 src/frd.c      	  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 229 src/frd.c      	  fwrite(&co[3*i],sizeof(double),1,f1);
i                 230 src/frd.c      	  fwrite(&co[3*i+1],sizeof(double),1,f1);
i                 231 src/frd.c      	  fwrite(&co[3*i+2],sizeof(double),1,f1);
i                 235 src/frd.c            for(i=0;i<*nk;i++){
i                 237 src/frd.c      	  fprintf(f1,"%3s%10" ITGFORMAT "%12.5E%12.5E%12.5E\n",m1,i+1,(float)co[3*i],
i                 238 src/frd.c      		(float)co[3*i+1],(float)co[3*i+2]);
i                 240 src/frd.c      	  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 241 src/frd.c      	  fwrite(&co[3*i],sizeof(double),1,f1);
i                 242 src/frd.c      	  fwrite(&co[3*i+1],sizeof(double),1,f1);
i                 243 src/frd.c      	  fwrite(&co[3*i+2],sizeof(double),1,f1);
i                 258 src/frd.c      	for(i=0;i<*ne0;i++){
i                 259 src/frd.c      	    if(ipkon[i]<=-1){
i                 261 src/frd.c      	    }else if(strcmp1(&lakon[8*i],"ESPRNGC")==0){
i                 263 src/frd.c      	    }else if(strcmp1(&lakon[8*i],"ESPRNGF")==0){
i                 265 src/frd.c                  }else if(strcmp1(&lakon[8*i],"DCOUP3D")==0){
i                 283 src/frd.c          for(i=0;i<*ne0;i++){
i                 284 src/frd.c            if(ipkon[i]<=-1){
i                 286 src/frd.c            }else if(strcmp1(&lakon[8*i],"F")==0){
i                 288 src/frd.c            }else if(strcmp1(&lakon[8*i],"ESPRNGC")==0){
i                 290 src/frd.c            }else if(strcmp1(&lakon[8*i],"ESPRNGF")==0){
i                 292 src/frd.c            }else if(strcmp1(&lakon[8*i],"DCOUP3D")==0){
i                 295 src/frd.c      	  indexe=ipkon[i];
i                 297 src/frd.c            imat=ielmat[i*mi[2]];
i                 298 src/frd.c            if(strcmp1(&lakon[8*i+3],"2")==0){
i                 302 src/frd.c      	if(((strcmp1(&lakon[8*i+6]," ")==0)||
i                 304 src/frd.c      	    (strcmp1(&lakon[8*i+6],"I")==0))&&
i                 305 src/frd.c                 (strcmp2(&lakon[8*i+6],"LC",2)!=0)){
i                 308 src/frd.c      		    m1,i+1,p4,p0,imat,m2);
i                 317 src/frd.c      	    iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 332 src/frd.c      	}else if(strcmp2(&lakon[8*i+6],"LC",2)==0){
i                 340 src/frd.c      	    if(ielmat[i*mi[2]+k]==0) break;
i                 372 src/frd.c      	}else if(strcmp1(&lakon[8*i+6],"B")==0){
i                 377 src/frd.c      	    fprintf(f1,"%3s%10" ITGFORMAT "%5s%5s%5" ITGFORMAT "\n",m1,i+1,p12,p0,imat);
i                 381 src/frd.c      	    iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 393 src/frd.c      	  if(strcmp1(&lakon[8*i+6],"A")==0) iaxial=1;
i                 396 src/frd.c                          m1,i+1,p10,p0,imat,m2);
i                 400 src/frd.c      	    iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 408 src/frd.c            }else if(strcmp1(&lakon[8*i+3],"8")==0){
i                 409 src/frd.c      	  if((strcmp1(&lakon[8*i+6]," ")==0)||
i                 416 src/frd.c      			  m1,i+1,p1,p0,imat,m2);
i                 420 src/frd.c      		  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 427 src/frd.c      	  }else if(strcmp1(&lakon[8*i+6],"B")==0){
i                 431 src/frd.c      	      if(strcmp1(&lakon[8*i+4],"R")==0){
i                 433 src/frd.c      		      fprintf(f1,"%3s%10" ITGFORMAT "%5s%5s%5" ITGFORMAT "\n",m1,i+1,p11,p0,imat);
i                 437 src/frd.c      		      iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 444 src/frd.c      	      }else if(strcmp1(&lakon[8*i+4],"I")==0){
i                 446 src/frd.c      		      fprintf(f1,"%3s%10" ITGFORMAT "%5s%5s%5" ITGFORMAT "\n",m1,i+1,p11,p0,imat);
i                 450 src/frd.c      		      iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 462 src/frd.c      	      if(strcmp1(&lakon[8*i+6],"A")==0) iaxial=1;
i                 463 src/frd.c      	      if((strcmp1(&lakon[8*i+4],"R")==0)||
i                 464 src/frd.c      		 (strcmp1(&lakon[8*i+4]," ")==0)){
i                 467 src/frd.c      			      m1,i+1,p9,p0,imat,m2);
i                 471 src/frd.c      		      iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 478 src/frd.c      	      }else if(strcmp1(&lakon[8*i+4],"I")==0){
i                 481 src/frd.c      			      m1,i+1,p9,p0,imat,m2);
i                 485 src/frd.c      		      iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 494 src/frd.c            }else if((strcmp1(&lakon[8*i+3],"10")==0)||
i                 495 src/frd.c                     (strcmp1(&lakon[8*i+3],"14")==0)){
i                 501 src/frd.c      		  m1,i+1,p6,p0,imat,m2);
i                 505 src/frd.c      	  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 512 src/frd.c            }else if(strcmp1(&lakon[8*i+3],"4")==0){
i                 518 src/frd.c      		  m1,i+1,p3,p0,imat,m2);
i                 522 src/frd.c      	  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 529 src/frd.c            }else if(strcmp1(&lakon[8*i+3],"15")==0){
i                 530 src/frd.c      	if((strcmp1(&lakon[8*i+6]," ")==0)||
i                 537 src/frd.c      		    m1,i+1,p5,p0,imat,m2);
i                 545 src/frd.c      	    iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 562 src/frd.c      	  if(strcmp1(&lakon[8*i+6],"A")==0) iaxial=1;
i                 565 src/frd.c      		    m1,i+1,p8,p0,imat,m2);
i                 569 src/frd.c      	    iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 577 src/frd.c            }else if(strcmp1(&lakon[8*i+3],"6")==0){
i                 578 src/frd.c      	  if((strcmp1(&lakon[8*i+6]," ")==0)||
i                 585 src/frd.c      			  m1,i+1,p2,p0,imat,m2);
i                 589 src/frd.c      		  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 600 src/frd.c      	      if(strcmp1(&lakon[8*i+6],"A")==0) iaxial=1;
i                 603 src/frd.c      			  m1,i+1,p7,p0,imat,m2);
i                 607 src/frd.c      		  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 615 src/frd.c            }else if((strcmp1(&lakon[8*i],"D")==0)&&(ithermal[1]>1)){
i                 621 src/frd.c      		  fprintf(f1,"%3s%10" ITGFORMAT "%5s%5s%5" ITGFORMAT "\n",m1,i+1,p11,p0,imat);
i                 625 src/frd.c      		  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 637 src/frd.c      		  fprintf(f1,"%3s%10" ITGFORMAT "%5s%5s%5" ITGFORMAT "\n",m1,i+1,p11,p0,imat);
i                 641 src/frd.c      		  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 653 src/frd.c      		  fprintf(f1,"%3s%10" ITGFORMAT "%5s%5s%5" ITGFORMAT "\n",m1,i+1,p12,p0,imat);
i                 657 src/frd.c      		  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 666 src/frd.c            }else if((strcmp1(&lakon[8*i],"E")==0)&&
i                 667 src/frd.c                     (strcmp1(&lakon[8*i+6],"A")==0)){
i                 672 src/frd.c      	  fprintf(f1,"%3s%10" ITGFORMAT "%5s%5s%5" ITGFORMAT "\n",m1,i+1,p11,p0,imat);
i                 675 src/frd.c      	  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                1021 src/frd.c          if((iaxial==1)&&(strcmp1(&filab[352],"I")==0)){for(i=0;i<*nk;i++){fn[1+i*mt]*=180.;fn[2+i*mt]*=180.;fn[3+i*mt]*=180.;}}
i                1024 src/frd.c          if((iaxial==1)&&(strcmp1(&filab[352],"I")==0)){for(i=0;i<*nk;i++){fn[1+i*mt]/=180.;fn[2+i*mt]/=180.;fn[3+i*mt]/=180.;}}
i                1095 src/frd.c          for(i=*ne-1;i>=0;i--){
i                1096 src/frd.c            if((strcmp1(&lakon[8*i+1],"S")!=0)||(strcmp1(&lakon[8*i+6],"C")!=0))
i                1099 src/frd.c          noutloc=*ne-i-1;
i                1112 src/frd.c          for(i=*ne-1;i>=0;i--){
i                1113 src/frd.c            if((strcmp1(&lakon[8*i+1],"S")!=0)||(strcmp1(&lakon[8*i+6],"C")!=0))
i                1115 src/frd.c            strcpy1(text,&lakon[8*i+7],1);
i                1117 src/frd.c            nodes=kon[ipkon[i]+nope-1];
i                1120 src/frd.c      	  for(j=0;j<6;j++)fprintf(f1,"%12.5E",(float)stx[6*mi[0]*i+j]);
i                1124 src/frd.c      	      ifl=(float)stx[6*mi[0]*i+j];
i                1187 src/frd.c          for(i=*ne-1;i>=0;i--){
i                1188 src/frd.c            if((strcmp1(&lakon[8*i+1],"S")!=0)||(strcmp1(&lakon[8*i+6],"C")!=0))
i                1191 src/frd.c          noutloc=*ne-i-1;
i                1199 src/frd.c          for(i=*ne-1;i>=0;i--){
i                1200 src/frd.c            if((strcmp1(&lakon[8*i+1],"S")!=0)||(strcmp1(&lakon[8*i+6],"C")!=0))
i                1202 src/frd.c            nope=atoi(&lakon[8*i+7])+1;
i                1203 src/frd.c            nodes=kon[ipkon[i]+nope-1];
i                1205 src/frd.c      	  fprintf(f1,"%3s%10" ITGFORMAT "%12.5E\n",m1,nodes,(float)ener[i*mi[0]]);
i                1208 src/frd.c      	  ifl=(float)ener[i*mi[0]];
i                1233 src/frd.c          for(i=0;i<*nstate_;i++){
i                1234 src/frd.c            ifieldstate[i]=1;icompstate[i]=i;
i                1530 src/frd.c          if((iaxial==1)&&(strcmp1(&filab[1222],"I")==0)){for(i=0;i<*nk;i++)v[1+i*mt]*=180.;}
i                1534 src/frd.c          if((iaxial==1)&&(strcmp1(&filab[1222],"I")==0)){for(i=0;i<*nk;i++)v[1+i*mt]/=180.;}
i                  40 src/frdcyc.c     ITG nkt,icntrl,*kont=NULL,*ipkont=NULL,*inumt=NULL,*ielmatt=NULL,net,i,l,
i                  66 src/frdcyc.c       for(i=0;i<*nk;i++) inocs[i]=-1;
i                  67 src/frdcyc.c       for(i=0;i<*ne;i++) ielcs[i]=-1;
i                  70 src/frdcyc.c     for(i=0;i<*mcs;i++){
i                  71 src/frdcyc.c       is=cs[17*i+4];
i                  73 src/frdcyc.c       ielset=cs[17*i+12];
i                  79 src/frdcyc.c           ielcs[iel]=i;
i                  89 src/frdcyc.c             inocs[node]=i;
i                  98 src/frdcyc.c             ielcs[iel]=i;
i                 108 src/frdcyc.c               inocs[node]=i;
i                 183 src/frdcyc.c       for(i=1;i<is;i++){
i                 185 src/frdcyc.c         theta=i*2.*pi/cs[17*jj];
i                 189 src/frdcyc.c   	  cot[3*l+i*3**nk]=cot[3*l];
i                 190 src/frdcyc.c   	  cot[1+3*l+i*3**nk]=cot[1+3*l]+theta;
i                 191 src/frdcyc.c   	  cot[2+3*l+i*3**nk]=cot[2+3*l];
i                 198 src/frdcyc.c   		  inotrt[2*l+i*2**nk]=inotrt[2*l];
i                 205 src/frdcyc.c           for(l=0;l<*nkon;l++){kont[l+i**nkon]=kon[l]+i**nk;}
i                 209 src/frdcyc.c                 ipkont[l+i**ne]=ipkon[l]+i**nkon;
i                 210 src/frdcyc.c                 ielmatt[mi[2]*(l+i**ne)]=ielmat[mi[2]*l];
i                 213 src/frdcyc.c                   lakont[l2+i*8**ne]=lakon[l2];
i                 216 src/frdcyc.c               else ipkont[l+i**ne]=-1;
i                 263 src/frdcyc.c       for(i=1;i<is;i++){
i                 265 src/frdcyc.c         for(l=0;l<*nk;l++){inumt[l+i**nk]=inum[l];}
i                 273 src/frdcyc.c                 vt[l+mt**nk*i]=v[l];
i                 281 src/frdcyc.c             if(inocs[l]==jj) t1t[l+*nk*i]=t1[l];
i                 290 src/frdcyc.c                 stnt[l+6**nk*i]=stn[l];
i                 301 src/frdcyc.c                 eent[l+6**nk*i]=een[l];
i                 312 src/frdcyc.c                 fnt[l+mt**nk*i]=fn[l];
i                 320 src/frdcyc.c             if(inocs[l]==jj) epnt[l+*nk*i]=epn[l];
i                 326 src/frdcyc.c             if(inocs[l]==jj) enernt[l+*nk*i]=enern[l];
i                 335 src/frdcyc.c                 xstatent[l+*nstate_**nk*i]=xstaten[l];
i                 346 src/frdcyc.c                 qfnt[l+3**nk*i]=qfn[l];
i                 357 src/frdcyc.c                 emnt[l+6**nk*i]=emn[l];
i                  31 src/frdheader.c   ITG i,ncomma;
i                  37 src/frdheader.c   for(i=10;i<70;i++)text[i]=' ';text[70]='\0';
i                  48 src/frdheader.c     for(i=8;i<70;i++)text[i]=' ';text[70]='\0';
i                  53 src/frdheader.c     for(i=8;i<70;i++)text[i]=' ';text[70]='\0';
i                  58 src/frdheader.c     for(i=9;i<70;i++)text[i]=' ';text[70]='\0';
i                  66 src/frdheader.c       for(i=8;i<24;i++)text[i]=' ';
i                  77 src/frdheader.c     for(i=10;i<70;i++)text[i]=' ';text[70]='\0';
i                  82 src/frdheader.c     for(i=10;i<70;i++)text[i]=' ';text[70]='\0';
i                  93 src/frdheader.c   for(i=0;i<75;i++)text[i]=' ';
i                  43 src/frdselect.c   ITG i,j,k,l,m,n,nksegment;
i                  50 src/frdselect.c     for(i=0;i<*nkcoords;i++){
i                  56 src/frdselect.c 	if(inum[i]<=0) continue;
i                  58 src/frdselect.c 	if(inum[i]>=0) continue;
i                  60 src/frdselect.c 	if(inum[i]==0) continue;
i                  68 src/frdselect.c 	      fprintf(f1,"%3s%10" ITGFORMAT "",m1,i+1);
i                  70 src/frdselect.c 	      iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                  75 src/frdselect.c 		    fprintf(f1,"%12.5E",(float)field1[i*nfield[0]+icomp[j]]);
i                  77 src/frdselect.c 		  ifl=(float)field1[i*nfield[0]+icomp[j]];
i                  82 src/frdselect.c 		  fprintf(f1,"%12.5E",(float)field2[i*nfield[1]+icomp[j]]);
i                  84 src/frdselect.c 		  ifl=(float)field2[i*nfield[1]+icomp[j]];
i                  95 src/frdselect.c 		  fprintf(f1,"%12.5E",(float)field1[i*nfield[0]+icomp[j]]);
i                  97 src/frdselect.c 		  ifl=(float)field1[i*nfield[0]+icomp[j]];
i                 102 src/frdselect.c 		  fprintf(f1,"%12.5E",(float)field2[i*nfield[1]+icomp[j]]);
i                 104 src/frdselect.c 		  ifl=(float)field2[i*nfield[1]+icomp[j]];
i                 119 src/frdselect.c 	  i=ialset[k]+l*nksegment-1;
i                 125 src/frdselect.c 	    if(inum[i]<=0) continue;
i                 127 src/frdselect.c 	    if(inum[i]>=0) continue;
i                 129 src/frdselect.c 	    if(inum[i]==0) continue;
i                 137 src/frdselect.c 		fprintf(f1,"%3s%10" ITGFORMAT "",m1,i+1);
i                 139 src/frdselect.c 		iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 144 src/frdselect.c 		    fprintf(f1,"%12.5E",(float)field1[i*nfield[0]+icomp[j]]);
i                 146 src/frdselect.c 		    ifl=(float)field1[i*nfield[0]+icomp[j]];
i                 151 src/frdselect.c 		    fprintf(f1,"%12.5E",(float)field2[i*nfield[1]+icomp[j]]);
i                 153 src/frdselect.c 		    ifl=(float)field2[i*nfield[1]+icomp[j]];
i                 164 src/frdselect.c 		    fprintf(f1,"%12.5E",(float)field1[i*nfield[0]+j]);
i                 166 src/frdselect.c 		    ifl=(float)field1[i*nfield[0]+icomp[j]];
i                 171 src/frdselect.c 		    fprintf(f1,"%12.5E",(float)field2[i*nfield[1]+j]);
i                 173 src/frdselect.c 		    ifl=(float)field2[i*nfield[1]+icomp[j]];
i                 189 src/frdselect.c 	    i=l+m*nksegment-1;
i                 195 src/frdselect.c 	      if(inum[i]<=0) continue;
i                 197 src/frdselect.c 	      if(inum[i]>=0) continue;
i                 199 src/frdselect.c 	      if(inum[i]==0) continue;
i                 207 src/frdselect.c 		  fprintf(f1,"%3s%10" ITGFORMAT "",m1,i+1);
i                 209 src/frdselect.c 		  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 214 src/frdselect.c 		      fprintf(f1,"%12.5E",(float)field1[i*nfield[0]+icomp[j]]);
i                 216 src/frdselect.c 		      ifl=(float)field1[i*nfield[0]+icomp[j]];
i                 221 src/frdselect.c 		      fprintf(f1,"%12.5E",(float)field2[i*nfield[1]+icomp[j]]);
i                 223 src/frdselect.c 		      ifl=(float)field2[i*nfield[1]+icomp[j]];
i                 234 src/frdselect.c 		      fprintf(f1,"%12.5E",(float)field1[i*nfield[0]+j]);
i                 236 src/frdselect.c 		      ifl=(float)field1[i*nfield[0]+icomp[j]];
i                 241 src/frdselect.c 		      fprintf(f1,"%12.5E",(float)field2[i*nfield[1]+j]);
i                 243 src/frdselect.c 		      ifl=(float)field2[i*nfield[1]+icomp[j]];
i                  29 src/frdvector.c   ITG i,k,l,m,nksegment;
i                  39 src/frdvector.c       for(i=0;i<*nkcoords;i++){
i                  40 src/frdvector.c 	if(inum[i]<=0) continue;
i                  42 src/frdvector.c 	  fprintf(f1,"%3s%10" ITGFORMAT "%12.5E%12.5E%12.5E\n",m1,i+1,
i                  43 src/frdvector.c                   (float)v[(mi[1]+1)*i+1],
i                  44 src/frdvector.c 		  (float)v[(mi[1]+1)*i+2],(float)v[(mi[1]+1)*i+3]);
i                  46 src/frdvector.c 	  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                  47 src/frdvector.c 	  ifl=(float)v[(mi[1]+1)*i+1];fwrite(&ifl,sizeof(float),1,f1);
i                  48 src/frdvector.c 	  ifl=(float)v[(mi[1]+1)*i+2];fwrite(&ifl,sizeof(float),1,f1);
i                  49 src/frdvector.c 	  ifl=(float)v[(mi[1]+1)*i+3];fwrite(&ifl,sizeof(float),1,f1);
i                  53 src/frdvector.c       for(i=0;i<*nkcoords;i++){
i                  54 src/frdvector.c 	if(inum[i]<=0) continue;
i                  55 src/frdvector.c 	if(inotr[2*i]==0){
i                  57 src/frdvector.c 	    fprintf(f1,"%3s%10" ITGFORMAT "%12.5E%12.5E%12.5E\n",m1,i+1,
i                  58 src/frdvector.c                     (float)v[(mi[1]+1)*i+1],
i                  59 src/frdvector.c 		    (float)v[(mi[1]+1)*i+2],(float)v[(mi[1]+1)*i+3]);
i                  61 src/frdvector.c 	    iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                  62 src/frdvector.c 	    ifl=(float)v[(mi[1]+1)*i+1];fwrite(&ifl,sizeof(float),1,f1);
i                  63 src/frdvector.c 	    ifl=(float)v[(mi[1]+1)*i+2];fwrite(&ifl,sizeof(float),1,f1);
i                  64 src/frdvector.c 	    ifl=(float)v[(mi[1]+1)*i+3];fwrite(&ifl,sizeof(float),1,f1);
i                  67 src/frdvector.c 	  FORTRAN(transformatrix,(&trab[7*(inotr[2*i]-1)],&co[3*i],a));
i                  69 src/frdvector.c 	    fprintf(f1,"%3s%10" ITGFORMAT "%12.5E%12.5E%12.5E\n",m1,i+1,
i                  70 src/frdvector.c 		    (float)(v[(mi[1]+1)*i+1]*a[0]+v[(mi[1]+1)*i+2]*a[1]+v[(mi[1]+1)*i+3]*a[2]),
i                  71 src/frdvector.c 		    (float)(v[(mi[1]+1)*i+1]*a[3]+v[(mi[1]+1)*i+2]*a[4]+v[(mi[1]+1)*i+3]*a[5]),
i                  72 src/frdvector.c 		    (float)(v[(mi[1]+1)*i+1]*a[6]+v[(mi[1]+1)*i+2]*a[7]+v[(mi[1]+1)*i+3]*a[8]));
i                  74 src/frdvector.c 	    iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                  75 src/frdvector.c 	    ifl=(float)v[(mi[1]+1)*i+1]*a[0]+v[(mi[1]+1)*i+2]*a[1]+v[(mi[1]+1)*i+3]*a[2];
i                  77 src/frdvector.c 	    ifl=(float)v[(mi[1]+1)*i+1]*a[3]+v[(mi[1]+1)*i+2]*a[4]+v[(mi[1]+1)*i+3]*a[5];
i                  79 src/frdvector.c 	    ifl=(float)v[(mi[1]+1)*i+1]*a[6]+v[(mi[1]+1)*i+2]*a[7]+v[(mi[1]+1)*i+3]*a[8];
i                  90 src/frdvector.c 	  i=ialset[k]+l*nksegment-1;
i                  91 src/frdvector.c 	  if(inum[i]<=0) continue;
i                  92 src/frdvector.c 	  if((*ntrans==0)||(strcmp1(&filabl[5],"G")==0)||(inotr[2*i]==0)){
i                  94 src/frdvector.c 		  fprintf(f1,"%3s%10" ITGFORMAT "%12.5E%12.5E%12.5E\n",m1,i+1,(float)v[(mi[1]+1)*i+1],
i                  95 src/frdvector.c 			  (float)v[(mi[1]+1)*i+2],(float)v[(mi[1]+1)*i+3]);
i                  97 src/frdvector.c 		  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                  98 src/frdvector.c 		  ifl=(float)v[(mi[1]+1)*i+1];fwrite(&ifl,sizeof(float),1,f1);
i                  99 src/frdvector.c 		  ifl=(float)v[(mi[1]+1)*i+2];fwrite(&ifl,sizeof(float),1,f1);
i                 100 src/frdvector.c 		  ifl=(float)v[(mi[1]+1)*i+3];fwrite(&ifl,sizeof(float),1,f1);
i                 103 src/frdvector.c 	      FORTRAN(transformatrix,(&trab[7*(inotr[2*i]-1)],&co[3*i],a));
i                 105 src/frdvector.c 		  fprintf(f1,"%3s%10" ITGFORMAT "%12.5E%12.5E%12.5E\n",m1,i+1,   
i                 106 src/frdvector.c 			  (float)(v[(mi[1]+1)*i+1]*a[0]+v[(mi[1]+1)*i+2]*a[1]+v[(mi[1]+1)*i+3]*a[2]),
i                 107 src/frdvector.c 			  (float)(v[(mi[1]+1)*i+1]*a[3]+v[(mi[1]+1)*i+2]*a[4]+v[(mi[1]+1)*i+3]*a[5]),
i                 108 src/frdvector.c 			  (float)(v[(mi[1]+1)*i+1]*a[6]+v[(mi[1]+1)*i+2]*a[7]+v[(mi[1]+1)*i+3]*a[8]));
i                 110 src/frdvector.c 		  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 111 src/frdvector.c 		  ifl=(float)v[(mi[1]+1)*i+1]*a[0]+v[(mi[1]+1)*i+2]*a[1]+v[(mi[1]+1)*i+3]*a[2];
i                 113 src/frdvector.c 		  ifl=(float)v[(mi[1]+1)*i+1]*a[3]+v[(mi[1]+1)*i+2]*a[4]+v[(mi[1]+1)*i+3]*a[5];
i                 115 src/frdvector.c 		  ifl=(float)v[(mi[1]+1)*i+1]*a[6]+v[(mi[1]+1)*i+2]*a[7]+v[(mi[1]+1)*i+3]*a[8];
i                 126 src/frdvector.c 	    i=l+m*nksegment-1;
i                 127 src/frdvector.c 	    if(inum[i]<=0) continue;
i                 128 src/frdvector.c 	    if((*ntrans==0)||(strcmp1(&filabl[5],"G")==0)||(inotr[2*i]==0)){
i                 130 src/frdvector.c 		    fprintf(f1,"%3s%10" ITGFORMAT "%12.5E%12.5E%12.5E\n",m1,i+1,(float)v[(mi[1]+1)*i+1],
i                 131 src/frdvector.c 			    (float)v[(mi[1]+1)*i+2],(float)v[(mi[1]+1)*i+3]);
i                 133 src/frdvector.c 		    iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 134 src/frdvector.c 		    ifl=(float)v[(mi[1]+1)*i+1];fwrite(&ifl,sizeof(float),1,f1);
i                 135 src/frdvector.c 		    ifl=(float)v[(mi[1]+1)*i+2];fwrite(&ifl,sizeof(float),1,f1);
i                 136 src/frdvector.c 		    ifl=(float)v[(mi[1]+1)*i+3];fwrite(&ifl,sizeof(float),1,f1);
i                 139 src/frdvector.c 	      FORTRAN(transformatrix,(&trab[7*(inotr[2*i]-1)],&co[3*i],a));
i                 141 src/frdvector.c 		  fprintf(f1,"%3s%10" ITGFORMAT "%12.5E%12.5E%12.5E\n",m1,i+1,   
i                 142 src/frdvector.c 			  (float)(v[(mi[1]+1)*i+1]*a[0]+v[(mi[1]+1)*i+2]*a[1]+
i                 143 src/frdvector.c 				  v[(mi[1]+1)*i+3]*a[2]),
i                 144 src/frdvector.c 			  (float)(v[(mi[1]+1)*i+1]*a[3]+v[(mi[1]+1)*i+2]*a[4]+
i                 145 src/frdvector.c 				  v[(mi[1]+1)*i+3]*a[5]),
i                 146 src/frdvector.c 			  (float)(v[(mi[1]+1)*i+1]*a[6]+v[(mi[1]+1)*i+2]*a[7]+
i                 147 src/frdvector.c 				  v[(mi[1]+1)*i+3]*a[8]));
i                 149 src/frdvector.c 		  iw=(int)(i+1);fwrite(&iw,sizeof(int),1,f1);
i                 150 src/frdvector.c 		  ifl=(float)v[(mi[1]+1)*i+1]*a[0]+v[(mi[1]+1)*i+2]*a[1]+v[(mi[1]+1)*i+3]*a[2];
i                 152 src/frdvector.c 		  ifl=(float)v[(mi[1]+1)*i+1]*a[3]+v[(mi[1]+1)*i+2]*a[4]+v[(mi[1]+1)*i+3]*a[5];
i                 154 src/frdvector.c 		  ifl=(float)v[(mi[1]+1)*i+1]*a[6]+v[(mi[1]+1)*i+2]*a[7]+v[(mi[1]+1)*i+3]*a[8];
i                  19 src/frecord.c    register ITG i, n, c;
i                  21 src/frecord.c    for (i=0; i<MAX_LINE_LENGTH-1; i++)
i                  23 src/frecord.c      string[i] = getc(handle1);
i                  24 src/frecord.c      if (string[i] == '\n')
i                  26 src/frecord.c        for (n=i+1; n<MAX_LINE_LENGTH; ++n) string[n] = '\0';
i                  27 src/frecord.c        return(i);
i                  29 src/frecord.c      if (string[i] == '\r')
i                  35 src/frecord.c        for (n=i+1; n<MAX_LINE_LENGTH; ++n) string[n] = '\0';
i                  36 src/frecord.c        return(i);
i                  38 src/frecord.c      else if (string[i] == (char)EOF)
i                  40 src/frecord.c        for (n=i+1; n<MAX_LINE_LENGTH; ++n) string[n] = '\0';
i                  41 src/frecord.c        return(i);
i                  45 src/getglobalresults.c 	indexe,istep,loadcase,nfaces,netet_,nktet=0,nfield,j,nodes[4],i,
i                  82 src/getglobalresults.c     for(i=0;i<MAX_LINE_LENGTH;i++){
i                  83 src/getglobalresults.c 	if(strcmp1(&datin[i]," ")==0){
i                  84 src/getglobalresults.c 	    datin[i]='\0';
i                  94 src/getglobalresults.c     for(i=0;i<*nboun;i++){
i                  95 src/getglobalresults.c 	if((xboun[i]<1.9232931375)&&(xboun[i]>1.9232931373)){
i                  96 src/getglobalresults.c 	    istep=iamboun[i];
i                 101 src/getglobalresults.c 	for(i=0;i<*nforc;i++){
i                 102 src/getglobalresults.c 	    if((xforc[i]<1.9232931375)&&(xforc[i]>1.9232931373)){
i                 103 src/getglobalresults.c 		istep=iamforc[i];
i                 109 src/getglobalresults.c 	for(i=0;i<*nload;i++){
i                 110 src/getglobalresults.c 	    if(strcmp1(&sideload[20*i+2],"SM")==0){
i                 111 src/getglobalresults.c 		istep=iamload[2*i];
i                 117 src/getglobalresults.c 	for(i=0;i<*nk;i++){
i                 118 src/getglobalresults.c 	    if((t1[i]<1.9232931375)&&(t1[i]>1.9232931373)){
i                 119 src/getglobalresults.c 		istep=iamt1[i];
i                 156 src/getglobalresults.c     for(i=0;i<anz[0].n;i++){
i                 157 src/getglobalresults.c 	if(node[i].nr>nktet) nktet=node[i].nr;
i                 163 src/getglobalresults.c     for (i=0;i<anz[0].n;i++){
i                 164 src/getglobalresults.c 	nodenr=node[i].nr;
i                 183 src/getglobalresults.c     for(i=0;i<anz[0].e;i++){
i                 184 src/getglobalresults.c 	ielemnr[i]=elem[i].nr;
i                 185 src/getglobalresults.c 	kontyp[i]=elem[i].type;
i                 186 src/getglobalresults.c 	ipkon[i]=indexe;
i                 187 src/getglobalresults.c 	if(kontyp[i]==1){
i                 189 src/getglobalresults.c 	}else if(kontyp[i]==2){
i                 191 src/getglobalresults.c 	}else if(kontyp[i]==3){
i                 193 src/getglobalresults.c 	}else if(kontyp[i]==4){
i                 195 src/getglobalresults.c 	}else if(kontyp[i]==5){
i                 197 src/getglobalresults.c 	}else if(kontyp[i]==6){
i                 201 src/getglobalresults.c 	    printf("         mesh not recognized; cgx element type=%" ITGFORMAT "\n",kontyp[i]);
i                 205 src/getglobalresults.c 	    kon[indexe++]=elem[i].nod[j];
i                 227 src/getglobalresults.c     for(i=0;i<ne;i++){
i                 228 src/getglobalresults.c 	type=kontyp[i];
i                 229 src/getglobalresults.c 	indexe=ipkon[i]-1;
i                 239 src/getglobalresults.c 		iparent[netet]=i+1;
i                 254 src/getglobalresults.c 		iparent[netet]=i+1;
i                 268 src/getglobalresults.c 	    iparent[netet]=i+1;
i                 282 src/getglobalresults.c 		iparent[netet]=i+1;
i                 297 src/getglobalresults.c 		iparent[netet]=i+1;
i                 312 src/getglobalresults.c 		iparent[netet]=i+1;
i                 334 src/getglobalresults.c     for(i=0;i<netet;i++){
i                 335 src/getglobalresults.c 	n1=kontet[4*i]-1;
i                 336 src/getglobalresults.c 	n2=kontet[4*i+1]-1;
i                 337 src/getglobalresults.c 	n3=kontet[4*i+2]-1;
i                 338 src/getglobalresults.c 	n4=kontet[4*i+3]-1;
i                 339 src/getglobalresults.c 	cgtet[3*i]=(cotet[3*n1]+cotet[3*n2]+cotet[3*n3]+cotet[3*n4])/4.;
i                 340 src/getglobalresults.c 	cgtet[3*i+1]=(cotet[3*n1+1]+cotet[3*n2+1]+cotet[3*n3+1]+cotet[3*n4+1])/4.;
i                 341 src/getglobalresults.c 	cgtet[3*i+2]=(cotet[3*n1+2]+cotet[3*n2+2]+cotet[3*n3+2]+cotet[3*n4+2])/4.;
i                 355 src/getglobalresults.c     for(i=0;i<netet;i++){
i                 356 src/getglobalresults.c 	nnx[i]=i+1;
i                 357 src/getglobalresults.c 	nny[i]=i+1;
i                 358 src/getglobalresults.c 	nnz[i]=i+1;
i                 359 src/getglobalresults.c 	x[i]=cgtet[3*i];
i                 360 src/getglobalresults.c 	y[i]=cgtet[3*i+1];
i                 361 src/getglobalresults.c 	z[i]=cgtet[3*i+2];
i                 362 src/getglobalresults.c 	xo[i]=x[i];
i                 363 src/getglobalresults.c 	yo[i]=y[i];
i                 364 src/getglobalresults.c 	zo[i]=z[i];
i                 380 src/getglobalresults.c     for(i=0;i<anz[0].l;i++){
i                 381 src/getglobalresults.c 	for(j=0;j<lcase[i].npheader;j++){
i                 382 src/getglobalresults.c 	    if(strcmp1(&lcase[i].pheader[j][5],"PSTEP")==0){
i                 383 src/getglobalresults.c 		strcpy1(text,&lcase[i].pheader[j][48],12);
i                 389 src/getglobalresults.c 	   (strcmp1(lcase[i].name,"NDTEMP")==0)){
i                 390 src/getglobalresults.c 	    loadcase=i;
i                 402 src/getglobalresults.c 	    printf("ERROR in getglobalresults: Could not read data for Dataset:%" ITGFORMAT "\n", i+1); 
i                 408 src/getglobalresults.c 	for(i=0;i<anz[0].n;i++){
i                 409 src/getglobalresults.c 	    nodenr=node[i].nr;
i                 418 src/getglobalresults.c     for(i=0;i<anz[0].l;i++){
i                 419 src/getglobalresults.c 	for(j=0;j<lcase[i].npheader;j++){
i                 420 src/getglobalresults.c 	    if(strcmp1(&lcase[i].pheader[j][5],"PSTEP")==0){
i                 421 src/getglobalresults.c 		strcpy1(text,&lcase[i].pheader[j][48],12);
i                 428 src/getglobalresults.c 	   (strcmp1(lcase[i].name,"DISP")==0)){
i                 429 src/getglobalresults.c 	    loadcase=i;
i                 441 src/getglobalresults.c 	    printf("ERROR in getglobalresults: Could not read data for Dataset:%" ITGFORMAT "\n", i+1); 
i                 447 src/getglobalresults.c 	for(i=0;i<anz[0].n;i++){
i                 448 src/getglobalresults.c 	    nodenr=node[i].nr;
i                 459 src/getglobalresults.c     for(i=0;i<anz[0].l;i++){
i                 460 src/getglobalresults.c 	for(j=0;j<lcase[i].npheader;j++){
i                 461 src/getglobalresults.c 	    if(strcmp1(&lcase[i].pheader[j][5],"PSTEP")==0){
i                 462 src/getglobalresults.c 		strcpy1(text,&lcase[i].pheader[j][48],12);
i                 468 src/getglobalresults.c 	   (strcmp1(lcase[i].name,"STRESS")==0)){
i                 469 src/getglobalresults.c 	    loadcase=i;
i                 481 src/getglobalresults.c 	    printf("ERROR in getglobalresults: Could not read data for Dataset:%" ITGFORMAT "\n", i+1); 
i                 487 src/getglobalresults.c 	for(i=0;i<anz[0].n;i++){
i                 488 src/getglobalresults.c 	    nodenr=node[i].nr;
i                 502 src/getglobalresults.c     for(i=0;i<anz[0].l;i++){
i                 503 src/getglobalresults.c 	for(j=0;j<lcase[i].npheader;j++){
i                 504 src/getglobalresults.c 	    if(strcmp1(&lcase[i].pheader[j][5],"PSTEP")==0){
i                 505 src/getglobalresults.c 		strcpy1(text,&lcase[i].pheader[j][48],12);
i                 511 src/getglobalresults.c 	   (strcmp1(lcase[i].name,"FORC")==0)){
i                 512 src/getglobalresults.c 	    loadcase=i;
i                 523 src/getglobalresults.c 	    printf("ERROR in getglobalresults: Could not read data for Dataset:%" ITGFORMAT "\n", i+1); 
i                 529 src/getglobalresults.c 	for(i=0;i<anz[0].n;i++){
i                 530 src/getglobalresults.c 	    nodenr=node[i].nr;
i                  45 src/inicont.c        *ilmpc=NULL,*nodempc=NULL,nmpc_,i,j,k,ncone;
i                  72 src/linstatic.c   ITG *inum=NULL,k,*icol=NULL,*irow=NULL,ielas,icmd=0,iinc=1,nasym=0,i,j,ic,ir,
i                 267 src/linstatic.c     for(i=0;i<*nboun;i++){
i                 268 src/linstatic.c 	if(strcmp1(&typeboun[i],"C")==0){
i                 269 src/linstatic.c 	    iretain[nretain]=i+1;
i                 270 src/linstatic.c 	    noderetain[nretain]=nodeboun[i];
i                 271 src/linstatic.c 	    ndirretain[nretain]=ndirboun[i];
i                 284 src/linstatic.c     for(i=0;i<nretain;i++){
i                 286 src/linstatic.c 	ic=*neq+iretain[i]-1;
i                 327 src/linstatic.c 	  xbounact[iretain[i]-1]=1.;
i                 346 src/linstatic.c 	  xbounact[iretain[i]-1]=0.;
i                 360 src/linstatic.c 	      submatrix[i*nretain+j]=fn[mt*(noderetain[j]-1)+ndirretain[j]];
i                  41 src/mastruct.c   ITG i,j,k,l,jj,ll,id,index,jdof1,jdof2,idof1,idof2,mpc1,mpc2,id1,id2,
i                  57 src/mastruct.c   for(i=0;i<mt**nk;++i){nactdof[i]=0;}
i                  61 src/mastruct.c       for(i=0;i<*ne;++i){
i                  63 src/mastruct.c 	  if(ipkon[i]<0) continue;
i                  64 src/mastruct.c 	  if(strcmp1(&lakon[8*i],"F")==0)continue;
i                  65 src/mastruct.c 	  indexe=ipkon[i];
i                  67 src/mastruct.c           if (strcmp1(&lakon[8*i+3],"8I")==0)nope=11;
i                  68 src/mastruct.c 	  else if(strcmp1(&lakon[8*i+3],"20")==0)nope=20;
i                  70 src/mastruct.c 	  else if(strcmp1(&lakon[8*i+3],"2")==0)nope=26;
i                  71 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"8")==0)nope=8;
i                  72 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"10")==0)nope=10;
i                  73 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"14")==0)nope=14;
i                  74 src/mastruct.c 	  else if ((strcmp1(&lakon[8*i+3],"4")==0)||
i                  75 src/mastruct.c                    (strcmp1(&lakon[8*i+2],"4")==0)) nope=4;
i                  76 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"15")==0)nope=15;
i                  77 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"6")==0)nope=6;
i                  78 src/mastruct.c 	  else if (strcmp1(&lakon[8*i],"E")==0){
i                  79 src/mastruct.c 	      if((strcmp1(&lakon[8*i+6],"C")==0)&&(*mortar==1)){
i                  80 src/mastruct.c 		  nope=kon[ipkon[i]-1];
i                  82 src/mastruct.c 		  lakonl[0]=lakon[8*i+7];
i                 101 src/mastruct.c       for(i=0;i<*ne;++i){
i                 103 src/mastruct.c 	  if(ipkon[i]<0) continue;
i                 104 src/mastruct.c 	  if(strcmp1(&lakon[8*i],"F")==0)continue;
i                 105 src/mastruct.c 	  indexe=ipkon[i];
i                 106 src/mastruct.c 	  if(strcmp1(&lakon[8*i+3],"20")==0)nope=20;
i                 107 src/mastruct.c 	  else if(strcmp1(&lakon[8*i+3],"2")==0)nope=26;
i                 108 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"8")==0)nope=8;
i                 109 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"10")==0)nope=10;
i                 110 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"14")==0)nope=14;
i                 111 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"4")==0)nope=4;
i                 112 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"15")==0)nope=15;
i                 113 src/mastruct.c 	  else if (strcmp1(&lakon[8*i+3],"6")==0)nope=6;
i                 114 src/mastruct.c 	  else if (strcmp1(&lakon[8*i],"E")==0){
i                 115 src/mastruct.c 	      if((strcmp1(&lakon[8*i+6],"C")==0)&&(*mortar==1)){
i                 116 src/mastruct.c 		  nope=kon[ipkon[i]-1];
i                 118 src/mastruct.c 		  lakonl[0]=lakon[8*i+7];
i                 121 src/mastruct.c 	  }else if (strcmp1(&lakon[8*i],"D ")==0){
i                 145 src/mastruct.c   for(i=0;i<*nmpc;++i){
i                 146 src/mastruct.c       if (strcmp1(&labmpc[20*i],"FLUID")==0) continue;
i                 147 src/mastruct.c       index=ipompc[i]-1;
i                 160 src/mastruct.c   for(i=0;i<*nboun;++i){
i                 161 src/mastruct.c       if(ndirboun[i]>mi[1]) continue;
i                 162 src/mastruct.c       if (strcmp1(&typeboun[i],"F")==0) continue;
i                 163 src/mastruct.c       nactdof[mt*(nodeboun[i]-1)+ndirboun[i]]=0;
i                 166 src/mastruct.c   for(i=0;i<*nmpc;++i){
i                 167 src/mastruct.c       if (strcmp1(&labmpc[20*i],"FLUID")==0) continue;
i                 168 src/mastruct.c       index=ipompc[i]-1;
i                 176 src/mastruct.c   for(i=0;i<*nk;++i){
i                 178 src/mastruct.c 	if(nactdof[mt*i+j]!=0){
i                 181 src/mastruct.c           nactdof[mt*i+j]=neq[0];
i                 184 src/mastruct.c 	    nactdof[mt*i+j]=0;
i                 190 src/mastruct.c   for(i=0;i<*nk;++i){
i                 191 src/mastruct.c       if(nactdof[mt*i]!=0){
i                 194 src/mastruct.c         nactdof[mt*i]=neq[1];
i                 197 src/mastruct.c 	  nactdof[mt*i]=0;
i                 215 src/mastruct.c     for(i=0;i<4**nk;++i){ipointer[i]=0;}
i                 221 src/mastruct.c     for(i=0;i<*ne;++i){
i                 223 src/mastruct.c       if(ipkon[i]<0) continue;
i                 224 src/mastruct.c       if(strcmp1(&lakon[8*i],"F")==0)continue;
i                 225 src/mastruct.c       indexe=ipkon[i];
i                 227 src/mastruct.c       if (strcmp1(&lakon[8*i+3],"8I")==0)nope=11;
i                 228 src/mastruct.c       else if(strcmp1(&lakon[8*i+3],"20")==0)nope=20;
i                 230 src/mastruct.c       else if(strcmp1(&lakon[8*i+3],"2")==0)nope=26;
i                 231 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"8")==0)nope=8;
i                 232 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"10")==0)nope=10;
i                 233 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"14")==0)nope=14;
i                 234 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"4")==0)nope=4;
i                 235 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"15")==0)nope=15;
i                 236 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"6")==0)nope=6;
i                 237 src/mastruct.c       else if (strcmp1(&lakon[8*i],"E")==0){
i                 238 src/mastruct.c 	  if((strcmp1(&lakon[8*i+6],"C")==0)&&(*mortar==1)){
i                 239 src/mastruct.c 	      nope=kon[ipkon[i]-1];
i                 241 src/mastruct.c 	      lakonl[0]=lakon[8*i+7];
i                 387 src/mastruct.c     for(i=0;i<*ne;++i){
i                 389 src/mastruct.c       if(ipkon[i]<0) continue;
i                 390 src/mastruct.c       if(strcmp1(&lakon[8*i],"F")==0)continue;
i                 391 src/mastruct.c       indexe=ipkon[i];
i                 392 src/mastruct.c       if(strcmp1(&lakon[8*i+3],"20")==0)nope=20;
i                 393 src/mastruct.c       else if(strcmp1(&lakon[8*i+3],"2")==0)nope=26;
i                 394 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"8")==0)nope=8;
i                 395 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"10")==0)nope=10;
i                 396 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"14")==0)nope=14;
i                 397 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"4")==0)nope=4;
i                 398 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"15")==0)nope=15;
i                 399 src/mastruct.c       else if (strcmp1(&lakon[8*i+3],"6")==0)nope=6;
i                 400 src/mastruct.c       else if (strcmp1(&lakon[8*i],"E")==0){
i                 401 src/mastruct.c 	  if((strcmp1(&lakon[8*i+6],"C")==0)&&(*mortar==1)){
i                 402 src/mastruct.c 	      nope=kon[ipkon[i]-1];
i                 404 src/mastruct.c 	      lakonl[0]=lakon[8*i+7];
i                 407 src/mastruct.c       }else if (strcmp1(&lakon[8*i],"D ")==0){
i                 555 src/mastruct.c     for(i=0;i<neq[2];++i){
i                 556 src/mastruct.c 	if(ipointer[i]==0){
i                 557 src/mastruct.c 	    if(i>=neq[1]) continue;
i                 561 src/mastruct.c 		    if(nactdof[mt*j+k]==i+1){
i                 573 src/mastruct.c 	istart=ipointer[i];
i                 577 src/mastruct.c 	    irow[istartold-1]=i+1;
i                 630 src/mastruct.c     for(i=0;i<neq[1];++i){icol[i]=0;}
i                 632 src/mastruct.c     for(i=0;i<nmast;++i){
i                 633 src/mastruct.c       if(mast1[i]==irow[i]){++isubtract;}
i                 635 src/mastruct.c 	mast1[i-isubtract]=mast1[i];
i                 636 src/mastruct.c 	irow[i-isubtract]=irow[i];
i                 637 src/mastruct.c 	if(k!=mast1[i]){
i                 638 src/mastruct.c 	  for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
i                 639 src/mastruct.c 	  k=mast1[i];
i                 649 src/mastruct.c     for(i=0;i<neq[1];++i){
i                 650 src/mastruct.c       if(jq[i+1]-jq[i]>0){
i                 651 src/mastruct.c 	isize=jq[i+1]-jq[i];
i                 652 src/mastruct.c 	FORTRAN(isortii,(&irow[jq[i]-1],&mast1[jq[i]-1],&isize,&kflag));
i                 664 src/mastruct.c 	for(i=neq[1];i<neq[2];++i){icol[i]=0;}
i                 665 src/mastruct.c 	for(i=nmast+isubtract;i<nmastboun;++i){
i                 669 src/mastruct.c 	    irow[i-isubtract]=mast1[i];
i                 670 src/mastruct.c 	    mast1[i-isubtract]=irow[i];
i                 674 src/mastruct.c 	    if(k!=irow[i]){
i                 675 src/mastruct.c 		for(l=k;l<irow[i];++l){jq[l]=i+1-isubtract;}
i                 676 src/mastruct.c 		k=irow[i];
i                 682 src/mastruct.c 	for(i=neq[1];i<neq[2];++i){
i                 683 src/mastruct.c 	    if(jq[i+1]-jq[i]>0){
i                 684 src/mastruct.c 		isize=jq[i+1]-jq[i];
i                 685 src/mastruct.c 		FORTRAN(isortii,(&irow[jq[i]-1],&mast1[jq[i]-1],&isize,&kflag));
i                  41 src/mastructcs.c   ITG i,j,k,l,jj,ll,id,index,jdof1,jdof2,idof1,idof2,mpc1,mpc2,id1,id2,
i                  58 src/mastructcs.c   for(i=0;i<mt**nk;++i){nactdof[i]=0;}
i                  62 src/mastructcs.c   for(i=0;i<*ne;++i){
i                  64 src/mastructcs.c     if(ipkon[i]<0) continue;
i                  65 src/mastructcs.c     indexe=ipkon[i];
i                  67 src/mastructcs.c     if (strcmp1(&lakon[8*i+3],"8I")==0)nope=11;
i                  68 src/mastructcs.c     else if(strcmp1(&lakon[8*i+3],"20")==0)nope=20;
i                  70 src/mastructcs.c     else if(strcmp1(&lakon[8*i+3],"2")==0)nope=26;
i                  71 src/mastructcs.c     else if (strcmp1(&lakon[8*i+3],"8")==0)nope=8;
i                  72 src/mastructcs.c     else if (strcmp1(&lakon[8*i+3],"10")==0)nope=10;
i                  73 src/mastructcs.c     else if ((strcmp1(&lakon[8*i+3],"4")==0)||
i                  74 src/mastructcs.c 	     (strcmp1(&lakon[8*i+2],"4")==0)) nope=4;
i                  75 src/mastructcs.c     else if (strcmp1(&lakon[8*i+3],"15")==0)nope=15;
i                  76 src/mastructcs.c     else if (strcmp1(&lakon[8*i+3],"6")==0)nope=6;
i                  77 src/mastructcs.c     else if (strcmp1(&lakon[8*i],"E")==0){
i                  78 src/mastructcs.c 	if((strcmp1(&lakon[8*i+6],"C")==0)&&(*mortar==1)){
i                  79 src/mastructcs.c 	    nope=kon[ipkon[i]-1];
i                  81 src/mastructcs.c 	    lakonl[0]=lakon[8*i+7];
i                 101 src/mastructcs.c   for(i=0;i<*nmpc;++i){
i                 102 src/mastructcs.c       index=ipompc[i]-1;
i                 114 src/mastructcs.c   for(i=0;i<*nboun;++i){
i                 115 src/mastructcs.c       if(ndirboun[i]>mi[1]) continue;
i                 116 src/mastructcs.c       nactdof[mt*(nodeboun[i]-1)+ndirboun[i]]=0;
i                 119 src/mastructcs.c   for(i=0;i<*nmpc;++i){
i                 120 src/mastructcs.c       index=ipompc[i]-1;
i                 128 src/mastructcs.c   for(i=0;i<*nk;++i){
i                 130 src/mastructcs.c 	if(nactdof[mt*i+j]!=0){
i                 132 src/mastructcs.c 	nactdof[mt*i+j]=neq[0];
i                 145 src/mastructcs.c   for(i=0;i<6**nk;++i){ipointer[i]=0;}
i                 147 src/mastructcs.c   for(i=0;i<*ne;++i){
i                 149 src/mastructcs.c     if(ipkon[i]<0) continue;
i                 150 src/mastructcs.c     indexe=ipkon[i];
i                 152 src/mastructcs.c     if(strcmp1(&lakon[8*i],"C3D8I")==0){nope=11;}
i                 153 src/mastructcs.c     else if(strcmp1(&lakon[8*i+3],"20")==0)nope=20;
i                 155 src/mastructcs.c     else if(strcmp1(&lakon[8*i+3],"2")==0)nope=26;
i                 156 src/mastructcs.c     else if (strcmp1(&lakon[8*i+3],"8")==0)nope=8;
i                 157 src/mastructcs.c     else if (strcmp1(&lakon[8*i+3],"10")==0)nope=10;
i                 158 src/mastructcs.c     else if (strcmp1(&lakon[8*i+3],"4")==0)nope=4;
i                 159 src/mastructcs.c     else if (strcmp1(&lakon[8*i+3],"15")==0)nope=15;
i                 160 src/mastructcs.c     else if (strcmp1(&lakon[8*i+3],"6")==0)nope=6;
i                 161 src/mastructcs.c     else if (strcmp1(&lakon[8*i],"E")==0){
i                 162 src/mastructcs.c 	if((strcmp1(&lakon[8*i+6],"C")==0)&&(*mortar==1)){
i                 163 src/mastructcs.c 	    nope=kon[ipkon[i]-1];
i                 165 src/mastructcs.c 	    lakonl[0]=lakon[8*i+7];
i                 457 src/mastructcs.c   for(i=0;i<neq[0];++i){
i                 458 src/mastructcs.c       if(ipointer[i]==0){
i                 459 src/mastructcs.c 	  if(i>=neq[1]) continue;
i                 463 src/mastructcs.c       istart=ipointer[i];
i                 467 src/mastructcs.c 	  irow[istartold-1]=i+1;
i                 502 src/mastructcs.c   for(i=0;i<neq[0];++i){icol[i]=0;}
i                 504 src/mastructcs.c   for(i=0;i<nmast;++i){
i                 505 src/mastructcs.c     if(mast1[i]==irow[i]){++isubtract;}
i                 507 src/mastructcs.c       mast1[i-isubtract]=mast1[i];
i                 508 src/mastructcs.c       irow[i-isubtract]=irow[i];
i                 509 src/mastructcs.c       if(k!=mast1[i]){
i                 510 src/mastructcs.c 	for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
i                 511 src/mastructcs.c 	k=mast1[i];
i                 519 src/mastructcs.c   for(i=0;i<neq[0];++i){
i                 520 src/mastructcs.c     if(jq[i+1]-jq[i]>0){
i                 521 src/mastructcs.c       isize=jq[i+1]-jq[i];
i                 522 src/mastructcs.c       FORTRAN(isortii,(&irow[jq[i]-1],&mast1[jq[i]-1],&isize,&kflag));
i                  40 src/mastructem.c   ITG i,j,k,l,jj,ll,id,index,jdof1,jdof2,idof1,idof2,mpc1,mpc2,id1,id2,
i                  56 src/mastructem.c   for(i=0;i<mt**nk;++i){nactdof[i]=0;}
i                  60 src/mastructem.c       for(i=0;i<*ne;++i){
i                  62 src/mastructem.c 	  if(ipkon[i]<0) continue;
i                  63 src/mastructem.c 	  indexe=ipkon[i];
i                  64 src/mastructem.c 	  if(strcmp1(&lakon[8*i+3],"20")==0)nope=20;
i                  65 src/mastructem.c 	  else if (strcmp1(&lakon[8*i+3],"8")==0)nope=8;
i                  66 src/mastructem.c 	  else if (strcmp1(&lakon[8*i+3],"10")==0)nope=10;
i                  67 src/mastructem.c 	  else if (strcmp1(&lakon[8*i+3],"4")==0) nope=4;
i                  68 src/mastructem.c 	  else if (strcmp1(&lakon[8*i+3],"15")==0)nope=15;
i                  69 src/mastructem.c 	  else if (strcmp1(&lakon[8*i+3],"6")==0)nope=6;
i                  77 src/mastructem.c 	  imat=ielmat[i*mi[2]];
i                 104 src/mastructem.c       for(i=0;i<*ne;++i){
i                 106 src/mastructem.c 	  if(ipkon[i]<0) continue;
i                 110 src/mastructem.c 	  imat=ielmat[i*mi[2]];
i                 114 src/mastructem.c 	  indexe=ipkon[i];
i                 115 src/mastructem.c 	  if(strcmp1(&lakon[8*i+3],"20")==0)nope=20;
i                 116 src/mastructem.c 	  else if (strcmp1(&lakon[8*i+3],"8")==0)nope=8;
i                 117 src/mastructem.c 	  else if (strcmp1(&lakon[8*i+3],"10")==0)nope=10;
i                 118 src/mastructem.c 	  else if (strcmp1(&lakon[8*i+3],"4")==0)nope=4;
i                 119 src/mastructem.c 	  else if (strcmp1(&lakon[8*i+3],"15")==0)nope=15;
i                 120 src/mastructem.c 	  else if (strcmp1(&lakon[8*i+3],"6")==0)nope=6;
i                 121 src/mastructem.c 	  else if (strcmp1(&lakon[8*i],"E")==0){
i                 122 src/mastructem.c 	      lakonl[0]=lakon[8*i+7];
i                 124 src/mastructem.c 	  else if (strcmp1(&lakon[8*i],"D ")==0){
i                 148 src/mastructem.c   for(i=0;i<*nmpc;++i){
i                 149 src/mastructem.c       index=ipompc[i]-1;
i                 163 src/mastructem.c   for(i=0;i<*nboun;++i){
i                 164 src/mastructem.c       if(ndirboun[i]>mi[1]){continue;}
i                 165 src/mastructem.c       nactdof[mt*(nodeboun[i]-1)+ndirboun[i]]=0;
i                 168 src/mastructem.c   for(i=0;i<*nmpc;++i){
i                 169 src/mastructem.c     index=ipompc[i]-1;
i                 177 src/mastructem.c   for(i=0;i<*nk;++i){
i                 179 src/mastructem.c 	if(nactdof[mt*i+j]!=0){
i                 182 src/mastructem.c           nactdof[mt*i+j]=neq[0];
i                 185 src/mastructem.c 	    nactdof[mt*i+j]=0;
i                 191 src/mastructem.c   for(i=0;i<*nk;++i){
i                 192 src/mastructem.c       if(nactdof[mt*i]!=0){
i                 195 src/mastructem.c         nactdof[mt*i]=neq[1];
i                 198 src/mastructem.c 	  nactdof[mt*i]=0;
i                 213 src/mastructem.c     for(i=0;i<4**nk;++i){ipointer[i]=0;}
i                 219 src/mastructem.c     for(i=0;i<*ne;++i){
i                 221 src/mastructem.c       if(ipkon[i]<0) continue;
i                 222 src/mastructem.c       indexe=ipkon[i];
i                 223 src/mastructem.c       if(strcmp1(&lakon[8*i+3],"20")==0)nope=20;
i                 224 src/mastructem.c       else if (strcmp1(&lakon[8*i+3],"8")==0)nope=8;
i                 225 src/mastructem.c       else if (strcmp1(&lakon[8*i+3],"10")==0)nope=10;
i                 226 src/mastructem.c       else if (strcmp1(&lakon[8*i+3],"4")==0)nope=4;
i                 227 src/mastructem.c       else if (strcmp1(&lakon[8*i+3],"15")==0)nope=15;
i                 228 src/mastructem.c       else if (strcmp1(&lakon[8*i+3],"6")==0)nope=6;
i                 363 src/mastructem.c     for(i=0;i<*ne;++i){
i                 365 src/mastructem.c       if(ipkon[i]<0) continue;
i                 369 src/mastructem.c       imat=ielmat[i*mi[2]];
i                 373 src/mastructem.c       indexe=ipkon[i];
i                 374 src/mastructem.c       if(strcmp1(&lakon[8*i+3],"20")==0)nope=20;
i                 375 src/mastructem.c       else if (strcmp1(&lakon[8*i+3],"8")==0)nope=8;
i                 376 src/mastructem.c       else if (strcmp1(&lakon[8*i+3],"10")==0)nope=10;
i                 377 src/mastructem.c       else if (strcmp1(&lakon[8*i+3],"4")==0)nope=4;
i                 378 src/mastructem.c       else if (strcmp1(&lakon[8*i+3],"15")==0)nope=15;
i                 379 src/mastructem.c       else if (strcmp1(&lakon[8*i+3],"6")==0)nope=6;
i                 380 src/mastructem.c       else if (strcmp1(&lakon[8*i],"E")==0){
i                 381 src/mastructem.c 	  lakonl[0]=lakon[8*i+7];
i                 383 src/mastructem.c       else if (strcmp1(&lakon[8*i],"D ")==0){
i                 521 src/mastructem.c     for(i=0;i<neq[2];++i){
i                 522 src/mastructem.c 	if(ipointer[i]==0){
i                 523 src/mastructem.c 	    if(i>=neq[1]) continue;
i                 527 src/mastructem.c 		    if(nactdof[mt*j+k]==i+1){
i                 539 src/mastructem.c 	istart=ipointer[i];
i                 543 src/mastructem.c 	    irow[istartold-1]=i+1;
i                 584 src/mastructem.c     for(i=0;i<neq[1];++i){icol[i]=0;}
i                 586 src/mastructem.c     for(i=0;i<nmast;++i){
i                 587 src/mastructem.c       if(mast1[i]==irow[i]){++isubtract;}
i                 589 src/mastructem.c 	mast1[i-isubtract]=mast1[i];
i                 590 src/mastructem.c 	irow[i-isubtract]=irow[i];
i                 591 src/mastructem.c 	if(k!=mast1[i]){
i                 592 src/mastructem.c 	  for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
i                 593 src/mastructem.c 	  k=mast1[i];
i                 603 src/mastructem.c     for(i=0;i<neq[1];++i){
i                 604 src/mastructem.c       if(jq[i+1]-jq[i]>0){
i                 605 src/mastructem.c 	isize=jq[i+1]-jq[i];
i                 606 src/mastructem.c 	FORTRAN(isortii,(&irow[jq[i]-1],&mast1[jq[i]-1],&isize,&kflag));
i                  32 src/mastructf.c   ITG i,j,k,l,index,idof1,idof2,node1,isubtract,nmast,ifree=0,istart,istartold,
i                  47 src/mastructf.c   for(i=0;i<*ne;i++){
i                  48 src/mastructf.c       idof1=i+1;
i                  49 src/mastructf.c       if(strcmp1(&lakon[8*i+3],"8")==0){
i                  51 src/mastructf.c       }else if(strcmp1(&lakon[8*i+3],"6")==0){
i                  57 src/mastructf.c       index=ipnei[i];
i                  72 src/mastructf.c   for(i=0;i<*neq;++i){
i                  73 src/mastructf.c       if(ipointer[i]==0){
i                  75 src/mastructf.c 	  printf("       element=%" ITGFORMAT "\n",i+1);
i                  78 src/mastructf.c       istart=ipointer[i];
i                  82 src/mastructf.c 	  irow[istartold-1]=i+1;
i                 123 src/mastructf.c   for(i=0;i<*neq;++i){icol[i]=0;}
i                 125 src/mastructf.c   for(i=0;i<nmast;++i){
i                 126 src/mastructf.c       if(mast1[i]==irow[i]){++isubtract;}
i                 128 src/mastructf.c 	  mast1[i-isubtract]=mast1[i];
i                 129 src/mastructf.c 	  irow[i-isubtract]=irow[i];
i                 130 src/mastructf.c 	  if(k!=mast1[i]){
i                 131 src/mastructf.c 	      for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
i                 132 src/mastructf.c 	      k=mast1[i];
i                 142 src/mastructf.c   for(i=0;i<*neq;++i){
i                 143 src/mastructf.c       if(jq[i+1]-jq[i]>0){
i                 144 src/mastructf.c 	  isize=jq[i+1]-jq[i];
i                 145 src/mastructf.c 	  FORTRAN(isortii,(&irow[jq[i]-1],&mast1[jq[i]-1],&isize,&kflag));
i                  35 src/mastructrad.c   ITG three=3,i,j,k,l,ii,jj,icav,jcav,*mast1rad=NULL,*irowrad=NULL,
i                  50 src/mastructrad.c       i=nloadtr[ii-1]-1;
i                  51 src/mastructrad.c       strcpy1(crcav,&sideload[20*i+17],three);
i                  68 src/mastructrad.c   for(i=0;i<*ntr;++i){
i                  69 src/mastructrad.c       if(ipointerrad[i]==0){
i                  71 src/mastructrad.c 	  printf("       DOF=%" ITGFORMAT "\n",i);
i                  74 src/mastructrad.c       istart=ipointerrad[i];
i                  78 src/mastructrad.c 	  irowrad[istartold-1]=i+1;
i                 115 src/mastructrad.c   for(i=0;i<*ntr;++i){icolrad[i]=0;}
i                 117 src/mastructrad.c   for(i=0;i<nmast;++i){
i                 118 src/mastructrad.c       if(mast1rad[i]==irowrad[i]){++isubtract;}
i                 120 src/mastructrad.c 	  mast1rad[i-isubtract]=mast1rad[i];
i                 121 src/mastructrad.c 	  irowrad[i-isubtract]=irowrad[i];
i                 122 src/mastructrad.c 	  if(k!=mast1rad[i]){
i                 123 src/mastructrad.c 	      for(l=k;l<mast1rad[i];++l){jqrad[l]=i+1-isubtract;}
i                 124 src/mastructrad.c 	      k=mast1rad[i];
i                 134 src/mastructrad.c   for(i=0;i<*ntr;++i){
i                 135 src/mastructrad.c       if(jqrad[i+1]-jqrad[i]>0){
i                 136 src/mastructrad.c 	  isize=jqrad[i+1]-jqrad[i];
i                 137 src/mastructrad.c 	  FORTRAN(isortii,(&irowrad[jqrad[i]-1],&mast1rad[jqrad[i]-1],
i                  34 src/matrixstorage.c   ITG i,j,l,*irow=NULL,*ai=NULL,*aj=NULL,kflag=2,ndim,jref,kstart,klen,
i                  73 src/matrixstorage.c     for(i=0;i<*nk;i++){
i                  74 src/matrixstorage.c       if(inotr[2*i]!=0){
i                  92 src/matrixstorage.c     for(i=0;i<*neq;i++){
i                  93 src/matrixstorage.c       ai[k]=i+1;
i                  94 src/matrixstorage.c       aj[k]=i+1;
i                  95 src/matrixstorage.c       aa[k]=ad[i];
i                  99 src/matrixstorage.c     for(i=0;i<*neq;i++){
i                 100 src/matrixstorage.c       for(j=0;j<icol[i];j++){
i                 101 src/matrixstorage.c         ai[k]=i+1;
i                 133 src/matrixstorage.c     for(i=0;i<*neq;i++){
i                 134 src/matrixstorage.c       ai[k]=i+1;
i                 135 src/matrixstorage.c       aj[k]=i+1;
i                 136 src/matrixstorage.c       aa[k]=ad[i];
i                 140 src/matrixstorage.c     for(i=0;i<*neq;i++){
i                 141 src/matrixstorage.c       for(j=0;j<icol[i];j++){
i                 142 src/matrixstorage.c         ai[k]=i+1;
i                 147 src/matrixstorage.c         aj[k]=i+1;
i                 157 src/matrixstorage.c     for(i=0;i<*neq;i++){
i                 185 src/matrixstorage.c     for(i=0;i<ndim2;i++){
i                 186 src/matrixstorage.c       il=ai[i];
i                 187 src/matrixstorage.c       ic=aj[i];
i                 213 src/matrixstorage.c       aa3[index+3*i3co+i3lo]=aa[i];
i                 220 src/matrixstorage.c     for (i=0;i<*nk;i++){
i                 221 src/matrixstorage.c       idof=nactdof[mt*i+1];
i                 223 src/matrixstorage.c       itrans=inotr[2*i];
i                 233 src/matrixstorage.c         FORTRAN(transformatrix,(&trab[7*itrans-7],&co[3*i],a));
i                 243 src/matrixstorage.c     for(i=0;i<npoint;i++){
i                 244 src/matrixstorage.c       i3l=ipoint[i]%neq3;
i                 245 src/matrixstorage.c       i3c=(ipoint[i]-i3l)/neq3;
i                 247 src/matrixstorage.c       FORTRAN(mult,(&aa3[9*i],&trans[9*i3c],&n));
i                 252 src/matrixstorage.c     for(i=0;i<npoint;i++){
i                 253 src/matrixstorage.c       i3l=ipoint[i]%neq3;
i                 255 src/matrixstorage.c       FORTRAN(mult,(&aa3[9*i],&trans[9*i3l],&n));
i                 261 src/matrixstorage.c     for(i=0;i<npoint;i++){
i                 262 src/matrixstorage.c       i3l=ipoint[i]%neq3;
i                 263 src/matrixstorage.c       i3c=(ipoint[i]-i3l)/neq3;
i                 273 src/matrixstorage.c         aa[k]=aa3[9*i+j];
i                 282 src/matrixstorage.c   for(i=0;i<*neq;i++){
i                 293 src/matrixstorage.c   for(i=0;i<ndim;i++){
i                 294 src/matrixstorage.c     fprintf(f2,"%" ITGFORMAT " %" ITGFORMAT " %20.13e\n",ai[i],aj[i],aa[i]);
i                 322 src/matrixstorage.c     for(i=0;i<*neq;i++){
i                 323 src/matrixstorage.c       ai[k]=i+1;
i                 324 src/matrixstorage.c       aj[k]=i+1;
i                 325 src/matrixstorage.c       aa[k]=adb[i];
i                 329 src/matrixstorage.c     for(i=0;i<*neq;i++){
i                 330 src/matrixstorage.c       for(j=0;j<icol[i];j++){
i                 331 src/matrixstorage.c         ai[k]=i+1;
i                 357 src/matrixstorage.c     for(i=0;i<*neq;i++){
i                 358 src/matrixstorage.c       ai[k]=i+1;
i                 359 src/matrixstorage.c       aj[k]=i+1;
i                 360 src/matrixstorage.c       aa[k]=adb[i];
i                 364 src/matrixstorage.c     for(i=0;i<*neq;i++){
i                 365 src/matrixstorage.c       for(j=0;j<icol[i];j++){
i                 366 src/matrixstorage.c         ai[k]=i+1;
i                 371 src/matrixstorage.c         aj[k]=i+1;
i                 381 src/matrixstorage.c     for(i=0;i<*neq;i++){
i                 409 src/matrixstorage.c     for(i=0;i<ndim2;i++){
i                 410 src/matrixstorage.c       il=ai[i];
i                 411 src/matrixstorage.c       ic=aj[i];
i                 437 src/matrixstorage.c       aa3[index+3*i3co+i3lo]=aa[i];
i                 444 src/matrixstorage.c     for (i=0;i<*nk;i++){
i                 445 src/matrixstorage.c       idof=nactdof[mt*i+1];
i                 447 src/matrixstorage.c       itrans=inotr[2*i];
i                 457 src/matrixstorage.c         FORTRAN(transformatrix,(&trab[7*itrans-7],&co[3*i],a));
i                 467 src/matrixstorage.c     for(i=0;i<npoint;i++){
i                 468 src/matrixstorage.c       i3l=ipoint[i]%neq3;
i                 469 src/matrixstorage.c       i3c=(ipoint[i]-i3l)/neq3;
i                 471 src/matrixstorage.c       FORTRAN(mult,(&aa3[9*i],&trans[9*i3c],&n));
i                 476 src/matrixstorage.c     for(i=0;i<npoint;i++){
i                 477 src/matrixstorage.c       i3l=ipoint[i]%neq3;
i                 479 src/matrixstorage.c       FORTRAN(mult,(&aa3[9*i],&trans[9*i3l],&n));
i                 485 src/matrixstorage.c     for(i=0;i<npoint;i++){
i                 486 src/matrixstorage.c       i3l=ipoint[i]%neq3;
i                 487 src/matrixstorage.c       i3c=(ipoint[i]-i3l)/neq3;
i                 497 src/matrixstorage.c         aa[k]=aa3[9*i+j];
i                 506 src/matrixstorage.c   for(i=0;i<*neq;i++){
i                 517 src/matrixstorage.c   for(i=0;i<ndim;i++){
i                 518 src/matrixstorage.c     fprintf(f3,"%" ITGFORMAT " %" ITGFORMAT " %20.13e\n",ai[i],aj[i],aa[i]);
i                 554 src/matrixstorage.c       for(i=0;i<*neq;i++){
i                 555 src/matrixstorage.c 	  inode=(ITG)((double)nactdofinv[(ITG)i]/mt)+1;
i                 556 src/matrixstorage.c 	  idir=nactdofinv[(ITG)i]-mt*(inode-1);
i                 560 src/matrixstorage.c       for(i=0;i<*neq/2;i++){
i                 561 src/matrixstorage.c 	  inode=(ITG)((double)nactdofinv[(ITG)i]/mt)+1;
i                 562 src/matrixstorage.c 	  idir=nactdofinv[(ITG)i]-mt*(inode-1);
i                 109 src/nonlingeo.c       nef=0,nface,nfreestream,nsolidsurf,i,indexe,icfd=0,id,*neij=NULL,
i                 216 src/nonlingeo.c 	      for(i=0;i<*ne;i++){
i                 217 src/nonlingeo.c 		  if(ipkon[i]<0) continue;
i                 218 src/nonlingeo.c 		  if(strcmp1(&lakon[i*8],"ED")==0){
i                 377 src/nonlingeo.c   for(i=0;i<*ne;++i){
i                 378 src/nonlingeo.c       if(ipkon[i]<0) continue;
i                 379 src/nonlingeo.c       indexe=ipkon[i];
i                 380 src/nonlingeo.c       if(strcmp1(&lakon[8*i],"F")==0){
i                 381 src/nonlingeo.c 	  icfd=1;nactdohinv[nef]=i+1;nef++;nactdoh[i]=nef;}
i                 952 src/nonlingeo.c       for(i=0;i<3;i++){cam[i]=0.;}for(i=3;i<5;i++){cam[i]=0.5;}
i                1076 src/nonlingeo.c 		      for(i=0;i<*ne;i++){
i                1077 src/nonlingeo.c 			  if(ipkon[i]<0) continue;
i                1079 src/nonlingeo.c 			      if(strcmp1(&lakon[i*8],"ES")==0){
i                1080 src/nonlingeo.c 				  if(strcmp1(&lakon[i*8+6],"C")==0){
i                1081 src/nonlingeo.c 				      imat=ielmat[i*mi[2]];
i                1346 src/nonlingeo.c 	  for(i=0;i<3;i++){cam[i]=0.;}for(i=3;i<5;i++){cam[i]=0.5;}
i                1390 src/nonlingeo.c 			  for(i=0;i<*ne;i++){
i                1391 src/nonlingeo.c 			      if(ipkon[i]<0) continue;
i                1393 src/nonlingeo.c 				  if(strcmp1(&lakon[i*8],"ES")==0){
i                1394 src/nonlingeo.c 				      if(strcmp1(&lakon[i*8+6],"C")==0){
i                1395 src/nonlingeo.c 					  imat=ielmat[i*mi[2]];
i                2061 src/nonlingeo.c 	for(i=0;i<*ntie;i++){	    
i                2062 src/nonlingeo.c 	    if(itietri[2*i+1]-itietri[2*i]+1>ntrimax)		
i                2063 src/nonlingeo.c 		ntrimax=itietri[2*i+1]-itietri[2*i]+1;  	
i                  40 src/pardiso.c    ITG i,j,k,l,maxfct=1,mnum=1,phase=12,nrhs=1,*perm=NULL,mtype,
i                  77 src/pardiso.c    for(i=0;i<64;i++){pt[i]=0;}
i                  98 src/pardiso.c  	  for(i=*neq-1;i>=0;--i){
i                  99 src/pardiso.c  	      for(j=0;j<icol[i];++j){
i                 103 src/pardiso.c  	      pointers[i]=k--;
i                 104 src/pardiso.c  	      icolpardiso[k]=i+1;
i                 105 src/pardiso.c  	      aupardiso[k]=ad[i];
i                 110 src/pardiso.c  	  for(i=*neq-1;i>=0;--i){
i                 111 src/pardiso.c  	      for(j=0;j<icol[i];++j){
i                 115 src/pardiso.c  	      pointers[i]=k--;
i                 116 src/pardiso.c  	      icolpardiso[k]=i+1;
i                 117 src/pardiso.c  	      aupardiso[k]=ad[i]-*sigma*adb[i];
i                 137 src/pardiso.c  	  for(i=0;i<*neq;i++){
i                 138 src/pardiso.c  	      for(j=0;j<icol[i];j++){
i                 140 src/pardiso.c  		      icolpardiso[k2]=i+1;
i                 149 src/pardiso.c  	  for(i=0;i<*neq;i++){
i                 150 src/pardiso.c  	      icolpardiso[k2]=i+1;
i                 151 src/pardiso.c  	      irowpardiso[k2]=i+1;
i                 152 src/pardiso.c  	      aupardiso[k2]=ad[i];
i                 168 src/pardiso.c  	  for(i=0;i<*neq;i++){
i                 169 src/pardiso.c  	      j=i+1;
i                 176 src/pardiso.c  		      pointers[i+1]=k+1;
i                 208 src/pardiso.c  	  for(i=0;i<*neq;i++){
i                 209 src/pardiso.c  	      for(j=0;j<icol[i];j++){
i                 210 src/pardiso.c  		  icolpardiso[k]=i+1;
i                 228 src/pardiso.c  	  for(i=0;i<*neq;i++){
i                 229 src/pardiso.c  	      j=i+1;
i                 239 src/pardiso.c  		      pointers[i+1]=k+1;
i                 268 src/pardiso.c  	  for(i=*neq-1;i>=0;i--){
i                 270 src/pardiso.c  	      for(j=jq[i+1]-1;j>=jq[i];j--){
i                 274 src/pardiso.c  	      icolpardiso[--k]=i+1;
i                 275 src/pardiso.c  	      aupardiso[k]=ad[i];
i                 276 src/pardiso.c  	      for(j=pointers[i+1]-1;j>=pointers[i];j--){
i                 280 src/pardiso.c  	      pointers[i+1]=l;
i                 299 src/pardiso.c        msglvl=0,i,error=0;
i                 324 src/pardiso.c    for(i=0;i<*neq;i++){b[i]=x[i];}
i                  39 src/pardiso_as.c   ITG i,j,k,maxfct=1,mnum=1,mtype=11,phase=12,nrhs=1,*perm=NULL,
i                  72 src/pardiso_as.c   for(i=0;i<64;i++){ptas[i]=0;}
i                  90 src/pardiso_as.c     for(i=0;i<*neq;i++){
i                  91 src/pardiso_as.c       pointersas[i]=k+1;
i                  95 src/pardiso_as.c       for(j=0;j<i;j++){
i                  96 src/pardiso_as.c 	ifortran=i+1;
i                 111 src/pardiso_as.c       irowpardisoas[k]=i+1;
i                 112 src/pardiso_as.c       aupardisoas[k]=ad[i];
i                 117 src/pardiso_as.c       for(j=jq[i];j<jq[i+1];j++){
i                 125 src/pardiso_as.c     for(i=0;i<*neq;i++){
i                 126 src/pardiso_as.c       pointersas[i]=k+1;
i                 130 src/pardiso_as.c       for(j=0;j<i;j++){
i                 131 src/pardiso_as.c 	ifortran=i+1;
i                 146 src/pardiso_as.c       irowpardisoas[k]=i+1;
i                 147 src/pardiso_as.c       aupardisoas[k]=ad[i]-*sigma*adb[i];
i                 152 src/pardiso_as.c       for(j=jq[i];j<jq[i+1];j++){
i                 172 src/pardiso_as.c       msglvl=0,i,error=0;
i                 188 src/pardiso_as.c   for(i=0;i<*neq;i++){b[i]=x[i];}
i                  74 src/pcgsolver.c   ITG i=0;
i                  79 src/pcgsolver.c   for (i=0; i<neq; i++)	--iz[i];
i                  80 src/pcgsolver.c   for (i=0; i<len; i++)	--ia[i];
i                 116 src/pcgsolver.c   for (i=0; i<neq; i++)	x[i] *= Factor[i];
i                 151 src/pcgsolver.c   ITG			i=0, k=0, ncg=0,iam,ier=0;
i                 160 src/pcgsolver.c   for (i=0; i<neq; i++)
i                 162 src/pcgsolver.c       x[i] = 0.0;	  
i                 163 src/pcgsolver.c       r[i] = b[i];
i                 164 src/pcgsolver.c       err=fabs(r[i]);
i                 212 src/pcgsolver.c 	  for (i=0; i<neq; i++)	g[i] = ekm1*g[i]-rho[i];
i                 219 src/pcgsolver.c 	  for (i=0; i<neq; i++)	g[i] = -rho[i];
i                 231 src/pcgsolver.c       for (i=0; i<neq; i++)
i                 233 src/pcgsolver.c 	  x[i] += qk*g[i];
i                 234 src/pcgsolver.c 	  r[i] += qk*z[i];
i                 235 src/pcgsolver.c 	  err=fabs(r[i]);
i                 283 src/pcgsolver.c   ITG			i=0, j=0, jlo=0, jup=0;
i                 287 src/pcgsolver.c   for (i=0; i<neq; i++)		d[i] =  1.0/sqrt(A[iz[i]]);
i                 291 src/pcgsolver.c   for (i=0; i<neq; i++)	b[i] *= -d[i];
i                 296 src/pcgsolver.c   for (i=1; i<neq; i++)
i                 298 src/pcgsolver.c       jlo = iz[i-1]+1;
i                 299 src/pcgsolver.c       jup = iz[i];
i                 300 src/pcgsolver.c       for (j=jlo; j<=jup; j++)	A[j] *= d[i]*d[ia[j]];
i                 326 src/pcgsolver.c   ITG				i=0, j=0, jlo=0, jup=0, k=0;
i                 331 src/pcgsolver.c   for (i=1; i<neq; i++)
i                 333 src/pcgsolver.c       z[i] = A[iz[i]]*p[i];
i                 337 src/pcgsolver.c       jlo = iz[i-1]+1;	
i                 341 src/pcgsolver.c       jup = iz[i]-1;				
i                 345 src/pcgsolver.c 	  z[i] += A[j]*p[k];
i                 346 src/pcgsolver.c 	  z[k] += A[j]*p[i];
i                 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;
i                 383 src/pcgsolver.c 	for (i=1; i<neq; i++)
i                 385 src/pcgsolver.c 		jlo = iz[i-1]+1;			/*..first non-zero element in current row......	*/
i                 386 src/pcgsolver.c 		jup = iz[i];				/*..diagonal element...........................	*/
i                 393 src/pcgsolver.c 	for (i=1; i<nILU; i++)
i                 395 src/pcgsolver.c 		jlo = iz[i-1]+1;			/*..first non-zero element in current row......	*/
i                 396 src/pcgsolver.c 		jup = iz[i];				/*..diagonal element...........................	*/
i                 416 src/pcgsolver.c 		id = iz[i];
i                 447 src/pcgsolver.c 	ITG				i=0, j=0, jlo=0, jup=0;
i                 451 src/pcgsolver.c 	for (i=1; i<neq; i++)
i                 454 src/pcgsolver.c 		jlo = iz[i-1]+1;			/*..first non-zero element in current row......	*/
i                 455 src/pcgsolver.c 		jup = iz[i];				/*..diagonal element in current row............	*/
i                 458 src/pcgsolver.c 		rho[i] = (r[i]-s)/C[jup];
i                 460 src/pcgsolver.c 	for (i=neq-1; i>0; i--)
i                 462 src/pcgsolver.c 		rho[i] /= C[iz[i]];
i                 463 src/pcgsolver.c 		jlo = iz[i-1]+1;			/*..first non-zero element in current row......	*/
i                 464 src/pcgsolver.c 		jup = iz[i]-1;				/*..diagonal element in current row............	*/
i                 466 src/pcgsolver.c 			rho[ia[j]] -= C[j]*rho[i];
i                 479 src/pcgsolver.c 	ITG 		i=0;
i                 482 src/pcgsolver.c 		for (i=0; i<n; i++)	       *Sum += a[i]*b[i];
i                 507 src/pcgsolver.c 	ITG			i=0, k=0, ncg=0,iam;
i                 518 src/pcgsolver.c 	for (i=0; i<neq; i++)
i                 520 src/pcgsolver.c 		x[i] =  0.0;					/*..start vector x0 = 0....................	*/
i                 521 src/pcgsolver.c 		r[i] =  b[i];					/*..residual vector r0 = Ax+b = b..........	*/
i                 522 src/pcgsolver.c 		p[i] = -r[i];					/*..initial search vector..................	*/
i                 523 src/pcgsolver.c 		err=fabs(r[i]);
i                 544 src/pcgsolver.c 			for (i=0; i<neq; i++)	p[i] = ekm1*p[i]-r[i];
i                 554 src/pcgsolver.c 		for (i=0; i<neq; i++)
i                 556 src/pcgsolver.c 			x[i] = x[i] + qk*p[i];
i                 557 src/pcgsolver.c 			r[i] = r[i] + qk*z[i];
i                 558 src/pcgsolver.c 			err=fabs(r[i]);
i                  37 src/precontact.c     ITG i,j,ntrimax,*nx=NULL,*ny=NULL,*nz=NULL,im,
i                  61 src/precontact.c     for(i=0;i<*ntie;i++){	    
i                  62 src/precontact.c 	if(itietri[2*i+1]-itietri[2*i]+1>ntrimax)		
i                  63 src/precontact.c 	    ntrimax=itietri[2*i+1]-itietri[2*i]+1;  	
i                 114 src/precontact.c     for(i=0;i<*ntie;i++){
i                 115 src/precontact.c 	ii=i+1;	   
i                 119 src/precontact.c 	if(tieset[i*(81*3)+80]=='C'){		
i                 120 src/precontact.c 	    nstart=itietri[2*i]-1;		
i                 121 src/precontact.c 	    ntri=itietri[2*i+1]-nstart;		
i                 140 src/precontact.c 	    for(l=itiefac[2*i];l<=itiefac[2*i+1];l++){
i                  26 src/preiter.c    ITG precFlg,niter=5000000,ndim,i,j,k,ier,*icol=NULL,*irow=NULL,
i                  45 src/preiter.c      for(i=0;i<*nzs;++i){
i                  46 src/preiter.c        irow_save[i]=irow[i];
i                  48 src/preiter.c      for(i=0;i<*neq;++i){
i                  49 src/preiter.c        icol_save[i]=icol[i];
i                  63 src/preiter.c    for(i=*neq-1;i>=0;--i){
i                  64 src/preiter.c      for(j=0;j<icol[i];++j){
i                  65 src/preiter.c        icol[--k]=i+1;
i                  71 src/preiter.c    for(i=0;i<*neq;++i){
i                  72 src/preiter.c      au[k]=ad[i];
i                  91 src/preiter.c    for(i=0;i<*neq;++i){
i                  92 src/preiter.c      b[i]=u[i];
i                 100 src/preiter.c      for(i=0;i<*nzs;++i){
i                 101 src/preiter.c        irow[i]=irow_save[i];
i                 103 src/preiter.c      for(i=0;i<*neq;++i){
i                 104 src/preiter.c        icol[i]=icol_save[i];
i                  34 src/radcyc.c     ITG i,is,nsegments,idtie,nkt,icntrl,imag=0,*kontri=NULL,mt=mi[1]+1,
i                  56 src/radcyc.c       for(i=0;i<*nk;i++) inocs[i]=-1;
i                  59 src/radcyc.c     for(i=0;i<*mcs;i++){
i                  60 src/radcyc.c       is=cs[17*i+4];
i                  62 src/radcyc.c       ielset=cs[17*i+12];
i                  77 src/radcyc.c             inocs[node]=i;
i                  95 src/radcyc.c               inocs[node]=i;
i                 109 src/radcyc.c     for(i=4**ntri;i<4**ntrit;i++) kontri[i]=0;
i                 111 src/radcyc.c     for(i=0;i<*ntri;i++){
i                 112 src/radcyc.c       node1=kontri[4*i];
i                 115 src/radcyc.c       node2=kontri[4*i+1];
i                 117 src/radcyc.c       node3=kontri[4*i+2];
i                 121 src/radcyc.c         j=i+k**ntri;
i                 125 src/radcyc.c         kontri[4*j+3]=kontri[4*i+3];
i                 141 src/radcyc.c       for(i=1;i<is;i++){
i                 143 src/radcyc.c         theta=i*2.*pi/cs[17*jj];
i                 147 src/radcyc.c   	  co[3*l+i*3**nk]=co[3*l];
i                 148 src/radcyc.c   	  co[1+3*l+i*3**nk]=co[1+3*l]-theta;
i                 149 src/radcyc.c   	  co[2+3*l+i*3**nk]=co[2+3*l];
i                  77 src/radflowload.c   ITG nhrs=1,info=0,i,j,iin=0,icntrl,icutb=0,iin_abs=0,mt=mi[1]+1,im,
i                 106 src/radflowload.c 	      for(i=0;i<mt**nk;i++) v[i]=vold[i];
i                 363 src/radflowload.c       for(i=0;i<*ntr;i++){
i                 364 src/radflowload.c 	  node=nelemload[2*nloadtr[i]-1];
i                 366 src/radflowload.c 	      tenv[i]=vold[mt*(node-1)]-physcon[0];
i                 368 src/radflowload.c 	      tenv[i]=xloadact[2*nloadtr[i]-1]-physcon[0];
i                 420 src/radflowload.c 	      for(i=0; i<num_cpus; i++)  {
i                 421 src/radflowload.c 		  ithread[i]=i;
i                 422 src/radflowload.c 		  pthread_create(&tid[i], NULL, (void *)calcviewmt, (void *)&ithread[i]);
i                 424 src/radflowload.c 	      for(i=0; i<num_cpus; i++)  pthread_join(tid[i], NULL);
i                 426 src/radflowload.c 	      for(i=0;i<*ntr;i++){
i                 428 src/radflowload.c 		      adview[i]+=adview[i+j**ntr];
i                 433 src/radflowload.c 	      for(i=0;i<2**nzsrad;i++){
i                 435 src/radflowload.c 		      auview[i]+=auview[i+j*2**nzsrad];
i                 549 src/radflowload.c void *calcviewmt(ITG *i){
i                 553 src/radflowload.c     indexad=*i**ntr1;
i                 554 src/radflowload.c     indexau=*i*2**nzsrad1;
i                 555 src/radflowload.c     indexdi=*i**ntrit1;
i                 558 src/radflowload.c     ntria=*i*nedelta+1;
i                 559 src/radflowload.c     ntrib=(*i+1)*nedelta;
i                  48 src/readfrd.c    register ITG i;
i                  53 src/readfrd.c      for(i=0; i<lcase[nr].ncomps; i++) SFREE(lcase[nr].dat[i]);
i                  57 src/readfrd.c      for(i=0; i<lcase[nr].npheader; i++) SFREE(lcase[nr].pheader[i]);
i                  60 src/readfrd.c    for(i=0; i<lcase[nr].ncomps; i++)
i                  62 src/readfrd.c      SFREE(lcase[nr].compName[i]);
i                  63 src/readfrd.c      SFREE(lcase[nr].icname[i]);
i                  90 src/readfrd.c    register ITG i=0, j=0;
i                 230 src/readfrd.c        for(i=0; i<node_field_size; i++) node[i].indx=-1;
i                 260 src/readfrd.c            for(i=anz->nmax+1; i<node_field_size; i++) node[i].indx=-1;
i                 294 src/readfrd.c         for(i=0; i<nod_data; i++)
i                 389 src/readfrd.c          for (i=0; i<elem_data; i++)
i                 493 src/readfrd.c                for (i=0; i<ipuf; i++)
i                 495 src/readfrd.c  		elem[anz->e].nod[i] = stoi(&rec_str[3+j*5], 1, 5);
i                 499 src/readfrd.c                    if (i<ipuf-1) length = frecord( handle, rec_str ); j=0;
i                 506 src/readfrd.c                for (i=0; i<ipuf; i++)
i                 508 src/readfrd.c                  elem[anz->e].nod[i] = stoi(&rec_str[3+j*10], 1, 10);
i                 512 src/readfrd.c                    if (i<ipuf-1) length = frecord( handle, rec_str ); j=0;
i                 592 src/readfrd.c          for(i=0;i<lcase[anz->l].npheader; i++)
i                 594 src/readfrd.c            if(compare(&lcase[anz->l].pheader[i][5],"PHID", 4)==4)
i                 596 src/readfrd.c              sscanf(lcase[anz->l].pheader[i],"%*s %" ITGFORMAT "", &ipuf);
i                 600 src/readfrd.c          if(ipuf!=-1) for(i=0;i<lcase[anz->l].npheader; i++)
i                 602 src/readfrd.c            if(compare(&lcase[anz->l].pheader[i][5],"PMODE", 5)==5)
i                 604 src/readfrd.c              sscanf(lcase[anz->l].pheader[i],"%*s %" ITGFORMAT "", &ipuf);
i                 643 src/readfrd.c                for(i=0; i<lcase[anz->l].ncomps; i++)
i                 646 src/readfrd.c                  lcase[anz->l].dat[i][nodenr]= value[i];
i                 692 src/readfrd.c                  i=6;
i                 714 src/readfrd.c  	        for(i=0; i<maxcomps; i++) lcase[anz->l].dat[i][nodenr]= stof(&rec_str[n+i*12], 1, 12);
i                 718 src/readfrd.c              else i=0;
i                 732 src/readfrd.c                lcase[anz->l].dat[i][nodenr]= stof(&rec_str[n+j*12], 1, 12);
i                 733 src/readfrd.c                i++;j++;
i                 734 src/readfrd.c              }while((j<6)&&(i<lcase[anz->l].ncomps));
i                 775 src/readfrd.c            for(i=0; i<(lcase[anz->l].ncomps); i++)
i                 777 src/readfrd.c              if ( (lcase[anz->l].compName[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
i                 779 src/readfrd.c              if ( (lcase[anz->l].icname[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
i                 781 src/readfrd.c              lcase[anz->l].max[i]=-MAX_FLOAT;
i                 782 src/readfrd.c              lcase[anz->l].min[i]=MAX_FLOAT;
i                 821 src/readfrd.c              for (i=0; i<ncomps; i++) if(lcase[anz->l].iexist[i]!=1) nvals++;
i                 841 src/readfrd.c                      for(i=0; i<nvals/6; i++) 
i                 854 src/readfrd.c                for(i=0; i<(lcase[anz->l].ncomps); i++)
i                 856 src/readfrd.c                  if ( (lcase[anz->l].dat[i] = (float *)malloc( (anz->nmax+1) * sizeof(float))) == NULL )
i                 858 src/readfrd.c                  for(j=0; j<=anz->nmax; j++) lcase[anz->l].dat[i][j]=0.;
i                 870 src/readfrd.c            for(i=0; i<ipuf; i++)
i                 872 src/readfrd.c              if ( (compName[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
i                 875 src/readfrd.c            for (i=0; i<ipuf; i++) strcpy(compName[i], lcase[anz->l].compName[ atoi(dat[i+3])-1 ]);
i                 876 src/readfrd.c            for (i=0; i<ipuf; i++) strcpy(lcase[anz->l].compName[i],compName[i]);
i                 887 src/readfrd.c            for (i=0; i<ipuf; i++) menu[i] = lcase[anz->l].menu[atoi(dat[i+3])-1];
i                 888 src/readfrd.c            for (i=0; i<ipuf; i++) lcase[anz->l].menu[i] =menu[i];
i                 889 src/readfrd.c            for (i=0; i<ipuf; i++) ictype[i] = lcase[anz->l].ictype[atoi(dat[i+3])-1];
i                 890 src/readfrd.c            for (i=0; i<ipuf; i++) lcase[anz->l].ictype[i] =ictype[i];
i                 891 src/readfrd.c            for (i=0; i<ipuf; i++) icind1[i] = lcase[anz->l].icind1[atoi(dat[i+3])-1];
i                 892 src/readfrd.c            for (i=0; i<ipuf; i++) lcase[anz->l].icind1[i] =icind1[i];
i                 893 src/readfrd.c            for (i=0; i<ipuf; i++) icind2[i] = lcase[anz->l].icind2[atoi(dat[i+3])-1];
i                 894 src/readfrd.c            for (i=0; i<ipuf; i++) lcase[anz->l].icind2[i] =icind2[i];
i                 895 src/readfrd.c            for (i=0; i<ipuf; i++) iexist[i] = lcase[anz->l].iexist[atoi(dat[i+3])-1];
i                 896 src/readfrd.c            for (i=0; i<ipuf; i++) lcase[anz->l].iexist[i] =iexist[i];
i                 898 src/readfrd.c            for(i=0; i<length; i++) SFREE(dat[i]);
i                 901 src/readfrd.c            for(i=0; i<ipuf; i++) SFREE(compName[i]);
i                 916 src/readfrd.c          for(i=0; i<lcase[anz->l].ncomps; i++)
i                 918 src/readfrd.c            if(lcase[anz->l].dat[i][node[j].nr] > lcase[anz->l].max[i])
i                 920 src/readfrd.c              lcase[anz->l].max[i]=lcase[anz->l].dat[i][node[j].nr];
i                 921 src/readfrd.c              lcase[anz->l].nmax[i]=node[j].nr;
i                 923 src/readfrd.c            if(lcase[anz->l].dat[i][node[j].nr] < lcase[anz->l].min[i])
i                 925 src/readfrd.c              lcase[anz->l].min[i]=lcase[anz->l].dat[i][node[j].nr];
i                 926 src/readfrd.c              lcase[anz->l].nmin[i]=node[j].nr;
i                 961 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
i                 963 src/readfrd.c          if ( (lcase[anz->l].compName[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
i                 965 src/readfrd.c          if ( (lcase[anz->l].icname[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
i                 967 src/readfrd.c          lcase[anz->l].max[i]=-MAX_FLOAT;
i                 968 src/readfrd.c          lcase[anz->l].min[i]=MAX_FLOAT;
i                 969 src/readfrd.c          if ( (lcase[anz->l].dat[i] = (float *)malloc( (anz->nmax+1) * sizeof(float))) == NULL )
i                 971 src/readfrd.c          for(j=0; j<=anz->nmax; j++) lcase[anz->l].dat[i][j]=0.;
i                 984 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
i                 986 src/readfrd.c          lcase[anz->l].menu[i] = 1;
i                 987 src/readfrd.c          lcase[anz->l].ictype[i] = 2;
i                 988 src/readfrd.c          lcase[anz->l].icind1[i] = i+1;
i                 989 src/readfrd.c          lcase[anz->l].icind2[i] = 0;
i                 990 src/readfrd.c          lcase[anz->l].iexist[i] = 0;
i                1012 src/readfrd.c              for(i=0; i<lcase[anz->l].ncomps; i++)
i                1013 src/readfrd.c                lcase[anz->l].dat[i][nodenr]= stof(&rec_str[8+i*12], 1, 12);
i                1017 src/readfrd.c              for(i=0; i<lcase[anz->l].ncomps; i++)
i                1018 src/readfrd.c                lcase[anz->l].dat[i][nodenr]= stof(&rec_str[13+i*12], 1, 12);
i                1025 src/readfrd.c          for(i=0; i<lcase[anz->l].ncomps; i++)
i                1027 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] > lcase[anz->l].max[i])
i                1029 src/readfrd.c              lcase[anz->l].max[i]=lcase[anz->l].dat[i][nodenr];
i                1030 src/readfrd.c              lcase[anz->l].nmax[i]=nodenr;
i                1032 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] < lcase[anz->l].min[i])
i                1034 src/readfrd.c              lcase[anz->l].min[i]=lcase[anz->l].dat[i][nodenr];
i                1035 src/readfrd.c              lcase[anz->l].nmin[i]=nodenr;
i                1069 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
i                1071 src/readfrd.c          if ( (lcase[anz->l].dat[i] = (float *)malloc( (anz->nmax+1) * sizeof(float))) == NULL )
i                1073 src/readfrd.c          if ( (lcase[anz->l].compName[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
i                1075 src/readfrd.c          if ( (lcase[anz->l].icname[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
i                1077 src/readfrd.c          lcase[anz->l].max[i]=-MAX_FLOAT;
i                1078 src/readfrd.c          lcase[anz->l].min[i]=MAX_FLOAT;
i                1079 src/readfrd.c          for(j=0; j<=anz->nmax; j++) lcase[anz->l].dat[i][j]=0.;
i                1092 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
i                1094 src/readfrd.c          lcase[anz->l].menu[i] = 1;
i                1095 src/readfrd.c          lcase[anz->l].ictype[i] = 4;
i                1096 src/readfrd.c          lcase[anz->l].iexist[i] = 0;
i                1138 src/readfrd.c                for(i=0; i<lcase[anz->l].ncomps; i++)
i                1139 src/readfrd.c                  lcase[anz->l].dat[i][nodenr]= stof(&rec_str[8+i*12], 1, 12);
i                1143 src/readfrd.c                for(i=0; i<lcase[anz->l].ncomps; i++)
i                1144 src/readfrd.c                  lcase[anz->l].dat[i][nodenr]= stof(&rec_str[13+i*12], 1, 12);
i                1152 src/readfrd.c          for(i=0; i<lcase[anz->l].ncomps; i++)
i                1154 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] > lcase[anz->l].max[i])
i                1156 src/readfrd.c              lcase[anz->l].max[i]=lcase[anz->l].dat[i][nodenr];
i                1157 src/readfrd.c              lcase[anz->l].nmax[i]=nodenr;
i                1159 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] < lcase[anz->l].min[i])
i                1161 src/readfrd.c              lcase[anz->l].min[i]=lcase[anz->l].dat[i][nodenr];
i                1162 src/readfrd.c              lcase[anz->l].nmin[i]=nodenr;
i                1196 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
i                1198 src/readfrd.c          if ( (lcase[anz->l].dat[i] = (float *)malloc( (anz->nmax+1) * sizeof(float))) == NULL )
i                1200 src/readfrd.c          if ( (lcase[anz->l].compName[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
i                1202 src/readfrd.c          if ( (lcase[anz->l].icname[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
i                1204 src/readfrd.c          lcase[anz->l].max[i]=-MAX_FLOAT;
i                1205 src/readfrd.c          lcase[anz->l].min[i]=MAX_FLOAT;
i                1206 src/readfrd.c          for(j=0; j<=anz->nmax; j++) lcase[anz->l].dat[i][j]=0.;
i                1219 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
i                1221 src/readfrd.c          lcase[anz->l].menu[i] = 1;
i                1222 src/readfrd.c          lcase[anz->l].ictype[i] = 1;
i                1223 src/readfrd.c          lcase[anz->l].icind1[i] = i+1;
i                1224 src/readfrd.c          lcase[anz->l].icind2[i] = 0;
i                1225 src/readfrd.c          lcase[anz->l].iexist[i] = 0;
i                1245 src/readfrd.c              for(i=0; i<lcase[anz->l].ncomps; i++)
i                1246 src/readfrd.c                lcase[anz->l].dat[i][nodenr]= stof(&rec_str[8+i*12], 1, 12);
i                1250 src/readfrd.c              for(i=0; i<lcase[anz->l].ncomps; i++)
i                1251 src/readfrd.c                lcase[anz->l].dat[i][nodenr]= stof(&rec_str[13+i*12], 1, 12);
i                1258 src/readfrd.c          for(i=0; i<lcase[anz->l].ncomps; i++)
i                1260 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] > lcase[anz->l].max[i])
i                1262 src/readfrd.c              lcase[anz->l].max[i]=lcase[anz->l].dat[i][nodenr];
i                1263 src/readfrd.c              lcase[anz->l].nmax[i]=nodenr;
i                1265 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] < lcase[anz->l].min[i])
i                1267 src/readfrd.c              lcase[anz->l].min[i]=lcase[anz->l].dat[i][nodenr];
i                1268 src/readfrd.c              lcase[anz->l].nmin[i]=nodenr;
i                1291 src/readfrd.c    for(i=0; i<anz_p; i++) SFREE(pheader[i]); SFREE(pheader);
i                1325 src/readfrd.c    register ITG i,j, n;
i                1359 src/readfrd.c    for(i=0; i<(lcase[lc].ncomps); i++)
i                1361 src/readfrd.c      if ( (lcase[lc].dat[i] = (float *)malloc( (anz->nmax+1) * sizeof(float))) == NULL )
i                1363 src/readfrd.c      for(j=0; j<=anz->nmax; j++) lcase[lc].dat[i][j]=0.;
i                1390 src/readfrd.c          for(i=0; i<ncomps; i++) length = frecord( handle, rec_str);
i                1399 src/readfrd.c            for(i=0; i<lcase[lc].ncomps; i++)
i                1402 src/readfrd.c              lcase[lc].dat[i][nodenr]= value[i];
i                1420 src/readfrd.c              i=6;
i                1440 src/readfrd.c            else for(i=0; i<maxcomps; i++) lcase[lc].dat[i][nodenr]= stof(&rec_str[n+i*12], 1, 12);
i                1442 src/readfrd.c          else i=0;
i                1451 src/readfrd.c            lcase[lc].dat[i][nodenr]= stof(&rec_str[n+j*12], 1, 12);
i                1452 src/readfrd.c            i++;j++;
i                1453 src/readfrd.c          }while((j<6)&&(i<lcase[lc].ncomps));
i                1459 src/readfrd.c        for(i=0; i<lcase[lc].ncomps; i++)
i                1461 src/readfrd.c          if(lcase[lc].dat[i][node[j].nr] > lcase[lc].max[i])
i                1463 src/readfrd.c            lcase[lc].max[i]=lcase[lc].dat[i][node[j].nr];
i                1464 src/readfrd.c            lcase[lc].nmax[i]=node[j].nr;
i                1466 src/readfrd.c          if(lcase[lc].dat[i][node[j].nr] < lcase[lc].min[i])
i                1468 src/readfrd.c            lcase[lc].min[i]=lcase[lc].dat[i][node[j].nr];
i                1469 src/readfrd.c            lcase[lc].nmin[i]=node[j].nr;
i                1510 src/readfrd.c    register ITG i=0,j, n;
i                1546 src/readfrd.c        for(i=0; i<ncomps; i++) *byte_offset+= frecord( handle, rec_str)+1;
i                1573 src/readfrd.c      for(i=0; i<lcase[lc].ncomps; i++)
i                1576 src/readfrd.c        dat[i]= value[i];
i                1615 src/readfrd.c          dat[i]= stof(&rec_str[n+j*12], 1, 12);
i                1616 src/readfrd.c          i++;j++;
i                1617 src/readfrd.c        }while((j<6)&&(i<lcase[lc].ncomps));
i                1626 src/readfrd.c        for (i=0; i<lcase[lc].ncomps; i++) if(lcase[lc].iexist[i]!=1) nvals++;
i                1633 src/readfrd.c          for(i=0; i<nvals/6; i++) 
i                1657 src/readfrd.c            i=6;
i                1677 src/readfrd.c          else for(i=0; i<maxcomps; i++) dat[i]= stof(&rec_str[n+i*12], 1, 12);
i                1679 src/readfrd.c        else i=0;
i                  36 src/readinput.c   ITG i,j,k,n,in=0,nlinemax=100000,irestartread,irestartstep,
i                  54 src/readinput.c   for(i=0;i<2*nentries;i++){ipoinp[i]=0;}
i                  91 src/readinput.c       i=-1;
i                  94 src/readinput.c 	      i++;
i                  95 src/readinput.c 	      if((buff[i]=='\0')||(buff[i]=='\n')||(buff[i]=='\r')||(k==1320)) break;
i                  96 src/readinput.c 	      if((buff[i]==' ')||(buff[i]=='\t')) continue;
i                  97 src/readinput.c 	      buff[k]=buff[i];
i                 102 src/readinput.c 	      i++;
i                 103 src/readinput.c 	      if((buff[i]=='\0')||(buff[i]=='\n')||(buff[i]=='\r')||(k==1320)) break;
i                 104 src/readinput.c 	      buff[k]=buff[i];
i                 198 src/readinput.c 	  for(i=0;i<n;i++){
i                 199 src/readinput.c 	      if(strcmp1(&textpart[(long long)132*i],"READ")==0){
i                 203 src/readinput.c 	      if(strcmp1(&textpart[(long long)132*i],"STEP")==0){
i                 204 src/readinput.c 		  irestartstep=atoi(&textpart[(long long)132*i+5]);
i                 362 src/readinput.c 	  for(i=0;i<n;i++){
i                 363 src/readinput.c 	      if(strcmp1(&textpart[(long long)132*i],"TYPE=TEMPERATURE")==0){
i                  41 src/remastruct.c         callfrommain,i,*irow=NULL,mt,im;
i                  70 src/remastruct.c     RENEW(irow,ITG,nzs[1]);for(i=0;i<nzs[1];i++) irow[i]=0;
i                 109 src/remastruct.c 	    RENEW(adb,double,neq[1]);for(i=0;i<neq[1];i++) adb[i]=0.;
i                 110 src/remastruct.c 	    RENEW(aub,double,nzs[1]);for(i=0;i<nzs[1];i++) aub[i]=0.;
i                  41 src/remastructar.c         callfrommain,i,*irow=NULL,mt;
i                  66 src/remastructar.c     RENEW(irow,ITG,nzs[1]);for(i=0;i<nzs[1];i++) irow[i]=0;
i                  41 src/remastructem.c         callfrommain,i,*irow=NULL,mt;
i                  70 src/remastructem.c     RENEW(irow,ITG,nzs[1]);for(i=0;i<nzs[1];i++) irow[i]=0;
i                  84 src/remastructem.c     RENEW(f,double,neq[1]);for(i=0;i<neq[1];i++) f[i]=0.;
i                  85 src/remastructem.c     RENEW(fext,double,neq[1]);for(i=0;i<neq[1];i++) fext[i]=0.;
i                  86 src/remastructem.c     RENEW(b,double,neq[1]);for(i=0;i<neq[1];i++) b[i]=0.;
i                  87 src/remastructem.c     RENEW(fini,double,neq[1]);for(i=0;i<neq[1];i++) fini[i]=0.;
i                  90 src/remastructem.c 	RENEW(aux2,double,neq[1]);for(i=0;i<neq[1];i++) aux2[i]=0.;
i                  91 src/remastructem.c 	RENEW(fextini,double,neq[1]);for(i=0;i<neq[1];i++) fextini[i]=0.;
i                  92 src/remastructem.c 	RENEW(adb,double,neq[1]);for(i=0;i<neq[1];i++) adb[i]=0.;
i                  93 src/remastructem.c 	RENEW(aub,double,nzs[1]);for(i=0;i<nzs[1];i++) aub[i]=0.;
i                  75 src/results.c          intpointvart,mt=mi[1]+1,i,j,*ithread=NULL;
i                 209 src/results.c  	for(i=0; i<num_cpus; i++)  {
i                 210 src/results.c  	    ithread[i]=i;
i                 211 src/results.c  	    pthread_create(&tid[i], NULL, (void *)resultsmechmt, (void *)&ithread[i]);
i                 213 src/results.c  	for(i=0; i<num_cpus; i++)  pthread_join(tid[i], NULL);
i                 215 src/results.c  	for(i=0;i<mt**nk;i++){
i                 216 src/results.c  	    fn[i]=fn1[i];
i                 218 src/results.c  	for(i=0;i<mt**nk;i++){
i                 220 src/results.c  		fn[i]+=fn1[i+j*mt**nk];
i                 293 src/results.c  	for(i=0; i<num_cpus; i++)  {
i                 294 src/results.c  	    ithread[i]=i;
i                 295 src/results.c  	    pthread_create(&tid[i], NULL, (void *)resultsthermmt, (void *)&ithread[i]);
i                 297 src/results.c  	for(i=0; i<num_cpus; i++)  pthread_join(tid[i], NULL);
i                 299 src/results.c  	for(i=0;i<*nk;i++){
i                 300 src/results.c  		fn[mt*i]=fn1[mt*i];
i                 302 src/results.c  	for(i=0;i<*nk;i++){
i                 304 src/results.c  		fn[mt*i]+=fn1[mt*i+j*mt**nk];
i                 356 src/results.c  void *resultsmechmt(ITG *i){
i                 360 src/results.c      indexfn=*i*mt1**nk1;
i                 361 src/results.c      indexqa=*i*3;
i                 362 src/results.c      indexnal=*i;
i                 367 src/results.c      nea=*i*nedelta+1;
i                 368 src/results.c      neb=(*i+1)*nedelta;
i                 370 src/results.c      if((*i==num_cpus-1)&&(neb<*ne1)) neb=*ne1;
i                 389 src/results.c  void *resultsthermmt(ITG *i){
i                 393 src/results.c      indexfn=*i*mt1**nk1;
i                 394 src/results.c      indexqa=*i*3;
i                 395 src/results.c      indexnal=*i;
i                 398 src/results.c      nea=*i*nedelta+1;
i                 399 src/results.c      neb=(*i+1)*nedelta;
i                 400 src/results.c      if((*i==num_cpus-1)&&(neb<*ne1)) neb=*ne1;
i                  78 src/resultsinduction.c         intpointvart,mt=mi[1]+1,i,j,*ithread=NULL,*islavsurf=NULL,
i                 200 src/resultsinduction.c 	for(i=0; i<num_cpus; i++)  {
i                 201 src/resultsinduction.c 	    ithread[i]=i;
i                 202 src/resultsinduction.c 	    pthread_create(&tid[i], NULL, (void *)resultsemmt, (void *)&ithread[i]);
i                 204 src/resultsinduction.c 	for(i=0; i<num_cpus; i++)  pthread_join(tid[i], NULL);
i                 244 src/resultsinduction.c 	for(i=0; i<num_cpus; i++)  {
i                 245 src/resultsinduction.c 	    ithread[i]=i;
i                 246 src/resultsinduction.c 	    pthread_create(&tid[i], NULL, (void *)resultsthermemmt, (void *)&ithread[i]);
i                 248 src/resultsinduction.c 	for(i=0; i<num_cpus; i++)  pthread_join(tid[i], NULL);
i                 250 src/resultsinduction.c 	for(i=0;i<*nk;i++){
i                 251 src/resultsinduction.c 		fn[mt*i]=fn1[mt*i];
i                 253 src/resultsinduction.c 	for(i=0;i<*nk;i++){
i                 255 src/resultsinduction.c 		fn[mt*i]+=fn1[mt*i+j*mt**nk];
i                 307 src/resultsinduction.c void *resultsemmt(ITG *i){
i                 312 src/resultsinduction.c     nea=*i*nedelta+1;
i                 313 src/resultsinduction.c     neb=(*i+1)*nedelta;
i                 315 src/resultsinduction.c     if((*i==num_cpus-1)&&(neb<*ne1)) neb=*ne1;
i                 326 src/resultsinduction.c void *resultsthermemmt(ITG *i){
i                 330 src/resultsinduction.c     indexfn=*i*mt1**nk1;
i                 331 src/resultsinduction.c     indexqa=*i*3;
i                 332 src/resultsinduction.c     indexnal=*i;
i                 335 src/resultsinduction.c     nea=*i*nedelta+1;
i                 336 src/resultsinduction.c     neb=(*i+1)*nedelta;
i                 337 src/resultsinduction.c     if((*i==num_cpus-1)&&(neb<*ne1)) neb=*ne1;
i                  34 src/sgi.c        ITG i,j,k,l,*pointers=NULL,method;
i                  54 src/sgi.c          for(i=*neq-1;i>=0;--i){
i                  55 src/sgi.c            for(j=0;j<icol[i];++j){
i                  59 src/sgi.c            pointers[i]=--k;
i                  60 src/sgi.c            irowsgi[k]=i;
i                  61 src/sgi.c            ausgi[k]=ad[i];
i                  66 src/sgi.c          for(i=*neq-1;i>=0;--i){
i                  67 src/sgi.c            for(j=0;j<icol[i];++j){
i                  71 src/sgi.c            pointers[i]=--k;
i                  72 src/sgi.c            irowsgi[k]=i;
i                  73 src/sgi.c            ausgi[k]=ad[i]-*sigma*adb[i];
i                 487 src/spooles.c  	    int nent,i,j;
i                 497 src/spooles.c  		for(i=0;i<*neq;i++){
i                 499 src/spooles.c  			if(fabs(ad[i*(int)*nzs+j])>1.e-20) nent++;
i                 576 src/spooles.c  		for(i=0;i<*neq;i++){
i                 578 src/spooles.c  			if(fabs(ad[i*(int)*nzs+j])>1.e-20){
i                 579 src/spooles.c  			    InpMtx_inputRealEntry(mtxA,j,i,
i                 580 src/spooles.c  						  ad[i*(int)*nzs+j]);
i                 720 src/spooles.c  		int i;
i                 724 src/spooles.c  		for (i = 0; i < size; i++) {
i                 725 src/spooles.c  			DenseMtx_setRealEntry(mtxB, i, 0, b[i]);
i                 752 src/spooles.c  		int i;
i                 753 src/spooles.c  		for (i = 0; i < size; i++) {
i                 754 src/spooles.c  			b[i] = DenseMtx_entries(mtxX)[i];
i                 807 src/spooles.c  	    int nent,i,j;
i                 817 src/spooles.c  		for(i=0;i<*neq;i++){
i                 819 src/spooles.c  			if(fabs(ad[i*(int)*nzs+j])>1.e-20) nent++;
i                 891 src/spooles.c  		for(i=0;i<*neq;i++){
i                 893 src/spooles.c  			if(fabs(ad[i*(int)*nzs+j])>1.e-20){
i                 894 src/spooles.c  			    InpMtx_inputRealEntry(mtxA,j,i,
i                 895 src/spooles.c  						  ad[i*(int)*nzs+j]);
i                1035 src/spooles.c  		int i;
i                1039 src/spooles.c  		for (i = 0; i < size; i++) {
i                1040 src/spooles.c  			DenseMtx_setRealEntry(mtxB, i, 0, b[i]);
i                1067 src/spooles.c  		int i;
i                1068 src/spooles.c  		for (i = 0; i < size; i++) {
i                1069 src/spooles.c  			b[i] = DenseMtx_entries(mtxX)[i];
i                  75 src/steadystate.c   ITG nev,i,j,k, *inum=NULL,*ipobody=NULL,inewton=0,nsectors,im,
i                 164 src/steadystate.c   for(i=neq[1];i>0;i--){
i                 165 src/steadystate.c       if(icol[i-1]>0){
i                 166 src/steadystate.c 	  *nzl=i;
i                 221 src/steadystate.c 	  for(i=0;i<*nload;i++){
i                 222 src/steadystate.c 	      iload=i+1;
i                 229 src/steadystate.c 	  for(i=0;i<*nforc;i++){
i                 230 src/steadystate.c 	      iforc=i+1;
i                 270 src/steadystate.c   for(i=0;i<nmddof;i++){imddof[i]-=1;}
i                 296 src/steadystate.c 	  for(i=0;i<nev;i++){
i                 297 src/steadystate.c 	      if(d[i]<0){d[i]=0.;}
i                 402 src/steadystate.c 	      for(i=nev;i<nev+nevd;i++){
i                 403 src/steadystate.c 		  if(d[i]<0){d[i]=0.;}
i                 412 src/steadystate.c 	  for(i=nev;i<nev+nevd;i++){nm[i]=nmd;}
i                 463 src/steadystate.c       for(i=0;i<*mcs;i++){
i                 464 src/steadystate.c 	  if(cs[17*i]>nsectors) nsectors=(ITG)(cs[17*i]+0.5);
i                 482 src/steadystate.c       for(i=3**nk;i<3**nk*nsectors;i++){co[i]=0.;}
i                 493 src/steadystate.c       for(i=*ne;i<*ne*nsectors;i++){ipkon[i]=-1;}
i                 514 src/steadystate.c       for(i=0;i<*nmpc;i++){ipompcold[i]=ipompc[i];}
i                 515 src/steadystate.c       for(i=0;i<3**mpcend;i++){nodempcold[i]=nodempc[i];}
i                 516 src/steadystate.c       for(i=0;i<*mpcend;i++){coefmpcold[i]=coefmpc[i];}
i                 517 src/steadystate.c       for(i=0;i<20**nmpc;i++){labmpcold[i]=labmpc[i];}
i                 518 src/steadystate.c       for(i=0;i<*nmpc;i++){ikmpcold[i]=ikmpc[i];}
i                 519 src/steadystate.c       for(i=0;i<*nmpc;i++){ilmpcold[i]=ilmpc[i];}
i                 586 src/steadystate.c   for(i=0;i<*ne;i++){
i                 587 src/steadystate.c       if(ipkon[i]==-1) continue;
i                 588 src/steadystate.c       if(strcmp1(&lakon[i*8],"ED")==0){
i                 638 src/steadystate.c       for(i=0;i<*nprint;i++){
i                 639 src/steadystate.c 	  if((strcmp1(&prlab[6*i],"S")==0)||
i                 640 src/steadystate.c 	     (strcmp1(&prlab[6*i],"E")==0)||
i                 641 src/steadystate.c 	     (strcmp1(&prlab[6*i],"PEEQ")==0)||
i                 642 src/steadystate.c 	     (strcmp1(&prlab[6*i],"ENER")==0)||
i                 643 src/steadystate.c 	     (strcmp1(&prlab[6*i],"SDV")==0)||
i                 644 src/steadystate.c 	     (strcmp1(&prlab[6*i],"RF")==0)) {intpointvar=1;break;}
i                 652 src/steadystate.c       for(i=0;i<*nprint;i++){
i                 653 src/steadystate.c 	  if((strcmp1(&prlab[6*i],"HFL")==0)||
i                 654 src/steadystate.c 	     (strcmp1(&prlab[6*i],"RFL")==0)) {intpointvar=1;break;}
i                 702 src/steadystate.c 	      for(i=0;i<nev;i++){e[i]=sqrt(d[i]);}
i                 707 src/steadystate.c 	      for(i=0;i<nev;i++){
i                 708 src/steadystate.c 		  e[i]=sqrt(sqrt(d[2*i]*d[2*i]+d[2*i+1]*d[2*i+1])+d[2*i])/sqrt(2.);
i                 713 src/steadystate.c 	  for(i=0;i<nev;i++){
i                 714 src/steadystate.c 	      if(i!=0){
i                 715 src/steadystate.c 		  if(fabs(e[i]-e[i-1])<1.e-5){continue;}
i                 717 src/steadystate.c 	      if(e[i]>=fmin){
i                 718 src/steadystate.c 		  if(e[i]<=fmax){
i                 722 src/steadystate.c                                  (freq[ndatatot]+e[i])/2.;}
i                 724 src/steadystate.c 			      freq[ndatatot+j]=(freq[ndatatot]+e[i])/2.+
i                 725 src/steadystate.c 				  (e[i]-freq[ndatatot])*pow(fabs(y),1./bias)*
i                 753 src/steadystate.c       for(i=0;i<*nboun;i++){
i                 754 src/steadystate.c 	  if(fabs(xboun[i])>1.e-10){
i                 769 src/steadystate.c 	  for(i=0;i<nev;i++){
i                 771 src/steadystate.c 		  if(sqrt(d[i])>(1.e-10)){
i                 772 src/steadystate.c 		      fric[i]=(alpham+betam*d[i]);
i                 777 src/steadystate.c 		      fric[i]=0.;
i                 780 src/steadystate.c 		  fric[2*i]=sqrt(sqrt(d[2*i]*d[2*i]+d[2*i+1]*d[2*i+1])+d[2*i])/
i                 782 src/steadystate.c 		  fric[2*i+1]=sqrt(sqrt(d[2*i]*d[2*i]+d[2*i+1]*d[2*i+1])-d[2*i])/
i                 784 src/steadystate.c 		  if(d[2*i+1]<0.) fric[2*i+1]=-fric[2*i+1];
i                 785 src/steadystate.c 		  fric[2*i]=alpham+betam*fric[2*i];
i                 786 src/steadystate.c 		  fric[2*i+1]=betam*fric[2*i+1];
i                 805 src/steadystate.c 	  for(i=0; i<imax; i++){
i                 807 src/steadystate.c 		  fric[i]=2.*sqrt(d[i])*xmodal[11+i];    
i                 809 src/steadystate.c 		  fric[2*i]=sqrt(sqrt(d[2*i]*d[2*i]+d[2*i+1]*d[2*i+1])+d[2*i])/
i                 811 src/steadystate.c 		  fric[2*i+1]=sqrt(sqrt(d[2*i]*d[2*i]+d[2*i+1]*d[2*i+1])-d[2*i])/
i                 813 src/steadystate.c 		  if(d[2*i+1]<0.) fric[2*i+1]=-fric[2*i+1];
i                 814 src/steadystate.c 		  fric[2*i]=2.*fric[2*i]*xmodal[11+i];
i                 815 src/steadystate.c 		  fric[2*i+1]=2.*fric[2*i+1]*xmodal[11+i];
i                 824 src/steadystate.c       for (i=0;i<*nforc;i++){
i                 825 src/steadystate.c 	  if(nodeforc[2*i+1]>=nsectors){
i                 826 src/steadystate.c 	      iphaseforc[i]=1;
i                 831 src/steadystate.c       for (i=0;i<*nload;i++){
i                 832 src/steadystate.c 	  if(nelemload[2*i+1]>=nsectors){
i                 833 src/steadystate.c 	      iphaseload[i]=1;
i                 839 src/steadystate.c 	  for (i=0;i<*nboun;i++){
i                 840 src/steadystate.c 	      if(nodeboun[i]>*nk){
i                 841 src/steadystate.c 		  iphaseboun[i]=1;
i                 842 src/steadystate.c 		  nodeboun[i]=nodeboun[i]-*nk;
i                 957 src/steadystate.c 	  for(i=0;i<6*mi[0]**ne;i++){eme[i]=0.;}
i                 984 src/steadystate.c 	      for(i=0;i<nev;i++){
i                 985 src/steadystate.c 		  FORTRAN(op,(&neq[1],&z[(long long)i*neq[1]],&zc[i*neq[1]],adc,auc,
i                 992 src/steadystate.c 	      for(i=0;i<nev*nev;i++){cc[i]=0.;}
i                 993 src/steadystate.c 	      for(i=0;i<nev;i++){
i                 994 src/steadystate.c 		  for(j=0;j<=i;j++){
i                 996 src/steadystate.c 			  cc[i*nev+j]+=z[(long long)j*neq[1]+k]*zc[i*neq[1]+k];
i                1003 src/steadystate.c 	      for(i=0;i<nev;i++){
i                1004 src/steadystate.c 		  for(j=i;j<nev;j++){
i                1005 src/steadystate.c 		      cc[i*nev+j]=cc[j*nev+i];
i                1026 src/steadystate.c 	  for (i=0;i<*nforc;i++){
i                1027 src/steadystate.c 	      xforcr[i]=xforcact[i]*(1-iphaseforc[i]);
i                1030 src/steadystate.c 	  for (i=0;i<*nload;i++){
i                1032 src/steadystate.c 		  xloadr[2*i+j]=xloadact[2*i+j]*(1-iphaseload[i]);
i                1036 src/steadystate.c 	  for(i=0;i<*nbody;i++){
i                1038 src/steadystate.c 		  xbodyr[7*i+j]=xbodyact[7*i+j];
i                1040 src/steadystate.c 	      if(ibody[3*i+2]==2){
i                1041 src/steadystate.c 		  xbodyr[7*i]=0.;
i                1065 src/steadystate.c 		      for(i=0;i<nev;i++){
i                1066 src/steadystate.c 			  i2=(long long)i*neq[1];
i                1067 src/steadystate.c 			  aa[i]=0.;
i                1070 src/steadystate.c 				  aa[i]+=z[i2+ikactmechr[j]]*br[ikactmechr[j]];
i                1074 src/steadystate.c 				  aa[i]+=z[i2+j]*br[j];
i                1079 src/steadystate.c 		      for(i=0;i<nev;i++){aa[i]=0.;}
i                1081 src/steadystate.c 			  for(i=0;i<nev;i++){
i                1085 src/steadystate.c 				      aa[i]+=z[(long long)i*nzdof+id-1]*br[ikactmechr[j]];
i                1093 src/steadystate.c 		      for(i=0;i<nev;i++){
i                1094 src/steadystate.c 			  aa[2*i]=0.;aa[2*i+1]=0.;
i                1097 src/steadystate.c 				  aa[2*i]+=z[(long long)2*i*neq[1]+ikactmechr[j]]*br[ikactmechr[j]];
i                1098 src/steadystate.c 				  aa[2*i+1]+=z[(long long)(2*i+1)*neq[1]+ikactmechr[j]]*br[ikactmechr[j]];
i                1102 src/steadystate.c 				  aa[2*i]+=z[(long long)2*i*neq[1]+j]*br[j];
i                1103 src/steadystate.c 				  aa[2*i+1]+=z[(long long)(2*i+1)*neq[1]+j]*br[j];
i                1108 src/steadystate.c 		      for(i=0;i<nev;i++){aa[2*i]=0.;aa[2*i+1]=0.;}
i                1110 src/steadystate.c 			  for(i=0;i<nev;i++){
i                1114 src/steadystate.c 				      aa[i]+=z[(long long)2*i*nzdof+id-1]*br[ikactmechr[j]];
i                1115 src/steadystate.c 				      aa[i]+=z[(long long)(2*i+1)*nzdof+id-1]*br[ikactmechr[j]];
i                1133 src/steadystate.c 	      for (i=0;i<*nboun;i++){
i                1134 src/steadystate.c 		  xbounr[i]=xbounact[i]*(1-iphaseboun[i]);
i                1138 src/steadystate.c 	      for(i=0;i<*nboun;i++){
i                1139 src/steadystate.c 		  ic=neq[1]+i;
i                1142 src/steadystate.c 		      fr[ir]=fr[ir]-au[j]*xbounr[i];
i                1171 src/steadystate.c 	  for (i=0;i<*nforc;i++){
i                1172 src/steadystate.c 	      xforci[i]=xforcact[i]*iphaseforc[i];
i                1175 src/steadystate.c 	  for (i=0;i<*nload;i++){
i                1177 src/steadystate.c 		  xloadi[2*i+j]=xloadact[2*i+j]*iphaseload[i];
i                1181 src/steadystate.c 	  for(i=0;i<*nbody;i++){
i                1183 src/steadystate.c 		  xbodyi[7*i+j]=xbodyact[7*i+j];
i                1185 src/steadystate.c 	      if(ibody[3*i+2]==1){
i                1186 src/steadystate.c 		  xbodyi[7*i]=0.;
i                1210 src/steadystate.c 		      for(i=0;i<nev;i++){
i                1211 src/steadystate.c 			  i2=(long long)i*neq[1];
i                1212 src/steadystate.c 			  bb[i]=0.;
i                1215 src/steadystate.c 				  bb[i]+=z[i2+ikactmechi[j]]*bi[ikactmechi[j]];
i                1219 src/steadystate.c 				  bb[i]+=z[i2+j]*bi[j];
i                1224 src/steadystate.c 		      for(i=0;i<nev;i++){bb[i]=0.;}
i                1226 src/steadystate.c 			  for(i=0;i<nev;i++){
i                1230 src/steadystate.c 				      bb[i]+=z[(long long)i*nzdof+id-1]*bi[ikactmechi[j]];
i                1238 src/steadystate.c 		      for(i=0;i<nev;i++){
i                1239 src/steadystate.c 			  bb[2*i]=0.;bb[2*i+1]=0.;
i                1242 src/steadystate.c 				  bb[2*i]+=z[(long long)2*i*neq[1]+ikactmechi[j]]*bi[ikactmechi[j]];
i                1243 src/steadystate.c 				  bb[2*i+1]+=z[(long long)(2*i+1)*neq[1]+ikactmechi[j]]*bi[ikactmechi[j]];
i                1247 src/steadystate.c 				  bb[2*i]+=z[(long long)2*i*neq[1]+j]*bi[j];
i                1248 src/steadystate.c 				  bb[2*i+1]+=z[(long long)(2*i+1)*neq[1]+j]*bi[j];
i                1253 src/steadystate.c 		      for(i=0;i<nev;i++){bb[2*i]=0.;bb[2*i+1]=0.;}
i                1255 src/steadystate.c 			  for(i=0;i<nev;i++){
i                1259 src/steadystate.c 				      bb[2*i]+=z[(long long)2*i*nzdof+id-1]*bi[ikactmechi[j]];
i                1260 src/steadystate.c 				      bb[2*i+1]+=z[(long long)(2*i+1)*nzdof+id-1]*bi[ikactmechi[j]];
i                1279 src/steadystate.c 	      for (i=0;i<*nboun;i++){
i                1280 src/steadystate.c 		  xbouni[i]=xbounact[i]*iphaseboun[i];
i                1284 src/steadystate.c 	      for(i=0;i<*nboun;i++){
i                1285 src/steadystate.c 		  ic=neq[1]+i;
i                1288 src/steadystate.c 		      fi[ir]=fi[ir]-au[j]*xbouni[i];
i                1316 src/steadystate.c 	      for(i=0;i<neq[1];i++){
i                1317 src/steadystate.c 		  br[i]+=freq[l]*(freq[l]*mubr[i]+alpham*mubi[i]+betam*fi[i]);
i                1318 src/steadystate.c 		  bi[i]+=freq[l]*(freq[l]*mubi[i]-alpham*mubr[i]-betam*fr[i]);
i                1323 src/steadystate.c 	      for(i=0;i<nev;i++){
i                1324 src/steadystate.c 		  aa[i]=0.;
i                1326 src/steadystate.c 		      aa[i]+=z[(long long)i*neq[1]+j]*br[j];
i                1330 src/steadystate.c 	      for(i=0;i<nev;i++){
i                1331 src/steadystate.c 		  bb[i]=0.;
i                1333 src/steadystate.c 		      bb[i]+=z[(long long)i*neq[1]+j]*bi[j];
i                1343 src/steadystate.c 		  for(i=0;i<nev;i++){
i                1344 src/steadystate.c 		      dd=pow(d[i]-pow(freq[l],2),2)+
i                1345 src/steadystate.c 			  pow(fric[i],2)*pow(freq[l],2);
i                1346 src/steadystate.c 		      bjr[i]=(aa[i]*(d[i]-freq[l]*freq[l])+
i                1347 src/steadystate.c 			      bb[i]*fric[i]*freq[l])/dd;
i                1348 src/steadystate.c 		      bji[i]=(bb[i]*(d[i]-freq[l]*freq[l])-
i                1349 src/steadystate.c 			      aa[i]*fric[i]*freq[l])/dd;
i                1358 src/steadystate.c 		  for(i=0;i<nev2;i++){
i                1360 src/steadystate.c 			  am[i*nev2+j]=0.;
i                1362 src/steadystate.c 		      bm[i]=0.;
i                1364 src/steadystate.c 		  for(i=0;i<nev;i++){
i                1365 src/steadystate.c 		      am[i*nev2+i]=d[i]-freq[l]*freq[l];
i                1366 src/steadystate.c 		      am[(i+nev)*nev2+i]=-fric[i]*freq[l];
i                1367 src/steadystate.c 		      bm[i]=aa[i];
i                1368 src/steadystate.c 		      am[i*nev2+nev+i]=-am[(i+nev)*nev2+i];
i                1369 src/steadystate.c 		      am[(i+nev)*nev2+nev+i]=am[i*nev2+i];
i                1370 src/steadystate.c 		      bm[nev+i]=bb[i];
i                1372 src/steadystate.c 			  am[(j+nev)*nev2+i]=am[(j+nev)*nev2+i]
i                1373 src/steadystate.c 			      -cc[i*nev+j]*freq[l];
i                1374 src/steadystate.c 			  am[j*nev2+nev+i]=am[j*nev2+nev+i]
i                1375 src/steadystate.c 			      +cc[i*nev+j]*freq[l];
i                1390 src/steadystate.c 		  for(i=0;i<nev;i++){
i                1391 src/steadystate.c 		      bjr[i]=bm[i];
i                1392 src/steadystate.c 		      bji[i]=bm[nev+i];
i                1405 src/steadystate.c 	      for(i=0;i<nev;i++){
i                1406 src/steadystate.c 		  ax[i]=-pow(freq[l],2)-freq[l]*fric[2*i+1]+d[2*i];
i                1407 src/steadystate.c 		  bx[i]=-freq[l]*fric[2*i]-d[2*i+1];
i                1409 src/steadystate.c 	      for(i=0;i<nev;i++){
i                1411 src/steadystate.c 		      am[j*nev2+i]=xmr[j*nev+i]*ax[j]+xmi[j*nev+i]*bx[j];
i                1412 src/steadystate.c 		      am[(j+nev)*nev2+i]=xmr[j*nev+i]*bx[j]-xmi[j*nev+i]*bx[j];
i                1413 src/steadystate.c 		      am[j*nev2+nev+i]=xmi[j*nev+i]*ax[j]-xmr[j*nev+i]*bx[j];
i                1414 src/steadystate.c 		      am[(j+nev)*nev2+nev+i]=xmi[j*nev+i]*bx[j]+xmr[j*nev+i]*ax[j];
i                1416 src/steadystate.c 		  bm[i]=aa[i]-bb[nev+i];
i                1417 src/steadystate.c 		  bm[nev+i]=bb[i]+aa[nev+i];
i                1432 src/steadystate.c 	      for(i=0;i<nev;i++){
i                1433 src/steadystate.c 		  bjr[i]=bm[i];
i                1434 src/steadystate.c 		  bji[i]=bm[nev+i];
i                1444 src/steadystate.c 	      for(i=0;i<nev;i++){
i                1445 src/steadystate.c 		  eig[i]=sqrt(d[i]);
i                1450 src/steadystate.c 	      for(i=0;i<nev;i++){
i                1451 src/steadystate.c 		  eig[2*i]=sqrt(sqrt(d[2*i]*d[2*i]+d[2*i+1]*d[2*i+1])+d[2*i])
i                1453 src/steadystate.c 		  eig[2*i+1]=sqrt(sqrt(d[2*i]*d[2*i]+d[2*i+1]*d[2*i+1])-d[2*i])
i                1455 src/steadystate.c 		  if(d[2*i+1]<0.) eig[2*i+1]=-eig[2*i+1];
i                1469 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                1470 src/steadystate.c 		      br[imddof[i]]=ubr[imddof[i]];
i                1478 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                1479 src/steadystate.c 		      br[imddof[i]]=0.;
i                1486 src/steadystate.c 		  for(i=0;i<neq[1];i++){
i                1488 src/steadystate.c 			  br[i]+=bjr[j]*z[(long long)j*neq[1]+i];
i                1492 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                1494 src/steadystate.c 			  br[imddof[i]]+=bjr[j]*z[(long long)j*neq[1]+imddof[i]];
i                1499 src/steadystate.c 	      for(i=0;i<nmddof;i++){
i                1500 src/steadystate.c 		  FORTRAN(nident,(izdof,&imddof[i],&nzdof,&id));
i                1502 src/steadystate.c 		      if(izdof[id-1]==imddof[i]){
i                1504 src/steadystate.c 			      br[imddof[i]]+=bjr[j]*z[(long long)j*nzdof+id-1];
i                1515 src/steadystate.c 		  i=imdnode[jj]-1;
i                1517 src/steadystate.c 		      vr[mt*i+j]=0.;
i                1528 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                1529 src/steadystate.c 		      bi[imddof[i]]=ubi[imddof[i]];
i                1537 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                1538 src/steadystate.c 		      bi[imddof[i]]=0.;
i                1545 src/steadystate.c 		  for(i=0;i<neq[1];i++){
i                1547 src/steadystate.c 			  bi[i]+=bji[j]*z[(long long)j*neq[1]+i];
i                1551 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                1553 src/steadystate.c 			  bi[imddof[i]]+=bji[j]*z[(long long)j*neq[1]+imddof[i]];
i                1558 src/steadystate.c 	      for(i=0;i<nmddof;i++){
i                1559 src/steadystate.c 		  FORTRAN(nident,(izdof,&imddof[i],&nzdof,&id));
i                1561 src/steadystate.c 		      if(izdof[id-1]==imddof[i]){
i                1563 src/steadystate.c 			      bi[imddof[i]]+=bji[j]*z[(long long)j*nzdof+id-1];
i                1575 src/steadystate.c 		  i=imdnode[jj]-1;
i                1577 src/steadystate.c 		      vi[mt*i+j]=0.;
i                1644 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                1645 src/steadystate.c 		      br[imddof[i]]=0.;
i                1730 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                1731 src/steadystate.c 		      bi[imddof[i]]=0.;
i                1746 src/steadystate.c 		      for(i=0;i<*nk;i++){
i                1748 src/steadystate.c 			      vreal=vr[mt*i+j];
i                1749 src/steadystate.c 			      va[mt*i+j]=sqrt(vr[mt*i+j]*vr[mt*i+j]+vi[mt*i+j]*vi[mt*i+j]);
i                1751 src/steadystate.c 				  if(vi[mt*i+j]>0.){vp[mt*i+j]=90.;}
i                1752 src/steadystate.c 				  else{vp[mt*i+j]=-90.;}
i                1755 src/steadystate.c 				  vp[mt*i+j]=atan(vi[mt*i+j]/vreal)*constant;
i                1756 src/steadystate.c 				  if(vreal<0.) vp[mt*i+j]+=180.;
i                1762 src/steadystate.c 			  i=imdnode[jj]-1;
i                1764 src/steadystate.c 			      vreal=vr[mt*i+j];
i                1765 src/steadystate.c 			      va[mt*i+j]=sqrt(vr[mt*i+j]*vr[mt*i+j]+vi[mt*i+j]*vi[mt*i+j]);
i                1767 src/steadystate.c 				  if(vi[mt*i+j]>0.){vp[mt*i+j]=90.;}
i                1768 src/steadystate.c 				  else{vp[mt*i+j]=-90.;}
i                1771 src/steadystate.c 				  vp[mt*i+j]=atan(vi[mt*i+j]/vreal)*constant;
i                1772 src/steadystate.c 				  if(vreal<0.) vp[mt*i+j]+=180.;
i                1780 src/steadystate.c 		      for(i=0;i<*nk;i++){
i                1781 src/steadystate.c 			  vreal=vr[mt*i];
i                1782 src/steadystate.c 			  va[mt*i]=sqrt(vr[mt*i]*vr[mt*i]+vi[mt*i]*vi[mt*i]);
i                1784 src/steadystate.c 			      if(vi[mt*i]>0){vp[mt*i]=90.;}
i                1785 src/steadystate.c 			      else{vp[mt*i]=-90.;}
i                1788 src/steadystate.c 			      vp[mt*i]=atan(vi[mt*i]/vreal)*constant;
i                1789 src/steadystate.c 			      if(vreal<0.) vp[mt*i]+=180.;
i                1794 src/steadystate.c 			  i=imdnode[jj]-1;
i                1795 src/steadystate.c 			  vreal=vr[mt*i];
i                1796 src/steadystate.c 			  va[mt*i]=sqrt(vr[mt*i]*vr[mt*i]+vi[mt*i]*vi[mt*i]);
i                1798 src/steadystate.c 			      if(vi[mt*i]>0){vp[mt*i]=90.;}
i                1799 src/steadystate.c 			      else{vp[mt*i]=-90.;}
i                1802 src/steadystate.c 			      vp[mt*i]=atan(vi[mt*i]/vreal)*constant;
i                1803 src/steadystate.c 			      if(vreal<0.) vp[mt*i]+=180.;
i                1818 src/steadystate.c 		      for(i=0;i<*nk;i++){
i                1820 src/steadystate.c 			      vreal=stnr[6*i+j];
i                1821 src/steadystate.c 			      stna[6*i+j]=sqrt(stnr[6*i+j]*stnr[6*i+j]+stni[6*i+j]*stni[6*i+j]);
i                1823 src/steadystate.c 				  if(stni[6*i+j]>0.){stnp[6*i+j]=90.;}
i                1824 src/steadystate.c 				  else{stnp[6*i+j]=-90.;}
i                1827 src/steadystate.c 				  stnp[6*i+j]=atan(stni[6*i+j]/vreal)*constant;
i                1828 src/steadystate.c 				  if(vreal<0.) stnp[6*i+j]+=180.;
i                1834 src/steadystate.c 			  i=imdnode[jj]-1;
i                1836 src/steadystate.c 			      vreal=stnr[6*i+j];
i                1837 src/steadystate.c 			      stna[6*i+j]=sqrt(stnr[6*i+j]*stnr[6*i+j]+stni[6*i+j]*stni[6*i+j]);
i                1839 src/steadystate.c 				  if(stni[6*i+j]>0.){stnp[6*i+j]=90.;}
i                1840 src/steadystate.c 				  else{stnp[6*i+j]=-90.;}
i                1843 src/steadystate.c 				  stnp[6*i+j]=atan(stni[6*i+j]/vreal)*constant;
i                1844 src/steadystate.c 				  if(vreal<0.) stnp[6*i+j]+=180.;
i                1883 src/steadystate.c 	  for (i=0;i<*nboun;i++){
i                1884 src/steadystate.c 	      if(iphaseboun[i]==1){
i                1885 src/steadystate.c 		  nodeboun[i]=nodeboun[i]+*nk;
i                1970 src/steadystate.c       for(i=0;i<*nboun;i++){
i                1971 src/steadystate.c 	  if(fabs(xboun[i])>1.e-10){
i                1985 src/steadystate.c 	  for(i=0;i<nev;i++){
i                1986 src/steadystate.c 	      if(sqrt(d[i])>(1.e-10)){
i                1987 src/steadystate.c 		  fric[i]=(alpham+betam*d[i]);
i                1992 src/steadystate.c 		  fric[i]=0.;
i                2011 src/steadystate.c 	  for(i=0; i<imax; i++){
i                2012 src/steadystate.c 	      fric[i]=2.*sqrt(d[i])*xmodal[11+i];     
i                2039 src/steadystate.c       for(i=0;i<l**nboun;i++){xbounacttime[i]=xbounact[i];}
i                2044 src/steadystate.c       for(i=0;i<*nforc;i++){
i                2046 src/steadystate.c 	      r[l]=xforcact[l**nforc+i];
i                2051 src/steadystate.c 	      xforcact[l**nforc+i]=r[l]/nfour*2.;
i                2053 src/steadystate.c 	  xforcact[i]=xforcact[i]/2.;
i                2056 src/steadystate.c       for(i=0;i<*nload;i++){
i                2058 src/steadystate.c 	      r[l]=xloadact[l*2**nload+2*i];
i                2063 src/steadystate.c 	      xloadact[l*2**nload+2*i]=r[l]/nfour*2.;
i                2065 src/steadystate.c 	  xloadact[2*i]=xloadact[2*i]/2.;
i                2068 src/steadystate.c       for(i=0;i<*nbody;i++){
i                2070 src/steadystate.c 	      r[l]=xbodyact[l**nbody+7*i];
i                2075 src/steadystate.c 	      xbodyact[l**nbody+7*i]=r[l]/nfour*2.;
i                2077 src/steadystate.c 	  xbodyact[7*i]=xbodyact[7*i]/2.;
i                2081 src/steadystate.c 	  for(i=0;i<*nboun;i++){
i                2083 src/steadystate.c 		  r[l]=xbounact[l**nboun+i];
i                2088 src/steadystate.c 		  xbounact[l**nboun+i]=r[l]/nfour*2.;
i                2090 src/steadystate.c 	      xbounact[i]=xbounact[i]/2.;
i                2149 src/steadystate.c 	  for(i=0;i<nev;i++){e[i]=sqrt(d[i]);}
i                2152 src/steadystate.c 	  for(i=0;i<nev;i++){
i                2153 src/steadystate.c 	      if(i!=0){
i                2154 src/steadystate.c 		  if(fabs(e[i]-e[i-1])<1.e-5){continue;}
i                2156 src/steadystate.c 	      if(e[i]>=fmin){
i                2157 src/steadystate.c 		  if(e[i]<=fmax){
i                2161 src/steadystate.c                                  (freqnh[ndatatot]+e[i])/2.;}
i                2163 src/steadystate.c 			      freqnh[ndatatot+j]=(freqnh[ndatatot]+e[i])/2.+
i                2164 src/steadystate.c 				  (e[i]-freqnh[ndatatot])*pow(fabs(y),1./bias)
i                2189 src/steadystate.c 	  for(i=0;i<6*mi[0]**ne;i++){eme[i]=0.;}
i                2260 src/steadystate.c 		  for(i=0;i<nev;i++){
i                2261 src/steadystate.c 		      FORTRAN(op,(&neq[1],&z[(long long)i*neq[1]],&zc[i*neq[1]],
i                2268 src/steadystate.c 		  for(i=0;i<nev*nev;i++){cc[i]=0.;}
i                2269 src/steadystate.c 		  for(i=0;i<nev;i++){
i                2270 src/steadystate.c 		      for(j=0;j<=i;j++){
i                2272 src/steadystate.c 			      cc[i*nev+j]+=z[(long long)j*neq[1]+k]*zc[i*neq[1]+k];
i                2279 src/steadystate.c 		  for(i=0;i<nev;i++){
i                2280 src/steadystate.c 		      for(j=i;j<nev;j++){
i                2281 src/steadystate.c 			  cc[i*nev+j]=cc[j*nev+i];
i                2289 src/steadystate.c 	      for(i=0;i<*nforc;i++){
i                2290 src/steadystate.c 		  xforcr[i]=xforcact[l**nforc+i];
i                2293 src/steadystate.c 	      for(i=0;i<*nload;i++){
i                2294 src/steadystate.c 		  xloadr[2*i]=xloadact[l*2**nload+2*i];
i                2297 src/steadystate.c 	      for(i=0;i<*nbody;i++){
i                2298 src/steadystate.c 		  xbodyr[7*i]=xbodyact[l**nbody+7*i];
i                2319 src/steadystate.c 		      for(i=0;i<nev;i++){
i                2320 src/steadystate.c 			  i2=(long long)i*neq[1];
i                2321 src/steadystate.c 			  aa[i]=0.;
i                2324 src/steadystate.c 				  aa[i]+=z[i2+ikactmech[j]]*br[ikactmech[j]];
i                2328 src/steadystate.c 				  aa[i]+=z[i2+j]*br[j];
i                2333 src/steadystate.c 		      for(i=0;i<nev;i++){aa[i]=0.;}
i                2335 src/steadystate.c 			  for(i=0;i<nev;i++){
i                2339 src/steadystate.c 				      aa[i]+=z[(long long)i*nzdof+id-1]*br[ikactmech[j]];
i                2348 src/steadystate.c 		  for(i=0;i<nev;i++){
i                2349 src/steadystate.c 		      bb[i]=0.;
i                2361 src/steadystate.c 		  for(i=0;i<neq[1];i++){bi[i]=0.;}
i                2363 src/steadystate.c 		  for(i=0;i<*nboun;i++){
i                2364 src/steadystate.c 		      xbounr[i]=xbounact[l**nboun+i];
i                2368 src/steadystate.c 		  for(i=0;i<*nboun;i++){
i                2369 src/steadystate.c 		      ic=neq[1]+i;
i                2372 src/steadystate.c 			  fr[ir]=fr[ir]-au[j]*xbounr[i];
i                2398 src/steadystate.c 		  for(i=0;i<neq[1];i++){
i                2399 src/steadystate.c 		      br[i]+=freq[l]*(freq[l]*mubr[i]);
i                2400 src/steadystate.c 		      bi[i]+=freq[l]*(-alpham*mubr[i]-betam*fr[i]);
i                2405 src/steadystate.c 		  for(i=0;i<nev;i++){
i                2406 src/steadystate.c 		      aa[i]=0.;
i                2408 src/steadystate.c 			  aa[i]+=z[(long long)i*neq[1]+j]*br[j];
i                2412 src/steadystate.c 		  for(i=0;i<nev;i++){
i                2413 src/steadystate.c 		      bb[i]=0.;
i                2415 src/steadystate.c 			  bb[i]+=z[(long long)i*neq[1]+j]*bi[j];
i                2423 src/steadystate.c 		  for(i=0;i<nev;i++){
i                2424 src/steadystate.c 		      dd=pow(d[i]-pow(freq[l],2),2)+
i                2425 src/steadystate.c 			  pow(fric[i],2)*pow(freq[l],2);
i                2426 src/steadystate.c 		      bjr[i]=(aa[i]*(d[i]-freq[l]*freq[l])+
i                2427 src/steadystate.c 			      bb[i]*fric[i]*freq[l])/dd;
i                2428 src/steadystate.c 		      bji[i]=(bb[i]*(d[i]-freq[l]*freq[l])-
i                2429 src/steadystate.c 			      aa[i]*fric[i]*freq[l])/dd;
i                2437 src/steadystate.c 		  for(i=0;i<nev2;i++){
i                2439 src/steadystate.c 			  am[i*nev2+j]=0.;
i                2441 src/steadystate.c 		      bm[i]=0.;
i                2443 src/steadystate.c 		  for(i=0;i<nev;i++){
i                2444 src/steadystate.c 		      am[i*nev2+i]=d[i]-freq[l]*freq[l];
i                2445 src/steadystate.c 		      am[(i+nev)*nev2+i]=-fric[i]*freq[l];
i                2446 src/steadystate.c 		      bm[i]=aa[i];
i                2447 src/steadystate.c 		      am[i*nev2+nev+i]=-am[(i+nev)*nev2+i];
i                2448 src/steadystate.c 		      am[(i+nev)*nev2+nev+i]=am[i*nev2+i];
i                2449 src/steadystate.c 		      bm[nev+i]=bb[i];
i                2451 src/steadystate.c 			  am[(j+nev)*nev2+i]=am[(j+nev)*nev2+i]
i                2452 src/steadystate.c 			      -cc[i*nev+j]*freq[l];
i                2453 src/steadystate.c 			  am[j*nev2+nev+i]=am[j*nev2+nev+i]
i                2454 src/steadystate.c 			      +cc[i*nev+j]*freq[l];
i                2469 src/steadystate.c 		  for(i=0;i<nev;i++){
i                2470 src/steadystate.c 		      bjr[i]=bm[i];
i                2471 src/steadystate.c 		      bji[i]=bm[nev+i];
i                2483 src/steadystate.c 		      for(i=0;i<nmddof;i++){
i                2484 src/steadystate.c 			  br[imddof[i]]=ubr[imddof[i]];
i                2492 src/steadystate.c 		      for(i=0;i<nmddof;i++){
i                2493 src/steadystate.c 			  br[imddof[i]]=0.;
i                2500 src/steadystate.c 		      for(i=0;i<neq[1];i++){
i                2502 src/steadystate.c 			      br[i]+=bjr[j]*z[(long long)j*neq[1]+i];
i                2506 src/steadystate.c 		      for(i=0;i<nmddof;i++){
i                2508 src/steadystate.c 			      br[imddof[i]]+=bjr[j]*z[(long long)j*neq[1]+imddof[i]];
i                2513 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                2514 src/steadystate.c 		      FORTRAN(nident,(izdof,&imddof[i],&nzdof,&id));
i                2516 src/steadystate.c 			  if(izdof[id-1]==imddof[i]){
i                2518 src/steadystate.c 				  br[imddof[i]]+=bjr[j]*z[(long long)j*nzdof+id-1];
i                2530 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                2531 src/steadystate.c 		      bi[imddof[i]]=0.;
i                2537 src/steadystate.c 		      for(i=0;i<neq[1];i++){
i                2539 src/steadystate.c 			      bi[i]+=bji[j]*z[(long long)j*neq[1]+i];
i                2543 src/steadystate.c 		      for(i=0;i<nmddof;i++){
i                2545 src/steadystate.c 			      bi[imddof[i]]+=bji[j]*z[(long long)j*neq[1]+imddof[i]];
i                2550 src/steadystate.c 		  for(i=0;i<nmddof;i++){
i                2551 src/steadystate.c 		      FORTRAN(nident,(izdof,&imddof[i],&nzdof,&id));
i                2553 src/steadystate.c 			  if(izdof[id-1]==imddof[i]){
i                2555 src/steadystate.c 				  bi[imddof[i]]+=bji[j]*z[(long long)j*nzdof+id-1];
i                2566 src/steadystate.c 		  for(i=0;i<neq[1];i++){
i                2567 src/steadystate.c 		      breal=br[i];
i                2568 src/steadystate.c 		      br[i]=sqrt(br[i]*br[i]+bi[i]*bi[i]);
i                2570 src/steadystate.c 			  if(bi[i]>0.){bi[i]=pi/2.;}
i                2571 src/steadystate.c 			  else{bi[i]=-pi/2.;}
i                2574 src/steadystate.c 			  bi[i]=atan(bi[i]/breal);
i                2575 src/steadystate.c 			  if(breal<0.){bi[i]+=pi;}
i                2582 src/steadystate.c 		      for(i=0;i<neq[1];i++){
i                2584 src/steadystate.c 			  bi[i]+=pi/2.;}
i                2591 src/steadystate.c 		      for(i=0;i<neq[1];i++){
i                2592 src/steadystate.c 			  btot[j*neq[1]+i]+=br[i]*cos(freq[l]*time+bi[i]);
i                2593 src/steadystate.c 			  bp[j*neq[1]+i]-=freq[l]*br[i]*sin(freq[l]*time+bi[i]);
i                2601 src/steadystate.c 		      i=imddof[jj];
i                2602 src/steadystate.c 		      breal=br[i];
i                2603 src/steadystate.c 		      br[i]=sqrt(br[i]*br[i]+bi[i]*bi[i]);
i                2605 src/steadystate.c 			  if(bi[i]>0.){bi[i]=pi/2.;}
i                2606 src/steadystate.c 			  else{bi[i]=-pi/2.;}
i                2609 src/steadystate.c 			  bi[i]=atan(bi[i]/breal);
i                2610 src/steadystate.c 			  if(breal<0.){bi[i]+=pi;}
i                2618 src/steadystate.c 			  i=imddof[jj];
i                2620 src/steadystate.c 			  bi[i]+=pi/2.;}
i                2628 src/steadystate.c 			  i=imddof[jj];
i                2629 src/steadystate.c 			  btot[j*neq[1]+i]+=br[i]*cos(freq[l]*time+bi[i]);
i                2630 src/steadystate.c 			  bp[j*neq[1]+i]-=freq[l]*br[i]*sin(freq[l]*time+bi[i]);
i                2642 src/steadystate.c 		      for(i=0;i<nmddof;i++){
i                2643 src/steadystate.c 			  br[imddof[i]]=0.;
i                2690 src/steadystate.c 		      i=imdnode[jj]-1;
i                2692 src/steadystate.c 			  vr[mt*i+j]=0.;
i                2837 src/steadystate.c       for(i=0;i<*nmpc;i++){ipompc[i]=ipompcold[i];}
i                2838 src/steadystate.c       for(i=0;i<3**mpcend;i++){nodempc[i]=nodempcold[i];}
i                2839 src/steadystate.c       for(i=0;i<*mpcend;i++){coefmpc[i]=coefmpcold[i];}
i                2840 src/steadystate.c       for(i=0;i<20**nmpc;i++){labmpc[i]=labmpcold[i];}
i                2841 src/steadystate.c       for(i=0;i<*nmpc;i++){ikmpc[i]=ikmpcold[i];}
i                2842 src/steadystate.c       for(i=0;i<*nmpc;i++){ilmpc[i]=ilmpcold[i];}
i                2853 src/steadystate.c       for(i=0;i<*nload;i++){
i                2854 src/steadystate.c 	  if(nelemload[2*i]<=*ne*nsectors){
i                2855 src/steadystate.c 	      nelemload[2*i]-=*ne*nelemload[2*i+1];
i                2857 src/steadystate.c 	      nelemload[2*i]-=*ne*(nsectors+nelemload[2*i+1]-1);
i                2873 src/steadystate.c       for(i=0;i<*nforc;i++){
i                2874 src/steadystate.c 	  if(nodeforc[2*i]<=*nk*nsectors){
i                2875 src/steadystate.c 	      nodeforc[2*i]-=*nk*nodeforc[2*i+1];
i                2877 src/steadystate.c 	      nodeforc[2*i]-=*nk*(nsectors+nodeforc[2*i+1]-1);
i                  17 src/stof.c       register ITG    n, i;
i                  21 src/stof.c       for (i=a-1; i<b; i++)
i                  24 src/stof.c         if ((i>=MAX_LINE_LENGTH)||(n>=MAX_LINE_LENGTH)) break;
i                  25 src/stof.c         puffer[n] = string[i];
i                  17 src/stoi.c       register ITG      n, i;
i                  21 src/stoi.c       for (i=a-1; i<b; i++)
i                  24 src/stoi.c         if ((i>=MAX_LINE_LENGTH)||(n>=MAX_LINE_LENGTH)) break;
i                  25 src/stoi.c         puffer[n] = string[i];
i                  18 src/stos.c       register ITG n, i;
i                  21 src/stos.c       for (i=a-1; i<b; i++)
i                  24 src/stos.c         if ((i>=MAX_LINE_LENGTH)||(n>=MAX_LINE_LENGTH)) break;
i                  25 src/stos.c         puffer[n] = string[i];
i                  33 src/stos.c       register ITG n, i;
i                  36 src/stos.c       for (i=a-1; i<b; i++)
i                  39 src/stos.c         if ((i>132)||(n>132)) break;
i                  40 src/stos.c         puffer[i] = string[n];
i                  30 src/strcmp2.c    ITG a,b,i;
i                  32 src/strcmp2.c    i=0;
i                  47 src/strcmp2.c      i++;
i                  48 src/strcmp2.c    }while((a==b)&&(i<length));
i                  26 src/strcpy1.c    ITG b,i,blank=0;
i                  28 src/strcpy1.c    for(i=0;i<length;i++) {
i                  15 src/strsplt.c    ITG i,j;
i                  27 src/strsplt.c    for(i=0; i<MAX_LINE_LENGTH; i++)
i                  29 src/strsplt.c      if(rec_str[i]==(char)EOF) {  break; } 
i                  30 src/strsplt.c      if(rec_str[i]=='\n') {  break; } 
i                  31 src/strsplt.c      if(rec_str[i]==0) {  break; } 
i                  32 src/strsplt.c      if((rec_str[i]==breakchar)&&(!skip_breakchar))
i                  48 src/strsplt.c        if(rec_str[i]=='"') skip_breakchar=!skip_breakchar;
i                  51 src/strsplt.c          dat[nextarg][letter]=rec_str[i];
i                  54 src/strsplt.c        else if(rec_str[i]!=' ')
i                  56 src/strsplt.c          dat[nextarg][letter]=rec_str[i];
i                  45 src/tau.c        ITG i,j,k,l,*irow=NULL;
i                  70 src/tau.c          for(i=*neq-1;i>=0;--i){
i                  71 src/tau.c            for(j=0;j<icol[i];++j){
i                  75 src/tau.c            pointtau[i]=--k;
i                  76 src/tau.c            irowtau[k]=i;
i                  77 src/tau.c            autau[k]=ad[i];
i                  82 src/tau.c          for(i=*neq-1;i>=0;--i){
i                  83 src/tau.c            for(j=0;j<icol[i];++j){
i                  87 src/tau.c            pointtau[i]=--k;
i                  88 src/tau.c            irowtau[k]=i;
i                  89 src/tau.c            autau[k]=ad[i]-*sigma*adb[i];
i                 128 src/tau.c        ITG i;
i                 148 src/tau.c        for(i=0;i<=*neq-1;++i){
i                 149 src/tau.c          b[i]=x[i];
i                  36 src/tiedcontact.c   ITG *itietri=NULL,*koncont=NULL,nconf,i,k,*nx=NULL,im,
i                 157 src/tiedcontact.c   for(i=0;i<*ntie;i++){
i                 158 src/tiedcontact.c     if(itietri[2*i+1]-itietri[2*i]+1>ntrimax)
i                 159 src/tiedcontact.c       ntrimax=itietri[2*i+1]-itietri[2*i]+1;
i                  33 src/writeheading.c     ITG i;
i                  49 src/writeheading.c     for(i=0;i<*nheading_;i++){
i                  50 src/writeheading.c 	strcpy1(text,&heading[66*i],66);