Сгладить двумерную матрицу, подститать сумму элементов ниже главной диагонали

#include <iostream.h>
#include <iomanip.h>
#include <conio.h>
#include <math.h>
 
//prototype functions
 
//functions vivoda matrix
void Vyvod (float **a, int nrow,int ncol);
void Vyvod1(float **a1,int nrow,int ncol);
 
//function sglazivaniya matrix
void Smooth_matrix(float **a,int nrow,int ncol,float **a1);
 
 
//function summ fabs elements down main diagonal
void Summ_fabs(float **a1,int nrow,int ncol,float summ);
 
int main()
{
int   nrow,// number row in starting matrix
      ncol,// number column in starting matrix
      i,j, // index of starting matrix
      m,n; // index of smooth matrix
float summ;// summ elements
 
clrscr();// clear screan
 
cout<<"Enter number row and column starting matrix: ";
cin>>nrow>>ncol;// enter number row and column of starting matrix
 
cout<<"Enter elements in matrix. It is size "<<nrow
    <<" on "<<ncol<<". "<<endl;
 
float **a=new float *[nrow];  // pick out memory for starting matrix
for (i=0;i<nrow;i++)
        a[i]=new float[ncol];
 
 
 
for (i=0;i<nrow;i++)      // enter matrix user
     for(j=0;j<ncol;j++)
        cin>>a[i][j];
 
float **a1=new float *[nrow];// pick out memory for smooth matrix
for (m=0;m<nrow;m++)
        a1[m]=new float[ncol];
 
cout<<"Starting matrix: "<<endl;
Vyvod(a,nrow,ncol);       // control print starting matrix
Smooth_matrix(a,nrow,ncol,a1);
cout<<"Smooth matrix is: "<<endl;
Vyvod1(a1,nrow,ncol);     // control print smooth matrix
Summ_fabs(a1,nrow,ncol,summ);
cout<<"Summ of elements down at main diagonal"
    <<"in smooth matrix is: "<<summ<<endl;
 
cout<<"For exit programm press any key...";
getch();
return 0;
}
 
void Vyvod(float **a,int nrow,int ncol){
        int i,j;
    cout.precision(3),cout.setf(ios::showpoint);
    cout.setf(ios::left,ios::adjustfield);
    cout.setf(ios::fixed,ios::floatfield);
    for(i=0;i<nrow;i++){
            for(j=0;j<ncol;j++)
                cout<<setw(4)<<a[i][j]<<" ";
        cout<<endl;
    }
}
 
void Vyvod1(float **a1,int nrow,int ncol){
        int m,n;
    cout.precision(3),cout.setf(ios::showpoint);
    cout.setf(ios::left,ios::adjustfield);
    cout.setf(ios::fixed,ios::floatfield);
    for(m=0;m<nrow;m++){
        for(n=0;n<ncol;n++)
        cout<<setw(4)<<a1[m][n]<<" ";
        cout<<endl;
    }
}
 
void Smooth_matrix (float **a,int nrow,int ncol,float **a1){
int i,j,m,n;
for (i=0,m=0;i<nrow,m<nrow;i++,m++){
    for (j=0,n=0;j<ncol,n<ncol;j++,n++){
    if(i==0||j==0)           a[i][j]=(a[i+1][j]+a[i][j+1])/2;
    if(i==0||j==ncol)        a[i][j]=(a[i+1][j]+a[i][j-1])/2;
    if(i==nrow||j==0)        a[i][j]=(a[i-1][j]+a[i][j+1])/2;
        if(i==nrow||j==ncol)     a[i][j]=(a[i-1][j]+a[i][j-1])/2;
    if(i==0||j>0||j<ncol)    a[i][j]=(a[i][j-1]+a[i][j+1]+a[i+1][j])/3;
    if(i==nrow||j>0||j<ncol) a[i][j]=(a[i][j-1]+a[i][j+1]+a[i-1][j])/3;
    if(i>0||i<nrow||j==0)    a[i][j]=(a[i-1][j]+a[i+1][j]+a[i][j+1])/3;
    if(i>0||i<nrow||j==ncol) a[i][j]=(a[i-1][j]+a[i+1][j]+a[i][j-1])/3;
    else if(i>0||i<nrow||j>0||j<ncol) a[i][j]=(a[i-1][j]+a[i][j+1]+a[i+1][j]+a[i][j-1])/4;
   a1[m][n]=a[i][j];
   }
  }
  return;
}
 
void Summ_fabs(float **a1,int nrow,int ncol,float summ)
{
int m,n;
summ=0;
for(m=1;m<nrow;m++){
   for(n=0;n<m;n++)
     summ+=fabs(a1[m][n]);
   }
}