SageMath计算螺旋线坐标第二类Christoffel符号报错排查

问题现象 #

使用SageMath编写代码计算螺旋线坐标的第二类Christoffel符号时，代码执行无报错，但生成的Christoffel对象Ga为空，调用Ga(All,All,All)时触发错误：

AttributeError: 'sage.symbolic.expression.Expression' object has no attribute 'together'

错误根源 #

你犯了两个关键错误：

1. 混用SymPy与SageMath符号系统：导入了SymPy的init_printing，但SageMath自带独立的符号计算框架，其张量类（Coordinates/MetricTensor/Christoffel）仅兼容自身的符号对象，与SymPy对象存在兼容性问题。
2. 雅可比矩阵计算方式错误：你用SymPy的逻辑生成度规矩阵，导致矩阵元素为SymPy表达式，Sage的MetricTensor无法正确解析，最终导致Christoffel符号计算失败，生成空对象。

修正后的代码 #

完全使用SageMath自带工具实现，代码如下：

# 定义Sage原生符号变量
ρ, ψ, ς, a = var('ρ ψ ς a')
# 创建坐标对象
x = Coordinates('x', [ρ, ψ, ς])

# 定义螺旋线坐标到笛卡尔坐标的映射
f = [ρ*cos(ψ + a*ς), ρ*sin(ψ + a*ς), ς]
# 用Sage原生方法计算雅可比矩阵
jac = jacobian(f, [ρ, ψ, ς])
# 计算度规张量（雅可比矩阵转置 × 原矩阵）
Metric = jac.transpose() * jac

# 初始化度规张量对象
g = MetricTensor('g', x, Metric)
# 生成Christoffel符号对象
Ga = Christoffel('Ga', g)
# 输出所有第二类Christoffel符号
Ga(All, All, All)

对比代码正常运行的原因 #

你提供的史瓦西度规代码完全基于SageMath原生工具：

• 符号变量由Sage生成
• 度规矩阵用Sage的diag函数创建，所有元素均为Sage原生符号表达式
• 整个流程与Sage张量类完全兼容，因此能正常计算并输出Christoffel符号

内容的提问来源于stack exchange，提问作者Hans-Friedrich Pfeiffer