问题标签 [sage]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
767 浏览

python - sage 中缺少 python 模块

我安装了 Sage 4.7.1 并遇到了一个奇怪的问题。我的许多使用 deepcopy() 和 uniq() 等函数的旧脚本不再将它们识别为全局名称。我已经能够通过一个一个地导入 python 模块来解决这个问题,但这很乏味。但是当我启动命令行 Sage 界面时,我可以在不导入复制模块的情况下键入“list2=deepcopy(list1)”,而且效果很好。命令行 Sage 怎么可能识别全局名称“deepcopy”,但如果我加载使用相同名称的脚本,它就无法识别它?

哎呀,对不起,还不熟悉stackoverflow。我输入:“sage_4.7.1/sage”来启动命令行界面;然后,我输入“load jbom.py”来加载我在python脚本中定义的所有函数。当我使用脚本中的一个函数时,它运行了几秒钟(复杂函数),然后到达我使用 Sage 通常作为全局名称(deepcopy、uniq 等)的一些函数的位置,但由于某种原因,我加载的脚本不知道函数是什么. 重申一下,正如我所描述的,我的脚本 jbom.py 曾经在我上次从事这项特定研究时工作。

如果我使用“load jbom.py”或“import jbom”也没有区别。这两种方法都获得了我在脚本中定义的函数(但在第二种情况下我必须使用 jbom。)并且都得到关于“deepcopy”不是全局名称的相同错误。

回复帝斯曼:我一直草率地描述了这个问题,对此我深表歉意。我创建了一个新脚本“experiment.py”,它的第一行是“import jbom”。执行experiment.py 中的函数可以识别jbom.py 中的函数,但无法识别deepcopy。我尝试将 jbom.py 加载为“加载 jbom.py”,我可以像几个月前一样使用这些功能。那么,这只是在没有正确使用导入/加载等的情况下分层脚本的问题吗?

已解决:我在 jbom.py 的开头添加了“from sage.all import *”,现在我可以加载 Experiment.py 并执行调用 jbom.py 函数的函数而没有任何问题。从关于导入/加载的 Sage 文档中,我无法真正说出我到底做错了什么。

0 投票
0 回答
204 浏览

sage - 使用 Sage 计算大量字段中的最大阶数

我正在 SAGE 中编写一个程序,除其他外,它必须计算“大”数字字段的最大阶数。我正在处理的数字字段不仅具有相当大的度数(到目前为止我必须处理的数字字段约为 40),而且判别式也非常大。
不幸的是,这使得无法使用 SAGE 用于计算最大订单的标准内置函数 - 该命令K.maximal_order()太费时了。让我具体说明以下代码中使用的符号:
Qa12 指定一个带有整数环的数字字段(即最大顺序)OO。需要计算的是OOK的扩展K的最大阶数Qa12。现在,在 MAGMA 中,这可以使用以下代码完成:

这在相当短的时间内运行,即通常不到一分钟。我没有找到直接翻译它的方法,因为 SAGE 中的前两行代码似乎没有类似的东西。很可能是我遗漏了一些东西,所以如果你认为你确实知道翻译前两行的方法,那么请告诉我。如果计算字段的扩展名,则代码应如下所示:

现在,这段代码有两个问题,第一个是所述数字字段的判别式很大,分解它是不可能的。换句话说,我还没有测试过for循环。你有什么方法可以避免处理这些巨大的判别器吗?for 循环是我不安全的第二个地方:因为(再次 - 据我所知)没有直接模拟pMaximalOrder(subOrderK,p)上面的 MAGMA 命令,我希望 SAGE 会“记住”之前的subOrderK值for 循环,因此重新创建了 MAGMA 命令的效果。真的是这样吗?如果注意,您是否找到避免此问题的方法?

注意:
几天前我在 ask.sagemath 上发布了基本相同的问题,但到目前为止无济于事。

0 投票
1 回答
196 浏览

python - 在 Sage 中使用 gnuplot 时出错,但在常规 Python 中工作正常

如果这不是提出这个问题的最佳场所,我深表歉意,但希望这里有人可以提供帮助。我想直接从 Sage 脚本中运行一些 gnuplot 命令,但我收到以下错误消息:

dyld:库未加载:/opt/local/lib/libfreetype.6.dylib 引用自:/opt/local/bin/gnuplot 原因:库版本不兼容:gnuplot 需要 14.0.0 或更高版本,但 libfreetype.6.dylib提供版本 10.0.0

如果我尝试在 Sage 中使用 gnuplotpy 界面,或者我只是使用 Sage 中的 os.system("gnuplot -e \"plot('sin(x)')\"") 之类的东西,则会出现此消息。但是,相同的 os.system(...) 命令在常规 python 中工作得很好。非常感谢。

0 投票
1 回答
547 浏览

sage - 在 Sage 中创建图形

我想在圣人中定义新图。令 V 为有限域 GF(q) 上的向量空间。图的顶点是来自 V 的 i 维子空间和来自 V 的 ni 维子空间,并且当且仅当两个子空间的直接和为 V 时,两个顶点是相邻的。

我无法在 sage 中定义此图。有什么建议吗?

0 投票
1 回答
775 浏览

list - sagemath:将列表中字典中的分数转换为 matplotlib 的小数

我正在尝试自动化圣人数学中一些效率图的过程,并且我有这个公式。Matplotlib.pyplot想要一个值列表。我要么得到一个列表中的字典,如下所述,要么得到一个如下形式的列表:

我想取这些分数,将它们转换为它们的数值近似值,并在 sage math[python] 中创建一个包含这些值的列表。该n()函数会将分数转换为小数,但我不确定如何隔离值以将其提供给n().

0 投票
0 回答
425 浏览

python - 求解包含大量重积分和导数的耦合迭代方程组

我正在尝试解决一个耦合迭代方程系统,每个方程都包含许多积分和导数。

首先,我使用最大值(嵌入 Sage 中)对其进行解析求解,但该解决方案过于依赖于我对未知函数所做的初始猜测,恒定的初始猜测几乎立即产生了响应,而符号函数在用作初始猜测时产生了系统深入计算,有时似乎永无止境。

然而,我用 Sage 尝试的实际上是我原始方程的简化版本,所以我认为可能是我别无选择,只能以数字方式处理积分和导数,但是,我有一些不可忽视的问题:

  1. 积分只允许有数值限制,并且不允许变量作为它们的上限(我认为也许数值方法算法比解析算法更快 - 尽管我在其计算中留下了变量或参数,但它只是没有'不能这样工作)。

  2. 被积函数也不能承认没有被积分的额外变量和参数。

  3. 导数函数本身就是一个很大的障碍,因为我无法计算偏导数或在积分的被积函数中使用导数。

为了摆脱数值导数的所有问题,我用符号 diff() 函数代替它,速度提高仍然有希望,但数值积分的问题仍然存在。

现在我有三个问题:

a-是否可以得出结论,除了离散方程并进行完整的数值处理而不是混合处理之外,我别无他法?

b-如果是这样,那么有什么方法可以自动执行此操作吗?我的方程不是使用 ODEint 的 DE 方程,否则它们是迭代方程,我有积分和导数,只是为了在每一步将我的未知数更新为更新的值。

c-如果我的计算量如此之大,是否有任何关于从 python 切换到 fortran 或类似的东西的建议?

此致

0 投票
1 回答
395 浏览

python - Python/在 Sage Notebook 中创建图表

我正在尝试在 Sage Notebook 中创建散点图。我正在读取 DataA.csv 中的数据,特别是第 1 行和第 3 行。我想根据这些数据创建散点图。一切正常,除了为每个单独的元组创建一个新的散点图,而不是将它们全部放在一个散点图上。关于如何解决此问题的任何想法。提前致谢

0 投票
1 回答
3247 浏览

python - 将 CSV 文件中特定列中的值相加(在 Sage Notebook 中)

我有一个 csv 文件,其中包含 4 个不同列中的大量值。使用 python,有没有办法将一个特定列中的值相加(比如第 1 列中的值)。我想找到一列中所有值的平均值

我可以使用这种方法隔离特定列,但是有没有办法修改它以便能够添加值并找到特定列中的平均值

提前致谢

0 投票
2 回答
492 浏览

math - 消除样本数据中的尖峰

我怎样才能摆脱离散数据集中的火花数据,但以“更顺畅”的方式?

举个例子

在此处输入图像描述

有两个火花,在 20000,但下一个在 600 也被认为是火花。

我已经设法将非常高的值归零,通过

不知道如何处理较低的。最大值应该达到 100,也许 beta 分布的参数需要更多的调整?

目前,它看起来像这样: 在此处输入图像描述

如果可能,请使用鼠尾草作为解释的参考。

0 投票
2 回答
9998 浏览

python - 如何在 sagemath 中添加标题、轴标签和图例?

我目前正在使用 sagemath 托管的在线工作簿制作一些图表。

这是我正在尝试生成图表的一些代码示例:

这是非常基本的——它本质上是两个并列的散点图。

有没有办法轻松添加轴标签、图例和可选的标题?

我设法破解了一个可以让我添加轴标签的解决方案,但它看起来非常丑陋和愚蠢。

有没有更简单的方法来添加轴标签、图例和标题?