Разобрать код программы и выявить ошибки - C (СИ)

Узнай цену своей работы

Формулировка задачи:

Начал работать с С , немогу разобрать код, почему не работает? Заранее спасибо, за помощь + к карме
Листинг программы 
  1. int nlfm (float *uzxt, float *uml, float *a, float *uzx, float *wzx, float *tzx, float *gs, float *gu, float *gw, float *gt, int n, int lp)
  2.  {
  3.     /* u*df/dx  + w*df/dz =sum {m,l} uzxt1{m,l}*fi{m,l}t;
  4.        lp=0 -- f=u; lp=1 -- f=w; lp=2 -- f= C(T);
  5.        uzxt1{m,l,2}=-1 :
  6.        fi (t)=exp(-uzxt1{m,l,4}*t)*(uzxt1{m,l,0}*sin(uzxt1{m,l,5}*t)+uzxt1{m,l,1}*cos(uzxt1{m,l,5}*t)
  7.        uzxt1{m,l,2}=1 :
  8.        fi (t)=(uzxt1{m,l,0}*exp(*(uzxt1{m,l,4}*t)+uzxt1{m,l,1}*exp(uzxt1{m,l,5}*t). */
  9.  
  10.     int i, j, k, l, m,  ml, mli, lnj, n1, n2, n3, n4, nl, nl1, k2, ki, k1, k3, j2, j3, l1, l2, lm, li, i2, lp1;
  11.     float *ay, *c, *d, *az,  *nxt, *u0, *tzx0, *e, *f, r1, r2, r3, g1, g2, g3, g4, g5, ggu, c0; 
  12.     ay=(float *)calloc(20,sizeof(float));
  13.     az=(float *)calloc(80,sizeof(float));
  14.     c=(float *)calloc(20,sizeof(float));
  15.     d=(float *)calloc(20,sizeof(float));
  16.     e=(float *)calloc(20,sizeof(float));
  17.     f=(float *)calloc(20,sizeof(float));
  18.     nxt=(float *)calloc(80,sizeof(float));
  19.     u0=(float *)calloc(80,sizeof(float));
  20.     tzx0=(float *)calloc(80,sizeof(float));
  21.     if ( tzx0!=NULL && u0!=NULL && c!=NULL && d!=NULL && e!=NULL && f!=NULL && nxt!=NULL && az!=NULL && ay!=NULL )
  22.      { 
  23.        n3=n*n; 
  24.        lp1=lp*n*n;
  25.        n4=n*n*n;
  26.        nl=4;
  27.        for (m=0;m<n;m++)
  28.          { mli=m*n*n;
  29.            ml=8*m*n;
  30.            lm=m*n;
  31.            for (l=0;l<n;l++)
  32.              { ggu=*(gs+lm+lp1+l);
  33.                lnj=l*n*n;
  34.                c0=0.; 
  35.                for (n1=0;n1<nl;n1++)
  36.                  { k1=n1*n;
  37.                    n2=6*n1;   
  38.                    for (l1=0;l1<nl;l1++)
  39.                      { l2=6*l1;
  40.                        nl1=8*(n*n1+l1);
  41.                        for (k=0;k<nl-1;k++)
  42.                          { k2=6*k;
  43.                            k3=8*(n*n1+k);
  44.                            ki=k*n;
  45.                            for (j=0;j<nl-1;j++)
  46.                              { j2=6*j;
  47.                                j3=8*(n*k+j);
  48.                                if ( lp==0 )
  49.                                 { r1=*(uml+k1+mli+l1)*(*(a+ki+lnj+j));
  50.                                   r2=*(uml+n4+k1+mli+l1)*(*(a+n4+ki+lnj+j));
  51.                                   for (i=0;i<8;i++) 
  52.                                     { *(c+i)=*(uzx+nl1+i); 
  53.                                       *(d+i)=*(uzx+j3+i); }
  54.  /* printf("\n\t m=%u  l=%u n1=%u l1=%u  k=%u j=%u  nl1=%u,  j3=%u;  r1=%8.3g, r2=%8.3g",m,l,n1,l1,k,j,nl1,j3,r1,r2);   */                              
  55.                                   pftgt (e,c,d,0);
  56.                                   for (i=0;i<8;i++) 
  57.                                     { *(c+i)=*(wzx+nl1+i); 
  58.                                       *(d+i)=*(uzx+j3+i); }
  59.                                   pftgt (ay,c,d,0);
  60.                                   for (i=0;i<3;i++)
  61.                                     { *(e+i)*=r1;
  62.                                       *(ay+i)*=r2; }
  63.                                   for (i=0;i<6;i++) *(e+6+i)=*(ay+i);
  64.                                   pslf1(ay,e,2,3,0);
  65.                                 }
  66.                                if ( lp==1 )
  67.                                 { r1=*(uml+k1+mli+l1)*(*(a+ki+lnj+j));
  68.                                   r2=*(uml+n4+k1+mli+l1)*(*(a+n4+ki+lnj+j));
  69.                                   for (i=0;i<8;i++) 
  70.                                     { *(c+i)=*(uzx+nl1+i); 
  71.                                       *(d+i)=*(wzx+j3+i); }
  72.                                  
  73.  /* printf("\n\t m=%u  l=%u n1=%u l1=%u  k=%u j=%u  nl1=%u,  j3=%u;  r1=%8.3g, r2=%8.3g",m,l,n1,l1,k,j,nl1,j3,r1,r2);  */                               
  74.                                   pftgt (e,c,d,0);
  75.                                   for (i=0;i<8;i++) 
  76.                                     { *(c+i)=*(wzx+nl1+i); 
  77.                                       *(d+i)=*(wzx+j3+i); }
  78.                                   pftgt (ay,c,d,0);
  79.                                   for (i=0;i<3;i++)
  80.                                     { *(e+i)*=r1;
  81.                                       *(ay+i)*=r2; }
  82.                                   for (i=0;i<6;i++) *(e+6+i)=*(ay+i);
  83.                                   pslf1(ay,e,2,3,0);
  84.                                 }
  85.                                if ( lp==2 )             
  86.                                 { r1=*(uml+k1+mli+l1)*(*(a+ki+lnj+j));
  87.                                   r2=*(uml+n4+k1+mli+l1)*(*(a+n4+ki+lnj+j));     
  88.                                   for (i=0;i<8;i++) 
  89.                                     { *(c+i)=*(uzx+nl1+i); 
  90.                                       *(d+i)=*(tzx+j3+i); }
  91. /*  printf("\n\t m=%u  l=%u n1=%u l1=%u  k=%u j=%u  nl1=%u,  j3=%u;  r1=%8.3g, r2=%8.3g",m,l,n1,l1,k,j,nl1,j3,r1,r2);   */                             
  92.                                   pftgt (e,c,d,0);
  93.                                   for (i=0;i<8;i++) 
  94.                                     { *(c+i)=*(wzx+nl1+i); 
  95.                                       *(d+i)=*(tzx+j3+i); }
  96.                                   pftgt (ay,c,d,0);
  97.                                   for (i=0;i<3;i++)
  98.                                     { *(e+i)*=r1;
  99.                                       *(ay+i)*=r2; }
  100.                                   for (i=0;i<6;i++) *(e+6+i)=*(ay+i);
  101.                                   pslf1(ay,e,2,3,0);
  102.                                 }
  103.                                for (i=0;i<6;i++) *(az+j2+i)=*(ay+i);
  104.                              }     /*  j  */
  105.                    /*   printf("\n\t m=%u  l=%u n1=%u l1=%u  k=%u",m,l,n1,l1,k);  */
  106.                            pslf1(ay,az,nl-1,3,0);
  107.                            for (i=0;i<6;i++) *(nxt+k2+i)=*(ay+i);
  108.                          }        /*  k  */
  109.                   /*   printf("\n\tm=%u  l=%u n1=%u  l1=%u ",m,l,n1,l1);  */
  110.                        pslf1(ay,nxt,nl-1,3,0);    
  111.                        for (i=0;i<6;i++) *(tzx0+l2+i)=*(ay+i); 
  112.                      }       /* l1  */
  113.               /*  printf("\n\tm=%u  l=%u n1=%u  n2=%u",m,l,n1,n2);  */           
  114.                    pslf1(ay,tzx0,nl,3,0);    
  115.                    for (i=0;i<6;i++) *(u0+n2+i)=*(ay+i);       
  116.                  }       /*  n1  */
  117.              /*  printf("\n\tm=%u  l=%u",m,l);   */             
  118.                pslf1(ay,u0,nl,3,0);    
  119.                *d=1;
  120.                *(d+1)=ggu;
  121.                yprv2 (ay,d,c,0,0);
  122.  
  123.                if ( lp==0 ) r2=*(gu+lm+l);                   
  124.                if ( lp==1 )
  125.                 { r2=*(gw+lm+l);
  126.                   c0=-r2; }
  127.                if ( lp==2 )
  128.                 { c0=*(gt+lm+l);
  129.                   r2=*(gt+lm+n3+l); }
  130.  
  131.                orig2 (d,c,0,0); 
  132.                l2=8*l;
  133.                *(uzxt+ml+l2)=r2;
  134.                *(uzxt+ml+l2+1)=ggu;
  135.                *(d+2)+=c0;
  136.                for (j=0;j<6;j++) *(uzxt+ml+l2+j+2)=*(d+j);
  137.              /*  orig3(d,ay,ggu,1);  */
  138.              }   /* l  */
  139.            printf("\n\t uzxt1 -- array uzxt(ml=%u;%u,%u)\n",ml,m,n);
  140.            for (i=0;i<n;i++) 
  141.              { l2=8*i; 
  142.                for (k=0;k<8;k++)  
  143.                    printf("%10.3g",*(uzxt+ml+l2+k)); PRN; }
  144.            PRN;  
  145.          }  /*  m  */
  146.        free(tzx0); free(u0);
  147.        free(nxt); 
  148.        free(d); free(c);
  149.        free(f); free(e);
  150.        free(az);free(ay);
  151.        return(1);
  152.      }
  153.     else
  154.      { printf("\n no memory for arrays in nlf1\n");
  155.        return(-1); }
  156.   }
  157.  
  158.  float tmlk1 (float *slz, float *slx, float *gt,  float t, float z, float x, int n)
  159.   { 
  160.     /* TT^(1)(x,z,t)=TT^(0)(x,z,t)+
  161.        tt1_{4,m,l}+e^(tt1_{6,m,l}*t)*[tt1_{2,m,l}*sh(tt1_{7,m,l}*t)+tt2_{3,m,l}*ch(tt1_{7,m,l}*t)] */    
  162.     int m, k, n1, mk;
  163.     float zz, xx, r1, r2, tts, ttx;
  164.     tts=0.;
  165.     n1=n+n;
  166.     for (m=0;m<n;m++)
  167.       { 
  168.     zz=cos(*(slz+n1+m)*z);
  169.         ttx=0.;      
  170.         for (k=0;k<n;k++)
  171.           { mk=8*(m*n+k);
  172.             xx=cos(*(slx+k+n1)*x); 
  173.             r2=*(gt+mk+4)+(*(gt+mk))*exp(-(*(gt+mk+1))*t);
  174.             if (*(gt+mk+5) == -1 )
  175.                r1=exp(*(gt+mk+6)*t)*(*(gt+mk+2)*sin(*(gt+mk+7)*t)+(*(gt+mk+3))*cos(*(gt+mk+7)*t));
  176.             else
  177.                r1=exp(*(gt+mk+6)*t)*(*(gt+mk+2)*sinh(*(gt+mk+7)*t)+(*(gt+mk+3))*cosh(*(gt+mk+7)*t));
  178.           /*     r1=*(gt+mk+2)*exp(*(gt+mk+6)*t)+(*(gt+mk+3))*exp(*(gt+mk+7)*t);  */
  179.             ttx+=(r1+r2)*xx;
  180.           }
  181.         tts+=ttx*zz;
  182.       }  
  183.     return(tts);
  184.   }
  185.  
  186.  float umlk1 (float *slz, float *slx, float *gt, float un, float t, float z, float x, int n)
  187.   { 
  188.    /*  UT^(1)(x,z,t)=ut_{0,m,l}+e^(ut_{3,m,l}*t)*[ut_{1,m,l}*sh(ut_{4,m,l}*t)+ut_{2,m,l}*ch(ut_{4,m,l}*t)]   */
  189.     int m, l, k1, mk;
  190.     float zz, xx, r1, r2, uts, utx;
  191.     uts=0.;
  192.     for (m=0;m<n;m++) 
  193.       { 
  194.         zz=(sin(*(slz+m)*z)-(*(slz+m))/un*cos(*(slz+m)*z)); 
  195.         utx=0.;
  196.         for (l=0;l<n;l++)
  197.           { 
  198.             mk=8*(l+m*n);
  199.             xx=sin(*(slx+l)*x); 
  200.             r2=*(gt+mk+4)+(*(gt+mk))*exp(-(*(gt+mk+1))*t);
  201.             if (*(gt+mk+5) == -1 )
  202.                r1=exp(*(gt+mk+6)*t)*(*(gt+mk+2)*sin(*(gt+mk+7)*t)+(*(gt+mk+3))*cos(*(gt+mk+7)*t));
  203.             else
  204.                r1=exp(*(gt+mk+6)*t)*(*(gt+mk+2)*sinh(*(gt+mk+7)*t)+(*(gt+mk+3))*cosh(*(gt+mk+7)*t));
  205.              /*  r1=*(gt+mk+2)*exp(*(gt+mk+6)*t)+(*(gt+mk+3))*exp(*(gt+mk+7)*t);  */
  206.             utx+=xx*(r1+r2);            
  207.           }
  208.         uts+=zz*utx;
  209.       }
  210.     return(uts);
  211.   }
  212.  
  213.  float wmlk1 (float *slz, float *slx, float *gt, float t, float z, float x, int n)
  214.   { 
  215.    /*  w^(1)(x,z,t) */
  216.     int m, k, l, mk;
  217.     float zz, xx, r1, r2, tts, ttx;
  218.     tts=0.;
  219.     for (m=0; m<n; m++)
  220.       { zz=cos(*(slz+n+m)*z);
  221.         ttx=0.;      
  222.         for (l=0;l<n;l++)
  223.           { 
  224.             mk=8*(l+m*n);
  225.             xx=cos(*(slx+l+n)*x); 
  226.             r2=*(gt+mk+4)+(*(gt+mk))*exp(-(*(gt+mk+1))*t);
  227.             if (*(gt+mk+5) == -1 )
  228.                r1=exp(*(gt+mk+6)*t)*(*(gt+mk+2)*sin(*(gt+mk+7)*t)+(*(gt+mk+3))*cos(*(gt+mk+7)*t));
  229.             else
  230.                r1=exp(*(gt+mk+6)*t)*(*(gt+mk+2)*sinh(*(gt+mk+7)*t)+(*(gt+mk+3))*cosh(*(gt+mk+7)*t));
  231.             /*   r1=*(gt+mk+2)*exp(*(gt+mk+6)*t)+(*(gt+mk+3))*exp(*(gt+mk+7)*t);  */
  232.             ttx+=xx*(r1+r2);            
  233.           }
  234.         tts+=ttx*zz;
  235.       }
  236.     return(tts);
  237.   }
  238.  
  239.  int grf0mlk (float *slz, float *slx, float *ut, float *gs, float *u, int n, int lp)
  240.   {
  241.     FILE  *fcl;
  242.     int i, k, l, m, l1, l2, n2, nz;  
  243.     float t, x, z, un, dt, dx, dz, *umlk;
  244.     nz=2*n;
  245.     n2=nz*nz;
  246.     un=*(u+10);
  247.     umlk=(float *)calloc(n2+1,sizeof(float));
  248.     if ( umlk!=NULL )
  249.      { 
  250.        dx=*u/nz;
  251.        dz=*(u+1)/nz;
  252.        t=1.;
  253.        z=0;
  254.        dt=0.2;  
  255.        for (l=0;l<nz;l++)
  256.          { l2=l*nz;  
  257.            x=0.;
  258.            for (k=0;k<nz;k++)  
  259.              { if ( lp==0 ) *(umlk+l2+k)=umlk0 (slz, slx, ut, gs, un, t, z, x, n);   
  260.                if ( lp==1 ) *(umlk+l2+k)=wmlk0 (slz, slx, ut, gs, t, z, x, n);   
  261.                if ( lp==2 ) *(umlk+l2+k)=tmlk0 (slz, slx, ut, gs, t, z, x, n);   
  262.                x=x+dx;
  263.              }
  264.            z=z+dz;
  265.          }
  266.        printf("\n\t umlk0(z,x,t=%5.2f) lp=%u; dx=%7.3g, dz=%7.3g dt=%5.2g -- array uzx1(%u,%u)\n",t,lp,dx,dz,dt,nz,nz);
  267.        for (l=0;l<nz;l++) 
  268.          { l2=l*nz; 
  269.            for (k=0;k<nz;k++)  
  270.                printf("%8.3g%c",*(umlk+l2+k), ( k%n==n-1 ) ? '\n' : ' '); 
  271.            PRN; }
  272.     
  273.        if ( lp==0 ) fcl = fopen("gduzx0.dat","w");
  274.        if ( lp==1 ) fcl = fopen("gdwzx0.dat","w");
  275.        if ( lp==2 ) fcl = fopen("gdtzx0.dat","w");
  276.        for (i=0; i<n2; i++) fprintf(fcl,"%g,",*(umlk+i));
  277.        fclose(fcl); 
  278.  
  279.        x=2.;
  280.        z=0.;
  281.        for (m=0;m<nz;m++) 
  282.          { t=0;
  283.            n2=m*n;
  284.            l2=m*nz;
  285.            for (l=0;l<nz;l++)
  286.              { if ( lp==0 ) *(umlk+l2+l)=umlk0 (slz, slx, ut, gs, un, t, z, x, n);   
  287.                if ( lp==1 ) *(umlk+l2+l)=wmlk0 (slz, slx, ut, gs, t, z, x, n);   
  288.                if ( lp==2 ) *(umlk+l2+l)=tmlk0 (slz, slx, ut, gs, t, z, x, n);   
  289.                t=t+dt;
  290.              }
  291.            z=z+dz;
  292.          }
  293.        printf("\n\t uml0(z,x=%5.2f,t) lp=%u  -- array uml0(%u,%u)\n",x,lp,nz,nz);
  294.        for (m=0; m<nz; m++) 
  295.          { n2=m*nz;
  296.            for (l=0; l<nz; l++)
  297.               printf("%8.3g%c",*(umlk+n2+l), ( l%n==n-1 ) ? '\n' : ' '); 
  298.            PRN;
  299.          }
  300.    
  301.        if ( lp==0 ) fcl = fopen("gduzt0.dat","w");
  302.        if ( lp==1 ) fcl = fopen("gdwzt0.dat","w");
  303.        if ( lp==2 ) fcl = fopen("gdtzt0.dat","w");
  304.        for (i=0; i<n2; i++) fprintf(fcl,"%g,",*(umlk+i));
  305.        fclose(fcl); 
  306.  
  307.        z=2.;
  308.        x=0.5;
  309.        for (m=0;m<nz;m++) 
  310.          { t=0;
  311.            n2=m*n;
  312.            l2=m*nz;
  313.            for (l=0;l<nz;l++)
  314.              { if ( lp==0 ) *(umlk+l2+l)=umlk0 (slz, slx, ut, gs, un, t, z, x, n);   
  315.                if ( lp==1 ) *(umlk+l2+l)=wmlk0 (slz, slx, ut, gs, t, z, x, n);   
  316.                if ( lp==2 ) *(umlk+l2+l)=tmlk0 (slz, slx, ut, gs, t, z, x, n);   
  317.                t=t+dt;
  318.              }
  319.            x=x+dx;
  320.          }
  321.        printf("\n\t uml0(x,t,z=%5.2f) lp=%u -- array uml0(%u,%u)\n",z,lp,nz,nz);
  322.        for (m=0; m<nz; m++) 
  323.          { n2=m*nz;
  324.            for (l=0; l<nz; l++)
  325.               printf("%8.3g%c",*(umlk+n2+l), ( l%n==n-1 ) ? '\n' : ' '); 
  326.            PRN;
  327.          }
  328.    
  329.        if ( lp==0 ) fcl = fopen("gduxt0.dat","w");
  330.        if ( lp==1 ) fcl = fopen("gdwxt0.dat","w");
  331.        if ( lp==2 ) fcl = fopen("gdtxt0.dat","w");
  332.        for (i=0; i<n2; i++) fprintf(fcl,"%g,",*(umlk+i));
  333.        fclose(fcl); 
  334.  
  335.        free(umlk);
  336.        return(1);
  337.      }
  338.     else
  339.      { printf("\n no memory for umlk in grf0mlk\n");
  340.        return(-1); }
  341.   }
  342.                   
  343.  int grfmlk (float *slz, float *slx, float *ut, float *u, int n, int lp)
  344.   {
  345.     FILE  *fcl;
  346.     int i, k, l, m, l1, l2, n2, nz;  
  347.     float t, x, z, un, dt, dx, dz, *umlk;
  348.     nz=2*n;
  349.     n2=nz*nz;
  350.     un=*(u+10);
  351.     umlk=(float *)calloc(n2+1,sizeof(float));
  352.     if ( umlk!=NULL )
  353.      { 
  354.   /*     printf("\n\t array ut(%u,%u,8)\n",n,n);
  355.        for (m=0; m<n; m++) 
  356.          { n2=8*m*n;
  357.            for (l=0; l<n; l++)
  358.              { l2=8*l;
  359.                for (i=0;i<8;i++)  printf("%10.3g",*(ut+n2+l2+i)); 
  360.                PRN; }
  361.          }
  362.  ... дополнения

Решение задачи: «Разобрать код программы и выявить ошибки»

textual
Листинг программы 
  1.  float tmlk1 (float *slz, float *slx, float *gt,  float t, float z, float x, int n);
  2.  float umlk1 (float *slz, float *slx, float *gt, float un, float t, float z, float x, int n);
  3.  float wmlk1 (float *slz, float *slx, float *gt, float t, float z, float x, int n);
  4.  int grf0mlk (float *slz, float *slx, float *ut, float *gs, float *u, int n, int lp);
  5.  int grfmlk (float *slz, float *slx, float *ut, float *u, int n, int lp);

ИИ поможет Вам:

  • решить любую задачу по программированию
  • объяснить код
  • расставить комментарии в коде
  • и т.д
Попробуйте бесплатно

Оцени полезность:

12   голосов , оценка 3.75 из 5

Нужна аналогичная работа?

Оформи быстрый заказ и узнай стоимость

Бесплатно
Оформите заказ и авторы начнут откликаться уже через 10 минут
Похожие ответы