我建议查看一个与时间相关的示例，因为如果您在该级别制定，PETSc 可以提供更多诊断。例如，您可以使用 Rosenbrock 方法、加法 Runge-Kutta IMEX 方法或其他方法。有些涉及牛顿迭代，但这不是必需的，并不总是最好的方法。要使用这些方法，您可以从提供残差函数开始，$f(t,u)$，如$\dot u = f(t,u)$，或者，$g$在$g(t,u,\dot u) = 0$. PETSc 可以通过着色自动计算雅可比行列式，如果您后来决定组装雅可比行列式以提高效率，PETSc 可以检查您组装的内容是否正确并显示错误。如果您推出自己的时间积分和自己的牛顿，除非您自己编写，否则这些算法变体和调试功能将不可用。

用户手册解释了接口。另请参阅特定 API 文档的手册页。手册页与示例交叉引用。一个很好的例子是Brusselator src/ts/examples/tutorials/ex25.c。您可能还会发现教程中的幻灯片很有用。如果您遇到困难，请随时在此处或petsc-users邮件列表中提问。

即使在串行中，您使用文件系统的建议也会成为一个巨大的瓶颈。它极大地使工作流程复杂化，并且会限制您可以使用的方法。但是，如果您只想测试一次矩阵，您可以PetscBinaryWrite()从 MATLAB 使用并加载系统并使用src/ksp/ksp/examples/tutorials/ex10.c.