fireshape的使用（基础篇）

已于 2023-04-12 19:39:03 修改
阅读量288 收藏
点赞数
分类专栏： 传统数值方法 PDE约束优化控制问题 文章标签： mesh
于 2022-05-17 22:59:00 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/forrestguang/article/details/124831700
版权
传统数值方法 同时被 2 个专栏收录
35 篇文章 18 订阅 ¥199.90 ¥299.90
订阅专栏
PDE约束优化控制问题
27 篇文章 3 订阅 ¥259.90 ¥399.90
订阅专栏

fireshape的使用安装以及基本算例参考链接添加链接描述，这里本人将针对这个介绍做一些中文的翻译以及在代码上做一些修改和解释
fireshape的安装参考前一篇博客添加链接描述
特殊说明，后面可能会用到fireshape.zoo库，这个库使用简单的source firedrake/activate/bin 无法调用，原因未知，这里的处理方式是，写一个pre.sh文件，文件内容是

source firedrake/bin/activate
cd fireshape-master
source activate.sh
cd

然后使用source pre.sh命令进入fireshape,firedrake环境即可。下面重点对代码进行解释和重新实现，结果都会用paraview展示出来。

无约束优化问题

了解本专栏 订阅专栏 解锁全文
Galerkin码农选手
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
专栏目录
订阅专栏
firedrake的安装（docker环境）
forrestguang的博客
04-19 1533
注意以下安装过程都是默认在docker环境下的服务器进行的，用户需要被服务器管理员添加锦docker组 1:百度搜索firedrake官网，进去，然后Ctrl + f，搜索docker，找到下面这个部分 2：点击进入上面提到的docker image，复制右边的命令：docker pull firedrakeproject/firedrake，把这段命令复制到服务器终端，然后服务器会开始下载跟firedrake相关的库，大概有两三个G. 3：下载完了以后，键入命令： docker run -it – p
firedrake创建mesh和求解具体问题
forrestguang的博客
05-06 363
界面问题 −∇⋅(σ∇u)+u=5,Ω-\nabla \cdot(\sigma \nabla u) + u = 5,\Omega−∇⋅(σ∇u)+u=5,Ω u∣∂Ω1=0u|_{\partial \Omega_1} = 0u∣∂Ω1​​=0 [σ∇u⋅n]=3,in∂Ω2[\sigma\nabla u \cdot n] = 3,in \partial \Omega_2[σ∇u⋅n]=3,in∂Ω2​ ∂Ω1\partial \Omega_1∂Ω1​为矩形边界 ∂Ω2\partial \Omega_2∂Ω2
有限元开源代码firedrake_Python编写的有限元
05-18
Firedrake is an automated system for the portable solution of partial differential equations using the finite element method (FEM). Firedrake enables users to employ a wide range of discretisations to an infinite variety of PDEs and employ either conventional CPUs or GPUs to obtain the solution.
完全空白的mac os 安装firdrake 和petsc-3.19.0
最新发布
weixin_44612116的博客
03-01 547
第六步，安装petsc，sh install.sh。第一步：建议应用商店安装xcode, 并终端安装。第三步：配置git和ssh 密钥。第四步： 安装homebrew。第二步： 开启终端代理。
seigen:基于有限元法的地震波弹性方程求解器
05-14
Seigen：通过代码生成进行地震学 概述 是基于有限元框架的弹性波方程求解器。 它是地震成像项目的一部分。 快速开始 Seigen需要安装Firedrake，并且必须在Firedrake虚拟环境中运行。 要首先安装Firedrake，请按照以下。 一旦安装了Firedrake并激活了虚拟环境，就可以使用以下命令安装Seigen： git clone https://github.com/opesci/seigen.git pip install -e seigen 执照 Seigen是开源软件，并根据发行。 有关更多信息，请参见文件LICENSE.md 。 接触 评论和反馈可以发送至 。
01.Linux基础篇
01-02
Linux基础讲解文档
Python入门——基础篇
06-27
Python语言入门，基础篇，思维导图格式，精简版（输入输出、数据类型、条件语句、循环、字符串等常用操作、文件操作、函数、类和对象，异常、模块和包），仅供参考
Git使用基础篇(一些常用命令和原理)
09-30
Git是一个分布式的版本控制工具，本篇文章从介绍Git开始，重点在于介绍Git的基本命令和使用技巧，让你尝试使用Git的同时，体验到原来一个版本控制工具可以对开发产生如此之多的影响
FPGA基础篇配套光盘
02-01
Altera FPGA/CPLD 基础篇 附带的光盘的资源，里面有书上的例子全套代码，欢迎下载
java程序设计基础篇
09-24
本书是java程序设计的基础篇。是英文版的。适合英语水平较好的java初学者观看。也适合开发人员的参考用书。中文版的可以进一步咨询QQ945745429
FireBird3的安装步骤(适用于切片工具和小猫统计)
追
05-24 2164
用于切片工具和小猫统计工具的FireBird 3安装非常简单。基本上按照默认设置就行。
有限元计算软件Fenics安装使用笔记（1）
塘朗晨光
11-03 9355
开篇写点个人思考： 来深圳快一年多了，越来越感觉圈子的重要性，我不是一个喜欢抱怨的人，但有时很多事情真的是身不由己，生活注定孤独，谁不想风花雪月，吟诗弹琴呢？找三五好友，有空得聚，没空得散，但很多人走着走着就散了，唯有惋惜而已，话题说的有点远了。回归正题。送可爱的你我一句话，一个很不错的小班在我情绪低落时送给我的。 猛兽都是独行的，只有牛羊才喜欢成群结队！ 最近接触《工业分析和数据建模》，老师推荐...
android shape的使用
hyrg_的专栏
11-29 445
shape用于设定形状，可以在selector，layout等里面使用，有6个子标签，各属性如下： xml version="1.0" encoding="utf-8"?> shape xmlns:android="http://schemas.android.com/apk/res/android" > 圆角 --> corners and
针对neumann边界条件的差分法代码
forrestguang的博客
10-11 4475
需要添加人工的Dirichlet边界 考虑的问题仍然是矩形区域的泊松方程。 如果整个区域都采用Neumann边界条件，即∂ΩN=∂Ω\partial\Omega_{N}=\partial\Omega∂ΩN​=∂Ω，则这个问题没有定解，因为对于某个精确解uuu而言，针对任意的常数C，u+cu+cu+c都是解 为此，我们常常需要添加一个人工提供的Dirichlet边界条件，比如说在这里，我们选取{−1}×[−1,1]\{-1\}\times[-1,1]{−1}×[−1,1]为Dirichlet边界，也就是说矩
五点差分法求解泊松方程-共轭梯度法和G-S迭代法
forrestguang的博客
12-30 2956
差分法求解矩形区域椭圆方程 Δu=f,x∈Ω=[−1,1]∗[−1,1]\Delta u = f, x \in \Omega = [-1,1]*[-1,1]Δu=f,x∈Ω=[−1,1]∗[−1,1] u=g,x∈∂Ωu = g, x \in \partial \Omegau=g,x∈∂Ω 在x轴和y轴上分别取分划长度为dx,dy，令M=int(2/dx)+1,N=int(2/dx)+1M = int(2/dx) + 1,N = int(2/dx) + 1M=int(2/dx)+1,N=int(2/dx)+
一维完全气体的Riemann问题求解
forrestguang的博客
05-04 2217
一维完全气体欧拉方程的引入 在这里我们将利用LF格式，MacCormack格式和Roe格式对一维完全气体欧拉方程进行数值求解。 {(ρρuE)t+(ρuρu2+pu(E+p))x=0,p=(γ−1)(E−12ρu2),γ=1.4.\left\{\begin{array}{c} \left(\begin{array}{c} \rho \\ \rho u \\ E \end{array}\right)_{t}+\left(\begin{array}{c} \rho u \\ \rho u^{2}+p \\ u(
有限元方法求解二维矩形区域椭圆方程
forrestguang的博客
09-06 1852
利用Galerkin方法 def UU(X, order): temp = 10*(X[:,0]+X[:,1])**2 + (X[:,0]-X[:,1])**2 + 0.5 if order[0]==0 and order[1]==0: return torch.log(temp) if order[0]==1 and order[1]==0: return temp**(-1) * (20*(X[:,0]+X[:,1]) + 2*(X[:,0]-X
区域分解与并行计算
forrestguang的博客
05-27 1809
区域分解算法 以下都以泊松方程来举例 问题引入 给定二维空间的泊松方程： −Δu=f,x∈Ω-\Delta u = f,x \in \Omega−Δu=f,x∈Ω u=g,x∈∂Ωu = g,x \in \partial \Omegau=g,x∈∂Ω 我们把区域分成若干部分分别求解，然后相邻区域做一个通信，这样就可以得到一个并行算法。下面我们先把区域分成两部分做一个介绍： 考虑区域Ω1,Ω2\Omega_1,\Omega_2Ω1​,Ω2​ 其中红色的是Ω1\Omega_1Ω1​，蓝色的是Ω2\Omega
二维空间的抛物型方程-含时间项的PDE差分法
forrestguang的博客
09-23 1662
使用迭代法求解 1：## 向前欧拉格式（注意，这种方法不稳定，不可取） ```python import numpy as np import matplotlib.pyplot as plt import torch import time import torch.nn as nn def UU(prob,X,t,ox,ot):#精确解，将空间坐标点排成（M*N，2）的形式计算 if prob == 1: if ox == [0,0] and ot == 0:
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

打赏作者

Galerkin码农选手

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

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

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

打赏作者

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

抵扣说明：

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

余额充值