回答情况：
提问时间：
问题标签：
1回答
我试图解决一个二维传热问题(我是一个MATLAB初学者)。我已经编写了下面的代码，当我运行它时，它不会停止。我试图使用error (err)合并一个收敛条件，但是每当我停止程序时，我总是发现我的错误是150。有人能帮助调整这段代码以获得收敛性检查吗？%===Solution of the problem of 2D in the bookmarks:Case2a=====
clear all;DY=3;
Lx= 60;
浏览 2修改于2017-11-27得票数 1
回答已采纳
我在Numpy/Scipy中四处寻找包含有限差分函数的模块。然而，我发现的最接近的东西是numpy.gradient()，它有利于二阶精度的一阶有限差分，但如果你想要更高的导数或更精确的方法，那就不多了。因此，我的问题是，是否有人知道任何模块( Numpy/Scipy或第三方模块的一部分)专门用于高阶(无论是精度还是导数)有限差分方法。注意，我说的是有限差<em
浏览 0修改于2013-09-24得票数 48
回答已采纳
1回答
我有一个复杂的有限差分模型，它用python编写，使用与下面示例代码相同的一般结构。它有两个for循环，一个用于每个迭代，然后在每个迭代中为x数组的每个位置设置一个循环。
浏览 0修改于2018-05-19得票数 1
回答已采纳
起初，这个问题根本不会优化，然后我意识到旧代码中有一些组件没有提供梯度，所以它们会自动地被有限差分。因此，我在这些组件中添加了fd‘’force_fd‘= True，但它仍然没有优化到正确的值。要使我的新代码优化到与旧的OpenMDAO代码相同的值，唯一的方法是将每个组件设置为有限差分，即使是在提供梯度的组件上。关于旧的和新的OpenMDAO之间的有限差分是如何工作的，我有几个问题：
当旧的OpenMDAO做自动有限</e
浏览 3提问于2016-02-22得票数 0
1回答
我正在调试我的约束随机梯度下降算法，论文建议用有限差分来检验梯度。我增加了一个惩罚函数，但模型不再收敛，所以我想检查我的梯度，如本文所建议的。
举个例子，z。
浏览 3提问于2016-01-12得票数 1
1回答
我想知道是否可以使用SSE (1,2,3,4，...)要优化以下循环：for (int j = l; j < size-1; ++j) for (int k = 1; k < size-1; ++k) v[j*size + k] = (u[j*size + k-1] + u[j*size + k+1]
+ u[(j-1)*size + k]+ u[(j+1)*size + k]) /
浏览 2提问于2010-11-19得票数 3
回答已采纳
3回答
我正在为一个二维PDE问题实现一个有限差分格式。我希望避免使用循环来产生有限的差异。例如，要生成u(x，y)_xx的二阶中心差，我可以将u(x，y)乘以如下：u_xy = (u_{i+1，j+1} + u_{i-1，j-1} - u_{i-1，j+1} - u_{i+1，j
浏览 9修改于2016-03-15得票数 1
回答已采纳
我可以手动编码有限差分格式，但我需要的不仅仅是2阶精度，可能是4阶甚至6阶。随着更高的精度顺序编码，它很快就变得乏味，特别是因为我需要它的数组的不同尺寸。
浏览 1提问于2018-04-06得票数 4
回答已采纳
2回答
我是库达的新手。我试图用Ricky动量的形式来解初值的波动方程。代码的性能为12 GFlops，尽管我的GPU性能为3900。为什么代码对我来说如此无效，我如何才能修复它呢？#include <iostream>#include "step.cu"#include "err.cu"using namespace std;
int main(int argc, char const *a
浏览 6提问于2021-12-23得票数 1
我想用一个隐式变量和一个非线性求解器来迭代它的值，在一个系统中使用有限差分进行导数计算。我知道牛顿解需要解析导数，但是非线性高斯-西德尔似乎也没有正确的迭代。
浏览 0提问于2017-11-09得票数 0
回答已采纳
1回答
嗨，下面我有一个解决非线性耦合PDE问题的代码，但是我需要实现周期性边界条件。周期性边界条件困扰着我，我应该在我的代码中添加什么来执行周期性边界条件？根据下面的模块算法建议进行更新。
其中a，b> 0。我的代码在下面
integer, parameter :: n = 10, A = 20
浏览 0修改于2016-03-10得票数 0
回答已采纳
1回答
下面是有限差分法的一个函数，这是一个非常标准的计算导数的方法，给定一些函数f(x)，一个网格(np.linspace)，以及网格的每一块之间的均匀距离(h)。
浏览 0提问于2018-07-05得票数 0
我正在尝试设计一种用于有限差分法的算法，但我有点困惑。所讨论的常微分方程是y''-5y'+10y = 10x，其中y(0)=0，y(1)=100。
浏览 16修改于2019-04-07得票数 0
2回答
让A是： std::vector<double> values_; A(const std::vector<double> &values) : values_(values){};
values_[i] += df;
virtual method2();}
overrides methods
浏览 3修改于2020-08-26得票数 1
腾讯云开发者
Copyright © 2013 - 2026 Tencent Cloud. All Rights Reserved. 腾讯云 版权所有
深圳市腾讯计算机系统有限公司 ICP备案/许可证号：粤B2-20090059 
粤公网安备44030502008569号
腾讯云计算（北京）有限责任公司 京ICP证150476号 | 京ICP备11018762号