使用此标记来解决有关在科学中使用计算机的问题其中许多问题也适用于scicomp.stackexchange.com。
我正在用厄米特（复）正半定矩阵进行一些计算。计算涉及取它们的指数。矩阵 A 的指数（以 2 为底），我将表示...
当 Const=False 和数字时，EXCEL LinEst 函数返回无意义 <1E-7
要重现我的问题，只需输入 A1: 9.1E-8 和 A2: 9.2E-8，然后标记这两个单元格并向下拖动（自动填充）到 A25。然后输入B1：“=1.1*A1”并将公式复制下来，...
我一直在阅读这篇文章《绘制正确分叉图的简单指南》，我想重现下图（图 10，第 2150011-7 页）： 我创建了这个程序，...
我现在正在研究两个多维数组 arr 和 cost。 arr 与尺寸 (width, height, m, n) 是密集的，而 cost 是与尺寸 (width, height, width, height) 稀疏的。 值：宽度和高度是...
V (cm3) T (°C) 1.032 10.0 1.094 29.5 1.156 50.0 1.215 69.5 1.273 90.0 编写一个 python 程序，通过假设求出系数 a 和 b 该液体体积随关系式 V = 1 + a...
我已经开始使用 valgrind，特别是用于代码分析的工具“callgrind”。 首先，我创建了一个用于数值积分的共享库（我已经实现了一个
我是 Rust 新手，我正在尝试使用 Rust 对多个矩阵执行顺序克罗内克积，即我想做 在Python中，我知道我能做到 将 numpy 导入为 np 从 functools 导入减少...
编码二维阵列中多个谐振子的四阶龙格-库塔方法。陷入单一振荡器
我之前在 Python 中使用四阶龙格库塔方法对一般（非线性）摆摆进行了数值求解。这是与此相关的问题，尝试求解第二个非线性摆......
我参考这个问题，已经有很好的答案了；但发现了不必要的操作（请参阅帖子中的讨论），我只是好奇我是否能成功消除这些...
尝试使用四阶龙格-库塔法求解非线性摆二阶微分方程，未得到预期结果
我是Python新手，所以我只将过程限制在matplotlib上，而不是进入NumPy。 我按照 Abhijit Kar Gupta 的书《Python 中的科学计算》，编写了 Python
我有一个多值函数，我想将其与 np.linspace 结合使用。每当函数中有 if/else 条件时，我都会收到值错误。 我期望应用一个函数...
我在使用 Pandas 对一组数据每行求和时遇到了很多麻烦。 我尝试更好地解释。 我有一个数据框（见下文），我想在其中计算“闪烁”事件的数量...
我发现自己经常遇到此类问题。我有一个类似的功能 def 计算（输入）： 结果=two_hour_computation（输入） 结果 = post_processing(结果) 返回结果...
我经常使用 python pandas 来分析实验室的测量数据。例如，如果我有一些“全局”参数“a = 5.2e-3”和“b = 3.21e-5”，这些参数对于...
我想用 sympy 求解具有周期解的三角方程。 例子： 2sin(1.5x) = 1 正弦(1.5x) = 0.5 1.5x = z 正弦(z) = 0.5 阿辛(0.5) = 0.5236 z1 = 0.5236 + 2k*pi (k 是 Z 的元素) (Z...
我必须反复评估以下函数作为 MCMC 链的一部分。它评估加权累积和并消耗大约 80% 的计算时间。 void method1(int n, double* x, double* ...
我正在尝试为以下泊松方程创建 PINN 模型： 这是自定义损失函数： def 自定义损失(y_t,y_p): 将 tf.GradientTape(persistent=True) 用作 g： 与 tf.Gradient...
[Mathematica Stack Exchange 的部分交叉帖子。在这里打开这个，这样我原则上可以在 numpy/scipy 生态系统中征求非数学解决方案。] 我正在尝试确定...
如何使用 OpenMP 并行化 COO 格式的稀疏矩阵-矩阵 (SpMM) 乘积？
我最近使用 COO 存储格式编写了稀疏矩阵向量 (SpMV) 产品的多线程版本： void dcoomv(SparseMatrixCOO *A, 双 *x, 双 *y) { 对于 (int i=0; im;...
Python3.9 中的 ModuleNotFoundError
我正在使用 python 在我大学的计算节点上进行一些数据处理。计算节点安装了多个版本的python。但是，每次我使用 $蟒蛇3