gemseo / core

base_formulation module

The base class for all formulations.

class gemseo.core.base_formulation.BaseFormulation(disciplines, objective_name, design_space, maximize_objective=False, grammar_type='JSONGrammar', **options)[source]

Bases: object

Base MDO formulation class to be extended in subclasses for use.

This class creates the MDOFunction instances computing the constraints, objective and observables from the disciplines and add them to the attached opt_problem.

It defines the multidisciplinary process, i.e. dataflow and workflow, implicitly.

By default,

  • the objective is minimized,

  • the type of a constraint is equality,

  • the activation value of a constraint is 0.

The link between the instances of MDODiscipline, the design variables and the names of the discipline outputs used as constraints, objective and observables is made with the MDOFunctionGenerator, which generates instances of MDOFunction from the disciplines.

Parameters:
  • disciplines (Sequence[MDODiscipline]) – The disciplines.

  • objective_name (str | Sequence[str]) – The name(s) of the discipline output(s) used as objective. If multiple names are passed, the objective will be a vector.

  • design_space (DesignSpace) – The design space.

  • maximize_objective (bool) –

    Whether to maximize the objective.

    By default it is set to False.

  • grammar_type (str) –

    The type of the input and output grammars, either MDODiscipline.JSON_GRAMMAR_TYPE or MDODiscipline.SIMPLE_GRAMMAR_TYPE.

    By default it is set to “JSONGrammar”.

  • **options (Any) – The options of the formulation.

add_constraint(output_name, constraint_type='eq', constraint_name=None, value=None, positive=False)[source]

Add a user constraint.

A user constraint is a design constraint in addition to the formulation specific constraints such as the targets (a.k.a. consistency constraints) in IDF.

The strategy of repartition of constraints is defined in the formulation class.

Parameters:
  • output_name (str) – The name of the output to be used as a constraint. For instance, if g_1 is given and constraint_type=”eq”, g_1=0 will be added as a constraint to the optimizer.

  • constraint_type (str) –

    The type of constraint, either “eq” for equality constraint or “ineq” for inequality constraint.

    By default it is set to “eq”.

  • constraint_name (str | None) – The name of the constraint to be stored, If None, the name is generated from the output name.

  • value (float | None) – The value of activation of the constraint. If None, the value is equal to 0.

  • positive (bool) –

    Whether to consider an inequality constraint as positive.

    By default it is set to False.

Return type:

None

add_observable(output_names, observable_name=None, discipline=None)[source]

Add an observable to the optimization problem.

The repartition strategy of the observable is defined in the formulation class.

Parameters:
  • output_names (str | Sequence[str]) – The name(s) of the output(s) to observe.

  • observable_name (str | None) – The name of the observable.

  • discipline (MDODiscipline | None) – The discipline computing the observed outputs. If None, the discipline is detected from inner disciplines.

Return type:

None

classmethod get_default_sub_options_values(**options)[source]

Get the default values of the sub-options of the formulation.

When some options of the formulation depend on higher level options, the default values of these sub-options may be obtained here, mainly for use in the API.

Parameters:

**options (str) – The options required to deduce the sub-options grammar.

Returns:

Either None or the sub-options default values.

Return type:

dict

get_expected_dataflow()[source]

Get the expected data exchange sequence.

This method is used for the XDSM representation and can be overloaded by subclasses.

Returns:

The expected sequence of data exchange where the i-th item is described by the starting discipline, the ending discipline and the coupling variables.

Return type:

list[tuple[gemseo.core.discipline.MDODiscipline, gemseo.core.discipline.MDODiscipline, list[str]]]

get_expected_workflow()[source]

Get the expected sequence of execution of the disciplines.

This method is used for the XDSM representation and can be overloaded by subclasses.

For instance:

  • [A, B] denotes the execution of A, then the execution of B

  • (A, B) denotes the concurrent execution of A and B

  • [A, (B, C), D] denotes the execution of A, then the concurrent execution of B and C, then the execution of D.

Returns:

A sequence of elements which are either an ExecutionSequence or a tuple of ExecutionSequence for concurrent execution.

Return type:

list[gemseo.core.execution_sequence.ExecutionSequence, tuple[gemseo.core.execution_sequence.ExecutionSequence]]

get_optim_variables_names()[source]

Get the optimization unknown names to be provided to the optimizer.

This is different from the design variable names provided by the user, since it depends on the formulation, and can include target values for coupling for instance in IDF.

Returns:

The optimization variable names.

Return type:

list[str]

get_sub_disciplines()[source]

Accessor to the sub-disciplines.

This method lists the sub scenarios’ disciplines.

Returns:

The sub-disciplines.

Return type:

list[gemseo.core.discipline.MDODiscipline]

classmethod get_sub_options_grammar(**options)[source]

Get the sub-options grammar.

When some options of the formulation depend on higher level options, the schema of the sub-options may be obtained here, mainly for use in the API.

Parameters:

**options (str) – The options required to deduce the sub-options grammar.

Returns:

Either None or the sub-options grammar.

Return type:

JSONGrammar

get_sub_scenarios()[source]

List the disciplines that are actually scenarios.

Returns:

The scenarios.

Return type:

list[Scenario]

get_top_level_disc()[source]

Return the disciplines which inputs are required to run the scenario.

A formulation seeks to compute the objective and constraints from the input variables. It structures the optimization problem into multiple levels of disciplines. The disciplines directly depending on these inputs are called top level disciplines.

By default, this method returns all disciplines. This method can be overloaded by subclasses.

Returns:

The top level disciplines.

Return type:

list[gemseo.core.discipline.MDODiscipline]

get_x_mask_x_swap_order(masking_data_names, all_data_names=None)[source]

Mask a vector from a subset of names, with respect to a set of names.

This method eventually swaps the order of the values if the order of the data names is inconsistent between these sets.

Parameters:
  • masking_data_names (Iterable[str]) – The names of the kept data.

  • all_data_names (Iterable[str] | None) – The set of all names. If None, use the design variables stored in the design space.

Returns:

The masked version of the input vector.

Raises:
  • IndexError – when the sizes of variables are inconsistent.

  • ValueError – when the names of variables are inconsistent.

Return type:

ndarray

get_x_names_of_disc(discipline)[source]

Get the design variables names of a given discipline.

Parameters:

discipline (MDODiscipline) – The discipline.

Returns:

The names of the design variables.

Return type:

list[str]

mask_x_swap_order(masking_data_names, x_vect, all_data_names=None)[source]

Mask a vector from a subset of names, with respect to a set of names.

This method eventually swaps the order of the values if the order of the data names is inconsistent between these sets.

Parameters:
  • masking_data_names (Iterable[str]) – The names of the kept data.

  • x_vect (ndarray) – The vector to mask.

  • all_data_names (Iterable[str] | None) – The set of all names. If None, use the design variables stored in the design space.

Returns:

The masked version of the input vector.

Raises:

IndexError – when the sizes of variables are inconsistent.

Return type:

ndarray

unmask_x_swap_order(masking_data_names, x_masked, all_data_names=None, x_full=None)[source]

Unmask a vector from a subset of names, with respect to a set of names.

This method eventually swaps the order of the values if the order of the data names is inconsistent between these sets.

Parameters:
  • masking_data_names (Iterable[str]) – The names of the kept data.

  • x_masked (ndarray) – The boolean vector to unmask.

  • all_data_names (Iterable[str] | None) – The set of all names. If None, use the design variables stored in the design space.

  • x_full (ndarray) – The default values for the full vector. If None, use the zero vector.

Returns:

The vector related to the input mask.

Raises:

IndexError – when the sizes of variables are inconsistent.

Return type:

ndarray

NAME: ClassVar[str] = 'MDOFormulation'

The name of the MDO formulation.

property design_space: DesignSpace

The design space on which the formulation is applied.

disciplines: Sequence[MDODiscipline]

The disciplines of the MDO process.

opt_problem: OptimizationProblem

The optimization problem generated by the formulation from the disciplines.

class gemseo.core.base_formulation.BaseFormulationsFactory(cls, module_names=None)[source]

Bases: object

A factory of BaseFormulation.

Parameters:
  • cls (type[Any]) –

  • module_names (Iterable[str] | None) –

create(formulation_name, disciplines, objective_name, design_space, maximize_objective=False, **options)[source]

Create a formulation.

Parameters:
  • formulation_name (str) – The name of a class implementing a formulation.

  • disciplines (Sequence[MDODiscipline]) – The disciplines.

  • objective_name (str) – The name(s) of the discipline output(s) used as objective. If multiple names are passed, the objective will be a vector.

  • design_space (DesignSpace) – The design space.

  • maximize_objective (bool) –

    Whether to maximize the objective.

    By default it is set to False.

  • **options – The options for the creation of the formulation.

Return type:

BaseFormulation

is_available(formulation_name)[source]

Check the availability of a formulation.

Parameters:

formulation_name (str) – The formulation name to check.

Returns:

Whether the formulation is available.

Return type:

bool

property formulations: list[str]

The available formulations.