Sympyで、分数が有理化されるのを回避したい

Question

前提

大学の卒業研究でアプリケーションを開発しています。解析解を求めるためにSympyを利用しているのですが、数値が簡単化されて表示されてしまい、困っています。

実現したいこと

ここに実現したいことを箇条書きで書いてください。

• 1/√2を表示したい場合にsqrt(2)/2と簡単化されてしまうのを回避したい

発生している問題・エラーメッセージ

エラーは表示されないが、計算結果が有理化されてしまう。

sqrt(2)/2

該当のソースコード

print(1/sympy.sqrt(2))

試したこと

インターネットで調べたところ、Sympyではこのような場合に必ず有理化されてしまうという書き込みが散見されました。（そんなはずはないと思うのですが......）
もしSympy以外の方法で、1/sqrt(2)を有利化せずに表示できるようでしたら教えていただきたいです。

補足情報（FW/ツールのバージョンなど）

Python 3.9.12

Answer 1

https://stackoverflow.com/questions/48839611/disabling-automatic-simplification-in-sympy#answer-48847102

A simpler way to disable automatic evaluation is to use context manager evaluate.

ただ

However, as the docstring of the context manager says, most of SymPy code expects automatic evaluation.

とあるので限定された用途しかなさそうですが。

結果例

>>> with sympy.evaluate(False):
...     print(1/sympy.sqrt(2))
1/(sqrt(2))

---

help(sympy.evaluate)の結果↓

Help on function evaluate in module sympy.core.parameters:

evaluate(x)
    Control automatic evaluation

    Explanation
    ===========

    This context manager controls whether or not all SymPy functions evaluate
    by default.

    Note that much of SymPy expects evaluated expressions.  This functionality
    is experimental and is unlikely to function as intended on large
    expressions.

    Examples
    ========

    >>> from sympy import evaluate
    >>> from sympy.abc import x
    >>> print(x + x)
    2*x
    >>> with evaluate(False):
    ...     print(x + x)
    x + x

---

>>> with sympy.evaluate(False):
...     x = 1/sympy.sqrt(2)
>>> x == 1/sympy.sqrt(2)
False

>>> x.equals(1/sympy.sqrt(2))
True

ですよねー、という感じがしますね。

Comments

[kenta_matsuda]

ご回答ありがとうございます。
なるほど、自動評価を回避するためにはevaluateをFalseにすればよかったのですね。
ただ、計算では自動評価された状態を用いるようなので、最終的な計算結果を1/sqrt(2)のような形で出力することは難しそうに思いました。（1/sqrt(2)*([1,0],[0,1])*([1],[0])=sqrt(2)/2*([1],[0])となり、1/sqrt(2)*([1],[0])にならないなど。）無難に自動評価の無効化を行わずに計算したいと思います。
お忙しいところご教示頂き、ありがとうございました。

Answer 2

N()あるいは.evalf()を使ってください。
参考：Numerical Evaluation

import sympy

expr = 1/sympy.sqrt(2)
print(expr) # sqrt(2)/2
print(sympy.N(expr))# 0.707106781186548
print(expr.evalf()) # 0.707106781186548

Comments

[kenta_matsuda]

ご回答頂きありがとうございます。参考ドキュメントの方も拝見しました。
この方法ですと計算結果が数値（0.707106781186548）で出力されてしまうのですが、これを有理化されていない分数の形（1/sqrt(2)）で出力する方法はありますか？