The invariant double-differential cross section, E1E2d6 delta/dp(1)3dp(2)3, and the double-spin asymmetry, A(LL), for inclusive multi-gamma pair production in which gamma-rays came from neutral mesons were measured with a 200 GeV/c longitudinally-polarized proton beam and a longitudinally-polarized proton target. Most of the multi-gamma pairs comes from two-jet type events which are sensitive to partonic interaction. The A(LL) values were found to be consistent with zero. The invariant double-differential cross section for inclusive pi0pi0 production was also measured. These measured cross sections are consistent with LUND Monte Carlo simulations. Using the LUND Monte Carlo simulation package with the Carlitz-Kaur model of spin dependent distribution functions of valence quarks, the A(LL) values have been compared with theoretical predictions of gluon polarization, DELTAG/G. The results put restrictions on the size of DELTAG/G in the region of 0.05 less-than-or-similar-to x less-than-or-similar-to 0.35.