C++ 中的 N 体模拟具有很大的动量守恒和巨大的能量偏差
作者:Malmel 提问时间:11/12/2023
我正在使用 Verlet 积分(更具体地说,Verlet 积分维基百科页面中“非恒定时差”部分中的最后一个方程)与 C++ 和 Eigen 库一起创建数值 N 体模拟。 我试过测量模拟的动量和能量守...
数值 问答列表
作者:Malmel 提问时间:11/12/2023
我正在使用 Verlet 积分(更具体地说,Verlet 积分维基百科页面中“非恒定时差”部分中的最后一个方程)与 C++ 和 Eigen 库一起创建数值 N 体模拟。 我试过测量模拟的动量和能量守...
作者:zzzbbx 提问时间:3/26/2011
是否可以在 Python 中将数组的元素设置为?NaN 此外,是否可以将变量设置为 +/- 无穷大?如果是这样,是否有任何函数可以检查一个数字是否是无穷大?...
作者:ocramz 提问时间:10/29/2015
介绍 我正在用 ;一些函数可以回调到步进例程、偏微分方程的优化或时间积分。inline-c 特别是在本机 C 语言中,回调可以对连续数组进行操作,通过指针修改它们,并将它们返回到一些不透明的(分布...
作者:Maroon Racoon 提问时间:11/7/2022
我的代码应该计算序列的第 100 个元素$x_0=1 ; x_i=\dfrac{x_{i-1}+1}{x_{i-1}+2}, i=1,2, \ldots$ 我写了迭代和递归函数,但结果不相等。是因为...
作者:tommsch 提问时间:5/22/2023
我正在寻找一种算法来生成随机无幂矩阵。 (我在这里有一个关于检查矩阵是否为无能的类似问题:确定矩阵是否为无能的数值算法) 我已经尝试了以下方法: 排名 1 矩阵 生成两个随机向量 v,w,其中 ,...
作者:KHTN 提问时间:10/28/2023
我有一个Matlab代码代码 model = createpde(); geometryFromEdges(model,@lshapeg); applyBoundaryCondition(mode...
作者:Matt Majic 提问时间:10/2/2023
说我想整合 f=@(x) [1,0]; integral2( @(x,y) f(x)*f(y)' ,0,1,0,1 ); 但这给出了错误 Integrand output size doe...
作者:فراز 提问时间:10/24/2023
我希望在 SciPy 中复制以下 MATLAB 代码的结果。 MATLAB版本: f = @(x, y) besselh(0, 2, x.^2 + y.^2); integral2(f, -0.1...
作者:Ferreira H. S. 提问时间:10/28/2023
我正在尝试评估(使用数值方法求解超越方程)当我改变初始猜测 M0 的值时,迭代次数会变化多少。我想我不应该期望它为如此不同的 M0 值返回相同数量的迭代,所以我认为这段代码可能不会使用我希望它使用的所...
作者:neuro_k 提问时间:10/31/2023
我需要解决以下问题来研究视觉感知。 我将在 4 个位置显示 4 张图像:[L1、L2、L3、L4]。图像将始终从一组 8 个图像中绘制:[1,2,3,4,5,6,7,8]。 有 8-choose-...