Source code for openmdao.components.demux_comp

"""Definition of the Demux Component."""

from six import iteritems

import numpy as np

from openmdao.core.explicitcomponent import ExplicitComponent

[docs]class DemuxComp(ExplicitComponent): """ Demux one or more inputs along a given axis. Attributes ---------- _vars : dict Container mapping name of variables to be demuxed with additional data. _output_names : dict Container mapping name of variables to be demuxed with associated outputs. """
[docs] def __init__(self, **kwargs): """ Instantiate DemuxComp and populate private members. Parameters ---------- **kwargs : dict Arguments to be passed to the component initialization method. """ super(DemuxComp, self).__init__(**kwargs) self._vars = {} self._output_names = {}
[docs] def initialize(self): """ Declare options. """ self.options.declare('vec_size', types=int, default=2, desc='The number of elements to be extracted from each input.')
[docs] def add_var(self, name, val=1.0, shape=None, units=None, desc='', axis=0): """ Add an input variable to be demuxed, and all associated output variables. Parameters ---------- name : str name of the variable in this component's namespace. val : float or list or tuple or ndarray or Iterable The initial value of the variable being added in user-defined units. Default is 1.0. shape : int or tuple or list or None Shape of this variable, only required if val is not an array. Default is None. units : str or None Units in which this input variable will be provided to the component during execution. Default is None, which means it is unitless. desc : str description of the variable axis : int The axis along which the elements will be selected. Note the axis must have length vec_size, otherwise a RuntimeError is raised at setup. """ self._vars[name] = {'val': val, 'shape': shape, 'units': units, 'desc': desc, 'axis': axis}
def _post_configure(self): """ Declare inputs, outputs, and derivatives for the demux component. """ # set static mode to False because we are doing things that would normally be done in setup self._static_mode = False opts = self.options vec_size = opts['vec_size'] for var, options in iteritems(self._vars): kwgs = dict(options) shape = options['shape'] size = axis = kwgs.pop('axis') if axis >= len(shape): raise RuntimeError("{}: Invalid axis ({}) for variable '{}' of " "shape {}".format(self.msginfo, axis, var, shape)) if shape[axis] != vec_size: raise RuntimeError("{}: Variable '{}' cannot be demuxed along axis {}. Axis size " "is {} but vec_size is {}.".format(self.msginfo, var, axis, shape[axis], vec_size)) self.add_input(var, **kwgs) template = np.reshape(np.arange(size), shape) self._output_names[var] = [] out_shape = list(shape) out_shape.pop(axis) if len(out_shape) == 0: out_shape = [1] for i in range(vec_size): out_name = '{0}_{1}'.format(var, i) self._output_names[var].append(out_name) self.add_output(name=out_name, val=options['val'], shape=out_shape, units=options['units'], desc=options['desc']) rs = np.arange( cs = np.atleast_1d(np.take(template, indices=i, axis=axis)) self.declare_partials(of=out_name, wrt=var, rows=rs, cols=cs, val=1.0) self._static_mode = True super(DemuxComp, self)._post_configure()
[docs] def compute(self, inputs, outputs): """ Demux the inputs into the appropriate outputs using numpy.take. Parameters ---------- inputs : Vector unscaled, dimensional input variables read via inputs[key] outputs : Vector unscaled, dimensional output variables read via outputs[key] """ opts = self.options vec_size = opts['vec_size'] for var in self._vars: ax = self._vars[var]['axis'] for i in range(vec_size): out_name = self._output_names[var][i] outputs[out_name] = np.reshape(np.take(inputs[var], indices=i, axis=ax), newshape=outputs[out_name].shape)