When I change optimizer to newton, it fails to run. Are there any other code that need to be changed?

Running Newton optimizer

epoch=00000
residual: eqn:5.6127, imp:0.011106
tmax=1
memory: 390 MiB
walltime: 1.863 s
walltime+callback: 1.864 s
u_00000.png
Traceback (most recent call last):
  File "e:\python\ODIL\odil-main\examples\heat_tmax\heat_tmax.py", line 222, in <module>
    main()
  File "e:\python\ODIL\odil-main\examples\heat_tmax\heat_tmax.py", line 200, in main
    odil.optimize(args, args.optimizer, problem, state, callback)
  File "D:\miniconda\Lib\site-packages\odil\util.py", line 316, in optimize
    return optimize_newton(args, problem, state, callback, **kwargs)
           ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "D:\miniconda\Lib\site-packages\odil\util.py", line 244, in optimize_newton
    vector, matrix = problem.linearize(state)
                     ^^^^^^^^^^^^^^^^^^^^^^^^
  File "D:\miniconda\Lib\site-packages\odil\core.py", line 1321, in linearize
    m = modsp.csr_array(garray)
        ^^^^^^^^^^^^^^^^^^^^^^^
  File "D:\miniconda\Lib\site-packages\scipy\sparse\_compressed.py", line 92, in __init__
    coo = self._coo_container(arg1, dtype=dtype)
          ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "D:\miniconda\Lib\site-packages\scipy\sparse\_coo.py", line 84, in __init__
    self._shape = check_shape(M.shape, allow_1d=is_array)
                  ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
  File "D:\miniconda\Lib\site-packages\scipy\sparse\_sputils.py", line 332, in check_shape
    raise ValueError('shape must be a 1- or 2-tuple of positive '
ValueError: shape must be a 1- or 2-tuple of positive integers

About tutorial of odil #3

Description

Metadata

Assignees

Labels

Type

Projects

Milestone

Relationships

Development

Issue actions