The discretization of Cahn-Hilliard equation with obstacle potential leads to a block 2 by 2 non-linear system, where the p1, 1q block has a non-linear and non-smooth term. Recently a globally convergent Newton Schur method was proposed for the non-linear Schur complement corresponding to this non-linear system. The solver may be seen as an inexact Uzawa method which has the falvour of an active set method in the sense that the active sets are first identified by solving a quadratic obstacle problem corresponding to the p1, 1q block of the block 2 by 2 nonlinear system, and a new decent direction is obtained after discarding the active set region. The problem becomes linear on nonactive set, and corresponds to solving a linear saddle point problem on truncated domains. For solving the quadratic obstacle problem, various optimal multigrid like methods have been proposed. In this paper solvers for the truncated saddle point problem is considered. Three preconditioners are considered, two of them have block diagonal structure, and the third one has block tridiagonal structure. One of the block diagonal preconditioners is obtained by adding certain scaling of stiffness and mass matrices, whereas, the remaining two involves Schur complement. Eigenvalue bound and condition number estimates are derived for the preconditioned untruncated problem. It is shown that the extreme eigenvalues of the preconditioned truncated system remain bounded by the extreme eigenvalues of the preconditioned untruncated system. Numerical experiments confirm the optimality of the solvers.


翻译:Cahn-Hilliard 等离散且有障碍潜力的离散方程式可导致第2轮2x2非线性系统,P1, 1q区块具有非线性和非线性术语。 最近为非线性Schur 补充非线性Schur 提议了一个全球趋同的牛顿 Schur 方法, 与这个非线性系统相对应。 解答器可能被视为一种不尽实际的Uzawa方法, 其特点是, 主动型组群首先通过解决与第2轮的p1、 1q 区块相对应的二次反线性障碍, 而在放弃活动性设定区域后, 获得了一个新的体面方向。 这个问题在非线性设置上成为线性牛顿 Schur, 与解决非线性 Schurnational Schur 相配对的线性点问题相对应。 各种最优的多格方法已被提出来研究。 三个前题被考虑过, 其中两个前题号已经解决了块性双面的底值, 第三组的极性极性前置值 和先验性前置的硬性前置系统 。

0
下载
关闭预览

相关内容

在数学中,鞍点或极大极小点是函数图形表面上的一点,其正交方向上的斜率(导数)都为零,但它不是函数的局部极值。鞍点是在某一轴向(峰值之间)有一个相对最小的临界点,在交叉轴上有一个相对最大的临界点。
【KDD2021】图神经网络,NUS- Xavier Bresson教授
专知会员服务
62+阅读 · 2021年8月20日
Stabilizing Transformers for Reinforcement Learning
专知会员服务
58+阅读 · 2019年10月17日
2019年机器学习框架回顾
专知会员服务
35+阅读 · 2019年10月11日
Hierarchically Structured Meta-learning
CreateAMind
26+阅读 · 2019年5月22日
Transferring Knowledge across Learning Processes
CreateAMind
27+阅读 · 2019年5月18日
深度自进化聚类:Deep Self-Evolution Clustering
我爱读PAMI
15+阅读 · 2019年4月13日
强化学习的Unsupervised Meta-Learning
CreateAMind
17+阅读 · 2019年1月7日
A Technical Overview of AI & ML in 2018 & Trends for 2019
待字闺中
16+阅读 · 2018年12月24日
推荐|深度强化学习聊天机器人(附论文)!
全球人工智能
4+阅读 · 2018年1月30日
【推荐】深度学习时序处理文献列表
机器学习研究会
7+阅读 · 2017年11月29日
【论文】变分推断(Variational inference)的总结
机器学习研究会
39+阅读 · 2017年11月16日
Arxiv
0+阅读 · 2021年11月11日
Arxiv
5+阅读 · 2018年10月4日
Paraphrase Generation with Deep Reinforcement Learning
Arxiv
3+阅读 · 2018年2月24日
VIP会员
相关资讯
Hierarchically Structured Meta-learning
CreateAMind
26+阅读 · 2019年5月22日
Transferring Knowledge across Learning Processes
CreateAMind
27+阅读 · 2019年5月18日
深度自进化聚类:Deep Self-Evolution Clustering
我爱读PAMI
15+阅读 · 2019年4月13日
强化学习的Unsupervised Meta-Learning
CreateAMind
17+阅读 · 2019年1月7日
A Technical Overview of AI & ML in 2018 & Trends for 2019
待字闺中
16+阅读 · 2018年12月24日
推荐|深度强化学习聊天机器人(附论文)!
全球人工智能
4+阅读 · 2018年1月30日
【推荐】深度学习时序处理文献列表
机器学习研究会
7+阅读 · 2017年11月29日
【论文】变分推断(Variational inference)的总结
机器学习研究会
39+阅读 · 2017年11月16日
Top
微信扫码咨询专知VIP会员