OpenMDAO两次生成着色文件

Question

我想知道这是否正确。我创建的大多数隐式和显式组件都使用以下行：

self.declare_coloring(wrt='*', method='cs', tol=1.0E-12, show_sparsity=True)

然后，当我找到运行我使用的驱动程序的文件时：

p.driver.declare_coloring()

在我的/coloring_files目录中，每个组件都有一个“col”和一个“磁盘”。

coloring_traj_phases_phase0_rhs_col_brakeThrottle.pkl    coloring_traj_phases_phase0_rhs_disc_implicitOutputs.pkl
coloring_traj_phases_phase0_rhs_col_implicitOutputs.pkl  coloring_traj_phases_phase0_rhs_disc_powerTrain.pkl
coloring_traj_phases_phase0_rhs_col_powerTrain.pkl       coloring_traj_phases_phase0_rhs_disc_spin.pkl
coloring_traj_phases_phase0_rhs_col_spin.pkl             coloring_traj_phases_phase0_rhs_disc_timeAdder.pkl
coloring_traj_phases_phase0_rhs_col_timeAdder.pkl        coloring_traj_phases_phase0_rhs_disc_timeSpace.pkl
coloring_traj_phases_phase0_rhs_col_timeSpace.pkl        coloring_traj_phases_phase0_rhs_disc_tracking.pkl
coloring_traj_phases_phase0_rhs_col_tracking.pkl         coloring_traj_phases_phase0_rhs_disc_tyreConstraint.pkl
coloring_traj_phases_phase0_rhs_col_tyreConstraint.pkl   coloring_traj_phases_phase0_rhs_disc_tyre.pkl
coloring_traj_phases_phase0_rhs_col_tyre.pkl             total_coloring.pkl
coloring_traj_phases_phase0_rhs_disc_brakeThrottle.pkl

是否需要这两组文件，还是重复两次操作？另外，我想知道用驱动程序声明着色是否使用CS以外的方法？我确实打算使用total_coloring.pkl进行静态着色。

Answer 1

Dymos可以使用两种转录方法之一: Radau伪谱法或高阶GaussLobatto法。

GaussLobatto方法是一个两步过程：

1. 在“离散化”节点上计算ODE。
2. 离散节点上的值和速率用于将状态率和状态速率插值到“配置”节点。
3. 在配置节点上使用步骤2中的插值状态值第二次对ODE进行计算。
4. 将插值率与ODE在配置节点上输出的速率(这些缺陷称为缺陷)进行比较--如果它们很小，则假设物理是精确的。

G29

Radau转录遵循类似的过程，但配置节点是离散节点的子集，因此不需要插值，ODE只需要计算一次。

如果将转录从dymos.GaussLobatto更改为dymos.Radau，那么每个ODE组件只有一个部分着色文件。否则，两者都需要把它们的颜色分开处理。