1. C++ / Говнокод #1272

    +76.3

    1. 001
    2. 002
    3. 003
    4. 004
    5. 005
    6. 006
    7. 007
    8. 008
    9. 009
    10. 010
    11. 011
    12. 012
    13. 013
    14. 014
    15. 015
    16. 016
    17. 017
    18. 018
    19. 019
    20. 020
    21. 021
    22. 022
    23. 023
    24. 024
    25. 025
    26. 026
    27. 027
    28. 028
    29. 029
    30. 030
    31. 031
    32. 032
    33. 033
    34. 034
    35. 035
    36. 036
    37. 037
    38. 038
    39. 039
    40. 040
    41. 041
    42. 042
    43. 043
    44. 044
    45. 045
    46. 046
    47. 047
    48. 048
    49. 049
    50. 050
    51. 051
    52. 052
    53. 053
    54. 054
    55. 055
    56. 056
    57. 057
    58. 058
    59. 059
    60. 060
    61. 061
    62. 062
    63. 063
    64. 064
    65. 065
    66. 066
    67. 067
    68. 068
    69. 069
    70. 070
    71. 071
    72. 072
    73. 073
    74. 074
    75. 075
    76. 076
    77. 077
    78. 078
    79. 079
    80. 080
    81. 081
    82. 082
    83. 083
    84. 084
    85. 085
    86. 086
    87. 087
    88. 088
    89. 089
    90. 090
    91. 091
    92. 092
    93. 093
    94. 094
    95. 095
    96. 096
    97. 097
    98. 098
    99. 099
    100. 100
    // Begin Init Filter
void initFilter(){
   int i,j;
   double Pgam,Ptet,Ppsi,
   sin_psi,cos_psi,sin_gam,cos_gam,cos_tet,
   Pl,Pb,uu[81],bb[81];  //  

     hfo=1.; Nsf=1;
     sugl[0]=1.5*M_PI/180.; sugl[1]=0.1*M_PI/180;sugl[2]=sugl[1];
     dbc=100.; dbc*=0.157236869549e-7; dlc=dbc; dac=sugl[0];
     Rk[0]=Omz/100000000.*hfo; Rk[1]=Omz/100000000.*hfo; Rk[2]=Rk[0];
     Rk[3]=Omz/100000000.*hfo; Rk[4]=Rk[3]; Rk[5]=Rk[3];
     for(i=0;i<3;i++) {
     DTet[i]=0.;
     Xs[i]=0.;
    }
    for(i=0;i<3;i++)  {
    som[i]=0.1*M_PI/(180.0*3600.0);          // som[i]=3.*M_PI/(180.0*3600.0);
   }
    som[0]=0.1*M_PI/(180.0*3600.0);          //som[0]=2.*M_PI/(180.0*3600.0);
    Tauom[0]=360.;                         //  Tauom[0]=300.;
    Tauom[1]=360.;                         //  Tauom[1]=100.;
    Tauom[2]=360.;                         //  Tauom[2]=300.;
 for(i=0;i<9;i++)  {
  for(j=0;j<9;j++) {
   Fko[i][j]=0.;
   Fo[i][j]=0.;
   Pko[i][j]=0.;
  }
  Xko[i]=0.;
  Fko[i][i]=1.;
  Fo[i][i]=1.;
 }

 for(i=0;i<6;i++)  {
  for(j=0;j<9;j++)
   Hizo[i][j]=0.;
 }

 for(i=0;i<9;i++)  {
  for(j=0;j<3;j++)
   Gko[i][j]=0.;
 }
//begin SKO for MNK
  Pgam=sugl[2]*sugl[2]; Ptet=sugl[1]*sugl[1]; Ppsi=sugl[0]*sugl[0];
  Pl=dlc*dlc; Pb=dbc*dbc;
  sin_psi=sin(Angl[0]); cos_psi=cos(Angl[0]);
  sin_gam=sin(Angl[2]); cos_gam=cos(Angl[2]);
  cos_tet=cos(Angl[1]);


 uu[0]=Ur[0][0]*Ur[0][0];           bb[0]=B3[0][0]*B3[0][0];
 uu[1]=Ur[0][0]*Ur[0][1];           bb[1]=B3[0][0]*B3[0][1];
 uu[2]=Ur[0][0]*Ur[0][2];           bb[2]=B3[0][0]*B3[0][2];
 uu[3]=Ur[0][0]*Ur[1][0];           bb[3]=B3[0][0]*B3[1][0];
 uu[4]=Ur[0][0]*Ur[1][1];           bb[4]=B3[0][0]*B3[1][1];
 uu[5]=Ur[0][0]*Ur[1][2];           bb[5]=B3[0][0]*B3[1][2];
 uu[6]=Ur[0][0]*Ur[2][0];           bb[6]=B3[0][0]*B3[2][0];
 uu[7]=Ur[0][0]*Ur[2][1];           bb[7]=B3[0][0]*B3[2][1];
 uu[8]=Ur[0][0]*Ur[2][2];           bb[8]=B3[0][0]*B3[2][2];

 uu[9]=Ur[0][1]*Ur[0][1];           bb[9]=B3[0][1]*B3[0][1];
 uu[10]=Ur[0][1]*Ur[0][2];          bb[10]=B3[0][1]*B3[0][2];
 uu[11]=Ur[0][1]*Ur[1][0];          bb[11]=B3[0][1]*B3[1][0];
 uu[12]=Ur[0][1]*Ur[1][1];          bb[12]=B3[0][1]*B3[1][1];
 uu[13]=Ur[0][1]*Ur[1][2];          bb[13]=B3[0][1]*B3[1][2];
 uu[14]=Ur[0][1]*Ur[2][0];          bb[14]=B3[0][1]*B3[2][0];
 uu[15]=Ur[0][1]*Ur[2][1];          bb[15]=B3[0][1]*B3[2][1];
 uu[16]=Ur[0][1]*Ur[2][2];          bb[16]=B3[0][1]*B3[2][2];

 uu[17]=Ur[0][2]*Ur[0][2];          bb[17]=B3[0][2]*B3[0][2];
 uu[18]=Ur[0][2]*Ur[1][0];          bb[18]=B3[0][2]*B3[1][0];
 uu[19]=Ur[0][2]*Ur[1][1];          bb[19]=B3[0][2]*B3[1][1];
 uu[20]=Ur[0][2]*Ur[1][2];          bb[20]=B3[0][2]*B3[1][2];
 uu[21]=Ur[0][2]*Ur[2][0];          bb[21]=B3[0][2]*B3[2][0];
 uu[22]=Ur[0][2]*Ur[2][1];          bb[22]=B3[0][2]*B3[2][1];
 uu[23]=Ur[0][2]*Ur[2][2];          bb[23]=B3[0][2]*B3[2][2];

 uu[24]=Ur[1][0]*Ur[1][0];          bb[24]=B3[1][0]*B3[1][0];
 uu[25]=Ur[1][0]*Ur[1][1];          bb[25]=B3[1][0]*B3[1][1];
 uu[26]=Ur[1][0]*Ur[1][2];          bb[26]=B3[1][0]*B3[1][2];
 uu[27]=Ur[1][0]*Ur[2][0];          bb[27]=B3[1][0]*B3[2][0];
 uu[28]=Ur[1][0]*Ur[2][1];          bb[28]=B3[1][0]*B3[2][1];
 uu[29]=Ur[1][0]*Ur[2][2];          bb[29]=B3[1][0]*B3[2][2];

 uu[30]=Ur[1][1]*Ur[1][1];          bb[30]=B3[1][1]*B3[1][1];
 uu[31]=Ur[1][1]*Ur[1][2];          bb[31]=B3[1][1]*B3[1][2];
 uu[32]=Ur[1][1]*Ur[2][0];          bb[32]=B3[1][1]*B3[2][0];
 uu[33]=Ur[1][1]*Ur[2][1];          bb[33]=B3[1][1]*B3[2][1];
 uu[34]=Ur[1][1]*Ur[2][2];          bb[34]=B3[1][1]*B3[2][2];

 uu[35]=Ur[1][2]*Ur[1][2];          bb[35]=B3[1][2]*B3[1][2];
 uu[36]=Ur[1][2]*Ur[2][0];          bb[36]=B3[1][2]*B3[2][0];
 uu[37]=Ur[1][2]*Ur[2][1];          bb[37]=B3[1][2]*B3[2][1];
 uu[38]=Ur[1][2]*Ur[2][2];          bb[38]=B3[1][2]*B3[2][2];

 uu[39]=Ur[2][0]*Ur[2][0];          bb[39]=B3[2][0]*B3[2][0];
 uu[40]=Ur[2][0]*Ur[2][1];          bb[40]=B3[2][0]*B3[2][1];
 uu[41]=Ur[2][0]*Ur[2][2];          bb[41]=B3[2][0]*B3[2][2];
 и так далее

    Не спрашивайте где это использовалось, спать лучше будете)))

    Запостил: guest, 28 Июня 2009

    Комментарии (11) RSS

    • ava guest 28.06.2009 01:06 # +1

      [quote=Говногость]Не спрашивайте где это использовалось, спать лучше будете)))[/quote]Да, да... Помню, помню. Кажется это использовалось в автоматической системе управления Московским метро.
      Ответить

    • ava guest 28.06.2009 12:40 # −2

      Валерий:
      Как бы расчёт физических свойств поверхности текстуры и перспективы фигур. Короче в любимой игрушке-стрелялки? :)
      Ответить

    • ava guest 28.06.2009 21:47 # 0

      запостивший:
      uu[42]=Ur[2][1]*Ur[2][1]; bb[42]=B3[2][1]*B3[2][1];
      uu[43]=Ur[2][1]*Ur[2][2]; bb[43]=B3[2][1]*B3[2][2];

      uu[44]=Ur[2][2]*Ur[2][2]; bb[44]=B3[2][2]*B3[2][2];

      Pur[0]=Ur[0][1]*Ur[0][1]*Pl; Pur[11]=Ur[0][0]*Ur[0][0]*Pl;
      Pur[1]=-Ur[0][1]*Ur[0][0]*Pl; Pur[12]=0;
      Pur[2]=0; Pur[13]=-Ur[0][0]*Ur[1][1]*Pl;
      Pur[3]=Ur[0][1]*Ur[1][1]*Pl; Pur[14]=Ur[0][0]*Ur[1][0]*Pl;
      Pur[4]=-Ur[0][1]*Ur[1][0]*Pl; Pur[16]=-Ur[0][0]*Ur[2][1]*Pl;
      Pur[5]=0; Pur[17]=Ur[0][0]*Ur[2][0]*Pl;
      Pur[6]=Ur[0][1]*Ur[2][1]*Pl; Pur[18]=0;
      Pur[7]=-Ur[0][1]*Ur[2][0]*Pl;
      Pur[8]=0;
      /* Pur[2][2]=Ur[0][2]*Ur[0][2]...=Pur[2][8]=...=Pur[8][2]=0 */

      Pur[33]=Ur[2][0]*Ur[2][0]*Pb+Ur[1][1]*Ur[1][1]*Pl;
      Pur[34]=Ur[2][0]*Ur[2][1]*Pb-Ur[1][1]*Ur[1][0]*Pl;
      Pur[35]=Ur[2][0]*Ur[2][2]*Pb;
      Pur[36]=Ur[1][1]*Ur[2][1]*Pl-Ur[2][0]*Ur[1][0]*Pb;
      Pur[37]=-Ur[2][0]*Ur[1][1]*Pb-Ur[1][1]*Ur[2][0]*Pl;
      Pur[38]=-Ur[2][0]*Ur[1][2]*Pb;

      Pur[44]=Ur[2][1]*Ur[2][1]*Pb+Ur[1][0]*Ur[1][0]*Pl;
      Pur[45]=Ur[2][1]*Ur[2][2]*Pb;
      Pur[46]=-Ur[2][1]*Ur[1][0]*Pb-Ur[1][0]*Ur[2][1]*Pl;
      Pur[47]=Ur[1][0]*Ur[2][0]*Pl-Ur[2][1]*Ur[1][1]*Pb;
      Pur[48]=-Ur[2][1]*Ur[1][2]*Pb;

      Pur[55]=Ur[2][2]*Ur[2][2]*Pb;
      Pur[56]=-Ur[2][2]*Ur[1][0]*Pb;
      Pur[57]=-Ur[2][2]*Ur[1][1]*Pb;
      Pur[58]=-Ur[2][2]*Ur[1][2]*Pb;

      Pur[66]=Ur[2][1]*Ur[2][1]*Pl+Ur[1][0]*Ur[1][0]*Pb;
      Pur[67]=Ur[1][0]*Ur[1][1]*Pb-Ur[2][1]*Ur[2][0]*Pl;
      Pur[68]=Ur[1][0]*Ur[1][2]*Pb;

      Pur[77]=Ur[1][1]*Ur[1][1]*Pb+Ur[2][0]*Ur[2][0]*Pl;
      Pur[78]=Ur[1][1]*Ur[1][2]*Pb;

      Pur[88]=Ur[1][2]*Ur[1][2]*Pb;


      Pb3[0]=B3[0][1]*B3[0][1]*Ppsi+B3[0][2]*B3[0][2]*
      sin_psi*sin_psi*Ptet+B3[2][0]*B3[2][0]*Pgam;
      Pb3[1]=-B3[0][1]*B3[0][0]*Ppsi-B3[0][2]*B3[0][2]*
      sin_psi*cos_psi*Ptet+B3[2][0]*B3[2][1]*Pgam;
      Pb3[2]=B3[0][2]*B3[1][2]*sin_psi*sin_gam*Ptet+
      B3[2][0]*B3[2][2]*Pgam;
      Pb3[3]=B3[0][2]*B3[1][2]*sin_psi*sin_psi*P
      Ответить

    • ava guest 28.06.2009 21:48 # 0

      запостивший:
      tet+
      B3[0][1]*B3[1][1]*Ppsi;
      Pb3[4]=-B3[0][1]*B3[1][0]*Ppsi-B3[0][2]*B3[1][2]*
      sin_psi*cos_psi*Ptet;
      Pb3[5]=sin_psi*cos_tet*Ptet;
      Pb3[6]=B3[0][1]*B3[2][1]*Ppsi-B3[2][0]*B3[0][0]*Pgam+
      B3[0][2]*B3[2][2]*sin_psi*sin_psi*Ptet;
      Pb3[7]=-B3[0][1]*B3[2][0]*Ppsi-B3[2][0]*B3[0][1]*Pgam-
      B3[0][2]*B3[2][2]*sin_psi*cos_psi*Ptet;
      Pb3[8]=-B3[0][2]*B3[1][2]*sin_psi*cos_gam*Ptet-
      B3[2][0]*B3[0][2]*Pgam;

      Pb3[11]=B3[0][0]*B3[0][0]*Ppsi+B3[0][2]*B3[0][2]*
      cos_psi*cos_psi*Ptet+B3[2][1]*B3[2][1]*Pgam;
      Pb3[12]=B3[2][1]*B3[2][2]*Pgam-B3[0][2]*B3[1][2]*
      sin_gam*cos_psi*Ptet;
      Pb3[13]=-B3[0][2]*B3[1][2]*cos_psi*sin_psi*Ptet-
      B3[0][0]*B3[1][1]*Ppsi;
      Pb3[14]=B3[0][2]*B3[1][2]*cos_psi*cos_psi*Ptet+
      B3[0][0]*B3[1][0]*Ppsi;
      Pb3[15]=-B3[0][2]*cos_psi*cos_tet*Ptet;

      Pb3[16]=-B3[0][0]*B3[2][1]*Ppsi-B3[2][1]*B3[0][0]*Pgam-
      B3[0][2]*B3[2][2]*sin_psi*cos_psi*Ptet;
      Pb3[17]=B3[0][0]*B3[2][0]*Ppsi-B3[2][1]*B3[0][1]*Pgam+
      B3[0][2]*B3[2][2]*cos_psi*cos_psi*Ptet;
      Pb3[18]=B3[0][2]*B3[1][2]*cos_psi*cos_gam*Ptet-
      B3[2][1]*B3[0][2]*Pgam;

      Pb3[22]=B3[1][2]*B3[1][2]*sin_gam*sin_gam*Ptet+B3[2][2]*B3[2][2]*Pgam;
      Pb3[23]=B3[1][2]*B3[1][2]*sin_psi*sin_gam*Ptet;
      Pb3[24]=-B3[1][2]*B3[1][2]*cos_psi*sin_gam*Ptet;
      Pb3[25]=B3[1][2]*cos_tet*sin_gam*Ptet;
      Pb3[26]=B3[1][2]*B3[2][2]*sin_gam*sin_psi*Ptet-B3[2][2]*B3[0][0]*Pgam;
      Pb3[27]=-B3[1][2]*B3[2][2]*sin_gam*cos_psi*Ptet-B3[2][2]*B3[0][1]*Pgam;
      Pb3[28]=-B3[1][2]*B3[1][2]*sin_gam*cos_gam*Ptet-B3[2][2]*B3[0][2]*Pgam;
      Ответить

    • ava guest 28.06.2009 21:49 # 0

      запостивший:
      Pb3[33]=B3[1][2]*B3[1][2]*sin_psi*sin_psi*Ptet+B3[1][1]*B3[1][1]*Ppsi;
      Pb3[34]=-B3[1][2]*B3[1][2]*sin_psi*cos_psi*Ptet-B3[1][1]*B3[1][0]*Ppsi;
      Pb3[35]=B3[1][2]*sin_psi*cos_tet*Ptet;
      Pb3[36]=B3[1][2]*B3[2][2]*sin_psi*sin_psi*Ptet+B3[1][1]*B3[2][1]*Ppsi;
      Pb3[37]=-B3[1][2]*B3[2][2]*sin_psi*cos_psi*Ptet-B3[1][1]*B3[2][0]*Ppsi;
      Pb3[38]=-B3[1][2]*B3[1][2]*sin_psi*cos_gam*Ptet;

      Pb3[44]=B3[1][2]*B3[1][2]*cos_psi*cos_psi*Ptet+B3[1][0]*B3[1][0]*Ppsi;
      Pb3[45]=-B3[1][2]*cos_psi*cos_tet*Ptet;
      Pb3[46]=-B3[1][2]*B3[2][2]*cos_psi*sin_psi*Ptet-B3[1][0]*B3[2][1]*Ppsi;
      Pb3[47]=B3[1][2]*B3[2][2]*cos_psi*cos_psi*Ptet+B3[1][0]*B3[2][0]*Ppsi;
      Pb3[48]=B3[1][2]*B3[1][2]*cos_psi*cos_gam*Ptet;

      Pb3[55]=cos_tet*cos_tet*Ptet;
      Pb3[56]=B3[2][2]*cos_tet*sin_psi*Ptet;
      Pb3[57]=-B3[2][2]*cos_tet*cos_psi*Ptet;
      Pb3[58]=-B3[1][2]*cos_tet*cos_gam*Ptet;

      Pb3[66]=B3[2][1]*B3[2][1]*Ppsi+B3[2][2]*B3[2][2]*
      sin_psi*sin_psi*Ptet+B3[0][0]*B3[2][1]*Pgam;
      Pb3[67]=-B3[2][1]*B3[2][0]*Ppsi-B3[2][2]*B3[2][2]*
      sin_psi*cos_psi*Ptet+B3[0][0]*B3[0][1]*Pgam;
      Pb3[68]=-B3[2][2]*B3[1][2]*sin_psi*cos_gam*Ptet+B3[0][0]*B3[0][2]*Pgam;

      Pb3[77]=B3[2][0]*B3[2][0]*Ppsi+B3[2][2]*B3[2][2]*
      cos_psi*cos_psi*Ptet+B3[0][1]*B3[0][1]*Pgam;
      Pb3[78]=B3[0][1]*B3[0][2]*Pgam+B3[2][2]*B3[1][2]*cos_psi*cos_gam*Ptet;

      Pb3[88]=B3[0][2]*B3[0][2]*Pgam+B3[1][2]*B3[1][2]*cos_gam*cos_gam*Ptet;

      Pko[0][0]=uu[9]*Pb3[0]+uu[30]*Pb3[11]+uu[42]*Pb3[22]+2*uu[12]*Pb3[1]+
      2*uu[15]*Pb3[2]+2*uu[33]*Pb3[12]+bb[0]*Pur[11]+bb[9]*Pur[44]+
      bb[17]*Pur[77]+2*bb[1]*Pur[14]+2*bb[2]*Pur[17]+
      2*bb[10]*Pur[47];
      Ответить

    • ava guest 28.06.2009 21:49 # 0

      запостивший:
      Pko[0][1]=uu[9]*Pb3[3]+uu[30]*Pb3[14]+uu[42]*Pb3[25]+uu[12]*(Pb3[4]+
      Pb3[13])+uu[33]*(Pb3[15]+Pb3[24])+uu[15]*(Pb3[5]+Pb3[23])+
      bb[3]*Pur[11]+bb[12]*Pur[44]+bb[20]*Pur[77]+(bb[4]+bb[11])*
      Pur[14]+(bb[5]+bb[18])*Pur[17]+(bb[13]+bb[19])*Pur[47];

      Pko[0][2]=uu[9]*Pb3[6]+uu[30]*Pb3[17]+uu[42]*Pb3[28]+uu[12]*(Pb3[7]+
      Pb3[16])+uu[15]*(Pb3[8]+Pb3[26])+uu[33]*(Pb3[18]+Pb3[27])+
      bb[6]*Pur[11]+bb[15]*Pur[44]+bb[23]*Pur[77]+(bb[7]+bb[14])*
      Pur[14]+(bb[8]+bb[21])*Pur[17]+(bb[16]+bb[22])*Pur[47];

      Pko[0][3]=uu[10]*Pb3[0]+uu[31]*Pb3[11]+uu[43]*Pb3[22]+(uu[13]+uu[19])*
      Pb3[1]+(uu[16]+uu[22])*Pb3[2]+(uu[34]+uu[37])*Pb3[12]+
      bb[9]*Pur[45]+bb[17]*Pur[78]+bb[10]*(Pur[48]+Pur[57]);

      Pko[0][4]=uu[10]*Pb3[3]+uu[13]*Pb3[4]+uu[16]*Pb3[5]+uu[19]*Pb3[13]+
      uu[34]*Pb3[15]+uu[22]*Pb3[23]+uu[37]*Pb3[24]+uu[43]*Pb3[25]+
      uu[31]*Pb3[14]+bb[12]*Pur[45]+bb[13]*Pur[48]+bb[19]*Pur[57]+
      bb[20]*Pur[78];

      Pko[0][5]=uu[10]*Pb3[6]+uu[13]*Pb3[7]+uu[16]*Pb3[8]+uu[19]*Pb3[16]+
      uu[31]*Pb3[17]+uu[34]*Pb3[18]+uu[22]*Pb3[26]+uu[37]*Pb3[27]+
      uu[43]*Pb3[28]+bb[15]*Pur[45]+bb[16]*Pur[48]+bb[22]*Pur[57]+
      bb[23]*Pur[78];

      Pko[1][1]=uu[9]*Pb3[33]+uu[30]*Pb3[44]+uu[42]*Pb3[55]+2*uu[12]*Pb3[34]+
      2*uu[15]*Pb3[35]+2*uu[33]*Pb3[45]+bb[24]*Pur[11]+bb[30]*
      Pur[44]+bb[35]*Pur[77]+2*bb[25]*Pur[14]+2*bb[26]*Pur[17]+
      2*bb[31]*Pur[47];

      Pko[1][2]=uu[9]*Pb3[36]+uu[30]*Pb3[47]+uu[42]*Pb3[58]+uu[12]*(Pb3[37]+
      Pb3[46])+uu[15]*(Pb3[38]+Pb3[56])+uu[33]*(Pb3[48]+Pb3[57])+
      bb[27]*Pur[11]+bb[33]*Pur[44]+bb[38]*Pur[77]+(bb[28]+bb[32])*
      Pur[14]+(bb[29]+bb[36])*Pur[17]+(bb[34]+bb[37])*Pur[47];

      Pko[1][3]=uu[10]*Pb3[3]+uu[19]*Pb3[4]+uu[22]*Pb3[5]+uu[13]*Pb3[13]+
      uu[31]*Pb3[14]+uu[37]*Pb3[15]+uu[16]*Pb3[23]+uu[34]*Pb3[24]+
      uu[43]*Pb3[25]+bb[12]*Pur[45]+bb[13]*Pur[57]+bb[19]*Pur[48]+
      Ответить

    • ava guest 28.06.2009 21:50 # 0

      запостивший:
      bb[20]*Pur[78];

      Pko[1][4]=uu[10]*Pb3[33]+uu[31]*Pb3[44]+uu[43]*Pb3[55]+(uu[13]+uu[19])*
      Pb3[34]+(uu[16]+uu[22])*Pb3[35]+(uu[34]+uu[37])*Pb3[45]+
      bb[30]*Pur[45]+bb[35]*Pur[78]+bb[31]*(Pur[48]+Pur[57]);

      Pko[1][5]=uu[19]*Pb3[46]+uu[31]*Pb3[47]+uu[34]*Pb3[48]+uu[22]*Pb3[56]+
      uu[37]*Pb3[57]+uu[43]*Pb3[58]+uu[10]*Pb3[36]+uu[13]*Pb3[37]+
      uu[16]*Pb3[38]+bb[33]*Pur[45]+bb[34]*Pur[48]+bb[37]*Pur[57]+
      bb[38]*Pur[78];

      Pko[2][2]=uu[9]*Pb3[66]+uu[30]*Pb3[77]+uu[42]*Pb3[88]+2*uu[12]*Pb3[67]+
      2*uu[15]*Pb3[68]+2*uu[33]*Pb3[78]+bb[39]*Pur[11]+bb[42]*
      Pur[44]+bb[44]*Pur[77]+2*bb[40]*Pur[14]+2*bb[41]*Pur[17]+
      2*bb[43]*Pur[47];

      Pko[2][3]=uu[10]*Pb3[6]+uu[19]*Pb3[7]+uu[22]*Pb3[8]+uu[13]*Pb3[16]+
      uu[31]*Pb3[17]+uu[37]*Pb3[18]+uu[16]*Pb3[26]+uu[34]*Pb3[27]+
      uu[43]*Pb3[28]+bb[15]*Pur[45]+bb[16]*Pur[57]+bb[22]*Pur[48]+
      bb[23]*Pur[78];

      Pko[2][4]=uu[10]*Pb3[36]+uu[19]*Pb3[37]+uu[21]*Pb3[38]+uu[13]*Pb3[46]+
      uu[31]*Pb3[47]+uu[37]*Pb3[48]+uu[16]*Pb3[56]+uu[34]*Pb3[57]+
      uu[43]*Pb3[58]+bb[33]*Pur[45]+bb[34]*Pur[57]+bb[37]*Pur[48]+
      bb[38]*Pur[78];
      Ответить

    • ava guest 28.06.2009 21:50 # 0

      запостивший:
      Pko[2][5]=uu[10]*Pb3[66]+uu[31]*Pb3[77]+uu[43]*Pb3[88]+(uu[13]+uu[19])*
      Pb3[67]+(uu[16]+uu[22])*Pb3[68]+(uu[34]+uu[37])*Pb3[78]+
      bb[42]*Pur[45]+bb[44]*Pur[78]+bb[43]*(Pur[48]+Pur[57]);

      Pko[3][3]=uu[17]*Pb3[0]+uu[35]*Pb3[11]+uu[44]*Pb3[22]+2*uu[20]*
      Pb3[1]+2*uu[23]*Pb3[2]+2*uu[38]*Pb3[12]+bb[9]*Pur[55]+
      bb[17]*Pur[88]+2*bb[10]*Pur[58];

      Pko[3][4]=uu[17]*Pb3[3]+uu[35]*Pb3[14]+uu[44]*Pb3[25]+uu[20]*(Pb3[4]+
      Pb3[13])+uu[23]*(Pb3[5]+Pb3[23])+uu[38]*(Pb3[15]+Pb3[24])+
      bb[12]*Pur[55]+bb[20]*Pur[88]+(bb[13]+bb[19])*Pur[58];

      Pko[3][5]=uu[17]*Pb3[6]+uu[35]*Pb3[17]+uu[44]*Pb3[28]+uu[20]*(Pb3[7]+
      Pb3[16])+uu[23]*(Pb3[8]+Pb3[26])+uu[38]*(Pb3[18]+Pb3[27])+
      bb[15]*Pur[55]+bb[23]*Pur[88]+(bb[16]+bb[22])*Pur[58];

      Pko[4][4]=uu[17]*Pb3[33]+uu[35]*Pb3[44]+uu[44]*Pb3[55]+2*uu[20]*
      Pb3[34]+2*uu[23]*Pb3[35]+2*uu[38]*Pb3[45]+bb[30]*Pur[55]+
      bb[35]*Pur[88]+2*bb[31]*Pur[58];

      Pko[4][5]=uu[17]*Pb3[36]+uu[35]*Pb3[47]+uu[44]*Pb3[58]+uu[20]*(Pb3[37]+
      Pb3[46])+uu[23]*(Pb3[38]+Pb3[56])+uu[38]*(Pb3[48]+Pb3[57])+
      bb[33]*Pur[55]+bb[38]*Pur[88]+(bb[34]+bb[37])*Pur[58];

      Pko[5][5]=uu[17]*Pb3[66]+uu[35]*Pb3[77]+uu[44]*Pb3[88]+2*uu[20]*
      Pb3[67]+2*uu[23]*Pb3[68]+2*uu[38]*Pb3[78]+bb[42]*Pur[55]+
      bb[44]*Pur[88]+2*bb[43]*Pur[58];
      Ответить

    • ava guest 28.06.2009 21:50 # 0

      запостивший:
      for(j=1;j<6;j++)
      Pko[j][0]=Pko[0][j];
      for(j=2;j<6;j++)
      Pko[j][1]=Pko[1][j];
      for(j=3;j<6;j++)
      Pko[j][2]=Pko[2][j];
      for(j=4;j<6;j++)
      Pko[j][3]=Pko[3][j];
      Pko[5][4]=Pko[4][5];

      Pko[6][6]=som[0]*som[0];
      Pko[7][7]=som[1]*som[1];
      Pko[8][8]=som[2]*som[2];

      Gko[6][0]=sqrt(hag*2./Tauom[0])*som[0];
      Gko[7][1]=sqrt(hag*2./Tauom[1])*som[1];
      Gko[8][2]=sqrt(hag*2./Tauom[2])*som[2];
      }
      Ответить

    • ava guest 28.06.2009 21:52 # −1

      запостивший:
      извините, что так всё засрал, но уж очень хотелось высказаться. в своё время мне надо было разобраться в этом коде... наверное не стоит говорить, что это как раз тот случай когда проще заново написать чем зарбираться
      Ответить

    • ava guest 29.06.2009 00:26 # +3

      Ну ты и наговнокодил. Наклал.
      Ответить

    Добавить комментарий