程序员经验分享-hwhale

辛积分器帮助

2024-01-07

我正在尝试开发物理模拟，并且想要实现四阶辛积分 http://en.wikipedia.org/wiki/Symplectic_integration方法。问题是我的数学一定是错误的，因为使用辛积分器时我的模拟根本不起作用（与模拟效果相当好的四阶龙格-库塔积分器相比）。我一直在谷歌上搜索这个问题，我能找到的都是关于这个主题的科学文章。我尝试采用文章中使用的方法，但没有运气。我想知道是否有人有使用辛积分器进行模拟的源代码，最好是模拟引力场，但任何辛积分器都可以。源代码使用什么语言并不重要，但我希望使用 C 风格语法的语言。谢谢！


正如您要求的源代码：来自HERE http://www.unige.ch/~hairer/software.html您可以下载 MATLAB 和 FORTRAN 代码，了解哈密顿系统的辛方法和可逆问题的对称方法。还有许多其他处理微分方程的方法。

And in THIS http://www.unige.ch/~hairer/preprints/gnicodes.ps.gz您可以在论文中找到算法的描述。

Edit

如果您使用数学this http://reference.wolfram.com/mathematica/tutorial/NDSolveSPRK.html也可能有帮助。

本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)

辛积分器帮助 的相关文章

  • AWS ElasticSearch Service - 从 CF 模板设置加密选项

    我正在创建一个云形成模板来在AWS中配置elasticsearch服务域 我想将加密下的此属性设置为 true 域的所有流量都需要 HTTPS 但我无法在 AWS 文档中找到执行此操作的方法 用于设置加密属性的其他选项 例如 启用静态数据加

  • ASP.NET MVC 6 (ASP.NET 5) 中的 Application_PreSendRequestHeaders 和 Application_BeginRequest

    如何在 ASP NET 5 MVC6 中使用这些方法 在 MVC5 中 我在 Global asax 中使用了它 现在呢 也许是入门班 protected void Application PreSendRequestHeaders obj

  • C 中的异或运算符

    在进行按位操作时 我在确定何时使用 XOR 运算符时遇到一些困难 按位与和或非常简单 当您想要屏蔽位时 请使用按位 AND 常见用例是 IP 寻址和子网掩码 当您想要打开位时 请使用包含或 然而 XOR 总是让我明白 我觉得如果在面试中被问

  • 如何在 C++ BOOST 中像图形一样加载 TIFF 图像

    我想要加载一个 tiff 图像 带有带有浮点值的像素的 GEOTIFF 例如 boost C 中的图形 我是 C 的新手 我的目标是使用从源 A 到目标 B 的双向 Dijkstra 来获得更高的性能 Boost GIL load tiif

  • 如何在 Jquery/Javascript 中绑定模糊和更改，但只触发一次函数？

    我试图在选择元素更改时触发函数 由于 Ipad 在 on change 方面遇到问题 我还想绑定到 blur 这在 Ipad 上工作得很好 但是我不希望两个事件都触发该函数两次 所以我需要某种挂钩来确保两个事件是否都触发change and

  • 使用 z = f(x, y) 形式的 B 样条方法来拟合 z = f(x)

    作为一个潜在的解决方案这个问题 https stackoverflow com questions 76476327 how to avoid creating many binary switching variables in gekk

  • Ruby 中的 url_encode

    I read 的文档url encode http rdoc info stdlib erb 1 9 3 ERB Util 3Aurl encode 是否有一个表可以准确地告诉我哪个字符被编码为什么 使用url encode ERB s u

  • jolt变换后json对象的排序

    Input The input json object 所需输出 Event1 Value1 Event2 collection of json objects Event3 The input json object 所以基本上输入 js

  • AWS DynamoDB 写后读一致性 - 理论上它是如何工作的？

    大多数nosql解决方案仅使用最终一致性 并且考虑到DynamoDB将数据复制到三个数据中心 如何保持写后读一致性 解决此类问题的通用方法是什么 我认为这很有趣 因为即使在 MySQL 复制中 数据也是异步复制的 我将详细告诉您 Dynam

  • 张量流中的复杂卷积

    我正在尝试运行一个简单的卷积 但包含复数 r np random random 1 10 10 10 i np random random 1 10 10 10 x tf complex r i conv layer tf layers c

  • CGImage/UIImage 在 UI 线程上延迟加载会导致卡顿

    我的程序显示一个水平滚动表面 从左到右平铺有 UIImageViews 代码在 UI 线程上运行 以确保新可见的 UIImageView 分配有新加载的 UIImage 加载发生在后台线程上 一切工作几乎都很好 除了每个图像变得可见时出现口

  • 使用 AppleScript 运行另一个应用程序而不将其显示在扩展坞上

    使用 AppleScript 您可以创建运行另一个应用程序的脚本 然后将该脚本本身另存为应用程序并将其放置在 Dock 中 问题 不是真正的问题 是 当您单击它时 它仍然会在扩展坞上显示其他应用程序 是否可以阻止其他应用程序在扩展坞中显示

  • 防止索引超出范围错误

    我想编写对某些条件的检查 而不必使用 try catch 并且我想避免出现 Index Out of Range 错误的可能性 if array Element 0 Object Length gt 0 array Element 1 Ob

  • 如何在 PHP 中从字符串类名实例化？ [关闭]

    Closed 这个问题需要细节或清晰度 help closed questions 目前不接受答案 如何创建返回方法名称的新实例 不幸的是我收到这个错误 错误 类名必须是有效的对象或字符串 这是我的代码 class Foo public f

  • 使用随机放置的 NaN 创建示例 numpy 数组

    出于测试目的 我想创建一个M by Nnumpy 数组与c随机放置的 NaN import numpy as np M 10 N 5 c 15 A np random randn M N A mask np nan 我在创建时遇到问题mas

  • 使用 libcurl 检查 SFTP 站点上是否存在文件

    我使用 C 和 libcurl 进行 SFTP FTPS 传输 在上传文件之前 我需要检查文件是否存在而不实际下载它 如果该文件不存在 我会遇到以下问题 set up curlhandle for the public private ke

  • 使用按位运算符相乘

    我想知道如何使用按位运算符将一系列二进制位相乘 但是 我有兴趣这样做来查找二进制值的十进制小数值 这是我正在尝试做的一个例子 假设 1010010 我想使用每个单独的位 以便将其计算为 1 2 1 0 2 2 1 2 3 0 2 4 虽然我

  • 恢复上传文件控制

    我确实阅读了以下帖子 C 暂停 恢复上传 https stackoverflow com questions 1048330 pause resume upload in c 使用 HTTP 恢复上传 https stackoverflow

  • 为什么使用 iPhone 或 iOS 设备在“iframe”中查看“position:fixed”时不起作用？

    我研究过 stackoverflow 似乎position fixed在 iOS 移动设备的 iframe 中 https stackoverflow com questions 15874910 position fixed and if

  • 两种情况或 if 哪个更快？ [关闭]

    Closed 这个问题不符合堆栈溢出指南 help closed questions 目前不接受答案 我必须制作一个 非常 轻的脚本 它将接受用户的选项并调用脚本中的函数来执行一些任务 现在我可以使用 IF 和 CASE 选项 但我想知道两

随机推荐

  • Log4Net停止记录，如何恢复？

    我遇到了 log4net c WinForms 的问题 如果出现短暂故障 例如由于 VMware 快照 网络共享在几秒钟内无法访问 它会停止记录 并且似乎没有记录的方法可以从此类故障中恢复中断 它只是停止记录 直到应用程序重新启动 版本 v

  • 验证地点所有者 Google Places API

    我正在计划一个应用程序 企业所有者可以在其中注册自己的位置 我想使用 Google Places API 来获取企业信息 这样 当新用户在应用程序中注册时 他会提供其位置的 ID 然后应用程序会从 API 检索该信息 我将需要某种方法来验证

  • 如何使用处理和 Twitter4j 从用户时间线获取 20 多个结果？

    我正在尝试制作一个草图 从特定用户的 Twitter 时间轴中获取最后 100 条推文 我正在使用 Twitter4j 草图工作正常 但我了解到 Twitter4J 默认情况下将时间线结果限制为 20 我已经看到此页面来了解如何在使用查询时

  • 从 WSL 选择 Python 解释器

    我使用的是 Windows 10 并且希望将 Python 的默认 VSCode 解释器设置为与我的 WSL 2 Ubuntu 中使用的解释器相同 以便我始终使用 一个 Python 当我单击 选择 Python 解释器 时 会出现查找解释

  • 有没有办法防止回调在破折号中触发？

    我有以下回调 app callback Output useridPicker value Output datePicker date Input url search def update form default url search

  • 使用 Google 位置 API

    请原谅我的无知 但经过几个小时的搜索 我运气不佳 无论如何 我正在尝试编写一个小型桌面应用程序 允许用户输入地址 然后以 GPS 坐标返回其大致位置 据我所知 Google 提供了一个地理编码 API 1 允许以下形式的请求 http ma

  • 在 pandas 数据框中添加一个具有 MultiIndex 列的字段

    我一直在寻找这个问题的答案 因为它看起来很简单 但还没有找到任何东西 如果我错过了什么 抱歉 我有 pandas 版本 0 10 0 并且我一直在尝试以下形式的数据 import pandas import numpy as np impo

  • 如何使用 jQuery 删除 HTML 和特殊字符？

    所以我想要做的是删除所有 HTML 表格标签元素 只留下链接标签 table border 0 cellpadding 1 cellspacing 0 tbody tr valign top td class headlines td td

  • 如何在 C++ 中显示浮点值的多个前导零？ [复制]

    这个问题在这里已经有答案了 在 C 程序中 我想显示一列浮点值 以便符号 数字和小数点全部对齐 必要时 多个前导零应填充每个值的整数部分 例如 A column of floating point values 000 0012 000 0

  • Amazon EC2 Linux 实例中的 GUI

    我想使用 Amazon EC2 但我需要知道是否存在带有 GUI 的 Linux 发行版的 AMI 这适用于 RHEL 6 3 64 位 注意 表示超级用户 可以替换为sudo yum groupinstall Desktop yum in

  • 当选择其他选择下拉菜单上的选项时，jQuery 显示/隐藏一个选择下拉菜单中的选项

    我需要在一个选择下拉列表上显示 隐藏选项 具体取决于另一个选择下拉选项 下面的代码显示了我想要实现的目标 如果 column select 选择菜单选项设置为 1 列 则 layout select 选择菜单必须仅显示 none 选项 如果

  • 使用 System.Uri 附加多个段

    var baseUri new Uri http localhost var uri1 new Uri baseUri 1 var uri2 new Uri uri1 2 不料 uri2 is http localhost 2 http l

  • 汇总空手道测试报告和持续测试

    使用 Karate 运行测试用例后 使用 Surefire 插件发布了一些 html 报告 特别是 我发现每个功能文件都有一个 html 报告 当从自动化管道运行测试时 这很不方便 就像在我的例子中 我使用 htmlpublish Jenk

  • 使用 Jenkins 和多配置复制工件时要聪明

    假设我有一组 虚构的 项目 FOO 和 BAR 这两个项目都有某种多配置选项 FOO 在轴上有一个矩阵X其取值为 x1 xn 所以 FOO 有 n 个版本 BAR 在轴上有一个矩阵Y其取值为 y1 ym 所以有 m 个 BAR 版本 然而

  • knitr mp4 电影嵌入在 Windows XP 上不起作用

    我用一个生成 mp4 电影的块编织了一个 Rmd 文件 到 html r clock fig width 7 fig height 6 fig show animate par mar rep 3 4 for i in seq pi 2 4

  • 使用 Forge 获取 BIM 360 文件的 Revit 版本

    有没有办法在 BIM 360 上查找文件的 Revit 版本而无需下载 rvt 文件 我发现了类似的问题here https stackoverflow com questions 49141340 are any of the forge

  • XGBoost 安装错误

    我在尝试编译 xgboost 时遇到错误 我没有 sudo 访问权限 这让事情变得更加困难 我运行了以下命令 git clone https github com dmlc xgboost git recursive cd xgboost

  • 在模态模式下显示来自其他进程的 WPF 窗口

    我有两个 WPF 应用程序和一个进程管理器 用于将数据从第一个 WPF 应用程序传递到第二个 WPF 应用程序 反之亦然 在一个用例中 我必须以模式模式在第二个应用程序的窗口 主窗口 上显示第一个应用程序的窗口 主窗口 因此 第二个 WPF

  • 魔法宝石：@user.activate！在临时服务器中引发“ActiveRecord::RecordNotSaved”异常，但在我的本地计算机中却没有引发异常

    即使正在发送激活邮件 但也会弹出此错误 完整的跟踪是 activerecord 3 2 2 lib active record persistence rb 104 in save activerecord 3 2 2 lib active

  • 辛积分器帮助

    我正在尝试开发物理模拟 并且想要实现四阶辛积分 http en wikipedia org wiki Symplectic integration方法 问题是我的数学一定是错误的 因为使用辛积分器时我的模拟根本不起作用 与模拟效果相当好的四

热门标签