insert_update_vinv

pyoptex.doe.cost_optimal.codex.formulas.insert_update_vinv(Vinv, Zs, pos, a, b, ratios)

Computes the update to Vinv based on the insertion of a row a into the design. b specifies possible additional group updates if necessary (e.g. when breaking up a group).

Parameters

Vinvnp.array(3d)

The current inverses of the observation covariance matrix.

Zslist(np.array(1d) or None)

The grouping matrices of all factors.

posint

The position where to insert the new run.

anp.array(1d)

The group of the newly inserted run.

blist(tuple(row_start, row_stop, group_from, group_to))

A list of operations to apply changing the groups of certain rows.

ratiosnp.array(2d)

The variance ratios of the factors in each row.

Returns

Zilist(np.array(1d) or None)

The new grouping matrices (non-ascending)

Vinvnp.array(3d)

The updated inverses observation covariance matrix.