MindSponge分子动力学模拟——安装与使用

技术背景

昇思MindSpore是由华为主导的一个,面向全场景构建最佳昇腾匹配、支持多处理器架构的开放AI框架。MindSpore不仅仅是软件层面的工具,更重要的是可以协同华为自研的昇腾Ascend平台,做到软硬件一体的行业解决方案。基于MindSpore的高通量、自动微分和自动并行等高级特性,以及对于python第三方库numpy的友好支持,我们开发了一款可以在MindSpore上进行分子动力学模拟的,模块化、高通量、端到端可微的下一代智能分子模拟程序库——MindSponge

本文我们的重点,是介绍如何安装与使用python语言和mindsponge分子动力学模拟框架。

免安装MindSponge

MindSponge是基于python和MindSpore框架来做的,因此更新迭代的频率会比较高。如果在本地工作空间用MindSponge源代码来构建一个项目,这样项目内的开发会相对更加的稳定,这是第一种建议的开发模式。首先,如果我们的项目只涉及到MindSponge的使用。然后clone到自己本地的一个工作空间:

$ git clone https://gitee.com/helloyesterday/mindsponge.git -b develop
正克隆到 'mindsponge'...
remote: Enumerating objects: 8691, done.
remote: Counting objects: 100% (5565/5565), done.
remote: Compressing objects: 100% (3164/3164), done.
remote: Total 8691 (delta 4144), reused 3107 (delta 2390), pack-reused 3126
接收对象中: 100% (8691/8691), 18.14 MiB | 1007.00 KiB/s, 完成.
处理 delta 中: 100% (6197/6197), 完成.

clone结束之后,会在workspace这个目录下生成一个mindsponge目录:

$ cd mindsponge/
$ ll
总用量 76
drwxrwxr-x  8 dechin dechin  4096 8月  16 10:20 ./
drwxrwxr-x  3 dechin dechin  4096 8月  16 10:20 ../
-rwxrwxr-x  1 dechin dechin  3644 8月  16 10:20 build.sh*
drwxrwxr-x  2 dechin dechin  4096 8月  16 10:20 docs/
drwxrwxr-x  8 dechin dechin  4096 8月  16 10:20 .git/
drwxrwxr-x  2 dechin dechin  4096 8月  16 10:20 .gitee/
-rw-rw-r--  1 dechin dechin  2208 8月  16 10:20 .gitignore
-rw-rw-r--  1 dechin dechin 11357 8月  16 10:20 LICENSE
drwxrwxr-x 14 dechin dechin  4096 8月  16 10:20 mindsponge/
-rw-rw-r--  1 dechin dechin    79 8月  16 10:20 pytest.ini
-rw-rw-r--  1 dechin dechin  2043 8月  16 10:20 README.en.md
-rw-rw-r--  1 dechin dechin  2486 8月  16 10:20 README.md
-rw-rw-r--  1 dechin dechin   128 8月  16 10:20 requirements.txt
-rw-rw-r--  1 dechin dechin  3323 8月  16 10:20 setup.py
drwxrwxr-x  3 dechin dechin  4096 8月  16 10:20 tests/
drwxrwxr-x  5 dechin dechin  4096 8月  16 10:20 tutorials/
-rw-rw-r--  1 dechin dechin     5 8月  16 10:20 version.txt

这样一来,我们就可以在这个路径下进行一些比较灵活的开发,比如新建一个workspace目录,然后在这个路径下开发:

$ mkdir workspace
$ touch water.py

例如,我们可以写一个这样的案例(因为是本地免安装的开发,因此需要用sys.path的方法,把mindsponge目录添加到系统路径中):

from mindspore import context
from mindspore.nn import Adam

import sys
sys.path.insert(0, '..')

from sponge import Sponge, Molecule, ForceField, set_global_units
from sponge.callback import WriteH5MD, RunInfo

if __name__ == "__main__":
    set_global_units('nm', 'kj/mol')
    context.set_context(mode=context.GRAPH_MODE, device_target='GPU')
    system = Molecule(template='water.spce.yaml')

    system.reduplicate([0.3, 0, 0])
    new_sys = system.copy([0, 0, -0.3])
    system.append(new_sys)

    potential = ForceField(system, parameters='SPCE')
    opt = Adam(system.trainable_params(), 1e-3)
    mini = Sponge(system, potential, opt)

    run_info = RunInfo(50)
    cb_h5md = WriteH5MD(system, 'tutorial_c01.h5md', save_freq=50, write_velocity=True, write_force=True)
    mini.run(500, callbacks=[run_info, cb_h5md])

这些代码的具体内容,我们都会在后面的文章中逐一讲解。这个案例运行的结果如下所示:

$ python3 water.py
[MindSPONGE] Started simulation at 2023-08-16 10:35:59
Warning! The optimizer "Adam<>" does not has the attribute "velocity".
[MindSPONGE] Step: 0, E_pot: 110.0423
[MindSPONGE] Step: 50, E_pot: 46.740437
[MindSPONGE] Step: 100, E_pot: -31.139118
[MindSPONGE] Step: 150, E_pot: -61.45553
[MindSPONGE] Step: 200, E_pot: -138.1261
[MindSPONGE] Step: 250, E_pot: -148.34193
[MindSPONGE] Step: 300, E_pot: -148.40369
[MindSPONGE] Step: 350, E_pot: -148.40425
[MindSPONGE] Step: 400, E_pot: -148.40425
[MindSPONGE] Step: 450, E_pot: -148.40456
[MindSPONGE] Finished simulation at 2023-08-16 10:36:05
[MindSPONGE] Simulation time: 5.86 seconds.
--------------------------------------------------------------------------------

简单来说,这个案例就是通过MindSponge构建了4个水分子这样的一个体系,然后用MindSpore内置的Adam优化器,做了一个能量极小化的操作。

编译安装MindSponge

编译构建版本的MindSponge,这里建议使用mindscience库上的内容:

$ git clone https://gitee.com/mindspore/mindscience.git -b dev-md
正克隆到 'mindscience'...
remote: Enumerating objects: 17524, done.
remote: Counting objects: 100% (1093/1093), done.
remote: Compressing objects: 100% (767/767), done.
remote: Total 17524 (delta 506), reused 645 (delta 316), pack-reused 16431
接收对象中: 100% (17524/17524), 322.82 MiB | 2.59 MiB/s, 完成.
处理 delta 中: 100% (10766/10766), 完成.
正在更新文件: 100% (2121/2121), 完成.
$ cd mindscience/MindSPONGE/
$ ll
总用量 100
drwxrwxr-x 6 dechin dechin  4096 8月  16 16:24 ./
drwxrwxr-x 9 dechin dechin  4096 8月  16 16:24 ../
drwxrwxr-x 5 dechin dechin  4096 8月  16 16:24 applications/
-rw-rw-r-- 1 dechin dechin  4949 8月  16 16:24 build.sh
-rw-rw-r-- 1 dechin dechin   793 8月  16 16:24 CMakeLists.txt
drwxrwxr-x 6 dechin dechin  4096 8月  16 16:24 docs/
-rw-rw-r-- 1 dechin dechin 11357 8月  16 16:24 LICENSE
-rw-rw-r-- 1 dechin dechin    69 8月  16 16:24 NOTICE
-rw-rw-r-- 1 dechin dechin   110 8月  16 16:24 OWNERS
-rw-rw-r-- 1 dechin dechin 10972 8月  16 16:24 README_CN.md
-rw-rw-r-- 1 dechin dechin 12128 8月  16 16:24 README.md
-rw-rw-r-- 1 dechin dechin  2384 8月  16 16:24 RELEASE_CN.md
-rw-rw-r-- 1 dechin dechin  2555 8月  16 16:24 RELEASE.md
-rw-rw-r-- 1 dechin dechin   171 8月  16 16:24 requirements.txt
-rw-rw-r-- 1 dechin dechin  3323 8月  16 16:24 setup.py
drwxrwxr-x 5 dechin dechin  4096 8月  16 16:24 src/
drwxrwxr-x 8 dechin dechin  4096 8月  16 16:24 tutorials/
-rw-rw-r-- 1 dechin dechin     6 8月  16 16:24 version.txt

关于具体的编译构建流程,都已经打包成了一个bash脚本,我们直接运行这个bash脚本即可完成编译:

$ bash build.sh -e gpu

构建完成后,会生成一个output目录,其中就有编译好的whl文件:

$ ll output/
总用量 1252
drwxrwxr-x 2 dechin dechin    4096 8月  17 09:10 ./
drwxrwxr-x 8 dechin dechin    4096 8月  17 09:10 ../
-rw-rw-r-- 1 dechin dechin 1267915 8月  17 09:10 mindsponge-1.0.0-py3-none-any.whl
-rw-rw-r-- 1 dechin dechin     100 8月  17 09:10 mindsponge-1.0.0-py3-none-any.whl.sha256

然后就可以进到目录,用pip进行安装:

$ cd output/
$ python3 -m pip install mindsponge-1.0.0-py3-none-any.whl

这里需要对应的调整版本号。安装完成后,可以查看到相应的软件版本信息:

$ python3 -m pip show mindsponge
Name: mindsponge
Version: 1.0.0
Summary: simulation package of next generation molecular modeling in mindspore
Home-page: https://www.mindspore.cn/
Author: The MindSpore Authors
Author-email: contact@mindspore.cn
License: Apache 2.0
Location: /home/dechin/anaconda3/envs/mindspore-latest/lib/python3.7/site-packages
Requires: descriptastorus, scipy, absl-py, pyyaml, dataclasses, biotite, h5py, biopython, numpy, glob2
Required-by: 

这样就基本完成安装了,我们就可以在同一个环境下的任何位置直接调用到mindsponge的包。不过需要注意的是,mindscience仓库和MindSponge仓库这两者的调用方式略有区别,在代码中只是需要调整引用的方式。其实使用编译安装的方法,引用也会更加快捷,不需要使用sys.path来添加路径,具体调整如下所示:

  • 免安装方法
import sys
sys.path.insert(0, 'your/mindsponge/path')
from sponge import Sponge
  • 编译构建方法
from sponge import Sponge

如下所示是一个使用MindSponge模拟一段时间的水溶液中的多肽:

总结概要

本文主要介绍MindSponge框架的免安装使用和编译构建使用这两种使用方法,MindSponge是一款可以在MindSpore上进行分子动力学模拟的,模块化、高通量、端到端可微的下一代智能分子模拟程序库。通过MindSponge,我们可以在GPU上更加便捷的开发分子动力学模拟算法和应用,并且原生的支持神经网络力场和增强采样,使得我们可以用一个框架就完成绝大部分的分子模拟工作。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:/a/597863.html

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

Gin 框架的使用

1、Gin 快速开发 1.1、环境准备 1.1.1、导入 gin 依赖 这里就叫 gin 依赖了&#xff0c;在 Goland 命令行中输入下面的命令&#xff1a; go get -u github.com/gin-gonic/gin 1.1.2、设置代理 如果下载失败&#xff0c;最好设置一下代理&#xff0c;在 cmd 命令行中输入下…

react【实用教程】 搭建开发环境(2024版)Vite+React (官方推荐)

以项目名 reactDemo为例 1. 下载脚手架 在目标文件夹中打开命令行 npm create vite2. 安装项目依赖 cd reactDemo npm i若安装失败&#xff0c;则修改下载源重试 npm config set registry https://registry.npmmirror.com3. 启动项目 npm run dev4. 预览项目 浏览器访问 http…

亚马逊FBA头程多少钱一公斤?FBA头程怎么收费?

在亚马逊的电商生态中&#xff0c;FBA服务已经成为许多卖家提升客户满意度和销售效率的重要工具&#xff0c;然而&#xff0c;对于使用FBA服务的卖家来说&#xff0c;选择一家合适的物流合作伙伴并了解其FBA头程的收费标准和计费方式同样至关重要&#xff0c;亚马逊FBA头程多少…

Elsevier——投稿系统遇到bug时的解决方法

重要&#xff1a;找期刊客服&#xff01;&#xff01;&#xff01; 一、方法&#xff1a; 1. 点击进入与官方客服的对话 2. 按要求输入个人信息 3. 输入遇到的问题 比如&#xff1a; 主题&#xff1a;The Current Status is jammed. 详细描述&#xff1a;The Current State o…

XSS-Labs 靶场通过解析(上)

前言 XSS-Labs靶场是一个专门用于学习和练习跨站脚本攻击&#xff08;XSS&#xff09;技术的在线平台。它提供了一系列的实验场景和演示&#xff0c;帮助安全研究人员、开发人员和安全爱好者深入了解XSS攻击的原理和防御方法。 XSS-Labs靶场的主要特点和功能包括&#xff1a;…

数据结构:线性表(详解)

线性表 线性表的知识框架&#xff1a; 线性表的定义&#xff1a; 线性表是具有相同数据类型的n(n > 0)个数据元素的有限序列&#xff0c;当n 0时线性表为一个空表。 若用L命名为线性表&#xff0c;则数据集合为L {a1,a2,…,an}&#xff0c;其中a1称为表头元素&#xff0c…

【方法】如何创建RAR格式压缩文件?

为了方便存储或者传输文件&#xff0c;我们经常会把文件打包成不同格式的压缩包&#xff0c;那如果想创建的是RAR格式的压缩包&#xff0c;要如何做呢&#xff1f; RAR是WinRAR软件独有的压缩格式&#xff0c;所以我们可以通过WinRAR软件来创建RAR格式压缩包。下面分享两种创建…

02_SpringBoot程序快速启动

目录 打包命令启动启动成功测试结果 打包 点击package打包命令&#xff0c;会生成target目录&#xff0c;目录下会有生成的jar包 命令启动 打开cmd命令窗口&#xff0c;进入子项目的target目录下,输入命令后&#xff0c;回车… java -jar .\note-boot-core-1.0-SNAPSHOT.j…

一起深度学习

CIFAR-10 卷积神经网络 下载数据集构建网络运行测试 下载数据集 batchsz 32cifar_train datasets.CIFAR10(data,trainTrue,transformtorchvision.transforms.Compose([torchvision.transforms.Resize((32,32)),torchvision.transforms.ToTensor()]),downloadTrue)cifar_train …

电脑录屏什么软件好?网友力荐的3款软件!

随着电脑的使用越来越广泛&#xff0c;电脑录屏软件也成为了人们日常生活中经常需要使用到的工具。无论是录制游戏画面、教程演示还是远程教育&#xff0c;一款优秀的电脑录屏软件都能为用户提供极大的帮助&#xff0c;可是电脑录屏什么软件好呢&#xff1f;本文将为大家介绍3款…

图形存储与处理在AI去衣技术中的关键角色

引言&#xff1a; 随着人工智能技术的不断进步&#xff0c;AI去衣技术作为一种颇具争议的应用&#xff0c;已经引起了广泛的关注。该技术依托于深度学习、计算机视觉等领域的先进成果&#xff0c;通过分析图像内容实现对人物衣物的识别和去除。在这一过程中&#xff0c;图形存储…

repo跟git的关系

关于repo 大都讲的太复杂了,大多是从定义角度跟命令角度去讲解,其实从现实项目使用角度而言repo很好理解. 我们都知道git是用来管理项目的,多人开发过程中git功能很好用.现在我们知道一个项目会用一个git仓库去管理,项目的开发过程中会使用git创建分支之类的来更好的维护项目代…

stateflow绝对时间逻辑实操

使用after运算符替换at运算符 如果将at运算符与绝对时间-时间逻辑一起使用,则在尝试模拟模型时会出现错误消息。请改用after运算符。 假设您想使用(5.33,秒)的转换来定义时间延迟。 将转换更改为after(5.33秒),如图所示。这样就不报错了。 使用带有后运算符的外部自循…

【源码+文档+安装教程】校园社团信息管理系统

摘要 随着信息技术在管理上越来越深入而广泛的应用&#xff0c;管理信息系统的实施在技术上已逐步成熟。本文介绍了校园社团信息管理系统的开发全过程。通过分析校园社团信息管理系统管理的不足&#xff0c;创建了一个计算机管理校园社团信息管理系统的方案。文章介绍了校园社团…

【EasySpider】EasySpider+mysql执行配置异常

问题 使用易采集工具操作时候&#xff0c;遇到一个执行异常&#xff0c;后来发现没有选择数据类型 Loading stealth.min.js MySQL config file path: ./mysql_config.json 成功连接到数据库。 Successfully connected to the database. Traceback (most recent call last):…

做了两年数仓,积累的12条SQL调优技巧

本文是作者本人做数仓调优时&#xff0c;所经常使用的SQL调优技巧&#xff0c;这些“技巧”也是经过日常不断摸索、问题排查以及网络检索并且经过本人在线上大规模使用过的&#xff0c;对于下面这12条&#xff08;不算多&#xff0c;但特别有用&#xff09;调优小“技巧”&…

[qnx] 通过zcu104 SD卡更新qnx镜像的步骤

0. 概述 本文演示如果给Xlinx zcu104开发板刷入自定义的qnx镜像 1.将拨码开关设置为SD卡启动 如下图所示&#xff0c;将1拨到On,2,3,4拨到Off&#xff0c;即为通过SD启动。 2.准备SD卡中的内容 首先需要将SD格式化为FAT32的&#xff08;如果已经是FAT32格式&#xff0c;则…

Python高级编程-DJango1

Python高级编程 灵感并不是在逻辑思考的延长线上产生 而是在破除逻辑或常识的地方才有灵感 目录 Python高级编程 1.python学习之前的准备 ​编辑 2.DJango 开发网站 3.创建项目 4.&#xff44;&#xff4a;&#xff41;&#xff4e;&#xff47;项目结构介绍 &#xff11;&…

转行HiL测试工程师

转行没方向&#xff1f;0基础也能转新能源汽车HiL测试岗位&#xff01; 都2024年了&#xff0c;不会还有同学想往软件测试、车载测试方向转吧&#xff01;996、卷经验、卷待遇… ❓❓❓❓想转行没有方向&#xff1f; 建议选择发展前景好的行业&#xff0c;转行前先找好行业&…

目标检测——打架视频数据集

引言 亲爱的读者们&#xff0c;您是否在寻找某个特定的数据集&#xff0c;用于研究或项目实践&#xff1f;欢迎您在评论区留言&#xff0c;或者通过公众号私信告诉我&#xff0c;您想要的数据集的类型主题。小编会竭尽全力为您寻找&#xff0c;并在找到后第一时间与您分享。 …