GammaLib

So far the GModelSpatial::flux method only works for circular regions and models. The method should be generalised so that it works for any kind of regions and models. This is probably best achieve by creating a model map for the spatial model and integrating the model map using the GSkyMap::flux() method, like in GModelSpatialDiffuseMap::flux():

    // Compute flux in sky region
    double flux = m_map.flux(region) * m_value.value();