Petsc求解poisson方程，两种方法的比较

最新推荐文章于 2023-12-13 20:02:30 发布
最新推荐文章于 2023-12-13 20:02:30 发布
阅读量1k 收藏 3
点赞数 1
分类专栏： 高性能计算
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/forrestguang/article/details/113800733
版权
高性能计算 专栏收录该内容
46 篇文章 29 订阅 ¥239.90 ¥99.00
订阅专栏

第一种方法只需要利用前面提到的矩阵和向量搭建

前面几篇博客提到了矩阵，向量的搭建，本人也介绍来求解任意大小的线性系统的方法，下面本人先介绍不调用Petsc任何有关的包的求解方法，为此先引入基本定义：
− Δ u = f , x ∈ Ω = ( 0 , 1 ) 2 -\Delta u = f,x \in \Omega=(0,1)^2

了解本专栏 订阅专栏 解锁全文
Galerkin码农选手
关注
  • 1
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
订阅专栏
python实现三角，矩形有限元求解二维泊松方程
forrestguang的博客
01-22 1513
重点强调：以下内容都是本人的上机作业，仅仅为了方便同僚学习交流，希望不会有抄袭我博客以及代码当成自己原创作业的现象发生 问题提出 {−Δu=f,x∈Ω,u=u0(x),x∈∂Ω0,∂u∂n=g−βu,x∈∂Ω1.\left \{ \begin{array}{r l} - \Delta u = f,x \in \Omega,\\ u = u_{0}(x),x \in \partial \Omega_{0},\\ \frac{\partial u}{\partial n} = g
基于PETSC的线性方程求解
forrestguang的博客
10-27 2627
使用KSP方法求解小型线性系统 需要注意的是，在PETSC中出现的语法和C语言基本相同，区别大部分都只存在于不同的表达方式 比如在C中，int表示整型，在PETSC中，PetscInt表示整型，这里大家百度用法，或者参考我上次的学习报告，我们这里重点介绍在PETSC中一些简单函数的定义调用，并且使用这种函数来求解实际问题（求解线性系统）。 我们慢慢来解释：第一句static char help[] = “Solve a 4x4 linear system using KSP.\n”;是为了对这个代码进行说明
SOLVERLAB：使用SALOMEMEDCoupling网格和PETSc矩阵的PDE求解器实验室
02-08
SOLVERLAB SOLVERLAB是一个几何和数值C ++ / Python库，为数值分析人员设计，他们致力于将一般形状和网格上的偏微分方程离散化，而更侧重于高级脚本编写。 目的是提供简单的MATLAB样式函数，以生成和操纵网格，字段和矩阵。 该库源自[CDMATH]（ ），这是一个同名的协作工作组。 它基于[SALOME]（ ）的[MEDcoupling]（ ）C ++ / python库 ）项目来处理网格和字段，并在C ++库[PETSC]（ ）上用于处理矩阵和线性求解器。 SOLVERLAB包括PDE系统，该系统是由核React堆堆芯建模产生的，涉及流体动力学，热和中子扩散以及固体弹性。 这是一个简单的环境，对学生和研究人员来说，是一种教学环境，可以在具有非结构化网格的一般几何形状上推广新的数值方法。 SOLVERLAB的主要研究目标是研究 一般网格上低马赫数可压缩流的
二维poisson方程c语言,二维Poisson方程并行求解算法.ppt
weixin_30587341的博客
05-24 515
二维Poisson方程并行求解算法第八讲二维 Poisson 方程的并行求解算法 主要内容 二维Poisson方程 二维Poisson方程 Jacobi 迭代 程序示例 并行算法 并行算法 并行算法 并行算法 并行算法 上机作业 * * 二维 Poisson 方程的差分离散 差分方程的 Jacobi 算法 串行算法 并行算法 红黑排序的 GS 算法 二维 Poisson 方程 其中 ?...
PETSc 库解线性方程组Ax = b
We need more……
04-16 3777
program main implicit none #include "finclude/petscsys.h" #include "finclude/petscvec.h" #include "finclude/petscmat.h" #include "finclude/petscksp.h" #include "finclude/petscpc.h" Vec x, b, u Mat A KSP ksp PC pc PetscReal no
Petsc求解一维非线性微分方程，SNES对象-2
forrestguang的博客
02-22 464
微分方程引入 −u′′+ρ(u)=f,x∈(0,1)-u^{''}+\rho\sqrt(u)=f,x \in(0,1)−u′′+ρ(​u)=f,x∈(0,1) u(0)=α,u(1)=βu(0)=\alpha,u(1)=\betau(0)=α,u(1)=β 我们假设剖分为N，那么就利用中心差分得到： h=1N−1h=\frac{1}{N-1}h=N−11​ −ui−1+2ui−ui+1−(R(ui)+fi)h2=0,1≤i≤N−2-u_{i-1}+2u_i-u_{i+1}-(R(u_i)+f_i)h^2=0
多重网格法解泊松方程（两步法）
程序员的博客
09-15 623
这里是一个一维泊松方程，采用两步法求解
含参泊松方程求解-初步探索
forrestguang的博客
03-16 961
问题引入 −Δu=f,x∈Ω=(0,1)2-\Delta u=f,x \in \Omega=(0,1)^2−Δu=f,x∈Ω=(0,1)2 u=g,x∈∂Ωu=g,x \in \partial \Omegau=g,x∈∂Ω u=exp⁡(−(x−a)2+(y−b)22l2),a,b,lu = \exp(-\frac{(x - a)^2 + (y - b)^2}{2l^2}),a,b,lu=exp(−2l2(x−a)2+(y−b)2​),a,b,l是参数，其中a∈[0,1],b∈[0,1],l∈[0.1,1]
poisson:使用有限元方法求解任意二维域上的泊松方程
06-12
任意二维域上的泊松 使用有限元方法求解具有 RHS f 和 Dirichlet 边界条件的任意二维域上的泊松方程。
python：二维有限元方法求解，泊松方程
走向CTO的路上...
12-13 1282
在求解泊松方程时，二维有限元方法将求解域分割成多个小区域（单元），并在每个单元内使用一个局部函数（有限元形函数）逼近原方程的解。《The Finite Element Method: Linear Static and Dynamic Finite Element Analysis》（作者：Thomas J. R. Hughes）- 这本书介绍了有限元方法在结构力学中的应用，包括泊松方程的求解。计算物理学：二维有限元方法在计算物理学中也有广泛的应用，如量子力学中的薛定谔方程、流体力学中的势流方程等。
中科院PETSc使用指南（中文）
09-15
对原版的petsc进行了翻译，特别适合于初级学习者使用。内容比较全面。
poisson-solver:泊松方程的一般求解器
04-04
泊松解算器 泊松方程的一般求解器 工作正在进行中
python求解偏微分方程
04-23
python求解偏微分方程，画图，容易修改。
基于PETSc 的有限元高性能求解方法 (2013年)
04-23
基于科学计算工具箱PETSc 提供的应用开发平台，针对单机多核计算机开发出一种线性并行求解器。用C++类的形式封装求解线性有限元问题必需的数据和操作，并行使用迭代法求解线性有限元问题。算例测试结果表明，在内存...
RSchr-5:用SLEPc和PETSc求解单电子Schrödinger方程-开源
04-29
RSchr-5计算单电子Schrödinger方程的最小特征值。 解决的问题是在有限域（即一个框）上定义的。 通过零Dirichlet边界条件解决了该问题。 所实现的算法使用以B样条为基础函数的有限元法。 离散化导致广义特征值问题...
偏微分方程组求解库PETSC源代码3.1
01-17
PETSc (Portable, Extensible Toolkit for Scientific Computation) 是美国能源部ODE2000支持开发的20多个ACTS工具箱之一，由Argonne国家实验室开发的可移植可扩展科学计算工具箱，主要用于在分布式存储环境高效求解...
多重网格算法求解二维泊松方程的python实现
forrestguang的博客
05-14 787
问题描述 二维泊松方程 {−Δu=f, in Ω,u=g, on ∂Ω,\left\{\begin{aligned} -\Delta u &=f, \text { in } \Omega, \\ u &= g ,\text { on }\partial \Omega,\\ \end{aligned}\right.{−Δuu​=f, in Ω,=g, on ∂Ω,​
泊松过程的模拟
weixin_42428226的博客
07-06 4365
1.
FFT解泊松方程
加菲猫的博客
07-26 6995
Pérez, Patrick, Gangnet M , Blake A . Poisson image editing[J]. ACM Transactions on Graphics, 2003, 22(3):313. 上一篇：泊松融合笔记 待求解的问题 Δf=∇⋅v(x,y)=div(v(x,y)),∀(x,y)∈Ω,withf∣R∖Ω=f∗∣R∖Ω\Delta f=\nabla\cd...
PETSc科学计算库
最新发布
12-28
PETSc（Portable, Extensible Toolkit for Scientific Computation）是一个并行科学计算库，它基于MPI（Message Passing Interface）实现了一套数据结构和库函数，用于解决基于偏微分方程的典型科学计算问题。PETSc支持多种编程语言，包括C、C++、Fortran和Python。 PETSc提供了丰富的功能和工具，可以用于求解各种科学计算问题，如线性和非线性方程组的求解、特征值问题的求解、优化问题的求解等。它具有高度的可扩展性和灵活性，可以在大规模并行计算环境中高效地进行计算。 以下是一个使用PETSc库求解线性方程组的Python示例： ```python from petsc4py import PETSc # 创建线性方程组的系数矩阵和右侧向量 A = PETSc.Mat().createAIJ(size=(3, 3)) A.setValues([[0, 0, 1], [0, 1, 0], [1, 0, 0]], [1, 2, 3]) A.assemble() b = PETSc.Vec().createSeq(size=3) b.setValues([4, 5, 6]) b.assemble() # 创建线性方程组的解向量 x = PETSc.Vec().createSeq(size=3) # 创建线性方程组的求解器 ksp = PETSc.KSP().create() ksp.setOperators(A) ksp.setFromOptions() # 求解线性方程组 ksp.solve(b, x) # 打印解向量 print("Solution vector:") x.view() # 释放资源 A.destroy() b.destroy() x.destroy() ksp.destroy() ``` 这个示例演示了如何使用PETSc库创建线性方程组的系数矩阵和右侧向量，然后使用PETSc提供的求解器求解线性方程组，并打印解向量。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

Galerkin码农选手

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值