Python +mongodb:如何快速计算大量向量近似度 - V2EX
V2EX = way to explore
V2EX 是一个关于分享和探索的地方
现在注册
已注册用户请  登录
推荐学习书目
Learn Python the Hard Way
Python Sites
PyPI - Python Package Index
http://diveintopython.org/toc/index.html
Pocoo
值得关注的项目
PyPy
Celery
Jinja2
Read the Docs
gevent
pyenv
virtualenv
Stackless Python
Beautiful Soup
结巴中文分词
Green Unicorn
Sentry
Shovel
Pyflakes
pytest
Python 编程
pep8 Checker
Styles
PEP 8
Google Python Style Guide
Code Style from The Hitchhiker's Guide
sunhk25
V2EX    Python

Python +mongodb:如何快速计算大量向量近似度

  •  
  •   sunhk25 2019-07-04 15:34:19 +08:00 3943 次点击
    这是一个创建于 2295 天前的主题,其中的信息可能已经有所发展或是发生改变。
    • 语料库 20 万句的向量( 250 维)值放到 DB 后,python 查询指定文章最相似的 TOP50。
    • 相似度的计算在用 numpy,计算时间是 2+秒,但是从 DB 把数据拿到内存就需要 7 秒。
    • 降维之类的方法不太懂,请教各位正常该如何实现相似检索?
    第 1 条附言    2019-07-05 18:08:17 +08:00
    Redis 的读取和 JSON.load 加起来的时间是 pickle 的几倍,不知道哪里姿势不对,为什么那么多人推荐 Redis。用的是 Python3.6。
    10 条回复    2021-12-23 23:30:49 +08:00
    chengxiao
        1
    chengxiao  
       2019-07-04 16:22:04 +08:00
    250 维....这不是数据库能处理的吧?
    lunaticus7
        2
    lunaticus7  
       2019-07-04 16:33:33 +08:00
    `但是从 DB 把数据拿到内存就需要 7 秒` 是指取 20000 * 250 向量?
    特征向量为什么要存 mangoDB 嘛,这点量也没多少,直接塞内存,精度可以转成 np.float32 甚至 np.float16,能省很多内存,这些精度做检索够用

    解决了了从数据库读特征向量的问题后,就可以优化检索计算了

    大方向就是分级检索:先用低运算量低精度的 metirc (句向量的话直接 cosine 就行)快速检索出一个较大的候选集,然后在候选集内部再用正常的高精度低速 metirc 得出最终结果


    懒得折腾可以直接上 fb 家的 faiss

    PS.你什么算法啊,20w 数据相似度居然要算 2s ?
    sunhk25
        3
    sunhk25  
    OP
       2019-07-04 16:59:09 +08:00
    @lunaticus7
    谢谢意见,7s 是指从数据库拿数据
    现在就是简单的计算内积,20 万的话就要 1G 内存
    250 * 20 * 10000 * 24 / (1024 * 1024) = 1144.4MB
    owenliang
        4
    owenliang  
       2019-07-04 17:58:20 +08:00
    不行就把计算摘出去,换个 c++/golang 做一下,1 秒内肯定没有问题。
    necomancer
        5
    necomancer  
       2019-07-04 22:56:28 +08:00   1
    numpy 的 dot 和 einsum 都比较慢,你可以考虑用 numba 的 guvectorize
    @guvectorize([(float64[:], float64[:], float64[:])],
    '(n),(n)->()', target='parallel')
    def my_inner_prod(a, b, ret):
    tmp1 = tmp2 = tmp3 = 0
    for i in range(a.shape[0]):
    tmp1 += a[i] * b[i]
    tmp2 += a[i] * a[i]
    tmp3 += b[i] * b[i]
    ret[0] = tmp1 / (tmp2 * tmp3) ** 0.5
    这个是 cos(theta),如果不除以模量则只要 tmp1 就可以了
    可以快很多。只要向量维度是对齐的,比如 (100000,250) . (100000,250) -> (100000,)
    或者(100000,250) . (1,250) -> (100000,)
    降维是不是可以考虑 PCA?
    necomancer
        6
    necomancer  
       2019-07-04 22:59:55 +08:00
    In [29]: a = np.random.random((200000,250))

    In [30]: my_inner_prod(a,a)
    Out[30]: array([1., 1., 1., ..., 1., 1., 1.])

    In [31]: %timeit my_inner_prod(a,a)
    49.3 ms ± 850 s per loop (mean ± std. dev. of 7 runs, 10 loops each)

    I7-3687u
    TobiahShaw
        7
    TobiahShaw  
       2019-07-05 12:00:39 +08:00
    想降维的话可以用下 PCA,分析一下主成分试一下,可以提前计算相似度存在 kv 服务,用的时候直接取
    ebingtel
        8
    ebingtel  
       2019-07-05 13:56:17 +08:00
    1、 预计算,不要求实时性的话
    2、用 redis 进行 bit 计算和存储
    sunhk25
        9
    sunhk25  
    OP
       2019-07-05 14:58:53 +08:00
    @ebingtel 正在调查用 memchached 还是 radis 呢。
    Redis 做缓存,mongodb 做持久化,两者结合的话会不会好一点?
    ldyisbes
        10
    ldyisbest  
       2021-12-23 23:30:49 +08:00
    试试向量数据库 milvus
    关于     帮助文档     自助推广系统     博客     API     FAQ     Solana     1169 人在线   最高记录 6679       Select Language
    创意工作者们的社区
    World is powered by solitude
    VERSION: 3.9.8.5 20ms UTC 17:39 PVG 01:39 LAX 10:39 JFK 13:39
    Do have faith in what you're doing.
    ubao msn snddm index pchome yahoo rakuten mypaper meadowduck bidyahoo youbao zxmzxm asda bnvcg cvbfg dfscv mmhjk xxddc yybgb zznbn ccubao uaitu acv GXCV ET GDG YH FG BCVB FJFH CBRE CBC GDG ET54 WRWR RWER WREW WRWER RWER SDG EW SF DSFSF fbbs ubao fhd dfg ewr dg df ewwr ewwr et ruyut utut dfg fgd gdfgt etg dfgt dfgd ert4 gd fgg wr 235 wer3 we vsdf sdf gdf ert xcv sdf rwer hfd dfg cvb rwf afb dfh jgh bmn lgh rty gfds cxv xcv xcs vdas fdf fgd cv sdf tert sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf sdf shasha9178 shasha9178 shasha9178 shasha9178 shasha9178 liflif2 liflif2 liflif2 liflif2 liflif2 liblib3 liblib3 liblib3 liblib3 liblib3 zhazha444 zhazha444 zhazha444 zhazha444 zhazha444 dende5 dende denden denden2 denden21 fenfen9 fenf619 fen619 fenfe9 fe619 sdf sdf sdf sdf sdf zhazh90 zhazh0 zhaa50 zha90 zh590 zho zhoz zhozh zhozho zhozho2 lislis lls95 lili95 lils5 liss9 sdf0ty987 sdft876 sdft9876 sdf09876 sd0t9876 sdf0ty98 sdf0976 sdf0ty986 sdf0ty96 sdf0t76 sdf0876 df0ty98 sf0t876 sd0ty76 sdy76 sdf76 sdf0t76 sdf0ty9 sdf0ty98 sdf0ty987 sdf0ty98 sdf6676 sdf876 sd876 sd876 sdf6 sdf6 sdf9876 sdf0t sdf06 sdf0ty9776 sdf0ty9776 sdf0ty76 sdf8876 sdf0t sd6 sdf06 s688876 sd688 sdf86