当前位置: 首页 >  技术分享 >  MindSponge分子动力学模拟——安装与使用(2023.08)

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

导读:技术背景.昇思MindSpore是由华为主导的一个,面向全场景构建最佳昇腾匹配、支持多处理器架构的开放AI框架。MindSpore不仅仅是软件层面的工具,更重要的是可以协同华为自研的昇腾Ascend平台,做到软硬件一体的行业解决方案。基于MindSpore的高通量、自动微分和自动

技术背景

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

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

免安装MindSponge

MindSponge是基于python和MindSpore框架来做的,因此更新迭代的频率会比较高。如果在本地工作空间用MindSponge源代码来构建一个项目,这样项目内的开发会相对更加的稳定,这是第一种建议的开发模式。首先,如果我们的项目只涉及到MindSponge的使用,用户可以在Gitee上先Fork一个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库上的内容,其主页地址为:https://gitee.com/mindspore/mindscience。具体的安装方法如下,首先我们也需要把这个代码仓clone到本地,这里我们一般使用dev- md分支:

$ 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: [email protected]
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上更加便捷的开发分子动力学模拟算法和应用,并且原生的支持神经网络力场和增强采样,使得我们可以用一个框架就完成绝大部分的分子模拟工作。

版权声明

本文首发链接为:https://www.cnblogs.com/dechinphy/p/ms-system.html

作者ID:DechinPhy

更多原著文章请参考:https://www.cnblogs.com/dechinphy/

打赏专用链接:https://www.cnblogs.com/dechinphy/gallery/image/379634.html

腾讯云专栏同步:https://cloud.tencent.com/developer/column/91958

CSDN同步链接:https://blog.csdn.net/baidu_37157624?spm=1008.2028.3001.5343

51CTO同步链接:https://blog.51cto.com/u_15561675

内容
  • 30个接口自动化测试面试题,赶紧收藏
    30个接口自动化测试面试题,赶紧
    2023-12-05
    1. 什么是接口自动化测试?.答:接口自动化测试是指使用自动化工具对接口进行测试,验证接口的正确性、稳定性和性能等方面的
  • 自动化横行的今天,手工测试员如何杀出一条血路?
    自动化横行的今天,手工测试员如何
    2023-12-05
    自动化测试是每个软件公司反复提及的,放眼望去,测试岗位的招聘要求里十有八九都会有一条“掌握自动化测试技巧”,甚至有的公司
  • OnionArch - 采用DDD+CQRS+.Net 7.0实现的洋葱架构
    OnionArch - 采用DD
    2023-12-04
    博主最近失业在家,找工作之余,看了一些关于洋葱(整洁)架构的资料和项目,有感而发,自己动手写了个洋葱架构解决方案,起名叫
  • 带你了解基于Ploto构建自动驾驶平台
    带你了解基于Ploto构建自动驾
    2023-12-04
    摘要: 华为云Solution as Code推出基于Ploto构建自动驾驶平台解决方案。.本文分享自华为云社区《基于P
  • Cola4.0 - DDD 设计思想
    Cola4.0 - DDD 设计
    2023-12-04
    cola前言.COLA提供了一整套代码架构,拿来即用。 其中包含了很多架构设计思想,包括讨论度很高的领域驱动设计DDD等
  • 手把手教你落地DDD
    手把手教你落地DDD
    2023-12-03
    一、前言 常见的DDD实现架构有很多种,如经典四层架构、六边形(适配器端口)架构、整洁架构(Clean.Architec
  • DDD架构为什么应该首选六边形架构?
    DDD架构为什么应该首选六边形架
    2023-12-03
    一、传统分层架构.分层架构的一个重要原则是:每层只能与位于其下方的层发生耦合。.分层架构分两种:一种是严格分层架构,规定
  • 【实践篇】手把手教你落地DDD
    【实践篇】手把手教你落地DDD
    2023-12-02
    1. 前言.常见的DDD实现架构有很多种,如经典四层架构、六边形(适配器端口)架构、整洁架构(Clean.Archite
  • 【实践篇】领域驱动设计:DDD工程参考架构
    【实践篇】领域驱动设计:DDD工
    2023-12-01
    背景.为什么要制定参考工程架构.不同团队落地DDD所采取的应用架构风格可能不同,并没有统一的、标准的DDD工程架构。有些
  • 一文揭秘DDD到底解决了什么问题
    一文揭秘DDD到底解决了什么问题
    2023-12-01
    DDD作为架构设计思想帮助微服务控制规模复杂度,那它是怎么做到的呢?.一、架构设计是为了解决系统复杂度.谈到架构,相信每
  • 智能车载设备
    智能车载设备
    2023-10-02
    智能车载设备.我们的智能车载设备是一款结合了最新科技和创新设计的汽车配件,旨在提升驾驶体验、提高安全性和为用户带来智能化
  • 智能医疗设备
    智能医疗设备
    2023-10-02
    智能医疗设备产品介绍.智能医疗设备是一款集成了先进技术的高科技产品,旨在提高医疗保健的效率和质量。该设备主要应用于医院、
  • 互联网金融服务平台
    互联网金融服务平台
    2023-10-01
    互联网金融服务平台.产品功能.个人理财:用户可以通过平台进行投资理财,选择适合自己的理财产品,实现资金增值。.贷款服务:
  • 电子元件电感
    电子元件电感
    2023-10-02
    电子元件电感.产品功能.电感是一种重要的电子元件,用于储存和释放电能,调节电路中的电流和电压。它在电子设备和通讯设备中起
  • 智能智能家居设备
    智能智能家居设备
    2023-10-05
    智能家居设备介绍.产品概述.我们当前运营的产品是智能家居设备,它是一款智能化的家居控制系*,旨在为用户提供更便捷、舒适的
  • 智能可穿戴设备
    智能可穿戴设备
    2023-10-03
    产品功能介绍.我们的智能可穿戴设备是一款集健康监测、运动追踪、通讯互动等多种功能于一体的产品。它采用先进的传感技术,可实
  • 智能手机
    智能手机
    2023-10-01
    产品功能介绍:智能手机.智能手机是一款集通讯、娱乐、办公等功能于一体的移动智能设备。首先,智能手机具有强大的通讯功能,支
  • 智能智能物流设备
    智能智能物流设备
    2023-10-04
    智能物流设备.1. 产品描述.智能物流设备是一款基于物联网技术的智能设备,主要用于运输、储存和**物流货物。其核心功能是
  • 智能智能娱乐设备
    智能智能娱乐设备
    2023-10-04
    产品功能介绍.1. 智能娱乐设备.我们的智能娱乐设备是一款结合了智能技术和娱乐功能的产品。它拥有丰富的娱乐资源,包括音乐
  • 智能电视
    智能电视
    2023-10-01
    产品功能介绍.智能电视是一款结合了传统电视和智能硬件的产品。它内置了智能操作系*,能够连接互联网并运行各种应用程序。智能