Ask Your Question
0

[CoCalc] [SageManifolds] Error : Black hole rendering with SageMath 2

asked 2023-11-14 16:47:40 +0100

wookhwi gravatar image

updated 2023-11-21 11:21:40 +0100

FrédéricC gravatar image

HI! I, Jaekwook Lee, am a graduate student at Chungbuk National university, Korea, and am currently working on particle trajectory around black holes.

The error source is https://nbviewer.org/github/sagemanif...

The error occur in Cell Number [38] and [39] as below


AttributeError                            Traceback (most recent call last)
Cell In [40], line 1
----> 1 get_ipython().run_cell_magic('cython', '', 'from libc.math cimport cos, sin\ncimport cython\n\n@cython.boundscheck(False)\n@cython.wraparound(False)\ncpdef intersection(list curve, float alpha, float beta, float dmin, float dmax):\n    """\n    Return True if the curve intersect the disk comprised between dmin and dmax\n    tilted of angles alpha and beta\n    """\n    cdef float x, y, z\n    cdef float x2, y2, z2\n    cdef float r, theta, phi\n    cdef int n, i\n    cdef float t\n    cdef float sinalpha, cosalpha\n    cdef float sinbeta, cosbeta\n    sinalpha = sin(alpha)\n    cosalpha = cos(alpha)\n    sinbeta = sin(beta)\n    cosbeta = cos(beta)\n    n = len(curve)\n    r, theta, phi = curve[0][2:5]\n    x, y, z = r*sin(theta)*cos(phi), r*sin(theta)*sin(phi), r*cos(theta)\n    x, y, z = x, y*cos(beta)-z*sin(beta), z*cos(beta)+y*sin(beta)\n    z = z*cos(alpha)+x*sin(alpha)\n    for i in range(1, n):        \n        r = curve[i][2]\n        theta = curve[i][3]\n        phi = curve[i][4]\n        x2, y2, z2 = r*sin(theta)*cos(phi), r*sin(theta)*sin(phi), r*cos(theta)\n        y2, z2 = y2*cosbeta-z2*sinbeta, z2*cosbeta+y2*sinbeta\n        x2, z2 = x2*cosalpha-z2*sinalpha, z2*cosalpha+x2*sinalpha\n        if z!=z2:\n            t = z/(z-z2)\n            if t>=0 and t<1 and curve[i][2]>dmin and curve[i][2]<dmax:\n                return True\n        x, y, z = x2, y2, z2\n    return False\n')
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/IPython/core/interactiveshell.py:2417, in InteractiveShell.run_cell_magic(self, magic_name, line, cell)
   2415 with self.builtin_trap:
   2416     args = (magic_arg_s, cell)
-> 2417     result = fn(*args, **kwargs)
   2418 return result
File /ext/sage/10.1/src/sage/repl/ipython_extension.py:364, in SageMagics.cython(self, line, cell)
    335 """
    336 Cython cell magic
    337 
   (...)
    361     test
    362 """
    363 from sage.misc.cython import cython_compile
--> 364 return cython_compile(cell)
File /ext/sage/10.1/src/sage/misc/cython.py:695, in cython_compile(code, **kwds)
    693 with open(tmpfile, 'w') as f:
    694     f.write(code)
--> 695 return cython_import_all(tmpfile, get_globals(), **kwds)
File /ext/sage/10.1/src/sage/misc/cython.py:585, in cython_import_all(filename, globals, **kwds)
    570 def cython_import_all(filename, globals, **kwds):
    571     """
    572     Imports all non-private (i.e., not beginning with an underscore)
    573     attributes of the specified Cython module into the given context.
   (...)
    583       code
    584     """
--> 585     m = cython_import(filename, **kwds)
    586     for k, x in m.__dict__.items():
    587         if k[0] != '_':
File /ext/sage/10.1/src/sage/misc/cython.py:556, in cython_import(filename, **kwds)
    539 def cython_import(filename, **kwds):
    540     """
    541     Compile a file containing Cython code, then import and return the
    542     module.  Raises an ``ImportError`` if anything goes wrong.
   (...)
    554     - the module that contains the compiled Cython code.
    555     """
--> 556     name, build_dir = cython(filename, **kwds)
    558     oldpath = sys.path
    559     try:
File /ext/sage/10.1/src/sage/misc/cython.py:380, in cython(filename, verbose, compile_message, use_cache, create_local_c_file, annotate, sage_namespace, create_local_so_file)
    378     from sage.misc.package_dir import cython_namespace_package_support
    379     with cython_namespace_package_support():
--> 380         ext, = cythonize([ext],
    381                          aliases=aliases,
    382                          include_path=includes,
    383                          compiler_directives=directives,
    384                          quiet=(verbose <= 0),
    385                          errors_to_stderr=False,
    386                          use_listing_file=True)
    387 finally:
    388     # Read the "listing file" which is the file containing
    389     # warning and error messages generated by Cython.
    390     try:
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/Cython/Build/Dependencies.py:968, in cythonize(module_list, exclude, nthreads, aliases, quiet, force, language, exclude_failures, **options)
    966 c_options = CompilationOptions(**options)
    967 cpp_options = CompilationOptions(**options); cpp_options.cplus = True
--> 968 ctx = c_options.create_context()
    969 options = c_options
    970 module_list, module_metadata = create_extension_list(
    971     module_list,
    972     exclude=exclude,
   (...)
    976     language=language,
    977     aliases=aliases)
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/Cython/Compiler/Main.py:596, in CompilationOptions.create_context(self)
    595 def create_context(self):
--> 596     return Context(self.include_path, self.compiler_directives,
    597                    self.cplus, self.language_level, options=self)
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/Cython/Compiler/Main.py:80, in Context.__init__(self, include_directories, compiler_directives, cpp, language_level, options)
     74 def __init__(self, include_directories, compiler_directives, cpp=False,
     75              language_level=None, options=None):
     76     # cython_scope is a hack, set to False by subclasses, in order to break
     77     # an infinite loop.
     78     # Better code organization would fix it.
---> 80     from . import Builtin, CythonScope
     81     self.modules = {"__builtin__" : Builtin.builtin_scope}
     82     self.cython_scope = CythonScope.create_cython_scope(self)
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/Cython/Compiler/CythonScope.py:5
      3 from .Symtab import ModuleScope
      4 from .PyrexTypes import *
----> 5 from .UtilityCode import CythonUtilityCode
      6 from .Errors import error
      7 from .Scanning import StringSourceDescriptor
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/Cython/Compiler/UtilityCode.py:3
      1 from __future__ import absolute_import
----> 3 from .TreeFragment import parse_from_strings, StringParseContext
      4 from . import Symtab
      5 from . import Naming
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/Cython/Compiler/TreeFragment.py:17
     15 from .Symtab import ModuleScope
     16 from . import PyrexTypes
---> 17 from .Visitor import VisitorTransform
     18 from .Nodes import Node, StatListNode
     19 from .ExprNodes import NameNode
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/Cython/Compiler/Visitor.py:16, in init Cython.Compiler.Visitor()
     14 from . import TypeSlots
     15 from . import Builtin
---> 16 from . import Nodes
     17 from . import ExprNodes
     18 from . import Errors
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/Cython/Compiler/Nodes.py:31
     29 from . import Options
     30 from . import DebugFlags
---> 31 from .Pythran import has_np_pythran, pythran_type, is_pythran_buffer
     32 from ..Utils import add_metaclass
     35 if sys.version_info[0] >= 3:
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/Cython/Compiler/Pythran.py:10
      7 import cython
      9 try:
---> 10     import pythran
     11     pythran_is_pre_0_9 = tuple(map(int, pythran.__version__.split('.')[0:2])) < (0, 9)
     12     pythran_is_pre_0_9_6 = tuple(map(int, pythran.__version__.split('.')[0:3])) < (0, 9, 6)
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/pythran/__init__.py:42
     40 import pythran.log
     41 from pythran.config import get_include
---> 42 from pythran.toolchain import (generate_cxx, compile_cxxfile, compile_cxxcode,
     43                                compile_pythrancode, compile_pythranfile,
     44                                test_compile)
     45 from pythran.spec import spec_parser
     46 from pythran.spec import load_specfile
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/pythran/toolchain.py:12
     10 from pythran.cxxgen import ReturnStatement
     11 from pythran.dist import PythranExtension, PythranBuildExt
---> 12 from pythran.middlend import refine, mark_unexported_functions
     13 from pythran.passmanager import PassManager
     14 from pythran.tables import pythran_ward
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/pythran/middlend.py:4
      1 """This module turns a python AST into an optimized, pythran compatible ast."""
      3 from pythran.analyses import ExtendedSyntaxCheck
----> 4 from pythran.optimizations import (ComprehensionPatterns, ListCompToGenexp,
      5                                    RemoveDeadFunctions)
      6 from pythran.transformations import (ExpandBuiltins, ExpandImports,
      7                                      ExpandImportAll, FalsePolymorphism,
      8                                      NormalizeCompare, NormalizeException,
   (...)
     15                                      NormalizeStaticIf, SplitStaticExpression,
     16                                      RemoveFStrings)
     17 import re
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/pythran/optimizations/__init__.py:21
     19 from .loop_full_unrolling import LoopFullUnrolling
     20 from .modindex import ModIndex
---> 21 from .pattern_transform import PatternTransform
     22 from .range_loop_unfolding import RangeLoopUnfolding
     23 from .range_based_simplify import RangeBasedSimplify
File /ext/sage/10.1/local/var/lib/sage/venv-python3.11.1/lib/python3.11/site-packages/pythran/optimizations/pattern_transform.py:346
    343 pattern_selector = type(known_pattern.pattern).__name__
    344 attr_name = "{}Patterns".format(pattern_selector)
    345 setattr(PatternTransform, attr_name,
--> 346         getattr(PatternTransform, attr_name) + (known_pattern,))
AttributeError: type object 'PatternTransform' has no attribute 'typePatterns'
edit retag flag offensive close merge delete

Comments


TypeError Traceback (most recent call last) File <timed eval="">:1 TypeError: intersection() takes 3 positional arguments but 5 were given

wookhwi gravatar imagewookhwi ( 2023-11-14 16:48:36 +0100 )edit

This looks like a Pythran error on your machine... What is your system ? How did you installed SageMath on it ? (it looks you are using SageMath 10.1).

eric_g gravatar imageeric_g ( 2023-11-14 17:30:34 +0100 )edit

I using SageMath on CoCalc's Jupyter Notebooks..,

https://cocalc.com/features/sage?utm_source=sagemath.org&utm_medium=landingpage (https://cocalc.com/features/sage?utm_...)

wookhwi gravatar imagewookhwi ( 2023-11-14 17:44:18 +0100 )edit

1 Answer

Sort by » oldest newest most voted
0

answered 2023-11-15 15:04:02 +0100

eric_g gravatar image

I've been in touch with the CoCalc staff about this error. It seems to be due to the implementation of Pythran in CoCalc's version of SageMath 10.1 and this might be fixed when SageMath 10.2 (to be released soon) will be installed in CoCalc.

Fortunately, there is an easy workaround: select the kernel "SageMath 10.0" to run the notebook, instead of the default "SageMath 10.1". To do this, it suffices to click on the blue "SageMath 10.1" in the top right menu and select "SageMath 10.0" at the bottom of the page. Then move to the "Kernel" button of the top menu and select "Restart kernel and run all". I've tried this solution on my CoCalc account and it works nicely for this notebook.

Note that the notebook is memory demanding and you might want to first try the lower resolution configuration

n_cpu = 4 # 4 Go Ram minimum
n_geod = 100
nx, ny = 180, 90

at the top of the notebook.

edit flag offensive delete link more

Comments

Thank you for your reply!!

As you answered, I implemented up to cell number [60] in that using Sagemath 10.0.

However, the following error occurs again from [61], is this also a version problem?


TimeoutError Traceback (most recent call last) File /ext/sage/10.0/local/var/lib/sage/venv-python3.11.1/lib/python3.11/urllib/request.py:1348, in AbstractHTTPHandler.do_open(self, http_class, req, **http_conn_args) 1347 try: -> 1348 h.request(req.get_method(), req.selector, req.data, headers, 1349 encode_chunked=req.has_header('Transfer-encoding')) 1350 except OSError as err: # timeout error File /ext/sage/10.0/local/var/lib/sage/venv-python3.11.1/lib/python3.11

wookhwi gravatar imagewookhwi ( 2023-11-18 07:38:26 +0100 )edit

It could be that you have no internet access from your CoCalc project. Have you enable the internet access to this project?

eric_g gravatar imageeric_g ( 2023-11-19 18:45:35 +0100 )edit

First of all, thank you for your kind response to my question!

I tried my CoCalc project with the Internet connected. However, errors continued in the same part. I have posted the entire contents of the error as a new post.

https://ask.sagemath.org/question/74413/cocalcsagemanifolds-error-black-hole-rendering-with-sagemath-3/ (https://ask.sagemath.org/question/744...)

wookhwi gravatar imagewookhwi ( 2023-11-21 07:12:00 +0100 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower

Stats

Asked: 2023-11-14 16:47:40 +0100

Seen: 277 times

Last updated: Nov 15 '23