Loader for Python classes implementing the Spectrum interface. More...

#include <GyotoPython.h>

Inheritance diagram for Gyoto::Spectrum::Python:

Public Types
typedef Gyoto::SmartPointer< Gyoto::SmartPointee >	Subcontractor_t(Gyoto::FactoryMessenger *, std::vector< std::string > const &)
	A subcontractor builds an object upon order from the Factory. More...

Public Member Functions
virtual Property const *	getProperties () const
	Get list of properties. More...

	Python (const Python &)

virtual Python *	clone () const
	Cloner.

virtual std::string	module () const
	Return module_.

virtual void	module (const std::string &)
	Set module_ and import the Python module. More...

virtual std::string	inlineModule () const
	Return inline_module_.

virtual void	inlineModule (const std::string &)
	Set inline_module_ and import the Python module. More...

virtual std::string	klass () const
	Retrieve class_.

virtual void	klass (const std::string &)
	Set class_ and instantiate the Python class. More...

virtual size_t	instance () const
	Retrieve pInstance_. More...

virtual void	instance (size_t address)

virtual void	detachInstance ()
	Detach pInstance_ and cached method pointers. More...

virtual void	attachInstance (PyObject *instance)
	Attach pInstance_ and cached method pointers. More...

virtual std::vector< double >	parameters () const
	Retrieve parameters_.

virtual void	parameters (const std::vector< double > &)
	Set parameters_ and send them to pInstance_. More...

virtual double	operator() (double nu) const
	I_nu = mySpectrum(nu), nu in Hz. Assumes optically thick regime.

virtual double	operator() (double nu, double opacity, double ds) const
	I_nu in optically thin regime. More...

virtual double	integrate (double nu1, double nu2)
	Integrate optically thick I_nu. More...

virtual bool	knowsProperty (const std::string &name) const

virtual void	set (std::string const &key, Value val)

virtual void	set (std::string const &key, Value val, std::string const &unit)

virtual void	set (Property const &p, Value val)

virtual void	set (Property const &p, Value val, std::string const &unit)

virtual Value	get (std::string const &key) const

virtual Value	get (std::string const &key, std::string const &unit) const

Value	get (Property const &p, std::string const &unit) const

Value	get (Property const &p) const

virtual int	setParameter (std::string name, std::string content, std::string unit)

virtual void	setParameter (Gyoto::Property const &p, std::string const &name, std::string const &content, std::string const &unit)
	Set parameter by Property (and name) More...

virtual void	fillProperty (Gyoto::FactoryMessenger *fmp, Property const &p) const

virtual void	fillElement (Gyoto::FactoryMessenger *fmp) const

void	setParameters (Gyoto::FactoryMessenger *fmp)

virtual double	integrate (double nu1, double nu2, const Spectrum::Generic *opacity, double ds)
	Integrate optically thin I_nu. More...

void	incRefCount ()
	Increment the reference counter. Warning: Don't mess with the counter.

int	decRefCount ()
	Decrement the reference counter and return current value. Warning: Don't mess with the counter.

int	getRefCount ()
	Get the current number of references.

virtual bool	isThreadSafe () const
	Whether this class is thread-safe. More...

Property const *	property (std::string const pname) const
	Find property by name. More...

std::string	describeProperty (Gyoto::Property const &p) const
	Format desrciption for a property. More...

void	help () const
	Print (to stdout) some help on this class. More...

virtual std::string	kind () const
	Get kind_.

virtual void	instance (PyObject *pinstance)
	Set pInstance_. More...

virtual bool	hasPythonProperty (std::string const &key) const

virtual void	setPythonProperty (std::string const &key, Value val)

virtual void	setPythonProperty (std::string const &key, Value val, std::string const &unit)

virtual Value	getPythonProperty (std::string const &key) const

virtual Value	getPythonProperty (std::string const &key, std::string const &unit) const

virtual Gyoto::Property::type_e	pythonPropertyType (std::string const &key) const

void	checkModuleForSingleClass ()
	Look for single class in pModule_. More...

PyObject *	instantiateClass (std::string &klass) const
	Creates an instance of class `klass` in module module_. More...

Public Attributes
	GYOTO_OBJECT_THREAD_SAFETY

Static Public Attributes
static GYOTO_OBJECT Property const	properties []

Protected Member Functions
virtual void	kind (const std::string)
	Set kind_. More...

Protected Attributes
PyObject *	pCall_
	Reference to ___call__. More...

PyObject *	pIntegrate_
	Reference to the (optional) integrate method.

bool	pCall_overloaded_
	Whether call is overloaded. More...

std::string	kind_
	The "kind" that is output in the XML entity. More...

std::vector< std::string >	plugins_
	The plug-ins that needs to be loaded to access this instance's class. More...

std::string	module_
	Name of the Python module that holds the class. More...

std::string	inline_module_
	Python source code for module that holds the class.

std::string	class_
	Name of the Python class that we want to expose. More...

std::vector< double >	parameters_
	Parameters that this class needs. More...

PyObject *	pModule_
	Reference to the python module once it has been loaded.

PyObject *	pInstance_
	Reference to the python instance once it has been instantiated.

PyObject *	pProperties_
	Reference to the properties member.

PyObject *	pSet_
	Reference to the (optional) Set method.

PyObject *	pGet_
	Reference to the (optional) Get method.

Friends
class	Gyoto::SmartPointer< Gyoto::Spectrum::Python >

Detailed Description

Loader for Python classes implementing the Spectrum interface.

It interfaces with a Python class which must implement at least the call method.

Sample XML file:

<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<Scenery>
  The goal of this example is to exhibit using a Metric and a Spectrum
  coded in Python. It is functionaly equivalent to
  example-fixed-star-minkowski-cartesian.xml in the main Gyoto example
  directory.
  The metric is implemented in gyoto_sample_metrics.py. The spectrum
  and opacity of the fixed star are implemented in
  gyoto_sample_spectra.py. Both files must be in the Python path or in
  the same directory as this XML file.
  An optically thin blob centered on the origin of the coordinate
  system, in flat space-time. Computation uses the Cartesian
  coordinates.
  <Metric kind = "Python" plugin="fallback:python*">
    We specify python3 as a "fallback" plug-in as the python plug-in
    can be compiled for several versions of Python. If loading this
    file from Python, the right python plug-in should have been loaded
    already. When loading this file from the gyoto command-line tool,
    gyoto will try to load all plug-ins starting with 'python' until
    it finds one that contains the metric class 'Python'.
    <Mass unit="sunmass"> 4e6 </Mass>
    <Cartesian/>
    <Module>gyoto_sample_metrics</Module>
    <Class>Minkowski</Class>
  </Metric>
  <Screen>
    <Distance unit="kpc"> 8 </Distance>
    <Time unit="yr"> 30e3 </Time>
    <FieldOfView unit="microas"> 150 </FieldOfView>
    In UTF-8 locales, "microas" may be written "µas".
    <Inclination unit="degree"> 90 </Inclination>
    In UTF-8 locales, "degree" may be written "°".
    <PALN> 0 </PALN>
    <Argument> 0 </Argument>
    <Resolution> 32 </Resolution>
  </Screen>
  <Astrobj kind = "FixedStar">
    <Radius> 12 </Radius>
    <Position> 0 0 0 </Position>
    <Spectrum kind="Python" plugin="fallback:python*">
      <Module>gyoto_sample_spectra</Module>
      <Class>PowerLaw</Class>
      <Constant>0.001</Constant>
      <Exponent>0.</Exponent>
    </Spectrum>
    <Opacity kind="Python" plugin="fallback:python*">
      <Module>gyoto_sample_spectra</Module>
      <Class>PowerLaw</Class>
      <Constant>0.001</Constant>
      <Exponent>0.</Exponent>
    </Opacity>
    <OpticallyThin/>
  </Astrobj>
  <Delta> 1e0 </Delta>
  <MinimumTime> 0. </MinimumTime>
  <Quantities>
    Intensity
  </Quantities>
  One can also specify a unit (if Gyoto was compiled with --with-udunits):
    Intensity[mJy/pix²]
    Intensity[mJy/µas²]
    Intensity[J.m-2.s-1.sr-1.Hz-1]
    Intensity[erg.cm-2.s-1.sr-1.Hz-1]
    Intensity[mJy.sr-1]
    Intensity[Jy.sr-1]
</Scenery>

Sample Python module:

'''Sample spectra for using with Gyoto Python plug-in
   Those classes demonstrate how to use Python classes as Gyoto
   Spectrum implementations using Gyoto's "python" plug-in. Note that
   this plug-in can be renamed to whatever matches the particular
   version of Python it has been built against (e.g. python3.4).
   The goal is to be able to instantiate these from XML, from
   Yorick... and even from Python using the gyoto extension...
   Synopsis:
   import gyoto.core
   gyoto.core.requirePlugin("python") # or python2.7 or python3.4...
   sp=gyoto.core.Spectrum("Python")
   sp.set("Module", "gyoto_sample_spectra")
   sp.set("Class", "PowerLaw") # or "BlackBody6000"
   sp.set("Parameters", (1., 2.))
   val=sp(3e8/2e-6)
   Classes that aim at implementing the Gyoto::Spectrum::Generic
   interface do so by providing the following members and methods:
   __call__(self, nu): mandatory;
   __setitem__: optional;
   integrate: optional;
   properties, set, get: optional.
'''
import math
import numpy
import gyoto.python
class BlackBody6000:
    '''Black-body spectrum at 6000K
    Parameters: none.
    This example is pretty minimal: it's a black-body with fixed
    temperature. We implement the only mandaroty function:
    __call__(self, nu).
    '''
    def __call__(self, nu):
        '''spectrum(frequency_in_Hz) = black-body distribution for T=6000K.
        This function implements only
        Gyoto::Spectrum::Python::operator()(double nu).
        It does so by not accepting the varargs argument. Contrast
        with the PowerLaw definition of __call__.
        '''
        temperature=6000.
        PLANCK_OVER_BOLTZMANN=4.7992373e-11
        return nu*nu*nu/(math.exp(PLANCK_OVER_BOLTZMANN*nu/temperature)-1.);
class PowerLaw(gyoto.python.SpectrumBase):
    '''Powerlaw spectrum
    Properties: (Constant, Exponent)
    This example is pretty complete. It implements everything usefull
    and some eye-candy.
    '''
    properties={"Constant": {"type": "double", "default": 0.},
                "Exponent": {"type": "double", "default": 0.}}
    def __call__(self, *args):
        '''spectrum(frequency_in_Hz) = Constant * nu**Exponent
        This function implements both
        Spectrum::Python::operator()(double nu).
        and
        Spectrum::Python::operator()(double nu, double opacity, double ds).
        This behavior is obtained by having the varargs *args as
        second argument instead of a normal variable.
        The second overloaded function is here exactly the same as the
        C++ generic implementation and therefore useless. It is here
        to illustrate the API.
        '''
        nu=args[0]
        if (len(args)==1):
            return self.Constant * math.pow(nu, self.Exponent)
        else:
            opacity=args[1]
            ds=args[2]
            thickness=(opacity*ds)
            if (thickness):
                return self(nu) * (1.-math.exp(-thickness))
            return 0.
    def integrate(self, nu1, nu2):
        '''
        If present, this function implements
        Gyoto::Spectrum::Python::integrate(double nu1, double nu2)
        If absent, the generic integrator is used.
        '''
        if (self.Exponent == -1.):
            return self.Constant * (math.log(nu2) -math.log(nu1))
        return self.Constant * (math.pow(nu2, self.Exponent+1)- math.pow(nu1, self.Exponent+1)) / (self.Exponent+1)

Member Typedef Documentation

◆ Subcontractor_t

typedef Gyoto::SmartPointer<Gyoto::SmartPointee> Gyoto::SmartPointee::Subcontractor_t(Gyoto::FactoryMessenger *, std::vector< std::string > const &)

inherited

A subcontractor builds an object upon order from the Factory.

Various classes need to provide a subcontractor to be able to instantiate themselves upon order from the Factory. A subcontractor is a function (often a static member function) which accepts a pointer to a FactoryMessenger as unique parameter, communicates with the Factory using this messenger to read an XML description of the object to build, and returns this objet. SmartPointee::Subcontractor_t* is just generic enough a typedef to cast to and from other subcontractor types: Astrobj::Subcontractor_t, Metric::Subcontractor_t, Spectrum::Subcontractor_t. A subcontractor needs to be registered using the relevant Register() function: Astrobj::Register(), Metric::Register(), Spectrum::Register().

Member Function Documentation

◆ attachInstance()

virtual void Gyoto::Spectrum::Python::attachInstance ( PyObject * instance )

virtual

Attach pInstance_ and cached method pointers.

Attaches pInstance_, #pProperties, pSet_ and pGet_. Increments their reference counters.

Derived classes should call Base::attachInstance and attach the other pointers.

Reimplemented from Gyoto::Python::Base.

◆ checkModuleForSingleClass()

void Gyoto::Python::Base::checkModuleForSingleClass ( )

inherited

Look for single class in pModule_.

If class_ is empty and pModule_ is not null, check whether there is a single class in pModule_. In this case set class_ to its name and return. If pModule_ contains no class, error out.

◆ describeProperty()

std::string Gyoto::Object::describeProperty ( Gyoto::Property const & p ) const

inherited

Format desrciption for a property.

Returns a string containing the name(s) and type of the property, as well as whether it supports unit.

◆ detachInstance()

virtual void Gyoto::Spectrum::Python::detachInstance ( )

virtual

Detach pInstance_ and cached method pointers.

Detaches (=calls Py_XDECREF and sets to NULL) #pProperties, pSet_ and pGet_.

Derived classes should detach the other pointers and call Base::detachInstance().

Reimplemented from Gyoto::Python::Base.

◆ getProperties()

virtual Property const* Gyoto::Spectrum::Python::getProperties ( ) const

virtual

Get list of properties.

This method is declared automatically by the GYOTO_OBJECT macro and defined automatically by the GYOTO_PROPERTY_END macro.

Reimplemented from Gyoto::Spectrum::Generic.

◆ help()

void Gyoto::Object::help ( ) const

inherited

Print (to stdout) some help on this class.

Describe all properties that this instance supports.

◆ inlineModule()

virtual void Gyoto::Spectrum::Python::inlineModule ( const std::string & )

virtual

Set inline_module_ and import the Python module.

Side effects:

sets module_ to "";
calls klass(class_) if class_ is already known, so module_ can be set (or reset) after class_.

Reimplemented from Gyoto::Python::Base.

◆ instance() [1/2]

virtual void Gyoto::Python::Base::instance ( PyObject * pinstance )

virtualinherited

Set pInstance_.

Detach the previously set instance, sets pModule_ to nullptr and class_ to "", then calls attachInstance() if instance is not null.

◆ instance() [2/2]

virtual size_t Gyoto::Spectrum::Python::instance ( ) const

virtual

Retrieve pInstance_.

Returns a borrowed reference to pInstance_, as an size_t integer.

◆ instantiateClass()

PyObject* Gyoto::Python::Base::instantiateClass ( std::string & klass ) const

inherited

Creates an instance of class klass in module module_.

Returns a new reference.

◆ integrate() [1/2]

virtual double Gyoto::Spectrum::Generic::integrate	(	double	nu1,
		double	nu2,
		const Spectrum::Generic *	opacity,
		double	ds
	)

virtualinherited

Integrate optically thin I_nu.

See operator()(double nu, double opacity, double ds) const

Parameters

nu1,nu2	boundaries for the integration
opacity	the frequency-dependent opacity law given as a pointer to a Gyoto::Spectrum::Generic sub-class instance
ds	the element length for spatial integration

Returns: I, the integral of I_nu between nu1 and nu2

◆ integrate() [2/2]

virtual double Gyoto::Spectrum::Python::integrate	(	double	nu1,
		double	nu2
	)

virtual

Integrate optically thick I_nu.

See operator()(double nu) const

Parameters

nu1,nu2 boundaries for the integration

Returns: I, the integral of I_nu between nu1 and nu2

Reimplemented from Gyoto::Spectrum::Generic.

◆ isThreadSafe()

virtual bool Gyoto::Object::isThreadSafe ( ) const

virtualinherited

Whether this class is thread-safe.

Return True if this object is thread-safe, i.e. if an instance and its clone can be used in parallel threads (in the context of Scenery::raytrace()). Known objects which are not thread-safe include Lorene metrics and everything from the Python plug-in.

The default implementation considers that the class itself is thread safe and recurses into the declared properties to check whether they are safe too. Classes that abide to the Object/Property paradigm and are themselves thread-safe have nothing special to do.

Objects that clone children in their copy constructor that are not declared as properties must take these children into account.

Classes that are never thread-safe must declare it. It acn be easily done using GYOTO_OBJECT_THREAD_SAFETY in the class declaration and GYOTO_PROPERTY_THREAD_UNSAFE in the class definition.

◆ kind()

virtual void Gyoto::Object::kind ( const std::string )

protectedvirtualinherited

Set kind_.

kind(const std::string) is protected because, for most Objects, it should not be changed in runtime.Set kind_

◆ klass()

virtual void Gyoto::Spectrum::Python::klass ( const std::string & c )

virtual

Set class_ and instantiate the Python class.

Sets pInstance_.

This generic implementation takes care of the common ground, but does not set 'this' or call parameters(parameters_). Therefore, all the derived classes should reimplement this method and at least call Python::Base::klass(c) and parameters(parameters_). Between the two is the right moment to check that the Python class implements the required API and to cache PyObject* pointers to class methods.

Reimplemented from Gyoto::Python::Base.

◆ module()

virtual void Gyoto::Spectrum::Python::module ( const std::string & )

virtual

Set module_ and import the Python module.

Side effects:

sets inline_module_ to "";
calls klass(class_) if class_ is already known, so module_ can be set (or reset) after class_.

Reimplemented from Gyoto::Python::Base.

◆ operator()()

virtual double Gyoto::Spectrum::Python::operator()	(	double	nu,
		double	opacity,
		double	ds
	)		const

virtual

I_nu in optically thin regime.

Generic implementation assumes emissivity = opacity.

Parameters

nu	frequency in Hz
opacity	such that opacity*ds=optical thickness.
ds	in geometrical units

Reimplemented from Gyoto::Spectrum::Generic.

◆ parameters()

virtual void Gyoto::Spectrum::Python::parameters ( const std::vector< double > & )

virtual

Set parameters_ and send them to pInstance_.

The parameters are sent to the class instance using the setitem method with numerical keys.

Reimplemented from Gyoto::Python::Base.

◆ property()

Property const* Gyoto::Object::property ( std::string const pname ) const

inherited

Find property by name.

Look into the Property list for a Property whose name (or name_false, for a boolean Property) is pname. Return a const pointer to the first such property found, or NULL if none is found.

◆ setParameter()

virtual void Gyoto::Object::setParameter	(	Gyoto::Property const &	p,
		std::string const &	name,
		std::string const &	content,
		std::string const &	unit
	)

virtualinherited

Set parameter by Property (and name)

This function is used when parsing an XML description, if Property (p) of this name is found (i.e. either p.name or p.name_false is equal to name). Implementation should fall-back on calling the direct's parent implementation:

class A: public Object {};
class B: public A { 
 using B::setParameter;
 virtual void setParameter(Gyoto::Property const &p,
                           std::string name,
                        std::string content,
                        std::string unit);
};
void B::setParameter(Gyoto::Property const &p,
                          std::string name,
                  std::string content,
                  std::string unit) {
  if (name=="Duff") doSomething(content, unit);
  else A::setParameter(p, name, content, unit);
}

Parameters

p	Property that matches name (p.name == name or p.name_false == name)
name	XML name of the parameter (XML entity)
content	string representation of the value
unit	string representation of the unit

Reimplemented in Gyoto::Astrobj::PolishDoughnut.

Member Data Documentation

◆ class_

std::string Gyoto::Python::Base::class_

protectedinherited

Name of the Python class that we want to expose.

Property name: Class.

◆ kind_

std::string Gyoto::Object::kind_

protectedinherited

The "kind" that is output in the XML entity.

E.g. for an Astrobj, fillElement() will ensure

is written.

◆ module_

std::string Gyoto::Python::Base::module_

protectedinherited

Name of the Python module that holds the class.

For instance, if the class is implemented in toto.py, the module name is "toto". Property name: Module.

◆ parameters_

std::vector<double> Gyoto::Python::Base::parameters_

protectedinherited

Parameters that this class needs.

A list of parameters (doubles) can be passed in the Property Parameters. They will be sent to the Python instance using setitem.

◆ pCall_

PyObject* Gyoto::Spectrum::Python::pCall_

protected

Reference to ___call__.

call is the method in the underlying Python class that implements Gyoto::Spectrum::Generic::operator()().

◆ pCall_overloaded_

bool Gyoto::Spectrum::Python::pCall_overloaded_

protected

Whether call is overloaded.

This is determined automatically by looking at the parameters accepted by call:

def __call__(self, nu)

In this case call is not overloaded and implements only virtual double operator()(double nu) const;

def __call__(self, *args)

In this case call is overloaded and implements both double operator()(double nu) const and virtual double operator()(double nu, double opacity, double ds) const.

◆ plugins_

std::vector<std::string> Gyoto::Object::plugins_

protectedinherited

The plug-ins that needs to be loaded to access this instance's class.

E.g. for an Astrobj, fillElement() will ensure

is written.

The documentation for this class was generated from the following file:

GyotoPython.h

Public Types

Public Member Functions

Public Attributes

Static Public Attributes

Protected Member Functions

Protected Attributes

Friends

Detailed Description

Member Typedef Documentation

◆ Subcontractor_t

Member Function Documentation

◆ attachInstance()

◆ checkModuleForSingleClass()

◆ describeProperty()

◆ detachInstance()

◆ getProperties()

◆ help()

◆ inlineModule()

◆ instance() [1/2]

◆ instance() [2/2]

◆ instantiateClass()

◆ integrate() [1/2]

◆ integrate() [2/2]

◆ isThreadSafe()

◆ kind()

◆ klass()

◆ module()

◆ operator()()

◆ parameters()

◆ property()

◆ setParameter()

Member Data Documentation

◆ class_

◆ kind_

◆ module_

◆ parameters_

◆ pCall_

◆ pCall_overloaded_

◆ plugins_