Eirik Hatlen / gammalib

…from HealPix version 3.20.

The bilinear interpolator class holds indices and weights for bilinear interpolation. It has an interpolation operator that operates on a linear array of double precision values.

… GSkymap object. Existing pixel values are kept, any new pixel values will be set to zero. Add an appropriate unit test.

The XML module now handles predefined entities. Special characters, such as " will now be transformed into their predefined entities, e.g. ", when writing the XML file. Upon reading, predefined entities will be converted into special chacaters, so that the user always deals with a human readable string.

A new parameter was introduce to specify the acceptable function decrease in an iteration. Acceptable means that the new parameters are kept for continuing the iterations. This behaviour existed already before, but a low threshold of 1e-6 did not really activate it. Now the threshold was raised to 1.0 and this help ctlike to converge quicker. Methods accept_dec() have been implemented to get and set the acceptable threshold.

Conflicts:
	ChangeLog
	inst/cta/src/GCTAResponse_helpers.cpp
	inst/cta/src/GCTASupport.hpp

Any elliptical model should be defined by semiminor and semimajor axes, these are in fact used in the CTA response computation. Hence they should live in the baseclass, not the derived class.

…ly the new Python 3.x division logic.

A constructor has been added to create a table with a specific number of rows and columns. This constructor can be used for table creation.

Set methods string(), real() and integer() have been added to set the table cells to specific values.

A save() method has been implemented that will save the CSV table into an ASCII file.

Hierarchical element access methods have been added to GXml and GXmlNode.Using a string in the format "params > param[1] > value > struct" one can in a single method call walk down the hierarchy of an XML document to access an element. The squared brackets allow specifying an element index in case that multiple elements with the same name exist. If squared brackets are omitted, an index of 0 is assumed. If the indicated element does not exist, the method returns a pointer to NULL.

The convert() and set() methods now support time references given in the UTC time system. Using convert() one can now extract for example MJD in the UTC time system. All this needs still to be documented.

This required implementation of a leap_seconds() computation method. The leap second table is statically implement in the method and needs manual update when leap seconds are added.

The new method stack_maps() now works on the object it is applied to instead of creating a new sky map. The reason is on the one hand semantic, as "stack_maps" is an action applied on the current object. For an object returning a skymap, the name should be something like "stacked" as the new map is a stacked version of the old one. But another reason is memory usage. The method now allocates a new memory space for the stacked map and releases the old memory once the stacking is done.

The new implementation is also more general as this method now works for both WCS and HEALPIX projections.

A unit test has been added to make sure that the maps were indeed stacked.

…ll GModelData derivatives

…or the observation Python module.

The classname() method returns in a human readable format the class name. This method is needed for full automatic typecasting in Python. Full automatic typecasting has for the moment been implemented for the observation module.

Automatic type casting was implemented in the GObservations.i interface using the runtime type casting information that is maintained by the SWIG interface. The function relies on a newly introduced type() method in GObservation that provides as a string the class type.

The former GCTAObservation type casting constructor has been removed as it now should be obsolete.