V2EX = way to explore
V2EX 是一个关于分享和探索的地方
Sign Up Now
For Existing Member  Sign In
V2EX  ›  necomancer  ›  全部回复第 9 页 / 共 32 页
回复总数  635
1 ... 5  6  7  8  9  10  11  12  13  14 ... 32  
2021 年 1 月 10 日
回复了 KIRAYOMATO 创建的主题 Python numpy 算逆矩阵出错是什么情况?
没看出这和摄动定理有啥关系,感觉只是数字太大的精度问题,原则上应该是

if np.linalg.cond(a) < 1/sys.float_info.epsilon:
....ainv = np.linalg.inv(a)
else:
....

你的体系 condition 已经是 1e26 了,float 没法做好的。
condition 巨大目前没有什么好办法,而我大概试了一下,好像转换成 np.longdouble 也不行……然后用 np.finfo 看了一下,好像 np.float 默认是 np.longdouble, np.longfloat....

建议方案:
-2. 如果是厄密矩阵,可以考虑 scipy.sparse.linalg.cg ,emm……最优先建议你重新考虑一下你的问题,重新设计解决方案和找观察量。而且矩阵元差这么大是所有矩阵元都是幂指数的形式么?看看有没有相关的数学性质能简化问题?
-1. 用 preconditioner 然后求解 Ax=I,一般试试 Jacobi preconditioner,我试了,你这个体系好像不太行……SOR 一类的我没试,不知道基于迭代的话效果会不会好一些
0. 试试 moore-penrose pseudoinverse 低精度先试试这个,说不定就好用了,比如你例子里如果
np.linalg.pinv(a, rcond=1e-20).dot(a)
看着效果还行~至少比你给的要好
1. 用别的库,支持更长的 double 的库,可以试试 mpmath,一个 python 库,据说只是速度略拉跨
2. 即便用了高精度的库,遇到大 condition 体系,一般来说也别直接就求逆,试试用 svd 求 u,s,vT,然后 a^-1=vs^{-1}uT,例如
u,s,vt=np.linalg.svd(a)
np.dot(vt.transpose(),np.dot(np.diag(s**-1),u.transpose()))
2021 年 1 月 6 日
回复了 liujing906qd 创建的主题 Android 请问年底了 3000-5000 这个预算推荐买什么手机呢
1+ 8t cyberpunk 2077 edition? 非常炫酷,最近考虑买一个
2020 年 12 月 23 日
回复了 lenkaren 创建的主题 Android 大法的 Xperia5II 是否值得考虑入手
最新用过的是 xperia xz2 港版,系统使用感觉还行,原生比较干净,玩游戏发热也少,视频音效的也都不错。总体还是不错的。
2020 年 12 月 22 日
回复了 James369 创建的主题 机器学习 机器学习的数学基础,包含哪些内容
linear algebra
vector calculus
tensor algebra
measure theory and probability (modern probability theory & random processes)
functional analysis
以及其必备基础知识。推荐 SUMS 系列,先找这些的书,然后看序,作者会说必备基础知识是什么。虽然不做具体研究的话不用太深入,但基础概念和一个很宏观的知识脉络必须有。
2020 年 12 月 7 日
回复了 Hlianbobo 创建的主题 Python 粗看 numpy 和 pandas 后的一点困惑
例如哪些数据处理是 excel 实现不了或者哪些海量数据的批量处理是 sql 或 vba 实现不了的?

完全没有,毕竟 python 和 excel+vba 都是图灵完备的。
2020 年 11 月 28 日
回复了 ericgui 创建的主题 程序员 迫于穷,是用 IDEA 社区版,还是 vscode ?
找个高校工作的对象?一次性解决好多问题……
@potato1949 忘了加狗头→_→
没那么夸张 @oneonesv 按照线性计算应该是~118w 就可以超越 3090 了。
2020 年 11 月 8 日
回复了 Dogod37 创建的主题 程序员 求推荐一个出差用的笔记本
alienware m15 系列
2020 年 11 月 7 日
回复了 villivateur 创建的主题 NAS 不要用普通硬盘做 NAS
@ungrown 是啊……我后来用了硬盘盒,独立供电,快 3 年了现在还是很稳。
2020 年 11 月 7 日
回复了 ericwood067 创建的主题 数学 请教一个复杂点的抛硬币的概率问题
3 和 1 其实很类似,比如你知道了 1 的概率 p,那么 2C50 p^2(1-p)^48
2020 年 11 月 6 日
回复了 sunshinev 创建的主题 程序员 Golang 开发的匿名互动在线聊天游戏,大家多提建议
灰常创意,期待能上线。但建议学习相关法律法规与政策,聊天类管理可严格了。
2020 年 11 月 6 日
回复了 necomancer 创建的主题 程序员 alienware a51m r2 i9 2080super 4k 细节
@20015jjw 移动的好像和台式不一样? 30xx 系列移动版是不是还没出。。。
@tcp 我没做。感觉免费的话折腾起来太麻烦,收费的话有违初衷。
google 的 blogspot 神器
2020 年 11 月 6 日
回复了 necomancer 创建的主题 程序员 alienware a51m r2 i9 2080super 4k 细节
@kile 不太到,光背个电脑电源加上包本身,10-12 公斤
@davidrogers 可以,耗电太丧病了……台式模式的时候我是这么干的
2020 年 11 月 5 日
回复了 necomancer 创建的主题 程序员 alienware a51m r2 i9 2080super 4k 细节
@nightwind 电流的滋滋滋声。
1 ... 5  6  7  8  9  10  11  12  13  14 ... 32  
About   ·   Help   ·   Advertise   ·   Blog   ·   API   ·   FAQ   ·   Solana   ·   2950 Online   Highest 6679   ·     Select Language
创意工作者们的社区
World is powered by solitude
VERSION: 3.9.8.5 · 35ms · UTC 05:55 · PVG 13:55 · LAX 22:55 · JFK 01:55
♥ Do have faith in what you're doing.