l                2054 src/CalculiX.h         ITG *i,ITG *l,ITG *ntri));
l                2314 src/CalculiX.h void FORTRAN(writesummary,(ITG *istep,ITG *j,ITG *icutb,ITG *l,double *ttime,
l                2317 src/CalculiX.h void FORTRAN(writesummarydiv,(ITG *istep,ITG *j,ITG *icutb,ITG *l,double *ttime,
l                  85 src/arpackcs.c     ielas=0,icmd=0,kk,l,nkt,icntrl,*kont=NULL,*ipkont=NULL,*inumt=NULL,
l                1079 src/arpackcs.c       for(l=0;l<*ne*ngraph;l++)ipkont[l]=-1;
l                1089 src/arpackcs.c       for(l=0;l<3**nk;l++){cot[l]=co[l];}
l                1090 src/arpackcs.c       if(*ntrans>0){for(l=0;l<*nk;l++){inotrt[2*l]=inotr[2*l];}}
l                1091 src/arpackcs.c       for(l=0;l<*nkon;l++){kont[l]=kon[l];}
l                1092 src/arpackcs.c       for(l=0;l<*ne;l++){ipkont[l]=ipkon[l];}
l                1093 src/arpackcs.c       for(l=0;l<8**ne;l++){lakont[l]=lakon[l];}
l                1094 src/arpackcs.c       for(l=0;l<*ne;l++){ielmatt[mi[2]*l]=ielmat[mi[2]*l];}
l                1108 src/arpackcs.c 	      for(l=0;l<*nk;l++){
l                1109 src/arpackcs.c 		  if(inocs[l]==jj){
l                1110 src/arpackcs.c 		      cot[3*l+i*3**nk]=cot[3*l];
l                1111 src/arpackcs.c 		      cot[1+3*l+i*3**nk]=cot[1+3*l]+theta;
l                1112 src/arpackcs.c 		      cot[2+3*l+i*3**nk]=cot[2+3*l];
l                1113 src/arpackcs.c 		      if(*ntrans>0){inotrt[2*l+i*2**nk]=inotrt[2*l];}
l                1116 src/arpackcs.c 	      for(l=0;l<*nkon;l++){kont[l+i**nkon]=kon[l]+i**nk;}
l                1117 src/arpackcs.c 	      for(l=0;l<*ne;l++){
l                1118 src/arpackcs.c 		  if(ielcs[l]==jj){
l                1119 src/arpackcs.c 		      if(ipkon[l]>=0){
l                1120 src/arpackcs.c 			  ipkont[l+i**ne]=ipkon[l]+i**nkon;
l                1121 src/arpackcs.c 			  ielmatt[mi[2]*(l+i**ne)]=ielmat[mi[2]*l];
l                1123 src/arpackcs.c 			      l2=8*l+l1;
l                1402 src/arpackcs.c 			  l=mt*l1+l2;
l                1403 src/arpackcs.c 			  vl[l2]=ctl*v[l]-stl*v[l+mt**nk];
l                1418 src/arpackcs.c 		      l=mt*l1+l2;
l                1419 src/arpackcs.c 		      vl[l2]=ctl*v[l]-stl*v[l+mt**nk];
l                1461 src/arpackcs.c 			      l=6*l1+l2;
l                1462 src/arpackcs.c 			      stnl[l2]=ctl*stn[l]-stl*stn[l+6**nk];
l                1524 src/arpackcs.c 				  l=6*l1+l2;
l                1525 src/arpackcs.c 				  stnl[l2]=ctl*stn[l]-stl*stn[l+6**nk];
l                1614 src/arpackcs.c 	      l=6*l1+l2;
l                1615 src/arpackcs.c 	      eenl[l2]=ctl*een[l]-stl*een[l+6**nk];
l                1677 src/arpackcs.c 		l=6*l1+l2;
l                1678 src/arpackcs.c 		eenl[l2]=ctl*een[l]-stl*een[l+6**nk];
l                1736 src/arpackcs.c     for(l=0;l<*nk;l++){inumt[l]=inum[l];}
l                1743 src/arpackcs.c       for(l=0;l<mt**nk;l++){vt[l]=v[l];}
l                1744 src/arpackcs.c       for(l=0;l<mt**nk;l++){vt[l+mt**nk*ngraph]=v[l+mt**nk];}}
l                1746 src/arpackcs.c       for(l=0;l<*nk;l++){t1t[l]=t1[l];};
l                1748 src/arpackcs.c 	for(l=0;l<6**nk;l++){stnt[l]=stn[l];}
l                1749 src/arpackcs.c 	for(l=0;l<6**nk;l++){stnt[l+6**nk*ngraph]=stn[l+6**nk];}}
l                1751 src/arpackcs.c 	for(l=0;l<6**nk;l++){eent[l]=een[l];};
l                1752 src/arpackcs.c 	for(l=0;l<6**nk;l++){eent[l+6**nk*ngraph]=een[l+6**nk];}}
l                1754 src/arpackcs.c       for(l=0;l<mt**nk;l++){fnt[l]=fn[l];}
l                1755 src/arpackcs.c       for(l=0;l<mt**nk;l++){fnt[l+mt**nk*ngraph]=fn[l+mt**nk];}}
l                1757 src/arpackcs.c       for(l=0;l<*nk;l++){enernt[l]=enern[l];};
l                1760 src/arpackcs.c       for(l=0;l<6*mi[0]**ne;l++){stxt[l]=stx[l];}
l                1761 src/arpackcs.c       for(l=0;l<6*mi[0]**ne;l++){stxt[l+6*mi[0]**ne*ngraph]=stx[l+6*mi[0]**ne];}}
l                1763 src/arpackcs.c 	for(l=0;l<6**nk;l++){emnt[l]=emn[l];};
l                1764 src/arpackcs.c 	for(l=0;l<6**nk;l++){emnt[l+6**nk*ngraph]=emn[l+6**nk];}}
l                1767 src/arpackcs.c 	for(l=0;l<6**nk;l++){cdnt[l]=cdn[l];}
l                1768 src/arpackcs.c 	for(l=0;l<6**nk;l++){cdnt[l+6**nk*ngraph]=cdn[l+6**nk];}}
l                1776 src/arpackcs.c         for(l=0;l<*nk;l++){inumt[l+i**nk]=inum[l];}
l                1793 src/arpackcs.c 			  l=mt*l1+l2;
l                1794 src/arpackcs.c 			  vt[l+mt**nk*i]=v[l];
l                1800 src/arpackcs.c                 l=mt*l1+l2;
l                1801 src/arpackcs.c                 vt[l+mt**nk*i]=ctl*v[l]-stl*v[l+mt**nk];
l                1821 src/arpackcs.c 			  l=mt*l1+l2;
l                1822 src/arpackcs.c 			  vt[l+mt**nk*(i+ngraph)]=v[l+mt**nk];
l                1828 src/arpackcs.c                 l=mt*l1+l2;
l                1829 src/arpackcs.c                 vt[l+mt**nk*(i+ngraph)]=stl*v[l]+ctl*v[l+mt**nk];
l                1836 src/arpackcs.c           for(l=0;l<*nk;l++){
l                1837 src/arpackcs.c 	      if(inocs[l]==jj) t1t[l+*nk*i]=t1[l];
l                1852 src/arpackcs.c 			  l=6*l1+l2;
l                1853 src/arpackcs.c 			  stnt[l+6**nk*i]=stn[l];
l                1859 src/arpackcs.c                 l=6*l1+l2;
l                1860 src/arpackcs.c                 stnt[l+6**nk*i]=ctl*stn[l]-stl*stn[l+6**nk];
l                1880 src/arpackcs.c 			  l=6*l1+l2;
l                1881 src/arpackcs.c 			  stnt[l+6**nk*(i+ngraph)]=stn[l+6**nk];
l                1887 src/arpackcs.c                 l=6*l1+l2;
l                1888 src/arpackcs.c                 stnt[l+6**nk*(i+ngraph)]=stl*stn[l]+ctl*stn[l+6**nk];
l                1905 src/arpackcs.c 			  l=6*l1+l2;
l                1906 src/arpackcs.c 			  eent[l+6**nk*i]=een[l];
l                1912 src/arpackcs.c                 l=6*l1+l2;
l                1913 src/arpackcs.c                 eent[l+6**nk*i]=ctl*een[l]-stl*een[l+6**nk];
l                1933 src/arpackcs.c 			  l=6*l1+l2;
l                1934 src/arpackcs.c 			  eent[l+6**nk*(i+ngraph)]=een[l+6**nk];
l                1940 src/arpackcs.c                 l=6*l1+l2;
l                1941 src/arpackcs.c                 eent[l+6**nk*(i+ngraph)]=stl*een[l]+ctl*een[l+6**nk];
l                1960 src/arpackcs.c 			  l=6*l1+l2;
l                1961 src/arpackcs.c 			  emnt[l+6**nk*i]=emn[l];
l                1967 src/arpackcs.c                 l=6*l1+l2;
l                1968 src/arpackcs.c                 emnt[l+6**nk*i]=ctl*emn[l]-stl*emn[l+6**nk];
l                1988 src/arpackcs.c 			  l=6*l1+l2;
l                1989 src/arpackcs.c 			  emnt[l+6**nk*(i+ngraph)]=emn[l+6**nk];
l                1995 src/arpackcs.c                 l=6*l1+l2;
l                1996 src/arpackcs.c                 emnt[l+6**nk*(i+ngraph)]=stl*emn[l]+ctl*emn[l+6**nk];
l                2016 src/arpackcs.c 			  l=6*l1+l2;
l                2017 src/arpackcs.c 			  cdnt[l+6**nk*i]=cdn[l];
l                2023 src/arpackcs.c                 l=6*l1+l2;
l                2024 src/arpackcs.c                 cdnt[l+6**nk*i]=ctl*cdn[l]-stl*cdn[l+6**nk];
l                2044 src/arpackcs.c 			  l=6*l1+l2;
l                2045 src/arpackcs.c 			  cdnt[l+6**nk*(i+ngraph)]=cdn[l+6**nk];
l                2051 src/arpackcs.c                 l=6*l1+l2;
l                2052 src/arpackcs.c                 cdnt[l+6**nk*(i+ngraph)]=stl*cdn[l]+ctl*cdn[l+6**nk];
l                2069 src/arpackcs.c 			  l=mt*l1+l2;
l                2070 src/arpackcs.c 			  fnt[l+mt**nk*i]=fn[l];
l                2076 src/arpackcs.c                 l=mt*l1+l2;
l                2077 src/arpackcs.c                 fnt[l+mt**nk*i]=ctl*fn[l]-stl*fn[l+mt**nk];
l                2097 src/arpackcs.c 		    l=mt*l1+l2;
l                2098 src/arpackcs.c 		    fnt[l+mt**nk*(i+ngraph)]=fn[l+mt**nk];
l                2104 src/arpackcs.c                 l=mt*l1+l2;
l                2105 src/arpackcs.c                 fnt[l+mt**nk*(i+ngraph)]=stl*fn[l]+ctl*fn[l+mt**nk];
l                2112 src/arpackcs.c           for(l=0;l<*nk;l++){
l                2113 src/arpackcs.c             if(inocs[l]==jj) enernt[l+*nk*i]=0.;
l                2133 src/arpackcs.c 	  l=mt*l1+l2;
l                2134 src/arpackcs.c 	  vreal=vt[l];
l                2135 src/arpackcs.c 	  vimag=vt[l+mt**nk*ngraph];
l                2136 src/arpackcs.c 	  vr[l]=sqrt(vreal*vreal+vimag*vimag);
l                2138 src/arpackcs.c 	    if(vimag>0){vi[l]=90.;}
l                2139 src/arpackcs.c 	    else{vi[l]=-90.;}
l                2142 src/arpackcs.c 	    vi[l]=atan(vimag/vreal)*constant;
l                2143 src/arpackcs.c 	    if(vreal<0) vi[l]+=180.;
l                2154 src/arpackcs.c 	  l=6*l1+l2;
l                2155 src/arpackcs.c 	  stnreal=stnt[l];
l                2156 src/arpackcs.c 	  stnimag=stnt[l+6**nk*ngraph];
l                2157 src/arpackcs.c 	  stnr[l]=sqrt(stnreal*stnreal+stnimag*stnimag);
l                2159 src/arpackcs.c 	    if(stnimag>0){stni[l]=90.;}
l                2160 src/arpackcs.c 	    else{stni[l]=-90.;}
l                2163 src/arpackcs.c 	    stni[l]=atan(stnimag/stnreal)*constant;
l                2164 src/arpackcs.c 	    if(stnreal<0) stni[l]+=180.;
l                2175 src/arpackcs.c 	  l=mt*l1+l2;
l                2176 src/arpackcs.c 	  fnreal=fnt[l];
l                2177 src/arpackcs.c 	  fnimag=fnt[l+mt**nk*ngraph];
l                2178 src/arpackcs.c 	  fnr[l]=sqrt(fnreal*fnreal+fnimag*fnimag);
l                2180 src/arpackcs.c 	    if(fnimag>0){fni[l]=90.;}
l                2181 src/arpackcs.c 	    else{fni[l]=-90.;}
l                2184 src/arpackcs.c 	    fni[l]=atan(fnimag/fnreal)*constant;
l                2185 src/arpackcs.c 	    if(fnreal<0) fni[l]+=180.;
l                2196 src/arpackcs.c 	  l=6*l1+l2;
l                2197 src/arpackcs.c 	  cdnreal=cdnt[l];
l                2198 src/arpackcs.c 	  cdnimag=cdnt[l+6**nk*ngraph];
l                2199 src/arpackcs.c 	  cdnr[l]=sqrt(cdnreal*cdnreal+cdnimag*cdnimag);
l                2201 src/arpackcs.c 	    if(cdnimag>0){cdni[l]=90.;}
l                2202 src/arpackcs.c 	    else{cdni[l]=-90.;}
l                2205 src/arpackcs.c 	    cdni[l]=atan(cdnimag/cdnreal)*constant;
l                2206 src/arpackcs.c 	    if(cdnreal<0) cdni[l]+=180.;
l                  78 src/complexfreq.c     iinc=0,l,iout=1,ielas,icmd,ifreebody,mode,m,nherm,
l                 620 src/complexfreq.c       for(l=0;l<nev;l++){
l                 621 src/complexfreq.c 	  if(d[l]<0.){
l                 622 src/complexfreq.c 	      eig[l]=0.;
l                 624 src/complexfreq.c 	      eig[l]=sqrt(d[l]);
l                 630 src/complexfreq.c       for(l=0;l<nev;l++){
l                 631 src/complexfreq.c 	  eig[2*l]=sqrt(sqrt(d[2*l]*d[2*l]+d[2*l+1]*d[2*l+1])+d[2*l])/sqrt(2.);
l                 632 src/complexfreq.c 	  eig[2*l+1]=sqrt(sqrt(d[2*l]*d[2*l]+d[2*l+1]*d[2*l+1])-d[2*l])/sqrt(2.);
l                 633 src/complexfreq.c 	  if(d[2*l+1]<0.) eig[2*l+1]=-eig[2*l+1];
l                 636 src/complexfreq.c   for(l=0;l<nev;l++){
l                 637 src/complexfreq.c       mode=l+1;
l                 639 src/complexfreq.c 	  eigxr[k]=xx[2*l*nev+2*k];
l                 640 src/complexfreq.c 	  eigxi[k]=xx[2*l*nev+2*k+1];
l                 654 src/complexfreq.c       for(l=0;l<nev;l++){
l                 660 src/complexfreq.c 		  zz[(long long)2*l*neqact+i]+=
l                 661 src/complexfreq.c 		      (xx[2*l*nev+2*k]*z[(long long)2*k*neqact+i]
l                 662 src/complexfreq.c 		       -xx[2*l*nev+2*k+1]*z[(long long)(2*k+1)*neqact+i]);
l                 666 src/complexfreq.c 		  zz[(long long)(2*l+1)*neqact+i]+=
l                 667 src/complexfreq.c 		      (xx[2*l*nev+2*k]*z[(long long)(2*k+1)*neqact+i]
l                 668 src/complexfreq.c 		       +xx[2*l*nev+2*k+1]*z[(long long)2*k*neqact+i]);
l                 690 src/complexfreq.c       for(l=0;l<nev;l++){
l                 693 src/complexfreq.c 	  FORTRAN(op,(&neq[1],&zz[l*neq[1]],z,adb,aub,jq,irow));
l                 695 src/complexfreq.c 	      sum+=zz[l*neq[1]+k]*z[k];
l                 700 src/complexfreq.c 	      zz[l*neq[1]+k]/=sum;
l                 717 src/complexfreq.c 	  for(l=0;l<nev;l++){
l                 719 src/complexfreq.c 	      FORTRAN(op,(&neq[1],&zz[l*neq[1]],z,adb,aub,jq,irow));
l                 720 src/complexfreq.c 	      for(m=l;m<nev;m++){
l                 722 src/complexfreq.c 		      xmr[l*nev+m]+=zz[m*neq[1]+k]*z[k];
l                 726 src/complexfreq.c 	      memcpy(&zi[0],&zz[(2*l+1)*neqact],sizeof(double)*neqact);
l                 727 src/complexfreq.c 	      for(k=0;k<neqact;k++){zi[neqact+k]=-zz[2*l*neqact+k];}
l                 730 src/complexfreq.c 	      for(m=l;m<nev;m++){
l                 732 src/complexfreq.c 		      xmi[l*nev+m]+=zz[m*neq[1]+k]*z[k];
l                 739 src/complexfreq.c 	  for(l=0;l<nev;l++){
l                 740 src/complexfreq.c 	      for(m=0;m<l;m++){
l                 741 src/complexfreq.c 		  xmr[l*nev+m]=xmr[m*nev+l];
l                 742 src/complexfreq.c 		  xmi[l*nev+m]=-xmi[m*nev+l];
l                 746 src/complexfreq.c 	  for(l=0;l<nev;l++){
l                 748 src/complexfreq.c 		  printf(" %f",xmr[m*nev+l]);
l                 753 src/complexfreq.c 	  for(l=0;l<nev;l++){
l                 755 src/complexfreq.c 		  printf(" %f",xmi[m*nev+l]);
l                 769 src/complexfreq.c       for(l=0;l<nev;l++){
l                 773 src/complexfreq.c 		      zz[(2*l+j)*neq[1]+i]+=xx[2*l*nev+2*k+j]*
l                 797 src/complexfreq.c       for(l=0;l<nev;l++){
l                 803 src/complexfreq.c 	  FORTRAN(op,(&neq[1],&zz[2*l*neq[1]],z,adb,aub,jq,irow));
l                 805 src/complexfreq.c 	      sum+=zz[2*l*neq[1]+k]*z[k];
l                 811 src/complexfreq.c 	  FORTRAN(op,(&neq[1],&zz[(2*l+1)*neq[1]],z,adb,aub,jq,irow));
l                 813 src/complexfreq.c 	      sum+=zz[(2*l+1)*neq[1]+k]*z[k];
l                 818 src/complexfreq.c 	      zz[2*l*neq[1]+k]/=sum;
l                 834 src/complexfreq.c 	  for(l=0;l<nev;l++){
l                 840 src/complexfreq.c 	      FORTRAN(op,(&neq[1],&zz[2*l*neq[1]],z,adb,aub,jq,irow));
l                 844 src/complexfreq.c 	      for(m=l;m<nev;m++){
l                 846 src/complexfreq.c 		      xmr[l*nev+m]+=zz[2*m*neq[1]+k]*z[k];
l                 849 src/complexfreq.c 		      xmi[l*nev+m]-=zz[(2*m+1)*neq[1]+k]*z[k];
l                 856 src/complexfreq.c 	      FORTRAN(op,(&neq[1],&zz[(2*l+1)*neq[1]],z,adb,aub,jq,irow));
l                 860 src/complexfreq.c 	      for(m=l;m<nev;m++){
l                 862 src/complexfreq.c 		      xmr[l*nev+m]+=zz[(2*m+1)*neq[1]+k]*z[k];
l                 865 src/complexfreq.c 		      xmi[l*nev+m]+=zz[2*m*neq[1]+k]*z[k];
l                 872 src/complexfreq.c 	  for(l=0;l<nev;l++){
l                 873 src/complexfreq.c 	      for(m=0;m<l;m++){
l                 874 src/complexfreq.c 		  xmr[l*nev+m]=xmr[m*nev+l];
l                 875 src/complexfreq.c 		  xmi[l*nev+m]=-xmi[m*nev+l];
l                 879 src/complexfreq.c 	  for(l=0;l<nev;l++){
l                 881 src/complexfreq.c 		  printf(" %f",xmr[m*nev+l]);
l                 886 src/complexfreq.c 	  for(l=0;l<nev;l++){
l                 888 src/complexfreq.c 		  printf(" %f",xmi[m*nev+l]);
l                1137 src/complexfreq.c   for(l=0;l<*ne*ngraph;l++)ipkont[l]=-1;
l                1147 src/complexfreq.c   for(l=0;l<3**nk;l++){cot[l]=co[l];}
l                1148 src/complexfreq.c   if(*ntrans>0){for(l=0;l<*nk;l++){inotrt[2*l]=inotr[2*l];}}
l                1149 src/complexfreq.c   for(l=0;l<*nkon;l++){kont[l]=kon[l];}
l                1150 src/complexfreq.c   for(l=0;l<*ne;l++){ipkont[l]=ipkon[l];}
l                1151 src/complexfreq.c   for(l=0;l<8**ne;l++){lakont[l]=lakon[l];}
l                1152 src/complexfreq.c   for(l=0;l<*ne;l++){ielmatt[mi[2]*l]=ielmat[mi[2]*l];}
l                1168 src/complexfreq.c 	for(l=0;l<*nk;l++){
l                1169 src/complexfreq.c 	  if(inocs[l]==jj){
l                1170 src/complexfreq.c 	    cot[3*l+i*3**nk]=cot[3*l];
l                1171 src/complexfreq.c 	    cot[1+3*l+i*3**nk]=cot[1+3*l]+theta;
l                1172 src/complexfreq.c 	    cot[2+3*l+i*3**nk]=cot[2+3*l];
l                1173 src/complexfreq.c 	    if(*ntrans>0){inotrt[2*l+i*2**nk]=inotrt[2*l];}
l                1176 src/complexfreq.c 	for(l=0;l<*nkon;l++){kont[l+i**nkon]=kon[l]+i**nk;}
l                1177 src/complexfreq.c 	for(l=0;l<*ne;l++){
l                1178 src/complexfreq.c 	  if(ielcs[l]==jj){
l                1179 src/complexfreq.c 	    if(ipkon[l]>=0){
l                1180 src/complexfreq.c 	      ipkont[l+i**ne]=ipkon[l]+i**nkon;
l                1181 src/complexfreq.c 	      ielmatt[mi[2]*(l+i**ne)]=ielmat[mi[2]*l];
l                1183 src/complexfreq.c 		l2=8*l+l1;
l                1388 src/complexfreq.c     for(l=0;l<*nk;l++){inumt[l]=inum[l];}
l                1399 src/complexfreq.c       for(l=0;l<mt**nk;l++){vt[l]=v[l];}
l                1400 src/complexfreq.c       for(l=0;l<mt**nk;l++){vt[l+mt**nk*ngraph]=v[l+mt**nk];}}
l                1402 src/complexfreq.c       for(l=0;l<*nk;l++){t1t[l]=t1[l];};
l                1404 src/complexfreq.c 	for(l=0;l<6**nk;l++){stnt[l]=stn[l];}
l                1405 src/complexfreq.c 	for(l=0;l<6**nk;l++){stnt[l+6**nk*ngraph]=stn[l+6**nk];}}
l                1407 src/complexfreq.c 	for(l=0;l<6**nk;l++){eent[l]=een[l];};
l                1408 src/complexfreq.c 	for(l=0;l<6**nk;l++){eent[l+6**nk*ngraph]=een[l+6**nk];}}
l                1410 src/complexfreq.c       for(l=0;l<mt**nk;l++){fnt[l]=fn[l];}
l                1411 src/complexfreq.c       for(l=0;l<mt**nk;l++){fnt[l+mt**nk*ngraph]=fn[l+mt**nk];}}
l                1413 src/complexfreq.c       for(l=0;l<*nk;l++){enernt[l]=enern[l];};
l                1415 src/complexfreq.c       for(l=0;l<6*mi[0]**ne;l++){stxt[l]=stx[l];}
l                1416 src/complexfreq.c       for(l=0;l<6*mi[0]**ne;l++){stxt[l+6*mi[0]**ne*ngraph]=stx[l+6*mi[0]**ne];}}
l                1418 src/complexfreq.c 	for(l=0;l<6**nk;l++){emnt[l]=emn[l];};
l                1419 src/complexfreq.c 	for(l=0;l<6**nk;l++){emnt[l+6**nk*ngraph]=emn[l+6**nk];}}
l                1432 src/complexfreq.c 	  for(l=0;l<*nk;l++){inumt[l+i**nk]=inum[l];}
l                1449 src/complexfreq.c 		      l=mt*l1+l2;
l                1450 src/complexfreq.c 		      vt[l+mt**nk*i]=v[l];
l                1456 src/complexfreq.c 		  l=mt*l1+l2;
l                1457 src/complexfreq.c 		  vt[l+mt**nk*i]=ctl*v[l]-stl*v[l+mt**nk];
l                1478 src/complexfreq.c 		      l=mt*l1+l2;
l                1479 src/complexfreq.c 		      vt[l+mt**nk*(i+ngraph)]=v[l+mt**nk];
l                1485 src/complexfreq.c 		  l=mt*l1+l2;
l                1486 src/complexfreq.c 		  vt[l+mt**nk*(i+ngraph)]=stl*v[l]+ctl*v[l+mt**nk];
l                1493 src/complexfreq.c 	    for(l=0;l<*nk;l++){
l                1494 src/complexfreq.c 	      if(inocs[l]==jj) t1t[l+*nk*i]=t1[l];
l                1509 src/complexfreq.c 		      l=6*l1+l2;
l                1510 src/complexfreq.c 		      stnt[l+6**nk*i]=stn[l];
l                1516 src/complexfreq.c 		  l=6*l1+l2;
l                1517 src/complexfreq.c 		  stnt[l+6**nk*i]=ctl*stn[l]-stl*stn[l+6**nk];
l                1537 src/complexfreq.c 		      l=6*l1+l2;
l                1538 src/complexfreq.c 		      stnt[l+6**nk*(i+ngraph)]=stn[l+6**nk];
l                1544 src/complexfreq.c 		  l=6*l1+l2;
l                1545 src/complexfreq.c 		  stnt[l+6**nk*(i+ngraph)]=stl*stn[l]+ctl*stn[l+6**nk];
l                1562 src/complexfreq.c 		      l=6*l1+l2;
l                1563 src/complexfreq.c 		      eent[l+6**nk*i]=een[l];
l                1569 src/complexfreq.c 		  l=6*l1+l2;
l                1570 src/complexfreq.c 		  eent[l+6**nk*i]=ctl*een[l]-stl*een[l+6**nk];
l                1590 src/complexfreq.c 				  l=6*l1+l2;
l                1591 src/complexfreq.c 				  eent[l+6**nk*(i+ngraph)]=een[l+6**nk];
l                1597 src/complexfreq.c 			  l=6*l1+l2;
l                1598 src/complexfreq.c 			  eent[l+6**nk*(i+ngraph)]=stl*een[l]+ctl*een[l+6**nk];
l                1615 src/complexfreq.c 		      l=6*l1+l2;
l                1616 src/complexfreq.c 		      emnt[l+6**nk*i]=emn[l];
l                1622 src/complexfreq.c 		  l=6*l1+l2;
l                1623 src/complexfreq.c 		  emnt[l+6**nk*i]=ctl*emn[l]-stl*emn[l+6**nk];
l                1643 src/complexfreq.c 				  l=6*l1+l2;
l                1644 src/complexfreq.c 				  emnt[l+6**nk*(i+ngraph)]=emn[l+6**nk];
l                1650 src/complexfreq.c 			  l=6*l1+l2;
l                1651 src/complexfreq.c 			  emnt[l+6**nk*(i+ngraph)]=stl*emn[l]+ctl*emn[l+6**nk];
l                1668 src/complexfreq.c 		      l=mt*l1+l2;
l                1669 src/complexfreq.c 		      fnt[l+mt**nk*i]=fn[l];
l                1675 src/complexfreq.c 		  l=mt*l1+l2;
l                1676 src/complexfreq.c 		  fnt[l+mt**nk*i]=ctl*fn[l]-stl*fn[l+mt**nk];
l                1696 src/complexfreq.c 				  l=mt*l1+l2;
l                1697 src/complexfreq.c 				  fnt[l+mt**nk*(i+ngraph)]=fn[l+mt**nk];
l                1703 src/complexfreq.c 			  l=mt*l1+l2;
l                1704 src/complexfreq.c 			  fnt[l+mt**nk*(i+ngraph)]=stl*fn[l]+ctl*fn[l+mt**nk];
l                1711 src/complexfreq.c 	    for(l=0;l<*nk;l++){
l                1712 src/complexfreq.c 	      if(inocs[l]==jj) enernt[l+*nk*i]=0.;
l                1734 src/complexfreq.c 	  l=mt*l1+l2;
l                1735 src/complexfreq.c 	  vreal=vt[l];
l                1736 src/complexfreq.c 	  vimag=vt[l+mt**nk*ngraph];
l                1737 src/complexfreq.c 	  vr[l]=sqrt(vreal*vreal+vimag*vimag);
l                1739 src/complexfreq.c 	    if(vimag>0){vi[l]=90.;}
l                1740 src/complexfreq.c 	    else{vi[l]=-90.;}
l                1743 src/complexfreq.c 	    vi[l]=atan(vimag/vreal)*constant;
l                1744 src/complexfreq.c 	    if(vreal<0) vi[l]+=180.;
l                1755 src/complexfreq.c 	  l=6*l1+l2;
l                1756 src/complexfreq.c 	  stnreal=stnt[l];
l                1757 src/complexfreq.c 	  stnimag=stnt[l+6**nk*ngraph];
l                1758 src/complexfreq.c 	  stnr[l]=sqrt(stnreal*stnreal+stnimag*stnimag);
l                1760 src/complexfreq.c 	    if(stnimag>0){stni[l]=90.;}
l                1761 src/complexfreq.c 	    else{stni[l]=-90.;}
l                1764 src/complexfreq.c 	    stni[l]=atan(stnimag/stnreal)*constant;
l                1765 src/complexfreq.c 	    if(stnreal<0) stni[l]+=180.;
l                  35 src/dfdbj.c      ITG j,j1,jdof,kdof,k,k1,l,id,index,ist,id1,ist1,index1,id2,ist2,index2,
l                  63 src/dfdbj.c    		  l=k-1;
l                  65 src/dfdbj.c    		      ikactcont[k]=ikactcont[l];
l                  66 src/dfdbj.c    		      ilactcont[k--]=ilactcont[l--];
l                  84 src/dfdbj.c    			  for(l=0; l<*nev; l++){
l                  85 src/dfdbj.c    			      dbcont[i4+l]-=sl*z[(long long)l**neq+kdof-1];
l                  92 src/dfdbj.c    			      for(l=0; l<*nev; l++){
l                  93 src/dfdbj.c    				dbcont[i4+l]-=sl*z[l**nzdof+id-1];
l                 116 src/dfdbj.c    					  for(l=0; l<*nev; l++){
l                 117 src/dfdbj.c    					    dbcont[i4+l]+=d1*z[(long long)l**neq+kdof-1];
l                 125 src/dfdbj.c    					    for(l=0; l<*nev; l++){
l                 126 src/dfdbj.c    					      dbcont[i4+l]+=d1*z[l**nzdof+id-1];
l                 172 src/dfdbj.c    				  l=k-1;
l                 174 src/dfdbj.c    				      ikactcont[k]=ikactcont[l];
l                 175 src/dfdbj.c    				      ilactcont[k--]=ilactcont[l--];
l                 194 src/dfdbj.c    					    for(l=0; l<*nev; l++){
l                 195 src/dfdbj.c    					      dbcont[i4+l]+=d1*z[(long long)l**neq+kdof-1];
l                 202 src/dfdbj.c    						for(l=0; l<*nev; l++){
l                 203 src/dfdbj.c    						  dbcont[i4+l]+=d1*z[l**nzdof+id-1];
l                 226 src/dfdbj.c    							    for(l=0; l<*nev; l++){
l                 227 src/dfdbj.c    							      dbcont[i4+l]-=d1*z[(long long)l**neq+kdof-1];
l                 234 src/dfdbj.c    								for(l=0; l<*nev; l++){
l                 235 src/dfdbj.c    								  dbcont[i4+l]-=d1*z[l**nzdof+id-1];
l                  78 src/dyna.c         iinc=0,jprint=0,l,iout,ielas=0,icmd,iprescribedboundary,init,ifreebody,
l                1530 src/dyna.c     	  for(l=0;l<nev;l++){
l                1531 src/dyna.c     	      zetaj=zeta[l];
l                1532 src/dyna.c     	      dj=d[l];
l                1536 src/dyna.c     	      if(fabs(d[l])<=1.e-10){
l                1537 src/dyna.c     		  aai=aa[l];
l                1538 src/dyna.c     		  bbi=bb[l];
l                1546 src/dyna.c     		  bj[l]=sum+cd[l]+dtime*cv[l];
l                1547 src/dyna.c     		  bjp[l]=sump+cv[l];
l                1560 src/dyna.c     		  FORTRAN(fsub,(&time,&dtime,&aa[l],&bb[l],&ddj,
l                1563 src/dyna.c     		  FORTRAN(fsub,(&time,&tstart,&aa[l],&bb[l],&ddj,
l                1570 src/dyna.c     		  bj[l]=sum/ddj+fexp*(fcos+zetaj/sqrt(1.-zetaj*zetaj)*fsin)*cd[l]+
l                1571 src/dyna.c     		      fexp*fsin*cv[l]/ddj;
l                1572 src/dyna.c     		  bjp[l]=sump/ddj+fexp*((-h1+ddj*h14)*fcos+(-ddj-h1*h14)*fsin)*cd[l]
l                1573 src/dyna.c     		    +fexp*(-h1*fsin+ddj*fcos)*cv[l]/ddj;
l                1588 src/dyna.c     		  FORTRAN(fsuper,(&time,&dtime,&aa[l],&bb[l],
l                1591 src/dyna.c     		  FORTRAN(fsuper,(&time,&tstart,&aa[l],&bb[l],
l                1598 src/dyna.c     		  bj[l]=sum/(2.*ddj)+(fexm+fexp)*cd[l]/2.+zetaj*(fexm-fexp)/(2.*
l                1599 src/dyna.c     		      sqrt(zetaj*zetaj-1.))*cd[l]+(fexm-fexp)*cv[l]/(2.*ddj);
l                1600 src/dyna.c     		  bjp[l]=sump/(2.*ddj)+(h1*fexm-h2*fexp)*cd[l]/2.
l                1601 src/dyna.c     		    +(h14*cd[l]+cv[l]/ddj)*(h1*fexm+h2*fexp)/2.;
l                1612 src/dyna.c     		  FORTRAN(fcrit,(&time,&dtime,&aa[l],&bb[l],&zetaj,&dj,
l                1615 src/dyna.c     		  FORTRAN(fcrit,(&time,&tstart,&aa[l],&bb[l],&zetaj,&dj,
l                1619 src/dyna.c     		  bj[l]=sum+fexp*((1.+h1*dtime)*cd[l]+dtime*cv[l]);
l                1620 src/dyna.c     		  bjp[l]=sump+fexp*(-h1*h1*dtime*cd[l]+
l                1621 src/dyna.c                                         (1.-h1*dtime)*cv[l]);
l                 100 src/dynacont.c   ITG i,j,k,l,init,*itg=NULL,ntg=0,maxlenmpc,icascade=0,loop,
l                 218 src/dynacont.c 	  l=k-1;
l                 220 src/dynacont.c 	      ikactcont[k--]=ikactcont[l--];
l                 418 src/dynacont.c     for(l=0;l<*nev;l++){
l                 419 src/dynacont.c       zetaj=zeta[l];
l                 420 src/dynacont.c       dj=d[l];
l                 424 src/dynacont.c       if(fabs(d[l])<=1.e-10){
l                 425 src/dynacont.c 	  aai=aa[l];
l                 426 src/dynacont.c 	  bbi=bb[l];
l                 434 src/dynacont.c 	  bjbas[l]=sum+cd[l]+*dtime*cv[l];
l                 435 src/dynacont.c 	  bjbasp[l]=sump+cv[l];
l                 448 src/dynacont.c 	FORTRAN(fsub,(time,dtime,&aa[l],&bb[l],&ddj,
l                 451 src/dynacont.c 	FORTRAN(fsub,(time,&tstart,&aa[l],&bb[l],&ddj,
l                 458 src/dynacont.c 	bjbas[l]=sum/ddj+fexp*(fcos+zetaj/sqrt(1.-zetaj*zetaj)*fsin)*cd[l]+
l                 459 src/dynacont.c 	  fexp*fsin*cv[l]/ddj;
l                 460 src/dynacont.c 	bjbasp[l]=sump/ddj+fexp*((-h1+ddj*h14)*fcos+(-ddj-h1*h14)*fsin)*cd[l]
l                 461 src/dynacont.c 	  +fexp*(-h1*fsin+ddj*fcos)*cv[l]/ddj;
l                 476 src/dynacont.c 	FORTRAN(fsuper,(time,dtime,&aa[l],&bb[l],
l                 479 src/dynacont.c 	FORTRAN(fsuper,(time,&tstart,&aa[l],&bb[l],
l                 486 src/dynacont.c 	bjbas[l]=sum/(2.*ddj)+(fexm+fexp)*cd[l]/2.+zetaj*(fexm-fexp)/(2.*sqrt(zetaj*zetaj-1.))*cd[l]+(fexm-fexp)*cv[l]/(2.*ddj);
l                 487 src/dynacont.c 	bjbasp[l]=sump/(2.*ddj)+(h1*fexm-h2*fexp)*cd[l]/2.+(h14*cd[l]+cv[l]/ddj)*(h1*fexm+h2*fexp)/2.;
l                 498 src/dynacont.c 	FORTRAN(fcrit,(time,dtime,&aa[l],&bb[l],&zetaj,&dj,
l                 501 src/dynacont.c 	FORTRAN(fcrit,(time,&tstart,&aa[l],&bb[l],&zetaj,&dj,
l                 505 src/dynacont.c 	bjbas[l]=sum+fexp*((1.+h1**dtime)*cd[l]+*dtime*cv[l]);
l                 506 src/dynacont.c 	bjbasp[l]=sump+fexp*(-h1*h1**dtime*cd[l]+(1.-h1**dtime)*cv[l]);
l                 517 src/dynacont.c     for(l=0;l<*nev;l++){
l                 518 src/dynacont.c       zetaj=zeta[l];
l                 519 src/dynacont.c       dj=d[l];
l                 523 src/dynacont.c       if(fabs(d[l])<=1.e-10){
l                 532 src/dynacont.c 	bjinc[l]=sum;
l                 533 src/dynacont.c 	bjincp[l]=sump;
l                 552 src/dynacont.c 	bjinc[l]=sum/ddj;
l                 553 src/dynacont.c 	bjincp[l]=sump/ddj;
l                 575 src/dynacont.c 	bjinc[l]=sum/(2.*ddj);
l                 576 src/dynacont.c 	bjincp[l]=sump/(2.*ddj);
l                 595 src/dynacont.c 	bjinc[l]=sum;
l                 596 src/dynacont.c 	bjincp[l]=sump;
l                 775 src/dynacont.c 		  for(l=0; l<*nev; l++){
l                 776 src/dynacont.c 		      bbb[l**nev+j]+=zl*dbcont[i2+l];
l                 790 src/dynacont.c 			  for(l=0; l<*nev; l++){
l                 791 src/dynacont.c 			      bbb[l**nev+j]+=zl*dbcont[i2+l];
l                 799 src/dynacont.c       for(l=0;l<*nev;l++){
l                 800 src/dynacont.c 	i1=l**nev;
l                 802 src/dynacont.c 	  if(j==l){lhs[i1+j]=1.;}else{lhs[i1+j]=0.;}
l                 805 src/dynacont.c 	dbj[l]=bjbas[l]+bjinc[l]*aaa[l]-bj[l];
l                  63 src/exo.c        ITG i,j,k,l,m,n,o,indexe,nemax,nlayer,noutloc,iset,iselect,ncomp,nope,
l                 241 src/exo.c          l=0;
l                 259 src/exo.c            elem_map[l] = i+1;
l                 269 src/exo.c      	  blkassign[l++]=1;
l                 279 src/exo.c      	    elem_map[l] = i+1;
l                 280 src/exo.c      	    blkassign[l++]=2;
l                 284 src/exo.c      	  blkassign[l++]=3;
l                 287 src/exo.c      	  blkassign[l++]=4;
l                 293 src/exo.c      	  blkassign[l++]=5;
l                 297 src/exo.c      	    blkassign[l++]=6;
l                 299 src/exo.c      	    blkassign[l++]=7;
l                 304 src/exo.c      	    blkassign[l++]=8;
l                 306 src/exo.c      	    blkassign[l++]=9;
l                 312 src/exo.c      	blkassign[l++]=10;
l                 315 src/exo.c      	blkassign[l++]=11;
l                 320 src/exo.c      	  blkassign[l++]=12;
l                 323 src/exo.c      	  blkassign[l++]=13;
l                 329 src/exo.c      	  blkassign[l++]=14;
l                 332 src/exo.c      	  blkassign[l++]=15;
l                 339 src/exo.c      	  blkassign[l++]=16;
l                 342 src/exo.c      	  blkassign[l++]=17;
l                 345 src/exo.c      	  blkassign[l++]=18;
l                 350 src/exo.c      	blkassign[l++]=19;
l                 389 src/exo.c          for(l=0;l<num_elem_blk;l++){
l                 395 src/exo.c      	if(blkassign[m]==l){
l                 408 src/exo.c            blksize[l]=j;
l                 409 src/exo.c            num_elem_in_blk=blksize[l];
l                 411 src/exo.c            connect = (ITG *) calloc (num_elem_in_blk*num_nodes_per_elem[l], sizeof(ITG));
l                 418 src/exo.c      	if (blkassign[o]==l){
l                 430 src/exo.c      	    for(l=0;l<mi[2];l++){
l                 431 src/exo.c      	      if(ielmat[i*mi[2]+l]==0) break;
l                 440 src/exo.c      	    for (j = 0; j <num_nodes_per_elem[l]; j++){
l                 444 src/exo.c      	    for (j = 0; j <num_nodes_per_elem[l]; j++){
l                 453 src/exo.c            switch (l)
l                 456 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "SPHERE", num_elem_in_blk, num_nodes_per_elem[l], num_attr);
l                 459 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "HEX", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 462 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "HEX", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 465 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "TRUSS", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 468 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "QUAD", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 471 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "HEX", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 474 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "TRUSS", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 477 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "TRUSS", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 480 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "SHELL", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 483 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "SHELL", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 486 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "TETRA", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 489 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "TETRA", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 492 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "WEDGE", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 495 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "HEX", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 498 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "WEDGE", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 501 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "WEDGE", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 508 src/exo.c      	  errr = ex_put_elem_block (exoid, l, "TRUSS", num_elem_in_blk, num_nodes_per_elem[l], num_attr);	  
l                 515 src/exo.c      	errr = ex_put_elem_conn (exoid, l, connect);
l                  47 src/exoselect.c   ITG i,j,k,l,m,o,nksegment;
l                  94 src/exoselect.c 	  for(l=0;l<*ngraph;l++){
l                  95 src/exoselect.c 	    i=ialset[k]+l*nksegment-1;
l                 117 src/exoselect.c 	  l=ialset[k-2];
l                 119 src/exoselect.c 	    l-=ialset[k];
l                 120 src/exoselect.c 	    if(l>=ialset[k-1]) break;
l                 123 src/exoselect.c 	      i=l+m*nksegment-1;
l                  32 src/exosetfind.c   ITG i,j,k,l,n,s,e,gen;
l                  53 src/exosetfind.c       gen=0; l=0;
l                  57 src/exosetfind.c 	  l-=3;
l                  62 src/exosetfind.c       l=e-s+1+gen+l;
l                  63 src/exosetfind.c       set_nums = (ITG *) calloc(l, sizeof(ITG));
l                  68 src/exosetfind.c       if (l>2){
l                  90 src/exosetfind.c       }else if(l>1){
l                  35 src/exovector.c   ITG i,j,k,l,m,n,ii,jj,kk;
l                  87 src/exovector.c 	  for(l=0;l<*ngraph;l++){
l                  88 src/exovector.c 	    i=ialset[k]+l*nksegment-1;
l                 101 src/exovector.c 	  l=ialset[k-2];
l                 103 src/exovector.c 	    l-=ialset[k];
l                 104 src/exovector.c 	    if(l>=ialset[k-1]) break;
l                 107 src/exovector.c 	      i=l+m*nksegment-1;      
l                  72 src/expand.c           ielas,icmd,kk,l,nkt,icntrl,imag=1,icomplex,kkv,kk6,iterm,
l                 251 src/expand.c   	    for(l=0;l<*nk;l++){
l                 252 src/expand.c   		if(inocs[l]==jj){
l                 253 src/expand.c   		    co[3*l+i*3**nk]=co[3*l];
l                 254 src/expand.c   		    co[1+3*l+i*3**nk]=co[1+3*l]+theta;
l                 255 src/expand.c   		    co[2+3*l+i*3**nk]=co[2+3*l];
l                 256 src/expand.c   		    if(*ntrans>0) inotr[2*l+i*2**nk]=inotr[2*l];
l                 259 src/expand.c   	    for(l=0;l<*nkon;l++){kon[l+i**nkon]=kon[l]+i**nk;}
l                 260 src/expand.c   	    for(l=0;l<*ne;l++){
l                 261 src/expand.c   		if(ielcs[l]==jj){
l                 262 src/expand.c   		    if(ipkon[l]>=0){
l                 263 src/expand.c   			ipkon[l+i**ne]=ipkon[l]+i**nkon;
l                 264 src/expand.c   			ielmat[mi[2]*(l+i**ne)]=ielmat[mi[2]*l];
l                 265 src/expand.c   			if(*norien>0) ielorien[l+i**ne]=ielorien[l];
l                 267 src/expand.c   			    l2=8*l+l1;
l                 271 src/expand.c   			ipkon[l+i**ne]=ipkon[l];
l                 554 src/expand.c   		l=mt*l1+l2;
l                 555 src/expand.c   		vt[l]=v[l];
l                 556 src/expand.c   		if(*nherm!=1)vti[l]=v[l+mt**nk];
l                 582 src/expand.c   				    l=mt*l1+l2;
l                 583 src/expand.c   				    vt[l+mt**nk*i]=v[l];
l                 584 src/expand.c   				    if(*nherm!=1)vti[l+mt**nk*i]=v[l+mt**nk];
l                 590 src/expand.c   			    l=mt*l1+l2;
l                 591 src/expand.c   			    vt[l+mt**nk*i]=ctl*v[l]-stl*v[l+mt**nk];
l                 592 src/expand.c   			    if(*nherm!=1)vti[l+mt**nk*i]=stl*v[l]+ctl*v[l+mt**nk];
l                 893 src/expand.c   	      for(l=1; l<*nsectors; l++){
l                 895 src/expand.c   		  ialset[tint-1]=(ialset[istartset[j]-1+k]!=-1)?ialset[istartset[j]-1+k]+*nk*l:-1;
l                 912 src/expand.c   	      for(l=1; l<*nsectors; l++){
l                 915 src/expand.c   		  ialset[tint-1]=(ialset[istartset[j]-1+k]!=-1)?(tint2+*ne*l)*10+(ialset[istartset[j]-1+k]-(tint2*10)):-1;
l                  40 src/frdcyc.c     ITG nkt,icntrl,*kont=NULL,*ipkont=NULL,*inumt=NULL,*ielmatt=NULL,net,i,l,
l                 163 src/frdcyc.c     for(l=0;l<3**nk;l++){cot[l]=co[l];}
l                 164 src/frdcyc.c     if(*ntrans>0){for(l=0;l<*nk;l++){inotrt[2*l]=inotr[2*l];}}
l                 169 src/frdcyc.c         for(l=0;l<*nkon;l++){kont[l]=kon[l];}
l                 170 src/frdcyc.c         for(l=0;l<*ne;l++){ipkont[l]=ipkon[l];}
l                 171 src/frdcyc.c         for(l=0;l<8**ne;l++){lakont[l]=lakon[l];}
l                 172 src/frdcyc.c         for(l=0;l<mi[2]**ne;l++){ielmatt[l]=ielmat[l];}
l                 187 src/frdcyc.c         for(l=0;l<*nk;l++){
l                 188 src/frdcyc.c           if(inocs[l]==jj){
l                 189 src/frdcyc.c   	  cot[3*l+i*3**nk]=cot[3*l];
l                 190 src/frdcyc.c   	  cot[1+3*l+i*3**nk]=cot[1+3*l]+theta;
l                 191 src/frdcyc.c   	  cot[2+3*l+i*3**nk]=cot[2+3*l];
l                 196 src/frdcyc.c   	  for(l=0;l<*nk;l++){
l                 197 src/frdcyc.c   	      if(inocs[l]==jj){
l                 198 src/frdcyc.c   		  inotrt[2*l+i*2**nk]=inotrt[2*l];
l                 205 src/frdcyc.c           for(l=0;l<*nkon;l++){kont[l+i**nkon]=kon[l]+i**nk;}
l                 206 src/frdcyc.c           for(l=0;l<*ne;l++){
l                 207 src/frdcyc.c             if(ielcs[l]==jj){
l                 208 src/frdcyc.c               if(ipkon[l]>=0){
l                 209 src/frdcyc.c                 ipkont[l+i**ne]=ipkon[l]+i**nkon;
l                 210 src/frdcyc.c                 ielmatt[mi[2]*(l+i**ne)]=ielmat[mi[2]*l];
l                 212 src/frdcyc.c                   l2=8*l+l1;
l                 216 src/frdcyc.c               else ipkont[l+i**ne]=-1;
l                 230 src/frdcyc.c     for(l=0;l<*nk;l++){inumt[l]=inum[l];}
l                 238 src/frdcyc.c       for(l=0;l<mt**nk;l++){vt[l]=v[l];};
l                 240 src/frdcyc.c       for(l=0;l<*nk;l++){t1t[l]=t1[l];};
l                 242 src/frdcyc.c       for(l=0;l<6**nk;l++){stnt[l]=stn[l];};
l                 244 src/frdcyc.c       for(l=0;l<6**nk;l++){eent[l]=een[l];};
l                 246 src/frdcyc.c       for(l=0;l<mt**nk;l++){fnt[l]=fn[l];};
l                 248 src/frdcyc.c       for(l=0;l<*nk;l++){epnt[l]=epn[l];};
l                 250 src/frdcyc.c       for(l=0;l<*nk;l++){enernt[l]=enern[l];};
l                 252 src/frdcyc.c       for(l=0;l<*nstate_**nk;l++){xstatent[l]=xstaten[l];};
l                 254 src/frdcyc.c       for(l=0;l<3**nk;l++){qfnt[l]=qfn[l];};
l                 257 src/frdcyc.c       for(l=0;l<6*mi[0]**ne;l++){stit[l]=sti[l];};
l                 259 src/frdcyc.c       for(l=0;l<6**nk;l++){emnt[l]=emn[l];};
l                 265 src/frdcyc.c         for(l=0;l<*nk;l++){inumt[l+i**nk]=inum[l];}
l                 272 src/frdcyc.c                 l=mt*l1+l2;
l                 273 src/frdcyc.c                 vt[l+mt**nk*i]=v[l];
l                 280 src/frdcyc.c           for(l=0;l<*nk;l++){
l                 281 src/frdcyc.c             if(inocs[l]==jj) t1t[l+*nk*i]=t1[l];
l                 289 src/frdcyc.c                 l=6*l1+l2;
l                 290 src/frdcyc.c                 stnt[l+6**nk*i]=stn[l];
l                 300 src/frdcyc.c                 l=6*l1+l2;
l                 301 src/frdcyc.c                 eent[l+6**nk*i]=een[l];
l                 311 src/frdcyc.c                 l=mt*l1+l2;
l                 312 src/frdcyc.c                 fnt[l+mt**nk*i]=fn[l];
l                 319 src/frdcyc.c           for(l=0;l<*nk;l++){
l                 320 src/frdcyc.c             if(inocs[l]==jj) epnt[l+*nk*i]=epn[l];
l                 325 src/frdcyc.c           for(l=0;l<*nk;l++){
l                 326 src/frdcyc.c             if(inocs[l]==jj) enernt[l+*nk*i]=enern[l];
l                 334 src/frdcyc.c                 l=*nstate_*l1+l2;
l                 335 src/frdcyc.c                 xstatent[l+*nstate_**nk*i]=xstaten[l];
l                 345 src/frdcyc.c                 l=3*l1+l2;
l                 346 src/frdcyc.c                 qfnt[l+3**nk*i]=qfn[l];
l                 356 src/frdcyc.c                 l=6*l1+l2;
l                 357 src/frdcyc.c                 emnt[l+6**nk*i]=emn[l];
l                  43 src/frdselect.c   ITG i,j,k,l,m,n,nksegment;
l                 118 src/frdselect.c 	for(l=0;l<*ngraph;l++){
l                 119 src/frdselect.c 	  i=ialset[k]+l*nksegment-1;
l                 184 src/frdselect.c 	l=ialset[k-2];
l                 186 src/frdselect.c 	  l-=ialset[k];
l                 187 src/frdselect.c 	  if(l>=ialset[k-1]) break;
l                 189 src/frdselect.c 	    i=l+m*nksegment-1;
l                  29 src/frdvector.c   ITG i,k,l,m,nksegment;
l                  89 src/frdvector.c 	for(l=0;l<*ngraph;l++){
l                  90 src/frdvector.c 	  i=ialset[k]+l*nksegment-1;
l                 121 src/frdvector.c 	l=ialset[k-2];
l                 123 src/frdvector.c 	  l-=ialset[k];
l                 124 src/frdvector.c 	  if(l>=ialset[k-1]) break;
l                 126 src/frdvector.c 	    i=l+m*nksegment-1;
l                 138 src/getglobalresults.c     anz->l=0;
l                 380 src/getglobalresults.c     for(i=0;i<anz[0].l;i++){
l                 418 src/getglobalresults.c     for(i=0;i<anz[0].l;i++){
l                 459 src/getglobalresults.c     for(i=0;i<anz[0].l;i++){
l                 502 src/getglobalresults.c     for(i=0;i<anz[0].l;i++){
l                 539 src/getglobalresults.c     for(j=0;j<anz->l;j++){
l                  41 src/mastruct.c   ITG i,j,k,l,jj,ll,id,index,jdof1,jdof2,idof1,idof2,mpc1,mpc2,id1,id2,
l                 256 src/mastruct.c 	  l=ll/mi[1];
l                 257 src/mastruct.c 	  m=ll-mi[1]*l;
l                 259 src/mastruct.c 	  node2=kon[indexe+l];
l                 424 src/mastruct.c 	  l=ll;
l                 426 src/mastruct.c 	  node2=kon[indexe+l];
l                 638 src/mastruct.c 	  for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
l                 645 src/mastruct.c     for(l=k;l<neq[1]+1;++l){jq[l]=nmast+1;}
l                 675 src/mastruct.c 		for(l=k;l<irow[i];++l){jq[l]=i+1-isubtract;}
l                 681 src/mastruct.c 	for(l=k;l<neq[2]+1;++l){jq[l]=nmastboun+1;}
l                  41 src/mastructcs.c   ITG i,j,k,l,jj,ll,id,index,jdof1,jdof2,idof1,idof2,mpc1,mpc2,id1,id2,
l                 185 src/mastructcs.c 	l=ll/3;
l                 186 src/mastructcs.c 	m=ll-3*l;
l                 188 src/mastructcs.c 	node2=kon[indexe+l];
l                 510 src/mastructcs.c 	for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
l                 517 src/mastructcs.c   for(l=k;l<neq[0]+1;++l){jq[l]=nmast+1;}
l                  40 src/mastructem.c   ITG i,j,k,l,jj,ll,id,index,jdof1,jdof2,idof1,idof2,mpc1,mpc2,id1,id2,
l                 241 src/mastructem.c 	  l=ll/mi[1];
l                 242 src/mastructem.c 	  m=ll-mi[1]*l;
l                 244 src/mastructem.c 	  node2=kon[indexe+l];
l                 400 src/mastructem.c 	  l=ll;
l                 402 src/mastructem.c 	  node2=kon[indexe+l];
l                 592 src/mastructem.c 	  for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
l                 599 src/mastructem.c     for(l=k;l<neq[1]+1;++l){jq[l]=nmast+1;}
l                  32 src/mastructf.c   ITG i,j,k,l,index,idof1,idof2,node1,isubtract,nmast,ifree=0,istart,istartold,
l                 131 src/mastructf.c 	      for(l=k;l<mast1[i];++l){jq[l]=i+1-isubtract;}
l                 138 src/mastructf.c   for(l=k;l<*neq+1;++l){jq[l]=nmast+1;}
l                  35 src/mastructrad.c   ITG three=3,i,j,k,l,ii,jj,icav,jcav,*mast1rad=NULL,*irowrad=NULL,
l                 123 src/mastructrad.c 	      for(l=k;l<mast1rad[i];++l){jqrad[l]=i+1-isubtract;}
l                 130 src/mastructrad.c   for(l=k;l<*ntr+1;++l){jqrad[l]=nmast+1;}
l                  34 src/matrixstorage.c   ITG i,j,l,*irow=NULL,*ai=NULL,*aj=NULL,kflag=2,ndim,jref,kstart,klen,
l                  98 src/matrixstorage.c     l=0;
l                 102 src/matrixstorage.c         aj[k]=irow[l];
l                 103 src/matrixstorage.c         aa[k]=au[l];
l                 104 src/matrixstorage.c         k++;l++;
l                 139 src/matrixstorage.c     l=0;
l                 143 src/matrixstorage.c         aj[k]=irow[l];
l                 144 src/matrixstorage.c         aa[k]=au[l];
l                 146 src/matrixstorage.c         ai[k]=irow[l];
l                 148 src/matrixstorage.c         aa[k]=au[l];
l                 150 src/matrixstorage.c         l++;
l                 328 src/matrixstorage.c     l=0;
l                 332 src/matrixstorage.c         aj[k]=irow[l];
l                 333 src/matrixstorage.c         aa[k]=aub[l];
l                 334 src/matrixstorage.c         k++;l++;
l                 363 src/matrixstorage.c     l=0;
l                 367 src/matrixstorage.c         aj[k]=irow[l];
l                 368 src/matrixstorage.c         aa[k]=aub[l];
l                 370 src/matrixstorage.c         ai[k]=irow[l];
l                 372 src/matrixstorage.c         aa[k]=aub[l];
l                 374 src/matrixstorage.c         l++;
l                  40 src/pardiso.c    ITG i,j,k,l,maxfct=1,mnum=1,phase=12,nrhs=1,*perm=NULL,mtype,
l                  94 src/pardiso.c        l=*nzs;
l                 100 src/pardiso.c  		  icolpardiso[--k]=irow[--l];
l                 101 src/pardiso.c  		  aupardiso[k]=au[l];
l                 112 src/pardiso.c  		  icolpardiso[--k]=irow[--l];
l                 113 src/pardiso.c  		  aupardiso[k]=au[l]-*sigma*aub[l];
l                 269 src/pardiso.c  	      l=k+1;
l                 280 src/pardiso.c  	      pointers[i+1]=l;
l                 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;
l                 407 src/pcgsolver.c 				for (l=llo; l<=lup; l++)
l                 409 src/pcgsolver.c 					if (ia[l]>ia[j])	break;
l                 410 src/pcgsolver.c 					if (ia[l]<ia[j])	continue;
l                 411 src/pcgsolver.c 					C[k] -= C[j]*C[l];
l                  38 src/precontact.c         l,nstart,kflag,ntri,ii;
l                 140 src/precontact.c 	    for(l=itiefac[2*i];l<=itiefac[2*i+1];l++){
l                 146 src/precontact.c 	        	mi,ncont,ipe,ime,pslavsurf,&ii,&l,&ntri));
l                  35 src/radcyc.c        node,i1,i2,nope,iel,indexe,j,k,ielset,node1,node2,node3,l,jj;
l                 145 src/radcyc.c         for(l=0;l<*nk;l++){
l                 146 src/radcyc.c           if(inocs[l]==jj){
l                 147 src/radcyc.c   	  co[3*l+i*3**nk]=co[3*l];
l                 148 src/radcyc.c   	  co[1+3*l+i*3**nk]=co[1+3*l]-theta;
l                 149 src/radcyc.c   	  co[2+3*l+i*3**nk]=co[2+3*l];
l                 119 src/readfrd.c    anz->u=anz->n=anz->e=anz->l=-1;
l                 534 src/readfrd.c        anz->l++;
l                 536 src/readfrd.c        printf ("reading Dataset No:%" ITGFORMAT "\n",anz->l+1);
l                 537 src/readfrd.c        if ( (lcase = (Datasets *)realloc((Datasets *)lcase, (anz->l+2) * sizeof(Datasets))) == NULL )
l                 540 src/readfrd.c        lcase[anz->l].handle=(FILE *)NULL;
l                 543 src/readfrd.c        lcase[anz->l].npheader=anz_p+1;
l                 544 src/readfrd.c        lcase[anz->l].pheader=pheader;
l                 545 src/readfrd.c        lcase[anz->l].fileptr=NULL;
l                 546 src/readfrd.c        lcase[anz->l].loaded=1;
l                 547 src/readfrd.c        lcase[anz->l].format_flag=format_flag;
l                 552 src/readfrd.c        stos( rec_str, 7, 12, lcase[anz->l].dataset_name);
l                 553 src/readfrd.c        lcase[anz->l].value=stof( rec_str, 13, 24 );
l                 571 src/readfrd.c        stos( rec_str, 37, 56, lcase[anz->l].dataset_text);
l                 572 src/readfrd.c        lcase[anz->l].analysis_type=stoi( rec_str, 57, 58 );
l                 573 src/readfrd.c        lcase[anz->l].step_number=stoi( rec_str, 59, 63 );
l                 576 src/readfrd.c          stos(rec_str,64,73,lcase[anz->l].analysis_name);
l                 580 src/readfrd.c          strcpy(lcase[anz->l].analysis_name,"");
l                 590 src/readfrd.c        if(lcase[anz->l].analysis_type>=2)
l                 592 src/readfrd.c          for(i=0;i<lcase[anz->l].npheader; i++)
l                 594 src/readfrd.c            if(compare(&lcase[anz->l].pheader[i][5],"PHID", 4)==4)
l                 596 src/readfrd.c              sscanf(lcase[anz->l].pheader[i],"%*s %" ITGFORMAT "", &ipuf);
l                 597 src/readfrd.c              if(ipuf>-1) sprintf(lcase[anz->l].dataset_text,"ND:%" ITGFORMAT "",ipuf);
l                 600 src/readfrd.c          if(ipuf!=-1) for(i=0;i<lcase[anz->l].npheader; i++)
l                 602 src/readfrd.c            if(compare(&lcase[anz->l].pheader[i][5],"PMODE", 5)==5)
l                 604 src/readfrd.c              sscanf(lcase[anz->l].pheader[i],"%*s %" ITGFORMAT "", &ipuf);
l                 605 src/readfrd.c              if(ipuf>-1) sprintf(&lcase[anz->l].dataset_text[strlen(lcase[anz->l].dataset_text)]," MODE:%" ITGFORMAT "",ipuf);
l                 613 src/readfrd.c          if (( format_flag==2)&&(lcase[anz->l].ncomps>0)&&(ncomps==lcase[anz->l].ncomps))
l                 620 src/readfrd.c              lcase[anz->l].loaded=0;
l                 626 src/readfrd.c                lcase[anz->l].fileptr=filepntr;
l                 628 src/readfrd.c                lcase[anz->l].handle=handle;
l                 629 src/readfrd.c                strcpy(lcase[anz->l].filename,datin);
l                 636 src/readfrd.c              if ( (value = (float *)realloc((float *)value, (lcase[anz->l].ncomps) * sizeof(float))) == NULL )
l                 641 src/readfrd.c                length=fread((float *)value,sizeof(float),lcase[anz->l].ncomps,handle);
l                 643 src/readfrd.c                for(i=0; i<lcase[anz->l].ncomps; i++)
l                 646 src/readfrd.c                  lcase[anz->l].dat[i][nodenr]= value[i];
l                 665 src/readfrd.c              lcase[anz->l].loaded=0;
l                 671 src/readfrd.c                lcase[anz->l].fileptr=filepntr;
l                 673 src/readfrd.c                lcase[anz->l].handle=handle;
l                 674 src/readfrd.c                strcpy(lcase[anz->l].filename,datin);
l                 688 src/readfrd.c              else if ( lcase[anz->l].irtype == 1 )
l                 695 src/readfrd.c                    lcase[anz->l].dat[0][nodenr]= stof(&rec_str[  13  ], 1, 12);
l                 696 src/readfrd.c                    lcase[anz->l].dat[1][nodenr]= stof(&rec_str[  25  ], 1, 12);
l                 697 src/readfrd.c                    lcase[anz->l].dat[2][nodenr]= stof(&rec_str[  37  ], 1, 12);
l                 698 src/readfrd.c                    lcase[anz->l].dat[3][nodenr]= stof(&rec_str[  49  ], 1, 12);
l                 699 src/readfrd.c                    lcase[anz->l].dat[4][nodenr]= stof(&rec_str[  61  ], 1, 12);
l                 700 src/readfrd.c                    lcase[anz->l].dat[5][nodenr]= stof(&rec_str[  73  ], 1, 12);
l                 704 src/readfrd.c                    lcase[anz->l].dat[0][nodenr]= stof(&rec_str[  8   ], 1, 12);
l                 705 src/readfrd.c                    lcase[anz->l].dat[1][nodenr]= stof(&rec_str[  20  ], 1, 12);
l                 706 src/readfrd.c                    lcase[anz->l].dat[2][nodenr]= stof(&rec_str[  32  ], 1, 12);
l                 707 src/readfrd.c                    lcase[anz->l].dat[3][nodenr]= stof(&rec_str[  44  ], 1, 12);
l                 708 src/readfrd.c                    lcase[anz->l].dat[4][nodenr]= stof(&rec_str[  56  ], 1, 12);
l                 709 src/readfrd.c                    lcase[anz->l].dat[5][nodenr]= stof(&rec_str[  68  ], 1, 12);
l                 714 src/readfrd.c  	        for(i=0; i<maxcomps; i++) lcase[anz->l].dat[i][nodenr]= stof(&rec_str[n+i*12], 1, 12);
l                 732 src/readfrd.c                lcase[anz->l].dat[i][nodenr]= stof(&rec_str[n+j*12], 1, 12);
l                 734 src/readfrd.c              }while((j<6)&&(i<lcase[anz->l].ncomps));
l                 739 src/readfrd.c            stos( rec_str, 6, 13, lcase[anz->l].name);
l                 740 src/readfrd.c            lcase[anz->l].ncomps = stoi(rec_str,14,18);
l                 741 src/readfrd.c            lcase[anz->l].irtype = stoi(rec_str,19,23);
l                 743 src/readfrd.c            if( lcase[anz->l].irtype > 2 )
l                 746 src/readfrd.c              anz->l--;
l                 750 src/readfrd.c            if ( (lcase[anz->l].nmax = (ITG *)malloc( (lcase[anz->l].ncomps) * sizeof(ITG))) == NULL )
l                 752 src/readfrd.c            if ( (lcase[anz->l].nmin = (ITG *)malloc( (lcase[anz->l].ncomps) * sizeof(ITG))) == NULL )
l                 754 src/readfrd.c            if ( (lcase[anz->l].max = (float *)malloc( (lcase[anz->l].ncomps) * sizeof(float))) == NULL )
l                 756 src/readfrd.c            if ( (lcase[anz->l].min = (float *)malloc( (lcase[anz->l].ncomps) * sizeof(float))) == NULL )
l                 758 src/readfrd.c            if ( (lcase[anz->l].compName = (char **)malloc( (lcase[anz->l].ncomps) * sizeof(char *))) == NULL )
l                 760 src/readfrd.c            if ( (lcase[anz->l].icname = (char **)malloc( (lcase[anz->l].ncomps) * sizeof(char *))) == NULL )
l                 762 src/readfrd.c            if ( (lcase[anz->l].menu = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 764 src/readfrd.c            if ( (lcase[anz->l].ictype = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 766 src/readfrd.c            if ( (lcase[anz->l].icind1 = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 768 src/readfrd.c            if ( (lcase[anz->l].icind2 = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 770 src/readfrd.c            if ( (lcase[anz->l].iexist = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 772 src/readfrd.c            if ( (lcase[anz->l].dat = (float **)malloc( (lcase[anz->l].ncomps) * sizeof(float *))) == NULL )
l                 774 src/readfrd.c    printf(" gen lc[%" ITGFORMAT "] ncomps:%" ITGFORMAT "\n",anz->l,lcase[anz->l].ncomps);
l                 775 src/readfrd.c            for(i=0; i<(lcase[anz->l].ncomps); i++)
l                 777 src/readfrd.c              if ( (lcase[anz->l].compName[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
l                 779 src/readfrd.c              if ( (lcase[anz->l].icname[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
l                 781 src/readfrd.c              lcase[anz->l].max[i]=-MAX_FLOAT;
l                 782 src/readfrd.c              lcase[anz->l].min[i]=MAX_FLOAT;
l                 787 src/readfrd.c            if(ncomps<lcase[anz->l].ncomps)
l                 789 src/readfrd.c              stos(rec_str, 6, 13, lcase[anz->l].compName[ncomps]);
l                 790 src/readfrd.c              lcase[anz->l].menu[ncomps] = stoi(rec_str,14,18);
l                 791 src/readfrd.c              lcase[anz->l].ictype[ncomps] = stoi(rec_str,19,23);
l                 792 src/readfrd.c              lcase[anz->l].icind1[ncomps] = stoi(rec_str,24,28);
l                 793 src/readfrd.c              lcase[anz->l].icind2[ncomps] = stoi(rec_str,29,33);
l                 794 src/readfrd.c              lcase[anz->l].iexist[ncomps] = stoi(rec_str,34,38);
l                 797 src/readfrd.c              if(!lcase[anz->l].iexist[ncomps]) ncomps++;
l                 800 src/readfrd.c                SFREE(lcase[anz->l].compName[ncomps]); lcase[anz->l].compName[ncomps]=NULL;
l                 801 src/readfrd.c                SFREE(lcase[anz->l].icname[ncomps]); lcase[anz->l].icname[ncomps]=NULL;
l                 807 src/readfrd.c  	    printf(" WARNING: unallocated component:%" ITGFORMAT " \"%s\" %" ITGFORMAT "\n", ncomps, &rec_str[5],lcase[anz->l].ncomps);
l                 814 src/readfrd.c            if(nentities==lcase[anz->l].ncomps)
l                 816 src/readfrd.c              lcase[anz->l].ncomps=ncomps;
l                 817 src/readfrd.c              if(lcase[anz->l].ncomps<6) maxcomps=lcase[anz->l].ncomps;
l                 821 src/readfrd.c              for (i=0; i<ncomps; i++) if(lcase[anz->l].iexist[i]!=1) nvals++;
l                 854 src/readfrd.c                for(i=0; i<(lcase[anz->l].ncomps); i++)
l                 856 src/readfrd.c                  if ( (lcase[anz->l].dat[i] = (float *)malloc( (anz->nmax+1) * sizeof(float))) == NULL )
l                 858 src/readfrd.c                  for(j=0; j<=anz->nmax; j++) lcase[anz->l].dat[i][j]=0.;
l                 875 src/readfrd.c            for (i=0; i<ipuf; i++) strcpy(compName[i], lcase[anz->l].compName[ atoi(dat[i+3])-1 ]);
l                 876 src/readfrd.c            for (i=0; i<ipuf; i++) strcpy(lcase[anz->l].compName[i],compName[i]);
l                 887 src/readfrd.c            for (i=0; i<ipuf; i++) menu[i] = lcase[anz->l].menu[atoi(dat[i+3])-1];
l                 888 src/readfrd.c            for (i=0; i<ipuf; i++) lcase[anz->l].menu[i] =menu[i];
l                 889 src/readfrd.c            for (i=0; i<ipuf; i++) ictype[i] = lcase[anz->l].ictype[atoi(dat[i+3])-1];
l                 890 src/readfrd.c            for (i=0; i<ipuf; i++) lcase[anz->l].ictype[i] =ictype[i];
l                 891 src/readfrd.c            for (i=0; i<ipuf; i++) icind1[i] = lcase[anz->l].icind1[atoi(dat[i+3])-1];
l                 892 src/readfrd.c            for (i=0; i<ipuf; i++) lcase[anz->l].icind1[i] =icind1[i];
l                 893 src/readfrd.c            for (i=0; i<ipuf; i++) icind2[i] = lcase[anz->l].icind2[atoi(dat[i+3])-1];
l                 894 src/readfrd.c            for (i=0; i<ipuf; i++) lcase[anz->l].icind2[i] =icind2[i];
l                 895 src/readfrd.c            for (i=0; i<ipuf; i++) iexist[i] = lcase[anz->l].iexist[atoi(dat[i+3])-1];
l                 896 src/readfrd.c            for (i=0; i<ipuf; i++) lcase[anz->l].iexist[i] =iexist[i];
l                 916 src/readfrd.c          for(i=0; i<lcase[anz->l].ncomps; i++)
l                 918 src/readfrd.c            if(lcase[anz->l].dat[i][node[j].nr] > lcase[anz->l].max[i])
l                 920 src/readfrd.c              lcase[anz->l].max[i]=lcase[anz->l].dat[i][node[j].nr];
l                 921 src/readfrd.c              lcase[anz->l].nmax[i]=node[j].nr;
l                 923 src/readfrd.c            if(lcase[anz->l].dat[i][node[j].nr] < lcase[anz->l].min[i])
l                 925 src/readfrd.c              lcase[anz->l].min[i]=lcase[anz->l].dat[i][node[j].nr];
l                 926 src/readfrd.c              lcase[anz->l].nmin[i]=node[j].nr;
l                 935 src/readfrd.c        anz->l++;
l                 936 src/readfrd.c        printf ("reading Dataset No:%" ITGFORMAT "\n",anz->l+1);
l                 937 src/readfrd.c        if ( (lcase = (Datasets *)realloc((Datasets *)lcase, (anz->l+2) * sizeof(Datasets))) == NULL )
l                 940 src/readfrd.c        lcase[anz->l].value=stof( rec_str, 13, 25 );
l                 941 src/readfrd.c        strcpy(lcase[anz->l].name,"DISP    ");
l                 942 src/readfrd.c        lcase[anz->l].ncomps = 3;
l                 943 src/readfrd.c        lcase[anz->l].irtype = 1;
l                 944 src/readfrd.c        lcase[anz->l].loaded = 1;
l                 945 src/readfrd.c        lcase[anz->l].format_flag=format_flag;
l                 947 src/readfrd.c        if ( (lcase[anz->l].nmax = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 949 src/readfrd.c        if ( (lcase[anz->l].nmin = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 951 src/readfrd.c        if ( (lcase[anz->l].max = (float *)malloc( lcase[anz->l].ncomps * sizeof(float))) == NULL )
l                 953 src/readfrd.c        if ( (lcase[anz->l].min = (float *)malloc( lcase[anz->l].ncomps * sizeof(float))) == NULL )
l                 955 src/readfrd.c        if ( (lcase[anz->l].dat = (float **)malloc( lcase[anz->l].ncomps * sizeof(float *))) == NULL )
l                 957 src/readfrd.c        if ( (lcase[anz->l].compName = (char **)malloc( lcase[anz->l].ncomps * sizeof(char *))) == NULL )
l                 959 src/readfrd.c        if ( (lcase[anz->l].icname = (char **)malloc( lcase[anz->l].ncomps * sizeof(char *))) == NULL )
l                 961 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
l                 963 src/readfrd.c          if ( (lcase[anz->l].compName[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
l                 965 src/readfrd.c          if ( (lcase[anz->l].icname[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
l                 967 src/readfrd.c          lcase[anz->l].max[i]=-MAX_FLOAT;
l                 968 src/readfrd.c          lcase[anz->l].min[i]=MAX_FLOAT;
l                 969 src/readfrd.c          if ( (lcase[anz->l].dat[i] = (float *)malloc( (anz->nmax+1) * sizeof(float))) == NULL )
l                 971 src/readfrd.c          for(j=0; j<=anz->nmax; j++) lcase[anz->l].dat[i][j]=0.;
l                 973 src/readfrd.c        if ( (lcase[anz->l].menu = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 975 src/readfrd.c        if ( (lcase[anz->l].ictype = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 977 src/readfrd.c        if ( (lcase[anz->l].icind1 = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 979 src/readfrd.c        if ( (lcase[anz->l].icind2 = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 981 src/readfrd.c        if ( (lcase[anz->l].iexist = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                 984 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
l                 986 src/readfrd.c          lcase[anz->l].menu[i] = 1;
l                 987 src/readfrd.c          lcase[anz->l].ictype[i] = 2;
l                 988 src/readfrd.c          lcase[anz->l].icind1[i] = i+1;
l                 989 src/readfrd.c          lcase[anz->l].icind2[i] = 0;
l                 990 src/readfrd.c          lcase[anz->l].iexist[i] = 0;
l                 993 src/readfrd.c        strcpy(lcase[anz->l].compName[0], "x       ");
l                 994 src/readfrd.c        strcpy(lcase[anz->l].compName[1], "y       ");
l                 995 src/readfrd.c        strcpy(lcase[anz->l].compName[2], "z       ");
l                1012 src/readfrd.c              for(i=0; i<lcase[anz->l].ncomps; i++)
l                1013 src/readfrd.c                lcase[anz->l].dat[i][nodenr]= stof(&rec_str[8+i*12], 1, 12);
l                1017 src/readfrd.c              for(i=0; i<lcase[anz->l].ncomps; i++)
l                1018 src/readfrd.c                lcase[anz->l].dat[i][nodenr]= stof(&rec_str[13+i*12], 1, 12);
l                1025 src/readfrd.c          for(i=0; i<lcase[anz->l].ncomps; i++)
l                1027 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] > lcase[anz->l].max[i])
l                1029 src/readfrd.c              lcase[anz->l].max[i]=lcase[anz->l].dat[i][nodenr];
l                1030 src/readfrd.c              lcase[anz->l].nmax[i]=nodenr;
l                1032 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] < lcase[anz->l].min[i])
l                1034 src/readfrd.c              lcase[anz->l].min[i]=lcase[anz->l].dat[i][nodenr];
l                1035 src/readfrd.c              lcase[anz->l].nmin[i]=nodenr;
l                1043 src/readfrd.c        anz->l++;
l                1044 src/readfrd.c        printf ("reading Dataset No:%" ITGFORMAT "\n",anz->l+1);
l                1045 src/readfrd.c        if ( (lcase = (Datasets *)realloc((Datasets *)lcase, (anz->l+2) * sizeof(Datasets))) == NULL )
l                1048 src/readfrd.c        lcase[anz->l].value=stof( rec_str, 13, 25 );
l                1049 src/readfrd.c        strcpy(lcase[anz->l].name,"STRESS  ");
l                1050 src/readfrd.c        lcase[anz->l].ncomps = 6;
l                1051 src/readfrd.c        lcase[anz->l].irtype = 1;
l                1052 src/readfrd.c        lcase[anz->l].loaded = 1;
l                1053 src/readfrd.c        lcase[anz->l].format_flag=format_flag;
l                1055 src/readfrd.c        if ( (lcase[anz->l].nmax = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1057 src/readfrd.c        if ( (lcase[anz->l].nmin = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1059 src/readfrd.c        if ( (lcase[anz->l].max = (float *)malloc( lcase[anz->l].ncomps * sizeof(float))) == NULL )
l                1061 src/readfrd.c        if ( (lcase[anz->l].min = (float *)malloc( lcase[anz->l].ncomps * sizeof(float))) == NULL )
l                1063 src/readfrd.c        if ( (lcase[anz->l].dat = (float **)malloc( lcase[anz->l].ncomps * sizeof(float *))) == NULL )
l                1065 src/readfrd.c        if ( (lcase[anz->l].compName = (char **)malloc( lcase[anz->l].ncomps * sizeof(char *))) == NULL )
l                1067 src/readfrd.c        if ( (lcase[anz->l].icname = (char **)malloc( lcase[anz->l].ncomps * sizeof(char *))) == NULL )
l                1069 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
l                1071 src/readfrd.c          if ( (lcase[anz->l].dat[i] = (float *)malloc( (anz->nmax+1) * sizeof(float))) == NULL )
l                1073 src/readfrd.c          if ( (lcase[anz->l].compName[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
l                1075 src/readfrd.c          if ( (lcase[anz->l].icname[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
l                1077 src/readfrd.c          lcase[anz->l].max[i]=-MAX_FLOAT;
l                1078 src/readfrd.c          lcase[anz->l].min[i]=MAX_FLOAT;
l                1079 src/readfrd.c          for(j=0; j<=anz->nmax; j++) lcase[anz->l].dat[i][j]=0.;
l                1081 src/readfrd.c        if ( (lcase[anz->l].menu = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1083 src/readfrd.c        if ( (lcase[anz->l].ictype = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1085 src/readfrd.c        if ( (lcase[anz->l].icind1 = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1087 src/readfrd.c        if ( (lcase[anz->l].icind2 = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1089 src/readfrd.c        if ( (lcase[anz->l].iexist = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1092 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
l                1094 src/readfrd.c          lcase[anz->l].menu[i] = 1;
l                1095 src/readfrd.c          lcase[anz->l].ictype[i] = 4;
l                1096 src/readfrd.c          lcase[anz->l].iexist[i] = 0;
l                1098 src/readfrd.c        lcase[anz->l].icind1[0] = 1;
l                1099 src/readfrd.c        lcase[anz->l].icind2[0] = 1;
l                1100 src/readfrd.c        lcase[anz->l].icind1[1] = 2;
l                1101 src/readfrd.c        lcase[anz->l].icind2[1] = 2;
l                1102 src/readfrd.c        lcase[anz->l].icind1[2] = 3;
l                1103 src/readfrd.c        lcase[anz->l].icind2[2] = 3;
l                1104 src/readfrd.c        lcase[anz->l].icind1[3] = 1;
l                1105 src/readfrd.c        lcase[anz->l].icind2[3] = 2;
l                1106 src/readfrd.c        lcase[anz->l].icind1[4] = 2;
l                1107 src/readfrd.c        lcase[anz->l].icind2[4] = 3;
l                1108 src/readfrd.c        lcase[anz->l].icind1[5] = 3;
l                1109 src/readfrd.c        lcase[anz->l].icind2[5] = 1;
l                1111 src/readfrd.c        strcpy(lcase[anz->l].compName[0], "xx      ");
l                1112 src/readfrd.c        strcpy(lcase[anz->l].compName[1], "yy      ");
l                1113 src/readfrd.c        strcpy(lcase[anz->l].compName[2], "zz      ");
l                1114 src/readfrd.c        strcpy(lcase[anz->l].compName[3], "xy      ");
l                1115 src/readfrd.c        strcpy(lcase[anz->l].compName[4], "yz      ");
l                1116 src/readfrd.c        strcpy(lcase[anz->l].compName[5], "zx      ");
l                1138 src/readfrd.c                for(i=0; i<lcase[anz->l].ncomps; i++)
l                1139 src/readfrd.c                  lcase[anz->l].dat[i][nodenr]= stof(&rec_str[8+i*12], 1, 12);
l                1143 src/readfrd.c                for(i=0; i<lcase[anz->l].ncomps; i++)
l                1144 src/readfrd.c                  lcase[anz->l].dat[i][nodenr]= stof(&rec_str[13+i*12], 1, 12);
l                1152 src/readfrd.c          for(i=0; i<lcase[anz->l].ncomps; i++)
l                1154 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] > lcase[anz->l].max[i])
l                1156 src/readfrd.c              lcase[anz->l].max[i]=lcase[anz->l].dat[i][nodenr];
l                1157 src/readfrd.c              lcase[anz->l].nmax[i]=nodenr;
l                1159 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] < lcase[anz->l].min[i])
l                1161 src/readfrd.c              lcase[anz->l].min[i]=lcase[anz->l].dat[i][nodenr];
l                1162 src/readfrd.c              lcase[anz->l].nmin[i]=nodenr;
l                1170 src/readfrd.c        anz->l++;
l                1171 src/readfrd.c        printf ("reading Dataset No:%" ITGFORMAT "\n",anz->l+1);
l                1172 src/readfrd.c        if ( (lcase = (Datasets *)realloc((Datasets *)lcase, (anz->l+2) * sizeof(Datasets))) == NULL )
l                1175 src/readfrd.c        lcase[anz->l].value=stof( rec_str, 13, 25 );
l                1176 src/readfrd.c        strcpy(lcase[anz->l].name,"TEMP    ");
l                1177 src/readfrd.c        lcase[anz->l].ncomps = 1;
l                1178 src/readfrd.c        lcase[anz->l].irtype = 1;
l                1179 src/readfrd.c        lcase[anz->l].loaded = 1;
l                1180 src/readfrd.c        lcase[anz->l].format_flag=format_flag;
l                1182 src/readfrd.c        if ( (lcase[anz->l].nmax = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1184 src/readfrd.c        if ( (lcase[anz->l].nmin = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1186 src/readfrd.c        if ( (lcase[anz->l].max = (float *)malloc( lcase[anz->l].ncomps * sizeof(float))) == NULL )
l                1188 src/readfrd.c        if ( (lcase[anz->l].min = (float *)malloc( lcase[anz->l].ncomps * sizeof(float))) == NULL )
l                1190 src/readfrd.c        if ( (lcase[anz->l].dat = (float **)malloc( lcase[anz->l].ncomps * sizeof(float *))) == NULL )
l                1192 src/readfrd.c        if ( (lcase[anz->l].compName = (char **)malloc( lcase[anz->l].ncomps * sizeof(char *))) == NULL )
l                1194 src/readfrd.c        if ( (lcase[anz->l].icname = (char **)malloc( lcase[anz->l].ncomps * sizeof(char *))) == NULL )
l                1196 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
l                1198 src/readfrd.c          if ( (lcase[anz->l].dat[i] = (float *)malloc( (anz->nmax+1) * sizeof(float))) == NULL )
l                1200 src/readfrd.c          if ( (lcase[anz->l].compName[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
l                1202 src/readfrd.c          if ( (lcase[anz->l].icname[i] = (char *)malloc( MAX_LINE_LENGTH * sizeof(char))) == NULL )
l                1204 src/readfrd.c          lcase[anz->l].max[i]=-MAX_FLOAT;
l                1205 src/readfrd.c          lcase[anz->l].min[i]=MAX_FLOAT;
l                1206 src/readfrd.c          for(j=0; j<=anz->nmax; j++) lcase[anz->l].dat[i][j]=0.;
l                1208 src/readfrd.c        if ( (lcase[anz->l].menu = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1210 src/readfrd.c        if ( (lcase[anz->l].ictype = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1212 src/readfrd.c        if ( (lcase[anz->l].icind1 = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1214 src/readfrd.c        if ( (lcase[anz->l].icind2 = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1216 src/readfrd.c        if ( (lcase[anz->l].iexist = (ITG *)malloc( lcase[anz->l].ncomps * sizeof(ITG))) == NULL )
l                1219 src/readfrd.c        for(i=0; i<lcase[anz->l].ncomps; i++)
l                1221 src/readfrd.c          lcase[anz->l].menu[i] = 1;
l                1222 src/readfrd.c          lcase[anz->l].ictype[i] = 1;
l                1223 src/readfrd.c          lcase[anz->l].icind1[i] = i+1;
l                1224 src/readfrd.c          lcase[anz->l].icind2[i] = 0;
l                1225 src/readfrd.c          lcase[anz->l].iexist[i] = 0;
l                1228 src/readfrd.c        strcpy(lcase[anz->l].compName[0], "Value   ");
l                1245 src/readfrd.c              for(i=0; i<lcase[anz->l].ncomps; i++)
l                1246 src/readfrd.c                lcase[anz->l].dat[i][nodenr]= stof(&rec_str[8+i*12], 1, 12);
l                1250 src/readfrd.c              for(i=0; i<lcase[anz->l].ncomps; i++)
l                1251 src/readfrd.c                lcase[anz->l].dat[i][nodenr]= stof(&rec_str[13+i*12], 1, 12);
l                1258 src/readfrd.c          for(i=0; i<lcase[anz->l].ncomps; i++)
l                1260 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] > lcase[anz->l].max[i])
l                1262 src/readfrd.c              lcase[anz->l].max[i]=lcase[anz->l].dat[i][nodenr];
l                1263 src/readfrd.c              lcase[anz->l].nmax[i]=nodenr;
l                1265 src/readfrd.c            if(lcase[anz->l].dat[i][nodenr] < lcase[anz->l].min[i])
l                1267 src/readfrd.c              lcase[anz->l].min[i]=lcase[anz->l].dat[i][nodenr];
l                1268 src/readfrd.c              lcase[anz->l].nmin[i]=nodenr;
l                1295 src/readfrd.c    anz->l++;
l                1313 src/readfrd.c    anz->olc = anz->l;
l                  28 src/readfrd.h    ITG   l;         /* number of loadcases (Datasets) */
l                  34 src/sgi.c        ITG i,j,k,l,*pointers=NULL,method;
l                  50 src/sgi.c        l=*nzs;
l                  56 src/sgi.c      	irowsgi[--k]=irow[--l]-1;
l                  57 src/sgi.c      	ausgi[k]=au[l];
l                  68 src/sgi.c      	irowsgi[--k]=irow[--l]-1;
l                  69 src/sgi.c      	ausgi[k]=au[l]-*sigma*aub[l];
l                  76 src/steadystate.c     iinc=0,l,iout,ielas,icmd,iprescribedboundary,ndata,nmd,nevd,
l                 956 src/steadystate.c       for(l=0;l<ndatatot;l=l+*jout){
l                 958 src/steadystate.c 	  time=freq[l]/(2.*pi);
l                1317 src/steadystate.c 		  br[i]+=freq[l]*(freq[l]*mubr[i]+alpham*mubi[i]+betam*fi[i]);
l                1318 src/steadystate.c 		  bi[i]+=freq[l]*(freq[l]*mubi[i]-alpham*mubr[i]-betam*fr[i]);
l                1344 src/steadystate.c 		      dd=pow(d[i]-pow(freq[l],2),2)+
l                1345 src/steadystate.c 			  pow(fric[i],2)*pow(freq[l],2);
l                1346 src/steadystate.c 		      bjr[i]=(aa[i]*(d[i]-freq[l]*freq[l])+
l                1347 src/steadystate.c 			      bb[i]*fric[i]*freq[l])/dd;
l                1348 src/steadystate.c 		      bji[i]=(bb[i]*(d[i]-freq[l]*freq[l])-
l                1349 src/steadystate.c 			      aa[i]*fric[i]*freq[l])/dd;
l                1365 src/steadystate.c 		      am[i*nev2+i]=d[i]-freq[l]*freq[l];
l                1366 src/steadystate.c 		      am[(i+nev)*nev2+i]=-fric[i]*freq[l];
l                1373 src/steadystate.c 			      -cc[i*nev+j]*freq[l];
l                1375 src/steadystate.c 			      +cc[i*nev+j]*freq[l];
l                1406 src/steadystate.c 		  ax[i]=-pow(freq[l],2)-freq[l]*fric[2*i+1]+d[2*i];
l                1407 src/steadystate.c 		  bx[i]=-freq[l]*fric[2*i]-d[2*i+1];
l                2021 src/steadystate.c       for(l=0;l<nfour;l++){
l                2023 src/steadystate.c 	  time=tmin+(tmax-tmin)*(double)l/(double)nfour;
l                2025 src/steadystate.c 	  FORTRAN(tempload,(xforcold,xforc,&xforcact[l**nforc],iamforc,nforc,
l                2026 src/steadystate.c 	    xloadold,xload,&xloadact[l*2**nload],iamload,nload,ibody,xbody,
l                2027 src/steadystate.c 	    nbody,xbodyold,&xbodyact[l*7**nbody],t1old,t1,t1act,
l                2029 src/steadystate.c             ithermal,nmethod,xbounold,xboun,&xbounact[l**nboun],iamboun,nboun,
l                2039 src/steadystate.c       for(i=0;i<l**nboun;i++){xbounacttime[i]=xbounact[i];}
l                2045 src/steadystate.c 	  for(l=0;l<nfour;l++){
l                2046 src/steadystate.c 	      r[l]=xforcact[l**nforc+i];
l                2050 src/steadystate.c 	  for(l=0;l<nfour;l++){
l                2051 src/steadystate.c 	      xforcact[l**nforc+i]=r[l]/nfour*2.;
l                2057 src/steadystate.c 	  for(l=0;l<nfour;l++){
l                2058 src/steadystate.c 	      r[l]=xloadact[l*2**nload+2*i];
l                2062 src/steadystate.c 	  for(l=0;l<nfour;l++){
l                2063 src/steadystate.c 	      xloadact[l*2**nload+2*i]=r[l]/nfour*2.;
l                2069 src/steadystate.c 	  for(l=0;l<nfour;l++){
l                2070 src/steadystate.c 	      r[l]=xbodyact[l**nbody+7*i];
l                2074 src/steadystate.c 	  for(l=0;l<nfour;l++){
l                2075 src/steadystate.c 	      xbodyact[l**nbody+7*i]=r[l]/nfour*2.;
l                2082 src/steadystate.c 	      for(l=0;l<nfour;l++){
l                2083 src/steadystate.c 		  r[l]=xbounact[l**nboun+i];
l                2087 src/steadystate.c 	      for(l=0;l<nfour;l++){
l                2088 src/steadystate.c 		  xbounact[l**nboun+i]=r[l]/nfour*2.;
l                2232 src/steadystate.c 	  for(l=0;l<nfour;l++){
l                2236 src/steadystate.c 	      freq[l]=freqnh[ii]*floor((l+1.)/2.+0.1);
l                2254 src/steadystate.c 	            ttime,&freq[l],istep,&iinc,ibody,clearini,&mortar,springarea,
l                2290 src/steadystate.c 		  xforcr[i]=xforcact[l**nforc+i];
l                2294 src/steadystate.c 		  xloadr[2*i]=xloadact[l*2**nload+2*i];
l                2298 src/steadystate.c 		  xbodyr[7*i]=xbodyact[l**nbody+7*i];
l                2364 src/steadystate.c 		      xbounr[i]=xbounact[l**nboun+i];
l                2399 src/steadystate.c 		      br[i]+=freq[l]*(freq[l]*mubr[i]);
l                2400 src/steadystate.c 		      bi[i]+=freq[l]*(-alpham*mubr[i]-betam*fr[i]);
l                2424 src/steadystate.c 		      dd=pow(d[i]-pow(freq[l],2),2)+
l                2425 src/steadystate.c 			  pow(fric[i],2)*pow(freq[l],2);
l                2426 src/steadystate.c 		      bjr[i]=(aa[i]*(d[i]-freq[l]*freq[l])+
l                2427 src/steadystate.c 			      bb[i]*fric[i]*freq[l])/dd;
l                2428 src/steadystate.c 		      bji[i]=(bb[i]*(d[i]-freq[l]*freq[l])-
l                2429 src/steadystate.c 			      aa[i]*fric[i]*freq[l])/dd;
l                2444 src/steadystate.c 		      am[i*nev2+i]=d[i]-freq[l]*freq[l];
l                2445 src/steadystate.c 		      am[(i+nev)*nev2+i]=-fric[i]*freq[l];
l                2452 src/steadystate.c 			      -cc[i*nev+j]*freq[l];
l                2454 src/steadystate.c 			      +cc[i*nev+j]*freq[l];
l                2581 src/steadystate.c 		  if((l!=0)&&(2*(ITG)floor(l/2.+0.1)==l)){
l                2592 src/steadystate.c 			  btot[j*neq[1]+i]+=br[i]*cos(freq[l]*time+bi[i]);
l                2593 src/steadystate.c 			  bp[j*neq[1]+i]-=freq[l]*br[i]*sin(freq[l]*time+bi[i]);
l                2616 src/steadystate.c 		  if((l!=0)&&(2*(ITG)floor(l/2.+0.1)==l)){
l                2629 src/steadystate.c 			  btot[j*neq[1]+i]+=br[i]*cos(freq[l]*time+bi[i]);
l                2630 src/steadystate.c 			  bp[j*neq[1]+i]-=freq[l]*br[i]*sin(freq[l]*time+bi[i]);
l                2682 src/steadystate.c 	  for(l=0;l<nfour;l++){
l                2683 src/steadystate.c 	      time=tmin+2.*pi/freqnh[ii]*(double)l/(double)nfour;
l                2701 src/steadystate.c 		    ndirboun,&xbounacttime[l**nboun],nboun,ipompc,nodempc,
l                2702 src/steadystate.c                     coefmpc,labmpc,nmpc,&btot[l*neq[1]],&bp[l*neq[1]],veold,&dtime,mi,
l                2711 src/steadystate.c 		      vold,&btot[l*neq[1]],nodeboun,ndirboun,
l                2712 src/steadystate.c                       &xbounacttime[l**nboun],nboun,
l                  41 src/storecontactdof.c   ITG j,j1,jdof,id,k,l,ist,index,node,ndir,*ikactcont=*ikactcontp;
l                  62 src/storecontactdof.c 		  l=k-1;
l                  64 src/storecontactdof.c 		      ikactcont[k--]=ikactcont[l--];
l                 100 src/storecontactdof.c 				  l=k-1;
l                 102 src/storecontactdof.c 				      ikactcont[k--]=ikactcont[l--];
l                  45 src/tau.c        ITG i,j,k,l,*irow=NULL;
l                  66 src/tau.c        l=*nzs;
l                  72 src/tau.c      	irowtau[--k]=irow[--l]-1;
l                  73 src/tau.c      	autau[k]=au[l];
l                  84 src/tau.c      	irowtau[--k]=irow[--l]-1;
l                  85 src/tau.c      	autau[k]=au[l]-*sigma*aub[l];