linear-algebra 相关问题

线性代数是数学的核心体,研究这些空间之间的向量空间和线性映射。主要主题包括:矩阵,向量,线性方程的同时系统,对角线形式,特征值,特征向量,投影,点,线,平面,超平面,旋转和反射。

如何实现点到面距离ICP

我正在实施2D ICP(迭代最近点)。我尝试使用点对点距离但损失很大。 事实上,两个相邻帧之间的丢失是可以接受的,但是,例如,当我

回答 1 投票 0

财富算法特征求解

是否存在 Fortunes eigensolve 算法的实现? 目前我正在尝试将其移植到 python,但到目前为止它没有提供正确的结果。 我唯一发现的是...

回答 1 投票 0

使用 numpy 计算非方阵的逆矩阵

有没有一种方法可以使用 numpy 计算 mxn 非方矩阵的逆矩阵?由于使用 la.inv(S) 似乎给了我一个 ValueError:预期方阵

回答 2 投票 0

如何使用 OpenMP 并行化 COO 格式的稀疏矩阵-矩阵 (SpMM) 乘积?

我最近使用 COO 存储格式编写了稀疏矩阵向量 (SpMV) 产品的多线程版本: void dcoomv(SparseMatrixCOO *A, 双 *x, 双 *y) { 对于 (int i=0; im;...

回答 1 投票 0

pytorch 或 numpy 中的稀疏和巨大矩阵乘法

我在大型稀疏矩阵上执行矩阵乘法时面临性能问题。我有一个场景,我需要将一个小尺寸向量 a 与一个巨大且高度稀疏的向量 b 相乘。

回答 1 投票 0

如何找到复数矩阵 A 使得 A * Transpose(A) = B,其中 B = Transpose(B) 并且是复数?

我想找到一个复值矩阵A,它满足属性A * Transpose(A) = B。矩阵B是复数且对称的(不是埃尔米特矩阵),即B = Transpose(B)。 我在朱莉娅工作。我是

回答 1 投票 0

在 3D 点阵中查找 2D 点阵的拓扑匹配算法

我有一个 3D 晶格,其单位单元(即最小重复单位)为 16 个点。由于它是 3D 晶格,因此它在所有 3 个维度(x、y、z)上都是周期性的。详细来说,晶胞看起来像这样:...

回答 1 投票 0

如何解决 OpenMP 内置向量缩减所遇到的这个问题?

我正在为 CSC 存储格式开发自己的稀疏 BLAS 函数实现。为此,我创建了以下数据结构: typedef 结构 SparseMatrixCSC { 整数米; // 数量...

回答 1 投票 0

如何在Python中高效计算lstsq 10k次?

我目前正在编写代码,将大量粒子(总共 10k)投影到 MLS(移动最小二乘法)表面上。这个任务需要调用 lstsq 函数 10k 次,每个矩阵...

回答 1 投票 0

单位向量旋转矩阵变换后的变形点云

我正在尝试将点云绕轴旋转一定角度而不扭曲其形状。每当我将其旋转非正交角度(任何不是 90 度的倍数)时,它都会变形......

回答 0 投票 0

高阶张量元素求和

我有一个维度为 $I imes I imes \dots imes I$ ($D$ 次)的 $D$ 路张量 $H$,它表示多项式的系数。为了更好的理解,请参考这篇文章:https:...

回答 1 投票 0

使用优化从一个高斯信号中减去另一个高斯信号

所以我有一个关于识别类似信号然后减去它们的问题,以便我们得到信号的基线值而不是某个零或负值。考虑示例图像 b...

回答 0 投票 0

随机选择的 m 个布尔 n 向量可能跨越整个空间的概率

假设 m = Ω(nlog 2n),一组随机选择的 m 个布尔 n 向量可能跨越整个空间的概率是多少? 我知道概率至少为 1/2,然后 m 随机选择布尔 n-vec...

回答 0 投票 0

不正确的矩阵4x4行列式计算

尝试用c++创建基础数学库进行学习,遇到了魔法问题... 通过这段代码做行列式: 行列式 A..D 之前已声明, // 矩阵结构: 结构{ ...

回答 0 投票 0

如何改进子矩阵-向量乘法c代码

我有用于计算通用子矩阵向量(gesubmv)乘法的 C 代码 y(rinds) = A(rinds, cinds)*x(cinds) + y(rinds),其中 rinds 和 cinds 是行和列的向量(不一定

回答 1 投票 0

按值传递列表列表不会在 Scheme 语言中更新

我正在尝试在 Scheme 中编写一个函数,它将使用高斯消除法计算给定矩阵(以列表列表的形式)的秩。 我会解释我想做什么......

回答 1 投票 0

为什么 numpy SVD 在给定一组恰好位于平面内的点时返回非法向量?

我正在尝试使用 numpy 的 SVD 方法 np.linalg.svd() 来拟合一架飞机。作为测试,我将在 R3 中使用两组点。在这两种情况下,所有点在第 3 维中的值为 100。由于设置...

回答 1 投票 0

使用 Fortran 中的 Strumpack 求解稀疏线性系统

我必须解决 CSR 格式的稀疏线性系统。我的代码使用库 Strumpack。 我不明白为什么它不返回解决方案。 CSR格式的矩阵由col_index构成,

回答 0 投票 0

Sympy 求解无法找到简单方程的零空间

我正在尝试使用 pythons 模块 sympy 查找各种函数集的零空间。它设法找到了某些集合的解决方案,例如 {(x - 1)!, x * (x - 2)!, (x - 2)!} 我的代码是 来自 sympy 导入

回答 0 投票 0

求解符号函数的零空间

我有符号函数(例如 x^2、x+1 等)并尝试找到等于零的线性组合。 我使用 Symbolic.jl 编写了一个 Julia 程序。在文档中,它说大多数 Julia m...

回答 1 投票 0

© www.soinside.com 2019 - 2024. All rights reserved.