ERROR: Running case-1.jl

[MohamedAbdrabou2025]

I am trying to run the scripts in this project and I got these error

Activating project at ~/Downloads/TL-FWPI.jl-0.0.3
┌ Warning: The active manifest file has dependencies that were resolved with a different julia version (1.8.5). Unexpected behavior may occur.
└ @ ~/Downloads/TL-FWPI.jl-0.0.3/Manifest.toml:0
┌ Warning: The project dependencies or compat requirements have changed since the manifest was last resolved.
│ It is recommended to Pkg.resolve() or consider Pkg.update() if necessary.
└ @ Pkg.API ~/.julia/juliaup/julia-1.11.6+0.aarch64.apple.darwin14/share/julia/stdlib/v1.11/Pkg/src/API.jl:1206
Precompiling project...
✗ JLD2
✗ PyCall
✗ GPUCompiler
✗ Zygote
✗ PyPlot
✗ AzureClusterlessHPC
✗ JOLI
✗ CUDA
✗ SetIntersectionProjection
✗ JUDI
✗ cuDNN
✗ MetaGraphs
✗ DrWatson
✗ Jutul
✗ NNlibCUDA
✗ Flux
✗ JutulDarcy
✗ InvertibleNetworks
✗ ImageSegmentation
✗ JutulDarcyRules
✗ Images
✗ FNO4CO2
0 dependencies successfully precompiled in 42 seconds. 324 already precompiled.
22 dependencies errored.
For a report of the errors see julia> err. To retry use pkg> precompile
Precompiling PyPlot...
Info Given PyPlot was explicitly requested, output will be shown live
ERROR: LoadError: ArgumentError: Source type for reinterpret must be isbits
Stacktrace:
[1] _reinterpret
@ ./reinterpretarray.jl:856 [inlined]
[2] reinterpret(::Type{UInt8}, x::MemoryRef{Int32})
@ Base ./essentials.jl:733
[3] _preserveas!(dest::Vector{UInt8}, ::Type{Cwstring}, x::String)
@ PyCall ~/.julia/packages/PyCall/twYvK/deps/depsutils.jl:67
[4] Py_SetPythonHome(libpy::Ptr{Nothing}, pyversion::VersionNumber, PYTHONHOME::String)
@ PyCall ~/.julia/packages/PyCall/twYvK/deps/depsutils.jl:36
[5] top-level scope
@ ~/.julia/packages/PyCall/twYvK/src/startup.jl:57
[6] include(mod::Module, _path::String)
@ Base ./Base.jl:562
[7] include(x::String)
@ PyCall ~/.julia/packages/PyCall/twYvK/src/PyCall.jl:1
[8] top-level scope
@ ~/.julia/packages/PyCall/twYvK/src/PyCall.jl:38
[9] include
@ ./Base.jl:562 [inlined]
[10] include_package_for_output(pkg::Base.PkgId, input::String, depot_path::Vector{String}, dl_load_path::Vector{String}, load_path::Vector{String}, concrete_deps::Vector{Pair{Base.PkgId, UInt128}}, source::String)
@ Base ./loading.jl:2881
[11] top-level scope
@ stdin:6
in expression starting at /Users/mohamedmousa/.julia/packages/PyCall/twYvK/src/startup.jl:41
in expression starting at /Users/mohamedmousa/.julia/packages/PyCall/twYvK/src/PyCall.jl:1
in expression starting at stdin:6
ERROR: LoadError: Failed to precompile PyCall [438e738f-606a-5dbb-bf0a-cddfbfd45ab0] to "/Users/mohamedmousa/.julia/compiled/v1.11/PyCall/jl_hO8SmW".
Stacktrace:
[1] error(s::String)
@ Base ./error.jl:35
[2] compilecache(pkg::Base.PkgId, path::String, internal_stderr::IO, internal_stdout::IO, keep_loaded_modules::Bool; flags::Cmd, cacheflags::Base.CacheFlags, reasons::Dict{String, Int64}, loadable_exts::Nothing)
@ Base ./loading.jl:3174
[3] (::Base.var"#1110#1111"{Base.PkgId})()
@ Base ./loading.jl:2579
[4] mkpidlock(f::Base.var"#1110#1111"{Base.PkgId}, at::String, pid::Int32; kwopts::@kwargs{stale_age::Int64, wait::Bool})
@ FileWatching.Pidfile ~/.julia/juliaup/julia-1.11.6+0.aarch64.apple.darwin14/share/julia/stdlib/v1.11/FileWatching/src/pidfile.jl:95
[5] #mkpidlock#6
@ ~/.julia/juliaup/julia-1.11.6+0.aarch64.apple.darwin14/share/julia/stdlib/v1.11/FileWatching/src/pidfile.jl:90 [inlined]
[6] trymkpidlock(::Function, ::Vararg{Any}; kwargs::@kwargs{stale_age::Int64})
@ FileWatching.Pidfile ~/.julia/juliaup/julia-1.11.6+0.aarch64.apple.darwin14/share/julia/stdlib/v1.11/FileWatching/src/pidfile.jl:116
[7] #invokelatest#2
@ ./essentials.jl:1057 [inlined]
[8] invokelatest
@ ./essentials.jl:1052 [inlined]
[9] maybe_cachefile_lock(f::Base.var"#1110#1111"{Base.PkgId}, pkg::Base.PkgId, srcpath::String; stale_age::Int64)
@ Base ./loading.jl:3698
[10] maybe_cachefile_lock
@ ./loading.jl:3695 [inlined]
[11] _require(pkg::Base.PkgId, env::String)
@ Base ./loading.jl:2565
[12] __require_prelocked(uuidkey::Base.PkgId, env::String)
@ Base ./loading.jl:2388
[13] #invoke_in_world#3
@ ./essentials.jl:1089 [inlined]
[14] invoke_in_world
@ ./essentials.jl:1086 [inlined]
[15] _require_prelocked(uuidkey::Base.PkgId, env::String)
@ Base ./loading.jl:2375
[16] macro expansion
@ ./loading.jl:2314 [inlined]
[17] macro expansion
@ ./lock.jl:273 [inlined]
[18] __require(into::Module, mod::Symbol)
@ Base ./loading.jl:2271
[19] #invoke_in_world#3
@ ./essentials.jl:1089 [inlined]
[20] invoke_in_world
@ ./essentials.jl:1086 [inlined]
[21] require(into::Module, mod::Symbol)
@ Base ./loading.jl:2260
[22] include
@ ./Base.jl:562 [inlined]
[23] include_package_for_output(pkg::Base.PkgId, input::String, depot_path::Vector{String}, dl_load_path::Vector{String}, load_path::Vector{String}, concrete_deps::Vector{Pair{Base.PkgId, UInt128}}, source::Nothing)
@ Base ./loading.jl:2881
[24] top-level scope
@ stdin:6
in expression starting at /Users/mohamedmousa/.julia/packages/PyPlot/NsAWt/src/PyPlot.jl:2
in expression starting at stdin:6
✗ PyCall
✗ PyPlot
4 dependencies successfully precompiled in 11 seconds. 31 already precompiled.

ERROR: LoadError: The following 1 direct dependency failed to precompile:

PyPlot

Failed to precompile PyPlot [d330b81b-6aea-500a-939a-2ce795aea3ee] to "/Users/mohamedmousa/.julia/compiled/v1.11/PyPlot/jl_1wQphk".
ERROR: LoadError: ArgumentError: Source type for reinterpret must be isbits
Stacktrace:
[1] _reinterpret
@ ./reinterpretarray.jl:856 [inlined]
[2] reinterpret(::Type{UInt8}, x::MemoryRef{Int32})
@ Base ./essentials.jl:733
[3] _preserveas!(dest::Vector{UInt8}, ::Type{Cwstring}, x::String)
@ PyCall ~/.julia/packages/PyCall/twYvK/deps/depsutils.jl:67
[4] Py_SetPythonHome(libpy::Ptr{Nothing}, pyversion::VersionNumber, PYTHONHOME::String)
@ PyCall ~/.julia/packages/PyCall/twYvK/deps/depsutils.jl:36
[5] top-level scope
@ ~/.julia/packages/PyCall/twYvK/src/startup.jl:57
[6] include(mod::Module, _path::String)
@ Base ./Base.jl:562
[7] include(x::String)
@ PyCall ~/.julia/packages/PyCall/twYvK/src/PyCall.jl:1
[8] top-level scope
@ ~/.julia/packages/PyCall/twYvK/src/PyCall.jl:38
[9] include
@ ./Base.jl:562 [inlined]
[10] include_package_for_output(pkg::Base.PkgId, input::String, depot_path::Vector{String}, dl_load_path::Vector{String}, load_path::Vector{String}, concrete_deps::Vector{Pair{Base.PkgId, UInt128}}, source::String)
@ Base ./loading.jl:2881
[11] top-level scope
@ stdin:6
in expression starting at /Users/mohamedmousa/.julia/packages/PyCall/twYvK/src/startup.jl:41
in expression starting at /Users/mohamedmousa/.julia/packages/PyCall/twYvK/src/PyCall.jl:1
in expression starting at stdin:6
ERROR: LoadError: Failed to precompile PyCall [438e738f-606a-5dbb-bf0a-cddfbfd45ab0] to "/Users/mohamedmousa/.julia/compiled/v1.11/PyCall/jl_hO8SmW".
Stacktrace:
[1] error(s::String)
@ Base ./error.jl:35
[2] compilecache(pkg::Base.PkgId, path::String, internal_stderr::IO, internal_stdout::IO, keep_loaded_modules::Bool; flags::Cmd, cacheflags::Base.CacheFlags, reasons::Dict{String, Int64}, loadable_exts::Nothing)
@ Base ./loading.jl:3174
[3] (::Base.var"#1110#1111"{Base.PkgId})()
@ Base ./loading.jl:2579
[4] mkpidlock(f::Base.var"#1110#1111"{Base.PkgId}, at::String, pid::Int32; kwopts::@kwargs{stale_age::Int64, wait::Bool})
@ FileWatching.Pidfile ~/.julia/juliaup/julia-1.11.6+0.aarch64.apple.darwin14/share/julia/stdlib/v1.11/FileWatching/src/pidfile.jl:95
[5] #mkpidlock#6
@ ~/.julia/juliaup/julia-1.11.6+0.aarch64.apple.darwin14/share/julia/stdlib/v1.11/FileWatching/src/pidfile.jl:90 [inlined]
[6] trymkpidlock(::Function, ::Vararg{Any}; kwargs::@kwargs{stale_age::Int64})
@ FileWatching.Pidfile ~/.julia/juliaup/julia-1.11.6+0.aarch64.apple.darwin14/share/julia/stdlib/v1.11/FileWatching/src/pidfile.jl:116
[7] #invokelatest#2
@ ./essentials.jl:1057 [inlined]
[8] invokelatest
@ ./essentials.jl:1052 [inlined]
[9] maybe_cachefile_lock(f::Base.var"#1110#1111"{Base.PkgId}, pkg::Base.PkgId, srcpath::String; stale_age::Int64)
@ Base ./loading.jl:3698
[10] maybe_cachefile_lock
@ ./loading.jl:3695 [inlined]
[11] _require(pkg::Base.PkgId, env::String)
@ Base ./loading.jl:2565
[12] __require_prelocked(uuidkey::Base.PkgId, env::String)
@ Base ./loading.jl:2388
[13] #invoke_in_world#3
@ ./essentials.jl:1089 [inlined]
[14] invoke_in_world
@ ./essentials.jl:1086 [inlined]
[15] _require_prelocked(uuidkey::Base.PkgId, env::String)
@ Base ./loading.jl:2375
[16] macro expansion
@ ./loading.jl:2314 [inlined]
[17] macro expansion
@ ./lock.jl:273 [inlined]
[18] __require(into::Module, mod::Symbol)
@ Base ./loading.jl:2271
[19] #invoke_in_world#3
@ ./essentials.jl:1089 [inlined]
[20] invoke_in_world
@ ./essentials.jl:1086 [inlined]
[21] require(into::Module, mod::Symbol)
@ Base ./loading.jl:2260
[22] include
@ ./Base.jl:562 [inlined]
[23] include_package_for_output(pkg::Base.PkgId, input::String, depot_path::Vector{String}, dl_load_path::Vector{String}, load_path::Vector{String}, concrete_deps::Vector{Pair{Base.PkgId, UInt128}}, source::Nothing)
@ Base ./loading.jl:2881
[24] top-level scope
@ stdin:6
in expression starting at /Users/mohamedmousa/.julia/packages/PyPlot/NsAWt/src/PyPlot.jl:2
in expression starting at stdin:
in expression starting at /Users/mohamedmousa/Downloads/TL-FWPI.jl-0.0.3/src/dummy_src_file.jl:171

julia>

[ziyiyin97]

Hi, thank you for raising this issue! I see that you are using a new julia version (1.11.6) and the codebase itself is using an older version of julia (1.8.5) as suggested in the log. Could you try to install julia 1.8.5 and see if the error persists? Also, once you reach a precompilation error, you could also try using Pkg; Pkg.resolve() and see if it resolves the dependency requirement errors.

[MohamedAbdrabou2025]

Ok. Will try that!
How about Python version?

[MohamedAbdrabou2025]

New ERROR:

ERROR: ArgumentError: ref of NULL PyObject
Stacktrace:
[1] _getproperty(o::PyObject, s::String)
@ PyCall ~/.julia/packages/PyCall/1gn3u/src/PyCall.jl:299
[2] __getproperty
@ ~/.julia/packages/PyCall/1gn3u/src/PyCall.jl:312 [inlined]
[3] getproperty
@ ~/.julia/packages/PyCall/1gn3u/src/PyCall.jl:317 [inlined]
[4] devito_interface(modelPy::PyObject, srcGeometry::GeometryIC{Float32}, srcData::Matrix{Float32}, recGeometry::GeometryIC{Float32}, recData::Nothing, dm::Nothing, options::JUDIOptions)
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/python_interface.jl:26
[5] time_modeling(model_full::Model, srcGeometry::GeometryIC{Float32}, srcData::Matrix{Float32}, recGeometry::GeometryIC{Float32}, recData::Nothing, dm::Nothing, op::Symbol, options::JUDIOptions)
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/time_modeling_serial.jl:39
[6] propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:9
[7] run_and_reduce(func::Function, #unused#::Nothing, nsrc::Int64, arg_func::JUDI.var"#212#213"{judiDataSourceModeling{Float32, :forward}, judiVector{Float32, Matrix{Float32}}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:37
[8] multi_src_propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:66
[9] *
@ ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/LinearOperators/operators.jl:172 [inlined]
[10] dummy_JUDI_operation()
@ Main ~/Downloads/TL-FWPI.jl-main/src/dummy_src_file.jl:144
[11] top-level scope
@ ~/Downloads/TL-FWPI.jl-main/scripts/case-1.jl:19

[MohamedAbdrabou2025]

New Error:

Building forward operator
ERROR: PyError ($(Expr(:escape, :(ccall(#= /Users/mohamedmousa/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:43 =# @pysym(:PyObject_Call), PyPtr, (PyPtr, PyPtr, PyPtr), o, pyargsptr, kw))))) <class 'AttributeError'>
AttributeError("'NoneType' object has no attribute 'shape'")
File "/Users/mohamedmousa/.julia/packages/JUDI/Bf5P3/src/pysource/interface.py", line 40, in forward_rec
rec, _, _ = forward(model, src_coords, rec_coords, wavelet, save=False,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.julia/packages/JUDI/Bf5P3/src/pysource/propagators.py", line 32, in forward
op = forward_op(model.physical_parameters, model.is_tti, model.is_viscoacoustic,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.julia/packages/JUDI/Bf5P3/src/pysource/operators.py", line 61, in call
value = self.func(*args, **kw)
^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.julia/packages/JUDI/Bf5P3/src/pysource/operators.py", line 114, in forward_op
op = Operator(pde + dft + g_expr + eq_save,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/operator/operator.py", line 172, in new
op = cls._build(expressions, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/operator/operator.py", line 205, in _build
irs, byproduct = cls._lower(expressions, profiler=profiler, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/operator/operator.py", line 284, in _lower
expressions = cls._lower_exprs(expressions, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/tools/timing.py", line 76, in call
retval = self.func(*args, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/operator/operator.py", line 350, in _lower_exprs
expressions = flatten([i._evaluate(expand=expand) for i in expressions])
^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/operations/interpolators.py", line 83, in _evaluate
return_value = self.operation(**kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/operations/interpolators.py", line 149, in operation
return self.interpolator._inject(expr=self.expr, field=self.field,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/operations/interpolators.py", line 441, in _inject
idx_subs, temps = self._interp_idx(fields, implicit_dims=implicit_dims,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/operations/interpolators.py", line 307, in _interp_idx
idx_subs.update({v: v.subs({k: p for (k, p) in zip(mapper, pos)})
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/sympy/core/basic.py", line 1171, in subs
rv = rv._subs(old, new, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/finite_differences/differentiable.py", line 169, in _subs
return self.func(*args, evaluate=False)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/tools/abc.py", line 160, in _rebuild
return cls(*args, **kwargs)
^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.local/lib/python3.12/site-packages/devito/types/basic.py", line 716, in new
args, kwargs = cls.args_setup(*args, **kwargs)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.julia/packages/JUDI/Bf5P3/src/pysource/sources.py", line 101, in args_setup
kwargs['nt'] = kwargs.get('time').shape[0]
^^^^^^^^^^^^^^^^^^^^^^^^

Stacktrace:
[1] pyerr_check
@ ~/.julia/packages/PyCall/1gn3u/src/exception.jl:75 [inlined]
[2] pyerr_check
@ ~/.julia/packages/PyCall/1gn3u/src/exception.jl:79 [inlined]
[3] _handle_error(msg::String)
@ PyCall ~/.julia/packages/PyCall/1gn3u/src/exception.jl:96
[4] macro expansion
@ ~/.julia/packages/PyCall/1gn3u/src/exception.jl:110 [inlined]
[5] #107
@ ~/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:43 [inlined]
[6] disable_sigint
@ ./c.jl:473 [inlined]
[7] __pycall!
@ ~/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:42 [inlined]
[8] _pycall!(ret::PyCall.PyObject, o::PyCall.PyObject, args::Tuple{PyCall.PyObject, Matrix{Float32}, Matrix{Float32}, Matrix{Float32}}, nargs::Int64, kw::PyCall.PyObject)
@ PyCall ~/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:29
[9] _pycall!(ret::PyCall.PyObject, o::PyCall.PyObject, args::Tuple{PyCall.PyObject, Matrix{Float32}, Matrix{Float32}, Matrix{Float32}}, kwargs::Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol}, NamedTuple{(:space_order, :f0), Tuple{Int64, Float32}}})
@ PyCall ~/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:11
[10] pycall(::PyCall.PyObject, ::Type{PyCall.PyArray}, ::PyCall.PyObject, ::Vararg{Any}; kwargs::Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol}, NamedTuple{(:space_order, :f0), Tuple{Int64, Float32}}})
@ PyCall ~/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:80
[11] wrapcall_data(::PyCall.PyObject, ::PyCall.PyObject, ::Vararg{Any}; kw::Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol}, NamedTuple{(:space_order, :f0), Tuple{Int64, Float32}}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/python_interface.jl:4
[12] devito_interface(modelPy::PyCall.PyObject, srcGeometry::GeometryIC{Float32}, srcData::Matrix{Float32}, recGeometry::GeometryIC{Float32}, recData::Nothing, dm::Nothing, options::JUDIOptions)
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/python_interface.jl:26
[13] time_modeling(model_full::Model, srcGeometry::GeometryIC{Float32}, srcData::Matrix{Float32}, recGeometry::GeometryIC{Float32}, recData::Nothing, dm::Nothing, op::Symbol, options::JUDIOptions)
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/time_modeling_serial.jl:39
[14] propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:9
[15] run_and_reduce(func::Function, #unused#::Nothing, nsrc::Int64, arg_func::JUDI.var"#212#213"{judiDataSourceModeling{Float32, :forward}, judiVector{Float32, Matrix{Float32}}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:37
[16] multi_src_propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:66
[17] *
@ ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/LinearOperators/operators.jl:172 [inlined]
[18] dummy_JUDI_operation()
@ Main ~/Downloads/TL-FWPI.jl-main/src/dummy_src_file.jl:144
[19] top-level scope
@ ~/Downloads/TL-FWPI.jl-main/scripts/case-1.jl:19

[MohamedAbdrabou2025]

Last Error:

ERROR: ArgumentError: ref of NULL PyObject
Stacktrace:
[1] _getproperty(o::PyObject, s::String)
@ PyCall ~/.julia/packages/PyCall/1gn3u/src/PyCall.jl:299
[2] __getproperty
@ ~/.julia/packages/PyCall/1gn3u/src/PyCall.jl:312 [inlined]
[3] getproperty
@ ~/.julia/packages/PyCall/1gn3u/src/PyCall.jl:317 [inlined]
[4] devito_interface(modelPy::PyObject, srcGeometry::GeometryIC{Float32}, srcData::Matrix{Float32}, recGeometry::GeometryIC{Float32}, recData::Nothing, dm::Nothing, options::JUDIOptions)
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/python_interface.jl:26
[5] time_modeling(model_full::Model, srcGeometry::GeometryIC{Float32}, srcData::Matrix{Float32}, recGeometry::GeometryIC{Float32}, recData::Nothing, dm::Nothing, op::Symbol, options::JUDIOptions)
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/time_modeling_serial.jl:39
[6] propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:9
[7] run_and_reduce(func::Function, #unused#::Nothing, nsrc::Int64, arg_func::JUDI.var"#212#213"{judiDataSourceModeling{Float32, :forward}, judiVector{Float32, Matrix{Float32}}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:37
[8] multi_src_propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:66
[9] *
@ ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/LinearOperators/operators.jl:172 [inlined]
[10] dummy_JUDI_operation()
@ Main ~/Downloads/TL-FWPI.jl-main/src/dummy_src_file.jl:144
[11] top-level scope
@ ~/Downloads/TL-FWPI.jl-main/scripts/case-1.jl:19

[mloubout]

This means JUDI failed to initialize and JUDI.ac is PyNull. So basically your JUDI install is broken for some reason. You can try to re-build it and test it individually to see what the actual error is

[MohamedAbdrabou2025]

Warning: The active manifest file has dependencies that were resolved with a different julia version (1.8.5). Unexpected behavior may occur.
└ @ ~/Downloads/TL-FWPI.jl-main/Manifest.toml:0
✗ GPUCompiler
✗ Zygote
✗ CUDA
✗ cuDNN
✗ NNlibCUDA
✗ Flux
✗ JutulDarcyRules
✗ InvertibleNetworks
✗ FNO4CO2
Precompiling packages finished.
0 dependencies successfully precompiled in 34 seconds. 311 already precompiled.
9 dependencies errored.
For a report of the errors see julia> err. To retry use pkg> precompile
Building forward operator
ERROR: PyCall.PyError("$(Expr(:escape, :(ccall(#= /Users/mohamedmousa/.julia/packages/PyCall/twYvK/src/pyfncall.jl:43 =# @pysym(:PyObject_Call), PyPtr, (PyPtr, PyPtr, PyPtr), o, pyargsptr, kw))))", PyCall.PyObject(Ptr{PyCall.PyObject_struct} @0x0000000139ffeb70), PyCall.PyObject(Ptr{PyCall.PyObject_struct} @0x000000015f73a080), PyCall.PyObject(Ptr{PyCall.PyObject_struct} @0x000000015f758500))
Stacktrace:
[1] pyerr_check
@ ~/.julia/packages/PyCall/twYvK/src/exception.jl:75 [inlined]
[2] pyerr_check
@ ~/.julia/packages/PyCall/twYvK/src/exception.jl:79 [inlined]
[3] _handle_error(msg::String)
@ PyCall ~/.julia/packages/PyCall/twYvK/src/exception.jl:96
[4] macro expansion
@ ~/.julia/packages/PyCall/twYvK/src/exception.jl:110 [inlined]
[5] #107
@ ~/.julia/packages/PyCall/twYvK/src/pyfncall.jl:43 [inlined]
[6] disable_sigint
@ ./c.jl:473 [inlined]
[7] __pycall!
@ ~/.julia/packages/PyCall/twYvK/src/pyfncall.jl:42 [inlined]
[8] _pycall!(ret::PyCall.PyObject, o::PyCall.PyObject, args::Tuple{…}, nargs::Int64, kw::PyCall.PyObject)
@ PyCall ~/.julia/packages/PyCall/twYvK/src/pyfncall.jl:29
[9] _pycall!(ret::PyCall.PyObject, o::PyCall.PyObject, args::Tuple{…}, kwargs::@kwargs{…})
@ PyCall ~/.julia/packages/PyCall/twYvK/src/pyfncall.jl:11
[10] pycall(::PyCall.PyObject, ::Type{…}, ::PyCall.PyObject, ::Vararg{…}; kwargs::@kwargs{…})
@ PyCall ~/.julia/packages/PyCall/twYvK/src/pyfncall.jl:80
[11] wrapcall_data(::PyCall.PyObject, ::PyCall.PyObject, ::Vararg{Any}; kw::@kwargs{space_order::Int64, f0::Float32})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/python_interface.jl:4
[12] devito_interface(modelPy::PyCall.PyObject, srcGeometry::GeometryIC{…}, srcData::Matrix{…}, recGeometry::GeometryIC{…}, recData::Nothing, dm::Nothing, options::JUDIOptions)
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/python_interface.jl:26
[13] time_modeling(model_full::Model, srcGeometry::GeometryIC{…}, srcData::Matrix{…}, recGeometry::GeometryIC{…}, recData::Nothing, dm::Nothing, op::Symbol, options::JUDIOptions)
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/time_modeling_serial.jl:39
[14] propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:9
[15] run_and_reduce(func::Function, ::Nothing, nsrc::Int64, arg_func::JUDI.var"#212#213"{…})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:37
[16] multi_src_propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:66
[17] *
@ ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/LinearOperators/operators.jl:172 [inlined]
[18] dummy_JUDI_operation()
@ Main ~/Downloads/TL-FWPI.jl-main/src/dummy_src_file.jl:144
[19] top-level scope
@ ~/Downloads/TL-FWPI.jl-main/scripts/baseline-FWI.jl:20
Some type information was truncated. Use show(err) to see complete types.

✗ GPUCompiler
✗ Zygote
✗ CUDA
✗ cuDNN
✗ NNlibCUDA
✗ Flux
✗ JutulDarcyRules
✗ InvertibleNetworks
✗ FNO4CO2
Precompiling packages finished.
0 dependencies successfully precompiled in 34 seconds. 311 already precompiled.
9 dependencies errored.
For a report of the errors see julia> err. To retry use pkg> precompile
Building forward operator
ERROR: PyError ($(Expr(:escape, :(ccall(#= /Users/mohamedmousa/.julia/packages/PyCall/twYvK/src/pyfncall.jl:43 =# @pysym(:PyObject_Call), PyPtr, (PyPtr, PyPtr, PyPtr), o, pyargsptr, kw))))) <class 'AttributeError'>
AttributeError("'NoneType' object has no attribute 'shape'")
File "/Users/mohamedmousa/.julia/packages/JUDI/Bf5P3/src/pysource/interface.py", line 40, in forward_rec
rec, _, _ = forward(model, src_coords, rec_coords, wavelet, save=False,
File "/Users/mohamedmousa/.julia/packages/JUDI/Bf5P3/src/pysource/propagators.py", line 32, in forward
op = forward_op(model.physical_parameters, model.is_tti, model.is_viscoacoustic,
File "/Users/mohamedmousa/.julia/packages/JUDI/Bf5P3/src/pysource/operators.py", line 61, in call
value = self.func(*args, **kw)
File "/Users/mohamedmousa/.julia/packages/JUDI/Bf5P3/src/pysource/operators.py", line 114, in forward_op
op = Operator(pde + dft + g_expr + eq_save,
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/operator/operator.py", line 161, in new
op = cls._build(expressions, **kwargs)
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/operator/operator.py", line 183, in _build
irs, byproduct = cls._lower(expressions, profiler=profiler, **kwargs)
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/operator/operator.py", line 257, in _lower
expressions = cls._lower_exprs(expressions, **kwargs)
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/tools/timing.py", line 76, in call
retval = self.func(*args, **kwargs)
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/operator/operator.py", line 325, in _lower_exprs
expressions = flatten([i._evaluate(expand=expand) for i in expressions])
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/operator/operator.py", line 325, in
expressions = flatten([i._evaluate(expand=expand) for i in expressions])
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/operations/interpolators.py", line 61, in _evaluate
return_value = self.operation(**kwargs)
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/operations/interpolators.py", line 122, in operation
return self.interpolator._inject(expr=self.expr, field=self.field,
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/operations/interpolators.py", line 375, in _inject
idx_subs, temps = self._interp_idx(fields, implicit_dims=implicit_dims,
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/operations/interpolators.py", line 247, in _interp_idx
idx_subs.update({v: v.subs({k: p for (k, p) in zip(mapper, pos)})
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/operations/interpolators.py", line 247, in
idx_subs.update({v: v.subs({k: p for (k, p) in zip(mapper, pos)})
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/sympy/core/basic.py", line 1040, in subs
rv = rv._subs(old, new, **kwargs)
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/finite_differences/differentiable.py", line 160, in _subs
return self.func(*args, evaluate=False)
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/tools/abc.py", line 152, in _rebuild
return cls(*args, **kwargs)
File "/Users/mohamedmousa/.julia/conda/3/aarch64/lib/python3.10/site-packages/devito/types/basic.py", line 850, in new
args, kwargs = cls.args_setup(*args, **kwargs)
File "/Users/mohamedmousa/.julia/packages/JUDI/Bf5P3/src/pysource/sources.py", line 101, in args_setup
kwargs['nt'] = kwargs.get('time').shape[0]

Stacktrace:
[1] pyerr_check
@ ~/.julia/packages/PyCall/twYvK/src/exception.jl:75 [inlined]
[2] pyerr_check
@ ~/.julia/packages/PyCall/twYvK/src/exception.jl:79 [inlined]
[3] _handle_error(msg::String)
@ PyCall ~/.julia/packages/PyCall/twYvK/src/exception.jl:96
[4] macro expansion
@ ~/.julia/packages/PyCall/twYvK/src/exception.jl:110 [inlined]
[5] #107
@ ~/.julia/packages/PyCall/twYvK/src/pyfncall.jl:43 [inlined]
[6] disable_sigint
@ ./c.jl:473 [inlined]
[7] __pycall!
@ ~/.julia/packages/PyCall/twYvK/src/pyfncall.jl:42 [inlined]
[8] _pycall!(ret::PyCall.PyObject, o::PyCall.PyObject, args::Tuple{…}, nargs::Int64, kw::PyCall.PyObject)
@ PyCall ~/.julia/packages/PyCall/twYvK/src/pyfncall.jl:29
[9] _pycall!(ret::PyCall.PyObject, o::PyCall.PyObject, args::Tuple{…}, kwargs::@kwargs{…})
@ PyCall ~/.julia/packages/PyCall/twYvK/src/pyfncall.jl:11
[10] pycall(::PyCall.PyObject, ::Type{…}, ::PyCall.PyObject, ::Vararg{…}; kwargs::@kwargs{…})
@ PyCall ~/.julia/packages/PyCall/twYvK/src/pyfncall.jl:80
[11] wrapcall_data(::PyCall.PyObject, ::PyCall.PyObject, ::Vararg{Any}; kw::@kwargs{space_order::Int64, f0::Float32})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/python_interface.jl:4
[12] devito_interface(modelPy::PyCall.PyObject, srcGeometry::GeometryIC{…}, srcData::Matrix{…}, recGeometry::GeometryIC{…}, recData::Nothing, dm::Nothing, options::JUDIOptions)
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/python_interface.jl:26
[13] time_modeling(model_full::Model, srcGeometry::GeometryIC{…}, srcData::Matrix{…}, recGeometry::GeometryIC{…}, recData::Nothing, dm::Nothing, op::Symbol, options::JUDIOptions)
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/time_modeling_serial.jl:39
[14] propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:9
[15] run_and_reduce(func::Function, ::Nothing, nsrc::Int64, arg_func::JUDI.var"#212#213"{…})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:37
[16] multi_src_propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/Modeling/propagation.jl:66
[17] *
@ ~/.julia/packages/JUDI/Bf5P3/src/TimeModeling/LinearOperators/operators.jl:172 [inlined]
[18] dummy_JUDI_operation()
@ Main ~/Downloads/TL-FWPI.jl-main/src/dummy_src_file.jl:144
[19] top-level scope
@ ~/Downloads/TL-FWPI.jl-main/scripts/baseline-FWI.jl:20
Some type information was truncated. Use show(err) to see complete types.

[mloubout]

As I mentionned above, please test JUDI by itself first as this is really hard to follow

[MohamedAbdrabou2025]

Building forward operator
/var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.c:11:10: fatal error: 'omp.h' file not found
11 | #include "omp.h"
| ^~~~~~~
1 error generated.
FAILED compiler invocation: gcc -march=native -O3 -g -fPIC -Wall -std=c99 -Wno-unused-result -Wno-unused-variable -Wno-unused-but-set-variable -ffast-math -shared /var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.c -lm -o /var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.dylib
Building forward operator
/var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.c:11:10: fatal error: 'omp.h' file not found
11 | #include "omp.h"
| ^~~~~~~
1 error generated.
FAILED compiler invocation: gcc -march=native -O3 -g -fPIC -Wall -std=c99 -Wno-unused-result -Wno-unused-variable -Wno-unused-but-set-variable -ffast-math -shared /var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.c -lm -o /var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.dylib
ERROR: PyError ($(Expr(:escape, :(ccall(#= /Users/mohamedmousa/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:43 =# @pysym(:PyObject_Call), PyPtr, (PyPtr, PyPtr, PyPtr), o, pyargsptr, kw))))) <class 'codepy.CompileError'>
CompileError('module compilation failed')
File "/Users/mohamedmousa/.julia/packages/JUDI/aCg1s/src/pysource/interface.py", line 44, in forward_rec
rec, _, I, _ = forward(model, src_coords, rec_coords, wavelet, save=False,
~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
f0=f0, illum=illum, fw=fw)
^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.julia/packages/JUDI/aCg1s/src/pysource/propagators.py", line 38, in forward
op = forward_op(model.physical_parameters, model.is_tti, model.is_viscoacoustic,
model.is_elastic, space_order, fw, model.spacing, save,
t_sub, model.fs, src_coords is not None, rcv_coords is not None,
nfreq(freq_list), dft_sub, ws is not None,
wr is not None, qwf is not None, nv_weights, illum)
File "/Users/mohamedmousa/.julia/packages/JUDI/aCg1s/src/pysource/operators.py", line 66, in call
value = self.func(*args, **kw)
File "/Users/mohamedmousa/.julia/packages/JUDI/aCg1s/src/pysource/operators.py", line 132, in forward_op
op.cfunction
File "/Users/mohamedmousa/.venvs/julia-py/lib/python3.13/site-packages/devito/operator/operator.py", line 829, in cfunction
self._jit_compile()
~~~~~~~~~~~~~~~~~^^
File "/Users/mohamedmousa/.venvs/julia-py/lib/python3.13/site-packages/devito/operator/operator.py", line 815, in _jit_compile
recompiled, src_file = self._compiler.jit_compile(self._soname, str(self))
~~~~~~~~~~~~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.venvs/julia-py/lib/python3.13/site-packages/devito/arch/compiler.py", line 399, in jit_compile
_, _, _, recompiled = compile_from_string(self, target, code, src_file,
~~~~~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
cache_dir=cache_dir, debug=debug,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
sleep_delay=sleep_delay)
^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.venvs/julia-py/lib/python3.13/site-packages/codepy/jit.py", line 439, in compile_from_string
toolchain.build_extension(ext_file, source_paths, debug=debug)
~~~~~~~~~~~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.venvs/julia-py/lib/python3.13/site-packages/codepy/toolchain.py", line 211, in build_extension
raise CompileError("module compilation failed")

Stacktrace:
[1] pyerr_check
@ ~/.julia/packages/PyCall/1gn3u/src/exception.jl:75 [inlined]
[2] pyerr_check
@ ~/.julia/packages/PyCall/1gn3u/src/exception.jl:79 [inlined]
[3] _handle_error(msg::String)
@ PyCall ~/.julia/packages/PyCall/1gn3u/src/exception.jl:96
[4] macro expansion
@ ~/.julia/packages/PyCall/1gn3u/src/exception.jl:110 [inlined]
[5] #107
@ ~/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:43 [inlined]
[6] disable_sigint
@ ./c.jl:473 [inlined]
[7] __pycall!
@ ~/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:42 [inlined]
[8] _pycall!(ret::PyObject, o::PyObject, args::Tuple{PyObject, Matrix{Float32}, Matrix{Float32}, Matrix{Float32}}, nargs::Int64, kw::PyObject)
@ PyCall ~/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:29
[9] _pycall!(ret::PyObject, o::PyObject, args::Tuple{PyObject, Matrix{Float32}, Matrix{Float32}, Matrix{Float32}}, kwargs::Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:fw, :f0, :illum), Tuple{Bool, Float32, Bool}}})
@ PyCall ~/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:11
[10] #pycall#112
@ ~/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:80 [inlined]
[11] #4
@ ~/.julia/packages/JUDI/aCg1s/src/utilities.jl:79 [inlined]
[12] (::JUDI.var"#1#2"{JUDI.var"#4#5"{Tuple{PyArray, PyObject}, Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:fw, :f0, :illum), Tuple{Bool, Float32, Bool}}}, PyObject, Tuple{PyObject, Matrix{Float32}, Matrix{Float32}, Matrix{Float32}}}})()
@ JUDI ~/.julia/packages/JUDI/aCg1s/src/utilities.jl:71
[13] lock(f::JUDI.var"#1#2"{JUDI.var"#4#5"{Tuple{PyArray, PyObject}, Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:fw, :f0, :illum), Tuple{Bool, Float32, Bool}}}, PyObject, Tuple{PyObject, Matrix{Float32}, Matrix{Float32}, Matrix{Float32}}}}, l::ReentrantLock)
@ Base ./lock.jl:185
[14] pylock(f::JUDI.var"#4#5"{Tuple{PyArray, PyObject}, Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:fw, :f0, :illum), Tuple{Bool, Float32, Bool}}}, PyObject, Tuple{PyObject, Matrix{Float32}, Matrix{Float32}, Matrix{Float32}}})
@ JUDI ~/.julia/packages/JUDI/aCg1s/src/utilities.jl:68
[15] rlock_pycall(::PyObject, ::Type{Tuple{PyArray, PyObject}}, ::PyObject, ::Vararg{Any}; kw::Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:fw, :f0, :illum), Tuple{Bool, Float32, Bool}}})
@ JUDI ~/.julia/packages/JUDI/aCg1s/src/utilities.jl:78
[16] wrapcall_data(::PyObject, ::PyObject, ::Vararg{Any}; kw::Base.Pairs{Symbol, Real, Tuple{Symbol, Symbol, Symbol}, NamedTuple{(:fw, :f0, :illum), Tuple{Bool, Float32, Bool}}})
@ JUDI ~/.julia/packages/JUDI/aCg1s/src/TimeModeling/Modeling/python_interface.jl:19
[17] devito_interface(modelPy::PyObject, srcGeometry::GeometryIC{Float32}, srcData::Matrix{Float32}, recGeometry::GeometryIC{Float32}, recData::Nothing, dm::Nothing, options::JUDIOptions, illum::Bool, fw::Bool)
@ JUDI ~/.julia/packages/JUDI/aCg1s/src/TimeModeling/Modeling/python_interface.jl:66
[18] macro expansion
@ ~/.julia/packages/JUDI/aCg1s/src/TimeModeling/Modeling/time_modeling_serial.jl:42 [inlined]
[19] macro expansion
@ ./timing.jl:382 [inlined]
[20] macro expansion
@ ~/.julia/packages/JUDI/aCg1s/src/utilities.jl:31 [inlined]
[21] _time_modeling(model_full::JUDI.IsoModel{Float32, 2}, srcGeometry::GeometryIC{Float32}, srcData::Matrix{Float32}, recGeometry::GeometryIC{Float32}, recData::Nothing, dm::Nothing, op::Symbol, options::JUDIOptions, fw::Bool, illum::Bool)
@ JUDI ~/.julia/packages/JUDI/aCg1s/src/TimeModeling/Modeling/time_modeling_serial.jl:41
[22] (::Base.var"#88#90"{Base.var"#88#89#91"{ExponentialBackOff, Nothing, typeof(JUDI._time_modeling)}})(::JUDI.IsoModel{Float32, 2}, ::Vararg{Any}; kwargs::Base.Pairs{Symbol, Union{}, Tuple{}, NamedTuple{(), Tuple{}}})
@ Base ./error.jl:309
[23] (::Base.var"#88#90"{Base.var"#88#89#91"{ExponentialBackOff, Nothing, typeof(JUDI._time_modeling)}})(::JUDI.IsoModel{Float32, 2}, ::Vararg{Any})
@ Base ./error.jl:291
[24] propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}}, illum::Bool)
@ JUDI ~/.julia/packages/JUDI/aCg1s/src/TimeModeling/Modeling/propagation.jl:9
[25] macro expansion
@ ~/.julia/packages/JUDI/aCg1s/src/TimeModeling/Modeling/propagation.jl:38 [inlined]
[26] macro expansion
@ ./timing.jl:382 [inlined]
[27] macro expansion
@ ~/.julia/packages/JUDI/aCg1s/src/utilities.jl:31 [inlined]
[28] run_and_reduce(func::typeof(JUDI.propagate), ::Nothing, nsrc::Int64, arg_func::JUDI.var"#324#325"{judiDataSourceModeling{Float32, :forward}, judiVector{Float32, Matrix{Float32}}, Bool}; kw::Nothing)
@ JUDI ~/.julia/packages/JUDI/aCg1s/src/TimeModeling/Modeling/propagation.jl:36
[29] run_and_reduce(func::Function, ::Nothing, nsrc::Int64, arg_func::Function)
@ JUDI ~/.julia/packages/JUDI/aCg1s/src/TimeModeling/Modeling/propagation.jl:35
[30] multi_src_propagate(F::judiDataSourceModeling{Float32, :forward}, q::judiVector{Float32, Matrix{Float32}})
@ JUDI ~/.julia/packages/JUDI/aCg1s/src/TimeModeling/Modeling/propagation.jl:79
[31] *
@ ~/.julia/packages/JUDI/aCg1s/src/TimeModeling/LinearOperators/operators.jl:174 [inlined]
[32] dummy_JUDI_operation()
@ Main ~/Downloads/TL-FWPI.jl-main/src/dummy_src_file.jl:144
[33] top-level scope
@ ~/Downloads/TL-FWPI.jl-main/scripts/case-1.jl:19

[mloubout]

Likely your JUDI is at an old version not compatible with latest devito

[MohamedAbdrabou2025]

What is your suggestion to solve this problem? Currently, I use Julia 1.8.5

[mloubout]

Upgrade JUDI.jl in the project. I. you need to install a more recent version

[MohamedAbdrabou2025]

(TL-FWPI) pkg> update JUDI
Updating registry at ~/.julia/registries/General.toml
No Changes to ~/Downloads/TL-FWPI.jl-main/Project.toml
No Changes to ~/Downloads/TL-FWPI.jl-main/Manifest.toml

julia> using Pkg

julia> Pkg.status("JUDI")
Status ~/Downloads/TL-FWPI.jl-main/Project.toml
⌃ [f3b833dc] JUDI v3.4.7
Info Packages marked with ⌃ have new versions available and may be upgradable.

Building forward operator
/var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.c:11:10: fatal error: 'omp.h' file not found
11 | #include "omp.h"
| ^~~~~~~
1 error generated.
FAILED compiler invocation: gcc -march=native -O3 -g -fPIC -Wall -std=c99 -Wno-unused-result -Wno-unused-variable -Wno-unused-but-set-variable -ffast-math -Wno-unused-result -shared -L/opt/homebrew/opt/llvm/lib /var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.c -lm -o /var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.dylib
Building forward operator
/var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.c:11:10: fatal error: 'omp.h' file not found
11 | #include "omp.h"
| ^~~~~~~
1 error generated.
FAILED compiler invocation: gcc -march=native -O3 -g -fPIC -Wall -std=c99 -Wno-unused-result -Wno-unused-variable -Wno-unused-but-set-variable -ffast-math -Wno-unused-result -shared -L/opt/homebrew/opt/llvm/lib /var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.c -lm -o /var/folders/4w/1d1cmn2n49g9x4szgzmqy2fh0000gp/T/devito-jitcache-uid502/ab5027f165a388695401019ccde7627a7b4573dd.dylib
ERROR: PyError ($(Expr(:escape, :(ccall(#= /Users/mohamedmousa/.julia/packages/PyCall/1gn3u/src/pyfncall.jl:43 =# @pysym(:PyObject_Call), PyPtr, (PyPtr, PyPtr, PyPtr), o, pyargsptr, kw))))) <class 'codepy.CompileError'>
CompileError('module compilation failed')
File "/Users/mohamedmousa/.julia/packages/JUDI/aCg1s/src/pysource/interface.py", line 44, in forward_rec
rec, _, I, _ = forward(model, src_coords, rec_coords, wavelet, save=False,
~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
f0=f0, illum=illum, fw=fw)
^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.julia/packages/JUDI/aCg1s/src/pysource/propagators.py", line 38, in forward
op = forward_op(model.physical_parameters, model.is_tti, model.is_viscoacoustic,
model.is_elastic, space_order, fw, model.spacing, save,
t_sub, model.fs, src_coords is not None, rcv_coords is not None,
nfreq(freq_list), dft_sub, ws is not None,
wr is not None, qwf is not None, nv_weights, illum)
File "/Users/mohamedmousa/.julia/packages/JUDI/aCg1s/src/pysource/operators.py", line 66, in call
value = self.func(*args, **kw)
File "/Users/mohamedmousa/.julia/packages/JUDI/aCg1s/src/pysource/operators.py", line 132, in forward_op
op.cfunction
File "/Users/mohamedmousa/.venvs/julia-py/lib/python3.13/site-packages/devito/operator/operator.py", line 829, in cfunction
self._jit_compile()
~~~~~~~~~~~~~~~~~^^
File "/Users/mohamedmousa/.venvs/julia-py/lib/python3.13/site-packages/devito/operator/operator.py", line 815, in _jit_compile
recompiled, src_file = self._compiler.jit_compile(self._soname, str(self))
~~~~~~~~~~~~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.venvs/julia-py/lib/python3.13/site-packages/devito/arch/compiler.py", line 399, in jit_compile
_, _, _, recompiled = compile_from_string(self, target, code, src_file,
~~~~~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
cache_dir=cache_dir, debug=debug,
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
sleep_delay=sleep_delay)
^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.venvs/julia-py/lib/python3.13/site-packages/codepy/jit.py", line 439, in compile_from_string
toolchain.build_extension(ext_file, source_paths, debug=debug)
~~~~~~~~~~~~~~~~~~~~~~~~~^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
File "/Users/mohamedmousa/.venvs/julia-py/lib/python3.13/site-packages/codepy/toolchain.py", line 211, in build_extension
raise CompileError("module compilation failed")

[MohamedAbdrabou2025]

I upgraded JUDI, but I thin there is another error:
a C compilation error from Devito/JUDI: the generated C kernel includes #include <omp.h>, but your toolchain can’t find OpenMP headers/libraries on macOS. On top of that, your Python/Julia stack is x86_64, while your Homebrew path in the compile line is /opt/homebrew/ (that’s arm64 Homebrew). Mixing arm64 Homebrew libs with an x86_64 Python process will cause trouble even after the header is found.

[mloubout]

Please refer to the devito documentation for this error on osx

• you need llvm and libomp (brew install llvm libomp_
• you need to set the path as per what brew says
• set the devito env (arch to clang and platform to m1)