谁帮我看看这个算法错在哪了，谢谢

shaohuan001

1. #include <math.h>
   
2. #include <iostream.h>
   
3. #include<fstream.h>
   
4. #include"iomanip.h"
   
5. 
   
6. int loc[200][4],e;
   
7. float a,b,g[12][12],h,anu;
   
8. MATD()
   
9. { int i,j,l,s,p,m;
   
10. float afh1[4],afh2[4],fh1,fh2,d,n,k1[3][3],ab;
    
11. 
    
12. afh1[0]=-1;afh1[1]=1;afh1[2]=1;afh1[3]=-1;
    
13. afh2[0]=-1;afh2[1]=-1;afh2[2]=1;afh2[3]=1;
    
14. d=e*h*h*h/(12*(1-anu*anu));
    
15. n=d/(60*a*b);
    
16. ab=b*b/(a*a);
    
17. 
    
18. for(i=0;i<4;i++)
    
19. for(j=0;j<4;j++)
    
20. {fh1=afh1[i]*afh1[j];
    
21. fh2=afh2[i]*afh2[j];
    
22. k1[0][0]=3*n*(15.0*(ab*fh1+fh2/ab)+(14-4*anu+5*ab+5.0/ab)*fh1*fh2);
    
23. k1[0][1]=-3*n*b*((2+3*anu+5.0/ab)*fh1*afh2[i]+15*afh2[i]/ab+5*anu*fh1*afh2[j]);
    
24. k1[0][2]=3*n*a*((2+3*anu+5*ab)*fh2*afh1[i]+15*afh1[i]*ab+5.0*anu*fh2*afh1[j]);
    
25. k1[1][0]=-3*n*b*((2+3*anu+5.0/ab)*fh1*afh2[j]+15*afh2[j]/ab+5*anu*fh1*afh2[i]);
    
26. k1[1][1]=n*b*b*(2*(1-anu)*fh1*(3+5*fh2)+5*(3+fh1)*(3+fh2)/ab);
    
27. k1[1][2]=-15*n*anu*a*b*(afh1[i]+afh1[j])*(afh2[i]+afh2[j]);
    
28. k1[2][0]=3*n*a*((2+3*anu+5*ab)*fh2*afh1[j]+15*afh1[j]*ab+5*anu*fh2*afh1[i]);
    
29. k1[2][1]=-15*n*anu*a*b*(afh1[i]+afh1[j])*(afh2[i]+afh2[j]);
    
30. k1[2][2]=n*a*a*(2*(1-anu)*fh2*(3+5*fh1)+5*(3+fh1)*(3+fh2)*ab);
    
31. 
    
32. for(m=0;m<3;m++)
    
33. for(l=0;l<3;l++)
    
34. {s=3*i+m;
    
35. p=3*j+l;
    
36. g[s][p]=k1[m][l]; } }
    
37. 
    
38. }
    
39. 
    
40. void main()
    
41. {  int i,j,l,v,maxnel,inode,nodei,idofn,nrows,nrowe,jnode,nodej,jdofn,
    
42. ncols,ncole,nel,maxnod,a1,b1;
    
43. float t[12][12],k2[12][12],sum[12][12],astif[100][100],ast[100][100],k[100][100]={0};
    
44. 
    
45. ifstream ifile("e:\\shaohuan001.txt");
    
46. ifile>>e>>h>>anu>>a>>b>>v>>maxnel>>maxnod;
    
47. 
    
48. for(i=0;i<12;i++)
    
49. for(j=0;j<12;j++)
    
50.   t[i][j]=0;
    
51. 
    
52. for(j=0;j<maxnel;j++)
    
53. {  
    
54. ifile>>loc[j][0]>>loc[j][1]>>loc[j][2]>>loc[j][3];
    
55. loc[j][0]--;loc[j][1]--;loc[j][2]--;loc[j][3]--;
    
56. 
    
57. if(loc[j][0]%6==0)
    
58. i=(loc[j][0]-1)%6;
    
59. else
    
60. i=loc[j][0]%6-1;
    
61. 
    
62. t[0][0]=cos(i*v);t[0][1]=-sin(i*v);t[1][0]=sin(i*v);t[1][1]=cos(i*v);  t[2][2]=1;
    
63.   t[3][3]=cos(i*v); t[3][4]=-sin(i*v);t[4][3]=sin(i*v);t[4][4]=cos(i*v); t[5][5]=1;
    
64.   t[6][6]=cos(i*v);t[6][7]=-sin(i*v); t[7][6]=sin(i*v);t[7][7]=cos(i*v);t[8][8]=1;
    
65.   t[9][9]=cos(i*v);t[9][10]=-sin(i*v);t[10][9]=sin(i*v);t[10][10]=cos(i*v);t[11][11]=1;
    
66. 
    
67.   MATD();
    
68. for(a1=0;a1<12;a1++)
    
69. for(b1=0;b1<12;b1++)
    
70. { sum[a1][b1]=0.0;
    
71. for(l=0;l<12;l++)
    
72. sum[a1][b1]=sum[a1][b1]+g[a1][l]*t[l][b1]; }
    
73. for(a1=0;a1<12;a1++)
    
74. for(b1=0;b1<12;b1++)
    
75. { k2[a1][b1]=0.0;
    
76. for(l=0;l<12;l++)
    
77.   k2[a1][b1]=k2[a1][b1]+t[l][a1]*sum[l][b1];
    
78. }
    
79. 
    
80. for(a1=0;a1<3*maxnod;a1++)
    
81. for(b1=0;b1<3*maxnod;b1++)
    
82.   astif[a1][b1]=0;
    
83. for(nel=0;nel<maxnel;nel++)
    
84. for(inode=0;inode<4;inode++)
    
85. {  nodei=loc[nel][inode];
    
86. for(idofn=0;idofn<3;idofn++)
    
87. {  nrows=nodei*3+idofn;
    
88.    nrowe=inode*3+idofn;
    
89. for(jnode=0;jnode<4;jnode++)
    
90. {  nodej=loc[nel][jnode];
    
91. for(jdofn=0;jdofn<3;jdofn++)
    
92. {   ncols=nodej*3+jdofn;
    
93.     ncole=jnode*3+jdofn;
    
94.     astif[nrows][ncols]=astif[nrows][ncols]+k2[nrowe][ncole];
    
95.     }    }    }    }
    
96. 
    
97. for(a1=0;a1<3*maxnod;a1++)
    
98. for(b1=0;b1<3*maxnod;b1++)
    
99. 
    
100. k[a1][b1]=k[a1][b1]+astif[a1][b1];}
     
101. 
     
102. 
     
103. for(a1=0;a1<3*maxnod;a1++)
     
104. {for(b1=0;b1<3*maxnod;b1++)
     
105. 
     
106. cout<<k[a1][b1]<<"\t";}
     
107. cout<<endl;
     
108. 
     
109. }

复制代码

[ 本帖最后由 风花雪月 于 2007-6-15 11:01 编辑 ]

风花雪月

是不能编译还是别的什么错误？请说明

cooller

自己写得程序，一句注释都没有，别人怎么看？建议加上注释。

yf2002043227

错误这么多?很多基本的错误,你是不是想让大家帮你调试啊,自己偷懒?:@@

appleseed05

没注释，看不懂啊