void Source_linear_Item_Calculation(double k_eff,double ****Cs,double ***q_sour,double ****source_linear,double ***cx,double ***cy)
{
    int g,g2,cell,reg,azimuth;
    double ****_fiss_Lsour=NULL,**** _Scatt_Lsour=NULL;
    double _XS_nf, _chi, _XS_s;
//    if ((q_scatter = (double *)calloc(NumGroups, sizeof(double))) == NULL || \
//        (q_fission = (double *)calloc(NumGroups, sizeof(double))) == NULL)
//        exit(120);

     for(cell=0;cell<NumCellsInCore;cell++){
        for(reg=0;reg<NumFSRsInCell[cell];reg++){
          for (azimuth = 0; azimuth < NumAzimuAng; azimuth++){
             for (g = 0; g < NumGroups; g++)
             {
                  _fiss_Lsour[g][cell][reg][azimuth]=0.0;                   这里出现problem recived signal SIGSEGV,Segmentation fault错误
                   _Scatt_Lsour[g][cell][reg][azimuth]=0.0;
             }
            }
        }
    }
    for(cell=0;cell<NumCellsInCore;cell++){
        for(reg=0;reg<NumFSRsInCell[cell];reg++){
          for (azimuth = 0; azimuth < NumAzimuAng; azimuth++){
             for (g = 0; g < NumGroups; g++)
             {

                 Cs[cell][reg][azimuth][g]=0.0;
          }
        }
     }
    }

      for(cell=0;cell<NumCellsInCore;cell++){
        for(reg=0;reg<NumFSRsInCell[cell];reg++){
          for (azimuth = 0; azimuth < NumAzimuAng; azimuth++){
             for (g = 0; g < NumGroups; g++)
             {

                 Cs[g][cell][reg][azimuth]=cx[cell][reg][g]*cos(EffAzimuAng[azimuth])+cy[cell][reg][g]*sin(EffAzimuAng[azimuth]);
             }
          }
        }
     }

    for(cell=0;cell<NumCellsInCore;cell++)
    {
        for(reg=0;reg<NumFSRsInCell[cell];reg++)
        {
            for(g=0;g<NumGroups;g++)
            {
                for (azimuth = 0; azimuth < NumAzimuAng; azimuth++)
                    {
                       for(g2=0;g2<NumGroups;g2++)
                          {
                             _XS_nf = GetXSs_tafnfchi_FSR(cell, reg, g2, MacroXS_nf, Inp_XSs_nf);
                             _chi = GetXSs_tafnfchi_FSR(cell, reg, g, Chi, Inp_Chi);
                             _XS_s = GetXSs_s_FSR(cell, reg, g2, g, MacroXS_s, Inp_XSs_p0);
                            _fiss_Lsour[g][cell][reg][azimuth] += (_chi /k_eff)*(_XS_nf*Cs[g2][cell][reg][azimuth])/(4.*3.14159265);
                            _Scatt_Lsour[g][cell][reg][azimuth] += (_XS_s*Cs[g2][cell][reg][azimuth])/(4.*3.14159265);
                           }
                          source_linear[g][cell][reg][azimuth]=_fiss_Lsour[g][cell][reg][azimuth]+_Scatt_Lsour[g][cell][reg][azimuth];
                          _fiss_Lsour[g][cell][reg][azimuth]=0.;
                          _Scatt_Lsour[g][cell][reg][azimuth]=0.;
                   }
            }
        }
    }
}

Segment fault一般因为下标越界或者指针乱飘。

学习了