GammaLib

The following test script

#!/usr/bin/env python
import gammalib

# Save map
map = gammalib.GSkyMap('CAR','GAL',0.0,0.0,1.0,1.0,360,180)
map.save('map.fits', True)

# Create diffuse map
diffuse_map = gammalib.GModelSpatialDiffuseMap('map.fits')
print('=== Diffuse map ===')
print(diffuse_map.value())

# Create diffuse cube
diffuse_cube = gammalib.GModelSpatialDiffuseCube('map.fits')
print('=== Diffuse cube ===')
print(diffuse_cube.value())
print(diffuse_cube.maps())
print(diffuse_cube.pixels())

=== Diffuse map ===
1.0
=== Diffuse cube ===
1.0
0
0

I decided to remove the GModelSpatialDiffuseCube::maps() and GModelSpatialDiffuseCube::pixels() methods since the same information can be accessed through the GModelSpatialDiffuseCube::cube().nmaps() and GModelSpatialDiffuseCube::cube().npix() methods.

I added a fetch_cube() call in the GModelSpatialDiffuseCube::cube() method so that the cube actually gets fetched before any access. This fixes also the GModelSpatialDiffuseCube::cube() method.

Two ctools unit tests needed to be adapted since they were using the GModelSpatialDiffuseCube::maps() and GModelSpatialDiffuseCube::pixels() methods.

Merged into devel.