跳到内容

JuliaInv/jInv.jl公司

生成状态 覆盖率状态 生成状态

jInv公司

jInv公司是Julia中PDE参数估计的灵活框架。它提供了易于扩展的用于PDE约束反问题的核心函数。我们的目标是高效并行地解决参数估计问题。有关更多详细信息,请参阅(http://arxiv.org/abs/1606.07399).

概述

jInv由五个子模块组成:

  1. ForwardShare公司-并行求解正问题的方法。
  2. 反向求解-反问题中常用的方法,如失配函数、正则化和数值优化。
  3. 网格-2D和3D中的规则网格和张量网格以及微分算子。
  4. 线性解算器-稀疏和(如果安装)直接线性解算器的接口,可用于求解离散化的PDE。
  5. 实用程序-实用程序函数

要求

jInv用于Julia版本0.7。并要求:

  1. Krylov方法.jl-求解(稀疏)线性系统的迭代方法。

附加(可选)软件包可实现更高的性能。ParSpMatVec.jl公司-用于稀疏矩阵向量乘积的共享存储器实现。jInv公司自动检测此包是否正确设置并默认使用。

附加附加软件包:

  1. MUMPSjInv.jl公司-MUMPS包装。用作直接线性解算器。
  2. “Pardiso.jl”-使用此包的扩展存在,但未针对Julia版本1.0进行测试。

jInv中的网格模块具有规则网格和张量网格,但也可以通过以下方式进行扩展

  1. JOcTree(作业目录树)-高效的空间自适应网格

安装

julia类型:

整件克隆(“https://github.com/JuliaInv/jInv.jl网站“,”jInv“)包装测试(“jInv”)

示例

一些反演示例可以在示例文件夹。

使用jInv的包

jInv为PDE参数估计问题提供了通用方法。为了将其用于应用,一些方法需要通过指定正向问题的求解器、灵敏度矩阵向量积等进行扩展。这应该组织在小模块中。目前实施的PDE模型如下:

  1. DivSigGrad.jl公司-静力学中的电导率反问题
  2. jInvSeismic.jl公司-地震反演软件包:全波形反演(jInvSeismic.FWI)和走时层析成像(jInvSeismic.EikonalInv)
  3. 麦克斯韦频率-麦克斯韦方程中电导率的反演

致谢

本材料部分基于国家科学基金会(批准号1522599)支持的工作。本材料中表达的任何观点、发现、结论或建议均为作者的观点,并不一定反映国家科学基金会的观点。