most recent 30 from stackoverflow.com 2026-04-24T12:40:38Z https://stackoverflow.com/feeds/tag?tagnames=mex&sort=newest https://creativecommons.org/licenses/by-sa/4.0/rdf https://stackoverflow.com/q/79750248 0 dogAwakeCat https://stackoverflow.com/users/2833518 2025-08-29T11:53:05Z 2025-08-29T13:28:12Z

<p>I can't compile the boost library on Windows 11, specifically serialization which is not header-only, for use with mex in MATLAB. I am compiling boost with <code>.\bootstrap.bat; .\b2.exe link=static runtime-link=static toolset=msvc address-model=64 architecture=x86 threading=multi --build-type=complete -a --reconfigure</code>, however, no library (<code>.lib</code>) files are created.</p> <p>Output (I've cut out bits that are not serialization related):</p> <pre><code>Building Boost.Build engine LOCALAPPDATA=C:\Users\user\AppData\Local Found with vswhere C:\Program Files\Microsoft Visual Studio\2022\Community Found with vswhere C:\Program Files\Microsoft Visual Studio\2022\Community Call_If_Exists &quot;C:\Program Files\Microsoft Visual Studio\2022\Community\Common7\Tools\..\..\VC\Auxiliary\Build\vcvarsall.bat&quot; AMD64 ********************************************************************** ** Visual Studio 2022 Developer Command Prompt v17.14.13 ** Copyright (c) 2025 Microsoft Corporation ********************************************************************** [vcvarsall.bat] Environment initialized for: 'x64' ### ### Using 'vc143' toolset. ### C:\Program Files\boost\boost_1_82_0\tools\build\src\engine&gt;&quot;cl&quot; /nologo -TP /wd4996 /wd4675 /EHs /GR /Zc:throwingNew /O2 /Ob2 /W3 /MD /Zc:forScope /Zc:wchar_t /Zc:inline /Gw /favor:blend /Feb2 -DNDEBUG builtins.cpp class.cpp command.cpp compile.cpp constants.cpp cwd.cpp debug.cpp debugger.cpp execcmd.cpp execnt.cpp execunix.cpp filent.cpp filesys.cpp fileunix.cpp frames.cpp function.cpp glob.cpp hash.cpp hcache.cpp hdrmacro.cpp headers.cpp jam.cpp jamgram.cpp lists.cpp make.cpp make1.cpp md5.cpp mem.cpp modules.cpp native.cpp object.cpp option.cpp output.cpp parse.cpp pathnt.cpp pathsys.cpp pathunix.cpp regexp.cpp rules.cpp scan.cpp search.cpp jam_strings.cpp startup.cpp subst.cpp sysinfo.cpp timestamp.cpp variable.cpp w32_getreg.cpp modules/order.cpp modules/path.cpp modules/property-set.cpp modules/regex.cpp modules/sequence.cpp modules/set.cpp /link kernel32.lib advapi32.lib user32.lib builtins.cpp class.cpp command.cpp compile.cpp constants.cpp cwd.cpp debug.cpp debugger.cpp execcmd.cpp execnt.cpp execunix.cpp filent.cpp filesys.cpp fileunix.cpp frames.cpp function.cpp glob.cpp hash.cpp hcache.cpp hdrmacro.cpp Generating Code... Compiling... headers.cpp jam.cpp jamgram.cpp lists.cpp make.cpp make1.cpp md5.cpp mem.cpp modules.cpp native.cpp object.cpp option.cpp output.cpp parse.cpp pathnt.cpp pathsys.cpp pathunix.cpp regexp.cpp rules.cpp scan.cpp Generating Code... Compiling... search.cpp jam_strings.cpp startup.cpp subst.cpp sysinfo.cpp timestamp.cpp variable.cpp w32_getreg.cpp order.cpp path.cpp property-set.cpp regex.cpp sequence.cpp set.cpp Generating Code... C:\Program Files\boost\boost_1_82_0\tools\build\src\engine&gt;dir *.exe Volume in drive C has no label. Volume Serial Number is 60A0-D31B Directory of C:\Program Files\boost\boost_1_82_0\tools\build\src\engine 08/29/2025 07:38 AM 263,168 b2.exe 1 File(s) 263,168 bytes 0 Dir(s) 1,340,648,009,728 bytes free Generating Boost.Build configuration in project-config.jam for msvc... Bootstrapping is done. To build, run: .\b2 To adjust configuration, edit 'project-config.jam'. Further information: - Command line help: .\b2 --help - Getting started guide: http://boost.org/more/getting_started/windows.html - Boost.Build documentation: http://www.boost.org/build/ Performing configuration checks - default address-model : 32-bit [1] - default architecture : arm [1] Building the Boost C++ Libraries. - compiler supports SSE2 : yes [2] - compiler supports SSE4.1 : yes [2] - has synchronization.lib : no [2] - has std::atomic_ref : yes [2] - has statx : yes [2] - has BCrypt API : no [2] - is Windows CE : yes [2] - has init_priority attribute : yes [2] - has stat::st_blksize : yes [2] - has stat::st_mtim : yes [2] - has stat::st_mtimensec : yes [2] - has stat::st_mtimespec : yes [2] - has stat::st_birthtim : yes [2] - has stat::st_birthtimensec : yes [2] - has stat::st_birthtimespec : yes [2] - has fdopendir(O_NOFOLLOW) : yes [2] - has POSIX *at APIs : yes [2] - cxx11_auto_declarations : yes [2] - cxx11_constexpr : yes [2] - cxx11_defaulted_functions : yes [2] - cxx11_final : yes [2] - cxx11_hdr_mutex : yes [2] - cxx11_hdr_tuple : yes [2] - cxx11_lambdas : yes [2] - cxx11_noexcept : yes [2] - cxx11_nullptr : yes [2] - cxx11_rvalue_references : yes [2] - cxx11_template_aliases : yes [2] - cxx11_thread_local : yes [2] - cxx11_variadic_templates : yes [2] - cxx11_auto_declarations : yes [3] - cxx11_constexpr : yes [3] - cxx11_defaulted_functions : yes [3] - cxx11_final : yes [3] - cxx11_hdr_mutex : yes [3] - cxx11_hdr_tuple : yes [3] - cxx11_lambdas : yes [3] - cxx11_noexcept : yes [3] - cxx11_nullptr : yes [3] - cxx11_rvalue_references : yes [3] - cxx11_template_aliases : yes [3] - cxx11_thread_local : yes [3] - cxx11_variadic_templates : yes [3] - has_icu builds : no [2] warning: Graph library does not contain MPI-based parallel components. note: to enable them, add &quot;using mpi ;&quot; to your user-config.jam. note: to suppress this message, pass &quot;--without-graph_parallel&quot; to bjam. - zlib : no [4] - bzip2 : no [4] - lzma : no [4] - zstd : no [4] - has_lzma_cputhreads builds : no [2] - cxx11_decltype : yes [2] - cxx11_basic_alignas : yes [2] - iconv (libc) : no [2] - iconv (separate) : no [2] - icu : no [2] - cxx11_defaulted_moves : yes [2] - cxx11_hdr_functional : yes [2] - cxx11_hdr_type_traits : yes [2] - cxx11_override : yes [2] - cxx11_range_based_for : yes [2] - cxx11_scoped_enums : yes [2] - cxx11_smart_ptr : yes [2] - cxx11_static_assert : yes [2] - native atomic int32 supported : yes [2] - has message compiler : no [2] - native syslog supported : yes [2] - pthread supports robust mutexes : yes [2] - compiler supports SSSE3 : yes [2] - compiler supports AVX2 : yes [2] - gcc visibility : yes [2] - sfinae_expr : yes [2] - cxx11_unified_initialization_syntax : yes [2] - cxx11_hdr_initializer_list : yes [2] - cxx11_hdr_chrono : yes [2] - cxx11_numeric_limits : yes [2] - cxx11_hdr_array : yes [2] - cxx11_hdr_atomic : yes [2] - cxx11_allocator : yes [2] - cxx11_explicit_conversion_operators : yes [2] - long double support : yes [2] warning: skipping optional Message Passing Interface (MPI) library. note: to enable MPI support, add &quot;using mpi ;&quot; to user-config.jam. note: to suppress this message, pass &quot;--without-mpi&quot; to bjam. note: otherwise, you can safely ignore this message. - std::fstream is moveable and swappable : yes [2] - Has Large File Support : no [2] - Has attribute init_priority : yes [2] - libbacktrace builds : no [2] - libbacktrace builds : no [3] - addr2line builds : no [2] - addr2line builds : no [3] - WinDbg builds : no [2] - WinDbg builds : no [3] - WinDbgCached builds : no [2] - WinDbgCached builds : no [3] - BOOST_COMP_GNUC &gt;= 4.3.0 : yes [2] - cxx11_hdr_thread : yes [2] - cxx11_hdr_regex : yes [2] - compiler supports SSE2 : yes [5] - compiler supports SSE4.1 : yes [5] - has synchronization.lib : no [5] - has std::atomic_ref : yes [5] - has statx : yes [5] - has BCrypt API : no [5] - is Windows CE : yes [5] - has init_priority attribute : yes [5] - has stat::st_blksize : yes [5] - has stat::st_mtim : yes [5] - has stat::st_mtimensec : yes [5] - has stat::st_mtimespec : yes [5] - has stat::st_birthtim : yes [5] - has stat::st_birthtimensec : yes [5] - has stat::st_birthtimespec : yes [5] - has fdopendir(O_NOFOLLOW) : yes [5] - has POSIX *at APIs : yes [5] - cxx11_auto_declarations : yes [5] - cxx11_constexpr : yes [5] - cxx11_defaulted_functions : yes [5] - cxx11_final : yes [5] - cxx11_hdr_mutex : yes [5] - cxx11_hdr_tuple : yes [5] - cxx11_lambdas : yes [5] - cxx11_noexcept : yes [5] - cxx11_nullptr : yes [5] - cxx11_rvalue_references : yes [5] - cxx11_template_aliases : yes [5] - cxx11_thread_local : yes [5] - cxx11_variadic_templates : yes [5] - cxx11_auto_declarations : yes [6] - cxx11_constexpr : yes [6] - cxx11_defaulted_functions : yes [6] - cxx11_final : yes [6] - cxx11_hdr_mutex : yes [6] - cxx11_hdr_tuple : yes [6] - cxx11_lambdas : yes [6] - cxx11_noexcept : yes [6] - cxx11_nullptr : yes [6] - cxx11_rvalue_references : yes [6] - cxx11_template_aliases : yes [6] - cxx11_thread_local : yes [6] - cxx11_variadic_templates : yes [6] - has_icu builds : no [5] - zlib : no (cached) [4] - bzip2 : no (cached) [4] - lzma : no (cached) [4] - zstd : no (cached) [4] - has_lzma_cputhreads builds : no [5] - cxx11_decltype : yes [5] - cxx11_basic_alignas : yes [5] - iconv (libc) : no [5] - iconv (separate) : no [5] - icu : no [5] - cxx11_defaulted_moves : yes [5] - cxx11_hdr_functional : yes [5] - cxx11_hdr_type_traits : yes [5] - cxx11_override : yes [5] - cxx11_range_based_for : yes [5] - cxx11_scoped_enums : yes [5] - cxx11_smart_ptr : yes [5] - cxx11_static_assert : yes [5] - native atomic int32 supported : yes [5] - has message compiler : no [5] - native syslog supported : yes [5] - pthread supports robust mutexes : yes [5] - compiler supports SSSE3 : yes [5] - compiler supports AVX2 : yes [5] - gcc visibility : yes [5] - sfinae_expr : yes [5] - cxx11_unified_initialization_syntax : yes [5] - cxx11_hdr_initializer_list : yes [5] - cxx11_hdr_chrono : yes [5] - cxx11_numeric_limits : yes [5] - cxx11_hdr_array : yes [5] - cxx11_hdr_atomic : yes [5] - cxx11_allocator : yes [5] - cxx11_explicit_conversion_operators : yes [5] - long double support : yes [5] - std::fstream is moveable and swappable : yes [5] - Has Large File Support : no [5] - Has attribute init_priority : yes [5] - libbacktrace builds : no [5] - libbacktrace builds : no [6] - addr2line builds : no [5] - addr2line builds : no [6] - WinDbg builds : no [5] - WinDbg builds : no [6] - WinDbgCached builds : no [5] - WinDbgCached builds : no [6] - BOOST_COMP_GNUC &gt;= 4.3.0 : yes [5] - cxx11_hdr_thread : yes [5] - cxx11_hdr_regex : yes [5] [1] msvc-6.0 [2] msvc-6.0/debug/address-model-64/architecture-x86/link-static/python-3.13/runtime-link-static/threadapi-win32/threading-multi/visibility-hidden [3] msvc-6.0/debug/address-model-64/architecture-x86/build-no/link-static/python-3.13/runtime-link-static/threadapi-win32/threading-multi/visibility-hidden [4] link-static [5] msvc-6.0/release/address-model-64/architecture-x86/link-static/python-3.13/runtime-link-static/threadapi-win32/threading-multi/visibility-hidden [6] msvc-6.0/release/address-model-64/architecture-x86/build-no/link-static/python-3.13/runtime-link-static/threadapi-win32/threading-multi/visibility-hidden Component configuration: - atomic : building - chrono : building - container : building - context : building - contract : building - coroutine : building - date_time : building - exception : building - fiber : building - filesystem : building - graph : building - graph_parallel : building - headers : building - iostreams : building - json : building - locale : building - log : building - math : building - mpi : building - nowide : building - program_options : building - python : building - random : building - regex : building - serialization : building - stacktrace : building - system : building - test : building - thread : building - timer : building - type_erasure : building - url : building - wave : building ...patience... ...patience... ...patience... ...patience... ...patience... ...patience... ...found 17967 targets... ...updating 1296 targets... boost-install.generate-cmake-config- bin.v2\libs\fiber\build\stage\boost_fiber-config.cmake boost-install.generate-cmake-config-version- bin.v2\libs\fiber\build\stage\boost_fiber-config-version.cmake boost-install.generate-cmake-config- bin.v2\libs\filesystem\build\stage\boost_filesystem-config.cmake boost-install.generate-cmake-config-version- bin.v2\libs\filesystem\build\stage\boost_filesystem-config-version.cmake compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_iserializer.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_oarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_xml_archive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_oserializer.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_pointer_iserializer.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_pointer_oserializer.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_text_iprimitive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_serializer_map.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_text_oprimitive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\binary_iarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\extended_type_info_typeid.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\extended_type_info.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\text_iarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\extended_type_info_no_rtti.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\text_oarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\polymorphic_oarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\polymorphic_iarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\binary_oarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\polymorphic_text_iarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\stl_port.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\polymorphic_text_oarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\polymorphic_binary_iarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\polymorphic_binary_oarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\polymorphic_xml_iarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\polymorphic_xml_oarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\xml_iarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\xml_grammar.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\xml_archive_exception.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\void_cast.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\xml_oarchive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_text_woprimitive.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\utf8_codecvt_facet.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\codecvt_null.obj compile-c-c++ bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_text_wiprimitive.obj msvc.archive bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\libboost_serialization-vc6-mt-s-x64-1_82.lib The filename, directory name, or volume label syntax is incorrect. if exist &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\libboost_serialization-vc6-mt-s-x64-1_82.lib&quot; DEL &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\libboost_serialization-vc6-mt-s-x64-1_82.lib&quot; /out:&quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\libboost_serialization-vc6-mt-s-x64-1_82.lib&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\archive_exception.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_archive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_iarchive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_iserializer.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_oarchive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_oserializer.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_pointer_iserializer.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_pointer_oserializer.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_serializer_map.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_text_iprimitive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_text_oprimitive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\basic_xml_archive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\binary_iarchive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\binary_oarchive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\extended_type_info.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\extended_type_info_typeid.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\extended_type_info_no_rtti.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\polymorphic_iarchive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\polymorphic_oarchive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\stl_port.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\text_iarchive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-static\runtime-link-static\threading-multi\text_oarchive.obj&quot; &quot;bin.v2\libs\serialization\build\msvc-6.0\release\address-model-64\architecture-x86\link-stati ... ...failed updating 76 targets... ...skipped 502 targets... ...updated 718 targets... </code></pre> <p>The issue I see is <code>The filename, directory name, or volume label syntax is incorrect.</code>, but I've no idea where to go from here.</p>

https://stackoverflow.com/q/79725398 5 dogAwakeCat https://stackoverflow.com/users/2833518 2025-08-04T21:12:43Z 2025-08-05T09:31:59Z

<p>I'm trying to use <code>std::format</code> in a mex file but get an <code>unresolved external symbol</code> error when I use <code>std</code> as a module.</p> <pre><code>// test_mex.cpp #include &quot;mex.h&quot; import std; void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { std::string temp{std::format(&quot;{:.2f}&quot;, 21.234)}; mexPrintf(&quot;%s\n&quot;, temp.c_str()); plhs[0] = mxCreateDoubleMatrix(10, 10, mxREAL); } </code></pre> <p>compiled with</p> <pre><code>!cl /std:c++latest /O2 /Oi /EHsc /MD /nologo /c -I&quot;C:\Program Files (x86)\Windows Kits\10\Include\10.0.26100.0\ucrt&quot; -I&quot;C:\Program Files\Microsoft Visual Studio\2022\Community\VC\Tools\MSVC\14.44.35207\include&quot; &quot;C:\Program Files\Microsoft Visual Studio\2022\Community\VC\Tools\MSVC\14.44.35207\modules\std.ixx&quot; mex -O COMPFLAGS=&quot;$COMPFLAGS -std:c++latest /O2 /Oi /EHsc /Wall /MD&quot; test_mex.cpp; </code></pre> <p>gives this error</p> <pre><code>Error using mex Creating library test_mex.lib and object test_mex.exp test_mex.obj : error LNK2019: unresolved external symbol &quot;public: static int const std::_General_precision_tables_2&lt;double&gt;::_Max_P&quot; (?_Max_P@?$_General_precision_tables_2@N@std@@2HB) referenced in function &quot;struct std::to_chars_result __cdecl std::_Floating_to_chars_general_precision&lt;double&gt;(char *,char * const,double,int)&quot; (??$_Floating_to_chars_general_precision@N@std@@YA?AUto_chars_result@0@PEADQEADNH@Z) test_mex.obj : error LNK2019: unresolved external symbol &quot;public: static int const std::_General_precision_tables_2&lt;float&gt;::_Max_P&quot; (?_Max_P@?$_General_precision_tables_2@M@std@@2HB) referenced in function &quot;struct std::to_chars_result __cdecl std::_Floating_to_chars_general_precision&lt;float&gt;(char *,char * const,float,int)&quot; (??$_Floating_to_chars_general_precision@M@std@@YA?AUto_chars_result@0@PEADQEADMH@Z) test_mex.mexw64 : fatal error LNK1120: 2 unresolved externals </code></pre> <p>If you don't use <code>std::format</code> and instead just specify the string, for example <code>std::string temp{&quot;test&quot;}</code> it compiles and runs. Furthermore, instead of using a module if you replace <code>import std;</code> with <code>#include &lt;format&gt;</code> again it compiles and runs. I'm using the Microsoft Visual C++ 2022 compiler.</p>

https://stackoverflow.com/q/79699706 6 Nuno https://stackoverflow.com/users/11848912 2025-07-13T02:27:36Z 2025-07-27T19:25:20Z

<p>This very likely is a bug in MATLAB R2025a on Ubuntu, and a bug report has been submitted to MathWorks. While they investigate, I would like to know: <strong>is there anything wrong with the Fortran MEX file (see below)?</strong></p> <p>In brief, I find that if the Fortran MEX file contains an internal subroutine used as an actual argument to another subroutine (which is allowed in Fortran 2008+, ref. <strong>Note 12.18 on page 290 of WD 1539-1 J3/10-007r1, F2008 Working Document</strong> and the discussions on <a href="https://fortran-lang.discourse.group/t/thread-safe-implementation-of-a-parametrized-objective-function-without-using-internal-subroutines/9919" rel="noreferrer">Fortran Discourse</a>), then MATLAB R2025a segfaults under Ubuntu 22.04 and 24.04. The problem did not occur until I upgraded to R2025a. Windows and macOS have not been tested.</p> <p>To reproduce the issue, execute the following in the MATLAB command window:</p> <pre class="lang-matlab prettyprint-override"><code>mex('-v', '-g', 'optimize.F90'); optimize(); </code></pre> <p>, where optimize.F90 is the following Fortran MEX file:</p> <pre class="lang-none prettyprint-override"><code>!----------------------------------------------------------------------! ! optimize.F90: a Fortran MEX file demonstrating a bug of MATLAB R2025a. ! To test, do the following in the MATLAB command window: ! ! mex -g optimize.F90; optimize(); ! ! In MATLAB R2025a under Ubuntu 24.04 with gfortran 13.3.0, this causes ! a segmentation fault. In earlier versions of MATLAB, it works fine. !----------------------------------------------------------------------! #include &quot;fintrf.h&quot; !----------------------------------------------------------------------! ! solver_mod: a module defining an optimization solver module solver_mod implicit none private public:: solver abstract interface ! OBJ: an abstract interface for an objective function subroutine OBJ(x, y) real, intent(in):: x real, intent(out):: y end subroutine OBJ end interface contains ! solver: a doing-nothing solver for demonstration. Its interface CANNOT ! be changed, e.g., to accept a hyper-parameter for `objective`. subroutine solver(objective) procedure(OBJ):: objective real f ! The following line causes a segmentation fault in MATLAB R2025a when ! `solver(objective)` is called from `mexFunction`. call objective(0.0, f) end subroutine solver end module solver_mod !----------------------------------------------------------------------! !----------------------------------------------------------------------! ! mexFunction: the entry point for the MEX function `optimize`, which ! calls `solver` from solver_mod to optimize `objective`. subroutine mexFunction(nlhs, plhs, nrhs, prhs) use solver_mod, only: solver implicit none integer nlhs, nrhs mwPointer plhs(*), prhs(*) ! `hyper_parameter`: a hyper-parameter for the objective function. In ! real applications, this parameter is not a constant, but a variable ! depending on the input arguments of `mexFunction`. real hyper_parameter, f hyper_parameter = 42 ! The following line works fine. call objective(0.0, f) ! The following line causes a segmentation fault in MATLAB R2025a. call solver(objective) contains !--------------------------------------------------------------! ! objective: a simple objective function for demonstration. ! F2008 allows to pass internal procedures as actual arguments. ! See Note 12.18 on page 290 of WD 1539-1 J3/10-007r1 (F2008 ! Working Document). We implement `objective` internally so that ! `hyper_parameter` is visible to it. Do NOT pass this parameter ! by a module variable, which is thread-unsafe. subroutine objective(x, y) real, intent(in):: x real, intent(out):: y y = (x+hyper_parameter)**2 end subroutine objective !--------------------------------------------------------------! end subroutine mexFunction !----------------------------------------------------------------------! </code></pre> <p>The MATLAB output, including crash log, is:</p> <pre class="lang-none prettyprint-override"><code>Verbose mode is on. No MEX options file identified; looking for an implicit selection. ... Looking for compiler 'gfortran' ... ... Executing command 'which gfortran' ...Yes ('/usr/bin/gfortran'). ... Looking for folder '/usr/bin' ...Yes. ... Executing command 'which gfortran' ...Yes ('/usr/bin/gfortran'). ... Executing command 'which gfortran' ...Yes ('/usr/bin/gfortran'). ... Executing command 'gfortran -print-file-name=libgfortran.so' ...Yes ('/usr/lib/gcc/x86_64-linux-gnu/10/libgfortran.so'). ... Looking for folder '/usr/lib/x86_64-linux-gnu' ...Yes. ... Executing command 'which gfortran' ...Yes ('/usr/bin/gfortran'). ... Executing command 'gfortran -dumpversion | awk -F '.' '$1&gt;=6' ' ...Yes ('10'). Found installed compiler 'gfortran'. ... Looking for compiler 'gfortran6-' ... ... Executing command 'which gfortran' ...Yes ('/usr/bin/gfortran'). ... Looking for folder '/usr/bin' ...Yes. ... Executing command 'which gfortran' ...Yes ('/usr/bin/gfortran'). ... Executing command 'gfortran -print-file-name=libgfortran.so' ...Yes ('/usr/lib/gcc/x86_64-linux-gnu/10/libgfortran.so'). ... Looking for folder '/usr/lib/x86_64-linux-gnu' ...Yes. ... Executing command 'which gfortran' ...Yes ('/usr/bin/gfortran'). ... Executing command 'gfortran -print-file-name=libgfortranbegin.a' ...Yes ('libgfortranbegin.a'). ... Looking for folder 'libgfortranbegin.a' ...No. Did not find installed compiler 'gfortran6-'. Options file details ------------------------------------------------------------------- Compiler location: /usr/bin Options file: /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/mexopts/gfortran6.xml CMDLINE2 : /usr/bin/gfortran -pthread -shared -g -Wl,--version-script,/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/lib/glnxa64/fortran_exportsmexfileversion.map /tmp/mex_1037477951999_9281/optimize.o /tmp/mex_1037477951999_9281/fortran_mexapi_version.o -Wl,-rpath-link,/opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64 -L&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/sys/os/glnxa64/orig&quot; -L&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64&quot; &quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/sys/os/glnxa64/orig/libstdc++.so.6&quot; -lmx -lmex -lmat -lm -L&quot;/usr/lib/x86_64-linux-gnu&quot; -lgfortran -z noexecstack -o optimize.mexa64 FC : /usr/bin/gfortran DEFINES : -DMATLAB_DEFAULT_RELEASE=R2017b -DUSE_MEX_CMD FFLAGS : -fexceptions -fbackslash -fPIC -fno-omit-frame-pointer -fallow-argument-mismatch -fdefault-integer-8 INCLUDE : -I&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/include&quot; -I&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/simulink/include&quot; FOPTIMFLAGS : -O2 FDEBUGFLAGS : -g LDF : /usr/bin/gfortran LDFLAGS : -pthread LDTYPE : -shared LINKEXPORT : -Wl,--version-script,/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/lib/glnxa64/fexport.map LINKEXPORTVER : -Wl,--version-script,/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/lib/glnxa64/fortran_exportsmexfileversion.map MWCPPLIB : &quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/sys/os/glnxa64/orig/libstdc++.so.6&quot; LINKLIBS : -Wl,-rpath-link,/opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64 -L&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/sys/os/glnxa64/orig&quot; -L&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64&quot; &quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/sys/os/glnxa64/orig/libstdc++.so.6&quot; -lmx -lmex -lmat -lm -L&quot;/usr/lib/x86_64-linux-gnu&quot; -lgfortran LDOPTIMFLAGS : -O LDDEBUGFLAGS : -g MW_GLIBC_SHIM : $MW_GLIBC_SHIM LDEXECSTACK : -z noexecstack OBJEXT : .o LDEXT : .mexa64 SETENV : FC=&quot;/usr/bin/gfortran&quot; FFLAGS=&quot;-fexceptions -fbackslash -fPIC -fno-omit-frame-pointer -fallow-argument-mismatch -fdefault-integer-8 -DMATLAB_DEFAULT_RELEASE=R2017b -DUSE_MEX_CMD &quot; FOPTIMFLAGS=&quot;-O2&quot; FDEBUGFLAGS=&quot;-g&quot; LD=&quot;/usr/bin/gfortran&quot; LDFLAGS=&quot;-pthread -shared -Wl,-rpath-link,/opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64 -L&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/sys/os/glnxa64/orig&quot; -L&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64&quot; &quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/sys/os/glnxa64/orig/libstdc++.so.6&quot; -lmx -lmex -lmat -lm -L&quot;/usr/lib/x86_64-linux-gnu&quot; -lgfortran -Wl,--version-script,/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/lib/glnxa64/fexport.map&quot; LDDEBUGFLAGS=&quot;-g&quot; GFORTRAN_INSTALLDIR : /usr/bin GFORTRAN_EXECUTABLE : /usr/bin/gfortran GFORTRAN_LIBDIR : /usr/lib/x86_64-linux-gnu GFORTRAN_VERSION : 10 MATLABROOT : /opt/hostedtoolcache/MATLAB/2025.1.999/x64 ARCH : glnxa64 SRC : &quot;/home/runner/work/test_matlab/test_matlab/segfault_2025a/optimize.F90&quot;;&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/version/fortran_mexapi_version.F&quot; OBJ : /tmp/mex_1037477951999_9281/optimize.o;/tmp/mex_1037477951999_9281/fortran_mexapi_version.o OBJS : /tmp/mex_1037477951999_9281/optimize.o /tmp/mex_1037477951999_9281/fortran_mexapi_version.o SRCROOT : /home/runner/work/test_matlab/test_matlab/segfault_2025a/optimize DEF : /tmp/mex_1037477951999_9281/optimize.def EXP : &quot;optimize.exp&quot; LIB : &quot;optimize.lib&quot; EXE : optimize.mexa64 ILK : &quot;optimize.ilk&quot; TEMPNAME : optimize EXEDIR : EXENAME : optimize MANIFEST : &quot;optimize.mexa64.manifest&quot; OPTIM : -g LINKOPTIM : -g CMDLINE1_0 : /usr/bin/gfortran -c -DMATLAB_DEFAULT_RELEASE=R2017b -DUSE_MEX_CMD -I&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/include&quot; -I&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/simulink/include&quot; -fexceptions -fbackslash -fPIC -fno-omit-frame-pointer -fallow-argument-mismatch -fdefault-integer-8 -g &quot;/home/runner/work/test_matlab/test_matlab/segfault_2025a/optimize.F90&quot; -o /tmp/mex_1037477951999_9281/optimize.o CMDLINE1_1 : /usr/bin/gfortran -c -DMATLAB_DEFAULT_RELEASE=R2017b -DUSE_MEX_CMD -I&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/include&quot; -I&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/simulink/include&quot; -fexceptions -fbackslash -fPIC -fno-omit-frame-pointer -fallow-argument-mismatch -fdefault-integer-8 -g &quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/version/fortran_mexapi_version.F&quot; -o /tmp/mex_1037477951999_9281/fortran_mexapi_version.o ------------------------------------------------------------------- Building with 'gfortran'. /usr/bin/gfortran -c -DMATLAB_DEFAULT_RELEASE=R2017b -DUSE_MEX_CMD -I&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/include&quot; -I&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/simulink/include&quot; -fexceptions -fbackslash -fPIC -fno-omit-frame-pointer -fallow-argument-mismatch -fdefault-integer-8 -g &quot;/home/runner/work/test_matlab/test_matlab/segfault_2025a/optimize.F90&quot; -o /tmp/mex_1037477951999_9281/optimize.o f951: Warning: Nonexistent include directory /opt/hostedtoolcache/MATLAB/2025.1.999/x64/simulink/include [-Wmissing-include-dirs] /usr/bin/gfortran -c -DMATLAB_DEFAULT_RELEASE=R2017b -DUSE_MEX_CMD -I&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/include&quot; -I&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/simulink/include&quot; -fexceptions -fbackslash -fPIC -fno-omit-frame-pointer -fallow-argument-mismatch -fdefault-integer-8 -g &quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/version/fortran_mexapi_version.F&quot; -o /tmp/mex_1037477951999_9281/fortran_mexapi_version.o f951: Warning: Nonexistent include directory /opt/hostedtoolcache/MATLAB/2025.1.999/x64/simulink/include [-Wmissing-include-dirs] { &quot;bundle.symbolic_name&quot; : &quot;be8e5504-e0f7-492e-b08a-03d11430c8f2&quot;, &quot;mw&quot; : { &quot;mex&quot; : { &quot;apiVersion&quot; : 0, &quot;release&quot; : &quot;R2025a&quot;, &quot;threadpoolSafe&quot; : 0 } } } /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/usResourceCompiler3 --manifest-add /tmp/mex_1037477951999_9281//mw_mex_tempmex_manifest.json --bundle-name mexVersioning --out-file /tmp/mex_1037477951999_9281//mw_mex_tempmex_bundle.zip /usr/bin/gfortran -pthread -shared -g -Wl,--version-script,/opt/hostedtoolcache/MATLAB/2025.1.999/x64/extern/lib/glnxa64/fortran_exportsmexfileversion.map /tmp/mex_1037477951999_9281/optimize.o /tmp/mex_1037477951999_9281/fortran_mexapi_version.o -Wl,-rpath-link,/opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64 -L&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/sys/os/glnxa64/orig&quot; -L&quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64&quot; &quot;/opt/hostedtoolcache/MATLAB/2025.1.999/x64/sys/os/glnxa64/orig/libstdc++.so.6&quot; -lmx -lmex -lmat -lm -L&quot;/usr/lib/x86_64-linux-gnu&quot; -lgfortran -z noexecstack -o optimize.mexa64 -Wl,--format=binary -Wl,/tmp/mex_1037477951999_9281//mw_mex_tempmex_bundle.zip -Wl,--format=default MEX completed successfully. -------------------------------------------------------------------------------- Segmentation violation detected at 2025-07-11 16:58:56 +0000 -------------------------------------------------------------------------------- Configuration: Crash Decoding : Disabled - No sandbox or build area path Crash Mode : continue (default) Default Encoding : UTF-8 Deployed : false GNU C Library : 2.39 stable Graphics Driver : Uninitialized software Graphics card 1 : 0x1414 ( 0x1414 ) 0x5353 Version 0.0.0.0 (0-0-0) MATLAB Architecture : glnxa64 MATLAB Entitlement ID : 11373431 MATLAB Root : /opt/hostedtoolcache/MATLAB/2025.1.999/x64 MATLAB Version : 25.1.0.2943329 (R2025a) OpenGL : software Operating System : Ubuntu 24.04.2 LTS Process ID : 9281 Processor ID : x86 Family 25 Model 1 Stepping 1, AuthenticAMD Session Key : 40xq29adxbsbdw85fpki4ykis Window System : No active display Fault Count: 1 Abnormal termination: Segmentation violation Current Thread: 'MCR 0 interpret' id 140534882498240 Register State (from fault): RAX = 00007fd0d3bf873c RBX = 00007fd0d3bf8910 RCX = 00007fd0ca49f11c RDX = 00007fd0d3bf8774 RSP = 00007fd0d3bf8718 RBP = 00007fd0d3bf8740 RSI = 00007fd0d3bf873c RDI = 00007fd0ca4a0000 R8 = 00007fd0d3bf8a00 R9 = 0000000000000000 R10 = 00007fd0d3bf8770 R11 = 0000000000000246 R12 = 00007fd0d3bf8900 R13 = 0000000000000000 R14 = 00007fd0d3bf8a00 R15 = 00007fd0d3bf8970 RIP = 00007fd0d3bf8774 EFL = 0000000000010202 CS = 0033 FS = 0000 GS = 0000 Stack Trace (from fault): [ 0] 0x00007fd0d3bf8774 &lt;unknown-module&gt;+00000000 [ 1] 0x00007fd0ca49f1d2 /home/runner/work/test_matlab/test_matlab/segfault_2025a/optimize.mexa64+00004562 mexfunction_+00000132 [ 2] 0x00007fd12df672d2 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmex.so+00975570 [ 3] 0x00007fd12df67358 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmex.so+00975704 [ 4] 0x00007fd12df67015 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmex.so+00974869 [ 5] 0x00007fd12df681a4 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmex.so+00979364 [ 6] 0x00007fd12df524fa /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmex.so+00890106 [ 7] 0x00007fd12e77a4fe /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_dispatcher.so+01549566 _ZN8Mfh_file20dispatch_file_commonEMS_FviPP11mxArray_tagiS2_EiS2_iS2_+00000142 [ 8] 0x00007fd12e77bc9a /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_dispatcher.so+01555610 [ 9] 0x00007fd12e77c04e /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_dispatcher.so+01556558 _ZN8Mfh_file8dispatchEiPSt10unique_ptrI11mxArray_tagN6matrix6detail17mxDestroy_deleterEEiPPS1_+00000030 [ 10] 0x00007fd12d8757a7 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwlxemainservices.so+02578343 [ 11] 0x00007fd120762fd0 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+11939792 [ 12] 0x00007fd1207566da /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+11888346 [ 13] 0x00007fd1206d9fd1 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+11378641 [ 14] 0x00007fd1203c86cb /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08160971 [ 15] 0x00007fd1203ca41c /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08168476 [ 16] 0x00007fd1203c71ff /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08155647 [ 17] 0x00007fd1203d8c55 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08227925 [ 18] 0x00007fd1203d94a9 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08230057 [ 19] 0x00007fd1203c6ff4 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08155124 [ 20] 0x00007fd1203c70ff /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08155391 [ 21] 0x00007fd12053583b /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+09656379 [ 22] 0x00007fd12053a7e4 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+09676772 [ 23] 0x00007fd12da0b9b4 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwlxemainservices.so+04241844 [ 24] 0x00007fd12d863683 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwlxemainservices.so+02504323 [ 25] 0x00007fd12d865c05 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwlxemainservices.so+02513925 [ 26] 0x00007fd12e77a4fe /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_dispatcher.so+01549566 _ZN8Mfh_file20dispatch_file_commonEMS_FviPP11mxArray_tagiS2_EiS2_iS2_+00000142 [ 27] 0x00007fd12e77bc9a /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_dispatcher.so+01555610 [ 28] 0x00007fd12e77c04e /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_dispatcher.so+01556558 _ZN8Mfh_file8dispatchEiPSt10unique_ptrI11mxArray_tagN6matrix6detail17mxDestroy_deleterEEiPPS1_+00000030 [ 29] 0x00007fd12d8757a7 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwlxemainservices.so+02578343 [ 30] 0x00007fd120762fd0 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+11939792 [ 31] 0x00007fd1207566da /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+11888346 [ 32] 0x00007fd1206d9fd1 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+11378641 [ 33] 0x00007fd1203c86cb /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08160971 [ 34] 0x00007fd1203ca41c /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08168476 [ 35] 0x00007fd1203c71ff /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08155647 [ 36] 0x00007fd1203d8c55 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08227925 [ 37] 0x00007fd1203d94a9 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08230057 [ 38] 0x00007fd1203c6ff4 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08155124 [ 39] 0x00007fd1203c70ff /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+08155391 [ 40] 0x00007fd12053583b /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+09656379 [ 41] 0x00007fd12053a7e4 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_lxe.so+09676772 [ 42] 0x00007fd12da0b9b4 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwlxemainservices.so+04241844 [ 43] 0x00007fd12d8df311 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwlxemainservices.so+03011345 [ 44] 0x00007fd12d8df50d /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwlxemainservices.so+03011853 [ 45] 0x00007fd12d9ac31d /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwlxemainservices.so+03851037 [ 46] 0x00007fd12d9ac54e /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwlxemainservices.so+03851598 [ 47] 0x00007fd12e356f0f /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwm_interpreter.so+01404687 _Z51inEvalCmdWithLocalReturnInDesiredWSAndPublishEventsRKNSt7__cxx1112basic_stringIDsSt11char_traitsIDsESaIDsEEEPibbP15inWorkSpace_tagN9MathWorks3lxe10EvalSourceE+00000063 [ 48] 0x00007fd12eadeecc /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00974540 _ZNK3iqm18InternalEvalPlugin24inEvalCmdWithLocalReturnERKNSt7__cxx1112basic_stringIDsSt11char_traitsIDsESaIDsEEEP15inWorkSpace_tag+00000108 [ 49] 0x00007fd12eadff35 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00978741 _ZN3iqm18InternalEvalPlugin7executeEP15inWorkSpace_tag+00000469 [ 50] 0x00007fd12eac7716 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00878358 [ 51] 0x00007fd12ea8c228 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00635432 [ 52] 0x00007fd12ea8c571 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00636273 [ 53] 0x00007fd13bf4954e /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmlutil.so+09737550 _ZNK14cmddistributor16IIPRunNowMessage7deliverERKN10foundation7msg_svc8exchange7RoutingE+00000030 [ 54] 0x00007fd13d3629cc /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwms.so+03549644 _ZN10foundation7msg_svc8exchange12MessageQueue7deliverERKN7mwboost10shared_ptrIKNS1_8EnvelopeEEE+00000252 [ 55] 0x00007fd13d36467f /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwms.so+03556991 [ 56] 0x00007fd13d349f6d /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwms.so+03448685 [ 57] 0x00007fd13d34ee74 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwms.so+03468916 [ 58] 0x00007fd13d34911f /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwms.so+03445023 [ 59] 0x00007fd13be86452 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmlutil.so+08938578 [ 60] 0x00007fd13be8d3d1 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmlutil.so+08967121 [ 61] 0x00007fd13ed2531e /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwrcf_framework.so+00312094 _ZN7mwboost6detail17shared_state_base13wait_internalERNS_11unique_lockINS_5mutexEEEb+00000222 [ 62] 0x00007fd12e93c612 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmcr.so+00706066 _ZN7mwboost6futureIvE3getEv+00000098 [ 63] 0x00007fd12e92ac0b /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmcr.so+00633867 [ 64] 0x00007fd13e347897 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmvm.so+03438743 _ZN14cmddistributor15PackagedTaskIIP10invokeFuncIN7mwboost8functionIFvvEEEEENS2_10shared_ptrINS2_6futureIDTclfp_EEEEEERKT_+00000071 [ 65] 0x00007fd13e347bc8 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmvm.so+03439560 _ZNSt17_Function_handlerIFN7mwboost3anyEvEZN14cmddistributor15PackagedTaskIIP10createFuncINS0_8functionIFvvEEEEESt8functionIS2_ET_EUlvE_E9_M_invokeERKSt9_Any_data+00000024 [ 66] 0x00007fd12eae9e3b /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+01019451 _ZN3iqm18PackagedTaskPlugin7executeEP15inWorkSpace_tag+00000091 [ 67] 0x00007fd12eac7716 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00878358 [ 68] 0x00007fd12ea934f9 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00664825 [ 69] 0x00007fd12e4ac0e9 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwbridge.so+00499945 [ 70] 0x00007fd12e4ac5c3 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwbridge.so+00501187 [ 71] 0x00007fd12e4c84ea /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwbridge.so+00615658 _Z22mnGetCommandLineBufferbRbN7mwboost8optionalIKP15inWorkSpace_tagEEbRKNS0_9function2IN6mlutil14cmddistributor17inExecutionStatusERKNSt7__cxx1112basic_stringIDsSt11char_traitsIDsESaIDsEEES4_EE+00000218 [ 72] 0x00007fd12e4c8881 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwbridge.so+00616577 _Z8mnParserv+00000513 [ 73] 0x00007fd12e964350 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmcr.so+00869200 [ 74] 0x00007fd13e347897 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmvm.so+03438743 _ZN14cmddistributor15PackagedTaskIIP10invokeFuncIN7mwboost8functionIFvvEEEEENS2_10shared_ptrINS2_6futureIDTclfp_EEEEEERKT_+00000071 [ 75] 0x00007fd13e347bc8 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmvm.so+03439560 _ZNSt17_Function_handlerIFN7mwboost3anyEvEZN14cmddistributor15PackagedTaskIIP10createFuncINS0_8functionIFvvEEEEESt8functionIS2_ET_EUlvE_E9_M_invokeERKSt9_Any_data+00000024 [ 76] 0x00007fd12eae9e3b /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+01019451 _ZN3iqm18PackagedTaskPlugin7executeEP15inWorkSpace_tag+00000091 [ 77] 0x00007fd12eac7716 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00878358 [ 78] 0x00007fd12ea918b2 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00657586 [ 79] 0x00007fd12ea9267d /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00661117 [ 80] 0x00007fd12ea92994 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwiqm.so+00661908 [ 81] 0x00007fd12e94bd07 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmcr.so+00769287 [ 82] 0x00007fd12e94c3ce /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmcr.so+00771022 [ 83] 0x00007fd12e94c69d /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwmcr.so+00771741 [ 84] 0x00007fd13ea0db17 /opt/hostedtoolcache/MATLAB/2025.1.999/x64/bin/glnxa64/libmwboost_thread.so.1.81.0+00043799 [ 85] 0x00007fd13f09caa4 /lib/x86_64-linux-gnu/libc.so.6+00641700 [ 86] 0x00007fd13f129c3c /lib/x86_64-linux-gnu/libc.so.6+01219644 This error was detected while a MEX-file was running. If the MEX-file is not an official MathWorks function, please examine its source code for errors. Please consult the External Interfaces Guide for information on debugging MEX-files. MATLAB is exiting because of fatal error exit status 255 </code></pre> <p>The same ran without error for me on earlier versions of MATLAB.</p> <p>For more details, see <a href="https://github.com/zequipe/test_matlab/tree/master/segfault_2025a" rel="noreferrer">https://github.com/zequipe/test_matlab/tree/master/segfault_2025a</a> .</p> <p>See <a href="https://github.com/zequipe/test_matlab/actions/runs/16096537922" rel="noreferrer">https://github.com/zequipe/test_matlab/actions/runs/16096537922</a> for a GitHub Actions run that reproduces the issue.</p>

https://stackoverflow.com/q/79532477 3 TheCubeDude https://stackoverflow.com/users/12590111 2025-03-25T01:04:37Z 2025-03-25T01:04:37Z

<p>I am trying to convert a matlab function to its equivalent CUDA MEX function. Below is the matlab function:</p> <pre><code>function [ShiftData_time] = simpledelayupdated_mat(RData, delay, fs) nfft = size(RData,1); binStart = floor(nfft/2); fftBin = (2*pi*ifftshift(((0:nfft-1)-binStart).'))/nfft; fftBin = fftBin'; RFFT = fft(RData, nfft, 1); ShiftData = RFFT .* ((exp(-1i * delay * fftBin)).'); ShiftData_time = ifft(ShiftData, nfft, 1); ShiftData_time = real(ShiftData_time); end </code></pre> <p>Here is the CUDA implementation using MEX :</p> <pre><code>#include &quot;mex.h&quot; #include &lt;cuda_runtime.h&gt; #include &lt;cufft.h&gt; #include &lt;math.h&gt; #define M_PI 3.14159 // Macro for cuFFT error checking. #define CUFFT_CHECK_ERROR(call) { cufftResult err = call; if(err != CUFFT_SUCCESS) { \ mexErrMsgIdAndTxt(&quot;CUDA:CUFFTError&quot;, &quot;CUFFT error: %d&quot;, err); } } // Macro for CUDA runtime error checking. #define CUDA_CHECK_ERROR(call) { cudaError_t err = call; if(err != cudaSuccess) { \ mexErrMsgIdAndTxt(&quot;CUDA:RuntimeError&quot;, &quot;CUDA error: %s&quot;, cudaGetErrorString(err)); } } // Updated CUDA kernel to apply the phase shift (delay) in the frequency domain. // Instead of launching one block per column with nfft threads, we now use a 2D grid. // grid.x corresponds to the column index, and grid.y is the number of blocks needed // to cover all frequency bins (with a fixed block size). __global__ void applyDelay(cufftDoubleComplex *data, const double *fftBin, double delay, int nfft, int nCols) { int col = blockIdx.x; // Column index. int blockId = blockIdx.y; // Block index within the column. int threadsPerBlock = blockDim.x; int threadId = threadIdx.x; // Compute global frequency bin index within this column. int index = blockId * threadsPerBlock + threadId; if (col &lt; nCols &amp;&amp; index &lt; nfft) { int idx = col * nfft + index; double phase = -delay * fftBin[index]; double cosPhase = cos(phase); double sinPhase = sin(phase); cufftDoubleComplex factor; factor.x = cosPhase; factor.y = sinPhase; cufftDoubleComplex orig = data[idx]; cufftDoubleComplex result; result.x = orig.x * factor.x - orig.y * factor.y; result.y = orig.x * factor.y + orig.y * factor.x; data[idx] = result; } } // The MEX gateway function. void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { // Check input count: RData, delay, fs. if(nrhs != 3) mexErrMsgIdAndTxt(&quot;MATLAB:simpledelayupdated:invalidNumInputs&quot;, &quot;Three inputs required: RData, delay, fs.&quot;); // Ensure RData is real and double precision. if(!mxIsDouble(prhs[0]) || mxIsComplex(prhs[0])) mexErrMsgIdAndTxt(&quot;MATLAB:simpledelayupdated:inputNotDouble&quot;, &quot;RData must be real double precision.&quot;); // Get input RData and its dimensions. double *RData = mxGetPr(prhs[0]); mwSize nfft = mxGetM(prhs[0]); // Number of rows. mwSize nCols = mxGetN(prhs[0]); // Number of columns. // Read scalar inputs: delay and fs (fs is provided for compatibility). double delay = mxGetScalar(prhs[1]); double fs = mxGetScalar(prhs[2]); // Not used directly. // Create output array (real, same size as input). plhs[0] = mxCreateDoubleMatrix(nfft, nCols, mxREAL); double *ShiftData_time = mxGetPr(plhs[0]); // Allocate host memory for complex input data. cufftDoubleComplex *h_complex = (cufftDoubleComplex*) mxMalloc(nfft * nCols * sizeof(cufftDoubleComplex)); for (mwSize col = 0; col &lt; nCols; col++) { for (mwSize row = 0; row &lt; nfft; row++) { int idx = col * nfft + row; h_complex[idx].x = RData[idx]; h_complex[idx].y = 0.0; } } // Precompute the FFT frequency bins. // This corresponds to the MATLAB expression: // fftBin = (2*pi*ifftshift(((0:nfft-1)-floor(nfft/2))))/nfft double *h_fftBin = (double*) mxMalloc(nfft * sizeof(double)); int binStart = (int)(nfft / 2); for (mwSize row = 0; row &lt; nfft; row++) { int shifted_idx = (row + binStart) % nfft; int centered = row - binStart; h_fftBin[shifted_idx] = (2.0 * M_PI * centered) / (double)nfft; } // Allocate device memory for the complex data and the fftBin vector. cufftDoubleComplex *d_data; CUDA_CHECK_ERROR(cudaMalloc((void**)&amp;d_data, nfft * nCols * sizeof(cufftDoubleComplex))); CUDA_CHECK_ERROR(cudaMemcpy(d_data, h_complex, nfft * nCols * sizeof(cufftDoubleComplex), cudaMemcpyHostToDevice)); double *d_fftBin; CUDA_CHECK_ERROR(cudaMalloc((void**)&amp;d_fftBin, nfft * sizeof(double))); CUDA_CHECK_ERROR(cudaMemcpy(d_fftBin, h_fftBin, nfft * sizeof(double), cudaMemcpyHostToDevice)); // Create a cuFFT plan for batched 1D FFTs (each column is one FFT of size nfft). cufftHandle plan; int rank = 1; int n[1] = {(int)nfft}; int istride = 1, idist = nfft; int ostride = 1, odist = nfft; CUFFT_CHECK_ERROR(cufftPlanMany(&amp;plan, rank, n, NULL, istride, idist, NULL, ostride, odist, CUFFT_Z2Z, nCols)); // Execute forward FFT (time to frequency domain). CUFFT_CHECK_ERROR(cufftExecZ2Z(plan, d_data, d_data, CUFFT_FORWARD)); // Launch the CUDA kernel to apply the phase shift. int threadsPerBlock = 256; // Use a safe block size. int blocksPerColumn = (nfft + threadsPerBlock - 1) / threadsPerBlock; dim3 grid(nCols, blocksPerColumn, 1); // grid.x: columns, grid.y: blocks per column. dim3 block(threadsPerBlock, 1, 1); applyDelay&lt;&lt;&lt;grid, block&gt;&gt;&gt;(d_data, d_fftBin, delay, nfft, nCols); CUDA_CHECK_ERROR(cudaDeviceSynchronize()); // Execute the inverse FFT to convert back to time domain. CUFFT_CHECK_ERROR(cufftExecZ2Z(plan, d_data, d_data, CUFFT_INVERSE)); // Allocate host memory for the inverse FFT result and copy from device. cufftDoubleComplex *h_result = (cufftDoubleComplex*) mxMalloc(nfft * nCols * sizeof(cufftDoubleComplex)); CUDA_CHECK_ERROR(cudaMemcpy(h_result, d_data, nfft * nCols * sizeof(cufftDoubleComplex), cudaMemcpyDeviceToHost)); // Normalize the result by dividing by nfft and take the real part. for(mwSize i = 0; i &lt; nfft * nCols; i++) { ShiftData_time[i] = h_result[i].x / ((double)nfft); } // Clean up resources. cufftDestroy(plan); cudaFree(d_data); cudaFree(d_fftBin); mxFree(h_complex); mxFree(h_fftBin); mxFree(h_result); } </code></pre> <p>This is how I am calling the Matlab function:</p> <pre><code>[temp] = simpledelayupdated(gather(RF_Arr_gpu), scalarDelay, sampling_Freq); </code></pre> <p>where RF_arr_gpu is a 2D MATLAB gpuArray, scalarDelay and sampling_Freq are scalars.</p> <p>The issue I'm facing is that the final output from both functions does not match when run on real data. However, I have tested individual components of the CUDA MEX function, and the differences between the MATLAB and CUDA MEX outputs are minimal.</p> <p>Can anyone help me identify where there could be issues in the code? I have attached the output from MATLAB and what we are currently getting from CUDA MEX. The MATLAB image is continuous while the CUDA MEX seems to be more discrete.</p> <p><img src="https://i.sstatic.net/ikgxXlj8.png" alt="MATLAB image" /> <img src="https://i.sstatic.net/Qs7B732n.png" alt="CUDA MEX image" /></p>

https://stackoverflow.com/q/79447789 4 Cal https://stackoverflow.com/users/29686189 2025-02-18T09:48:59Z 2025-02-20T21:01:44Z

<p>I have matrices <strong>A</strong> (m by v) and <strong>B</strong> (v by n). I also have a logical matrix <strong>L</strong> (m by n).</p> <p>I am interested in calculating only the values in <strong>A</strong> * <strong>B</strong> that correspond to logical values in <strong>L</strong> (values of 1s). Essentially I am interested in the quantity ( <strong>A</strong> * <strong>B</strong> ) .* <strong>L</strong> .</p> <p>For my problem, a typical L matrix has less than 0.1% percent of its values as 1s; the vast majority of the values are 0s. Thus, it makes no sense for me to literally perform ( <strong>A</strong> * <strong>B</strong> ) .* <strong>L</strong> , it would actually be faster to loop over each row of <strong>A</strong> * <strong>B</strong> that I want to compute, but even that is inefficient.</p> <hr /> <p><strong>Possible solution (need help vectorizing this code if possible)</strong></p> <p>My particular problem may have a nice solution given that the logical matrix <strong>L</strong> has a nice structure.</p> <p>Here's an example of L for a very small scale example (in most applications L is much much bigger and has much fewer 1-yellow entries, and many more 0-blue entries). <img src="https://i.imgur.com/P4TLMiG.png" alt="" /></p> <p>This L matrix is nice in that it can be represented as something like a permuted block matrix. This L in particular is composed of 9 &quot;blocks&quot; of 1s, where each block of 1s has its own set of row and column indices, defining a particular submatrix in <strong>A</strong> * <strong>B</strong>. For instance, the highlighted area here can be seen the values of 1 as a particular submatrix in L. <img src="https://i.imgur.com/JjsoeW0.png" alt="" /></p> <p>My solution was to do this: break the problem into submatrices over these blocks, and do matrix multiplications over each submatrix. I can get the row indices and column indices per each block's submatrix in L, organized in two cell lists &quot;rowidxs_list&quot; and &quot;colidxs_list&quot;, both with the number of cells equal to the number of blocks. For instance in the block example I gave, subblock 1, I could calculate those particular values in A * B by simply doing A( rowidxs_list{1} , : ) * B( : , colidxs_list{1} ) .</p> <p>That means that if I precomputed rowidxs_list and colidxs_list (ignore the costs of calculating these lists, they are negligable for my application), then my problem of calculating <strong>C</strong> = ( <strong>A</strong> * <strong>B</strong> ) .* <strong>L</strong> could effectively be done by:</p> <pre><code>C = sparse( m,n ) for i = 1:length( rowidxs_list ) C( rowidxs_list{i} , colidxs_list{i} ) = ... A( rowidxs_list{i} , : ) * B( : , colidxs_list{i} ); end </code></pre> <p>This seems like it would be the most efficient way to solve this problem <strong>if I knew how to vectorize this for loop</strong>. Does anyone see a way to vectorize this?</p> <p>There may be ways to vectorize if certain things hold, e.g. only if rowidxs_list and colidxs_list are matrix arrays instead of cell lists of lists (where each column in an array is an index list, thus replacing use of rowidxs_list{i} with rowidxs_list(i,:) ). I'd prefer to use cell lists here if possible since different lists can have different numbers of elements.</p> <hr /> <p><strong>other suggested solution (creating a mex file?)</strong></p> <p>I first posted this question on the /r/matlab subreddit, <a href="https://www.reddit.com/r/matlab/comments/1iro2m4/need_to_vectorize_efficiently_calculating_only/" rel="nofollow noreferrer">see here for the reddit thread</a>. The user &quot;qtac&quot; recommended that a C-MEX function linking to C programming language:</p> <blockquote> <p>My gut feeling is the only way to really optimize this is with a C-MEX solution; otherwise, you are going to get obliterated by overhead from subsref in these loops. With C you could loop over L until you find a nonzero element, and then do only the row-column dot product needed to populate that specific element. You will miss out on a lot of the BLAS optimizations but the computational savings may make up for it. Honestly I bet an LLM could write 90%+ of that MEX function for you; it's a well-formulated problem.</p> </blockquote> <p>I think this could be a good solution to pursue, but I'd like other opinions as well.</p>

https://stackoverflow.com/q/78876822 0 Reesh https://stackoverflow.com/users/26752151 2024-08-15T21:00:12Z 2024-08-16T21:21:20Z

<p>I am running this software called REShAPE on Windows 10 using Matlab 2019b and VS Studio 2017. although I get this error:</p> <pre><code> MATLAB: &gt; In vl_compilenn&gt;activate_nvcc (line 768) MATLAB: In vl_compilenn (line 300) MATLAB: In segmfunc (line 63)] MATLAB: vl_compilenn: CUDA: using NVCC 'C:\Program Files\NVIDIA GPU Computing Toolkit\CUDA\v10.1\bin\nvcc' (100205). MATLAB: vl_compilenn: CUDA: determining GPU compute capability (use the 'CudaArch' option to override) MATLAB: vl_compilenn: CUDA: cannot determine the capabilities of the installed GPU and/or CUDA; falling back to default MATLAB: vl_compilenn: CUDA: NVCC architecture string: '-gencode=arch=compute_20,code=\&quot;sm_20,compute_20\&quot; -gencode=arch=compute_30,code=\&quot;sm_30,compute_30\&quot;'. MATLAB: Warning: Support for toolboxdir('distcomp') will be removed in a future release. Use toolboxdir('parallel') instead.] MATLAB: &gt; In toolboxdir (line 49) MATLAB: In vl_compilenn (line 383) MATLAB: In segmfunc (line 63)] MATLAB: vl_compilenn: * Compiler and linker configurations * MATLAB: vl_compilenn: intermediate build products directory: C:\Users\gelme\REShAPE\Segmentation\matconvnet\matlab\mex\.build MATLAB: vl_compilenn: MEX files: C:\Users\gelme\REShAPE\Segmentation\matconvnet\matlab\mex/ MATLAB: vl_compilenn: Base options: -DENABLE_GPU -DENABLE_DOUBLE -O -DNDEBUG MATLAB: vl_compilenn: MEX CXX: -R2018a -D__SSSE3__ MATLAB: vl_compilenn: MEX CXXFLAGS: MATLAB: vl_compilenn: MEX CXXOPTIMFLAGS: MATLAB: vl_compilenn: MEX LINK: -R2018a -lmwblas -L&quot;C:\Program Files\NVIDIA GPU Computing Toolkit\CUDA\v10.1\lib\x64&quot; -lcudart -lcublas -lgpu MATLAB: vl_compilenn: MEX LINK LDFLAGS: MATLAB: vl_compilenn: MEX LINK LDOPTIMFLAGS: MATLAB: vl_compilenn: MEX LINK LINKLIBS: gdiplus.lib MATLAB: vl_compilenn: MEX CUDA: -R2018a MATLAB: vl_compilenn: MEX CUDA CXXFLAGS: MATLAB: vl_compilenn: MEX CUDA CXXOPTIMFLAGS: MATLAB: vl_compilenn: NVCC: -D_FORCE_INLINES --std=c++11 -I&quot;C:\Program Files\MATLAB\R2019b\extern\include&quot; -I&quot;C:\Program Files\MATLAB\R2019b\toolbox\parallel\gpu\extern\include&quot; -gencode=arch=compute_20,code=\&quot;sm_20,compute_20\&quot; -gencode=arch=compute_30,code=\&quot;sm_30,compute_30\&quot; --compiler-options=/MD MATLAB: vl_compilenn: * Reading images * MATLAB: vl_compilenn: vl_imreadjpeg enabled MATLAB: vl_compilenn: image library: gdiplus MATLAB: vl_compilenn: image library compile flags: MATLAB: vl_compilenn: image library link flags: gdiplus.lib MATLAB: vl_compilenn: NVCC CC: &quot;C:\Program Files\NVIDIA GPU Computing Toolkit\CUDA\v10.1\bin\nvcc&quot; -c -o -v &quot;C:\Users\gelme\REShAPE\Segmentation\matconvnet\matlab\mex\.build\bits\data.obj&quot; &quot;C:\Users\gelme\REShAPE\Segmentation\matconvnet\matlab\src\bits\data.cu&quot; -DENABLE_GPU -DENABLE_DOUBLE -O -DNDEBUG -D_FORCE_INLINES --std=c++11 -I&quot;C:\Program Files\MATLAB\R2019b\extern\include&quot; -I&quot;C:\Program Files\MATLAB\R2019b\toolbox\parallel\gpu\extern\include&quot; -gencode=arch=compute_20,code=\&quot;sm_20,compute_20\&quot; -gencode=arch=compute_30,code=\&quot;sm_30,compute_30\&quot; --compiler-options=/MD MATLAB: nvcc fatal : '-DNDEBUG': expected a number MATLAB: Command &quot;C:\Program Files\NVIDIA GPU Computing Toolkit\CUDA\v10.1\bin\nvcc&quot; -c -o -v &quot;C:\Users\gelme\REShAPE\Segmentation\matconvnet\matlab\mex\.build\bits\data.obj&quot; &quot;C:\Users\gelme\REShAPE\Segmentation\matconvnet\matlab\src\bits\data.cu&quot; -DENABLE_GPU -DENABLE_DOUBLE -O -DNDEBUG -D_FORCE_INLINES --std=c++11 -I&quot;C:\Program Files\MATLAB\R2019b\extern\include&quot; -I&quot;C:\Program Files\MATLAB\R2019b\toolbox\parallel\gpu\extern\include&quot; -gencode=arch=compute_20,code=\&quot;sm_20,compute_20\&quot; -gencode=arch=compute_30,code=\&quot;sm_30,compute_30\&quot; --compiler-options=/MD failed. </code></pre> <p>This seems to be an issue with the vl.compilenn.m showed in the image. It's still using nvcc when I am trying to use mex. How can I fix this issue that it's giving me in the image?</p> <p>This is the code I currently have for vl.compilenn</p> <pre><code>function vl_compilenn(varargin) %VL_COMPILENN Compile the MatConvNet toolbox. % opts.continue = false; opts.enableGpu = false; opts.enableImreadJpeg = true; opts.enableCudnn = false; opts.enableDouble = true; opts.imageLibrary = [] ; opts.imageLibraryCompileFlags = {} ; opts.imageLibraryLinkFlags = [] ; opts.verbose = 0; opts.debug = false; opts.cudaMethod = 'mex' ; opts.cudaRoot = [] ; opts.cudaArch = [] ; opts.defCudaArch = [... '-gencode=arch=compute_20,code=\&quot;sm_20,compute_20\&quot; '... '-gencode=arch=compute_30,code=\&quot;sm_30,compute_30\&quot;']; opts.mexConfig = '' ; opts.mexCudaConfig = '' ; opts.cudnnRoot = 'local/cudnn' ; opts.preCompileFn = [] ; % &lt;### Inserted by matconvnet_patch.py ###&gt; opts.mexFlags = {'-largeArrayDims'}; opts = vl_argparse(opts, varargin); % -------------------------------------------------------------------- % Files to compile % -------------------------------------------------------------------- arch = computer('arch') ; check_compability(arch); if isempty(opts.imageLibrary) switch arch case 'glnxa64', opts.imageLibrary = 'libjpeg' ; case 'maci64', opts.imageLibrary = 'quartz' ; case 'win64', opts.imageLibrary = 'gdiplus' ; end end if isempty(opts.imageLibraryLinkFlags) switch opts.imageLibrary case 'libjpeg', opts.imageLibraryLinkFlags = {'-ljpeg'} ; case 'quartz', opts.imageLibraryLinkFlags = {'-framework Cocoa -framework ImageIO'} ; case 'gdiplus', opts.imageLibraryLinkFlags = {'gdiplus.lib'} ; end end lib_src = {} ; mex_src = {} ; % Files that are compiled as CPP or CU depending on whether GPU support % is enabled. if opts.enableGpu, ext = 'cu' ; else ext='cpp' ; end lib_src{end+1} = fullfile(root,'matlab','src','bits',['data.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['datamex.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['nnconv.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['nnfullyconnected.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['nnsubsample.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['nnpooling.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['nnnormalize.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['nnnormalizelp.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['nnbnorm.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['nnbias.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['nnbilinearsampler.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src','bits',['nnroipooling.' ext]) ; mex_src{end+1} = fullfile(root,'matlab','src',['vl_nnconv.' ext]) ; mex_src{end+1} = fullfile(root,'matlab','src',['vl_nnconvt.' ext]) ; mex_src{end+1} = fullfile(root,'matlab','src',['vl_nnpool.' ext]) ; mex_src{end+1} = fullfile(root,'matlab','src',['vl_nnnormalize.' ext]) ; mex_src{end+1} = fullfile(root,'matlab','src',['vl_nnnormalizelp.' ext]) ; mex_src{end+1} = fullfile(root,'matlab','src',['vl_nnbnorm.' ext]) ; mex_src{end+1} = fullfile(root,'matlab','src',['vl_nnbilinearsampler.' ext]) ; mex_src{end+1} = fullfile(root,'matlab','src',['vl_nnroipool.' ext]) ; mex_src{end+1} = fullfile(root,'matlab','src',['vl_taccummex.' ext]) ; switch arch case {'glnxa64','maci64'} % not yet supported in windows mex_src{end+1} = fullfile(root,'matlab','src',['vl_tmove.' ext]) ; end % CPU-specific files lib_src{end+1} = fullfile(root,'matlab','src','bits','impl','im2row_cpu.cpp') ; lib_src{end+1} = fullfile(root,'matlab','src','bits','impl','copy_cpu.cpp') ; lib_src{end+1} = fullfile(root,'matlab','src','bits','impl','tinythread.cpp') ; lib_src{end+1} = fullfile(root,'matlab','src','bits','imread.cpp') ; % GPU-specific files if opts.enableGpu lib_src{end+1} = fullfile(root,'matlab','src','bits','impl','im2row_gpu.cu') ; lib_src{end+1} = fullfile(root,'matlab','src','bits','impl','copy_gpu.cu') ; lib_src{end+1} = fullfile(root,'matlab','src','bits','datacu.cu') ; mex_src{end+1} = fullfile(root,'matlab','src','vl_cudatool.cu') ; end % cuDNN-specific files if opts.enableCudnn end % Other files if opts.enableImreadJpeg mex_src{end+1} = fullfile(root,'matlab','src', ['vl_imreadjpeg.' ext]) ; mex_src{end+1} = fullfile(root,'matlab','src', ['vl_imreadjpeg_old.' ext]) ; lib_src{end+1} = fullfile(root,'matlab','src', 'bits', 'impl', ['imread_' opts.imageLibrary '.cpp']) ; end % -------------------------------------------------------------------- % Setup CUDA toolkit % -------------------------------------------------------------------- if opts.enableGpu opts.verbose &amp;&amp; fprintf('%s: * CUDA configuration *\n', mfilename) ; % Find the CUDA Devkit if isempty(opts.cudaRoot), opts.cudaRoot = search_cuda_devkit(opts) ; end opts.verbose &amp;&amp; fprintf('%s:\tCUDA: using CUDA Devkit ''%s''.\n', ... mfilename, opts.cudaRoot) ; opts.nvccPath = fullfile(opts.cudaRoot, 'bin', 'nvcc') ; switch arch case 'win64', opts.cudaLibDir = fullfile(opts.cudaRoot, 'lib', 'x64') ; case 'maci64', opts.cudaLibDir = fullfile(opts.cudaRoot, 'lib') ; case 'glnxa64', opts.cudaLibDir = fullfile(opts.cudaRoot, 'lib64') ; end % Set the nvcc method as default for Win platforms if strcmp(arch, 'win64') &amp;&amp; isempty(opts.cudaMethod) opts.cudaMethod = 'nvcc'; end % Activate the CUDA Devkit cuver = activate_nvcc(opts.nvccPath) ; opts.verbose &amp;&amp; fprintf('%s:\tCUDA: using NVCC ''%s'' (%d).\n', ... mfilename, opts.nvccPath, cuver) ; % Set the CUDA arch string (select GPU architecture) if isempty(opts.cudaArch), opts.cudaArch = get_cuda_arch(opts) ; end opts.verbose &amp;&amp; fprintf('%s:\tCUDA: NVCC architecture string: ''%s''.\n', ... mfilename, opts.cudaArch) ; end if opts.enableCudnn opts.cudnnIncludeDir = fullfile(opts.cudnnRoot, 'include') ; switch arch case 'win64', opts.cudnnLibDir = fullfile(opts.cudnnRoot, 'lib', 'x64') ; case 'maci64', opts.cudnnLibDir = fullfile(opts.cudnnRoot, 'lib') ; case 'glnxa64', opts.cudnnLibDir = fullfile(opts.cudnnRoot, 'lib64') ; end end % -------------------------------------------------------------------- % Compiler options % -------------------------------------------------------------------- % Build directories flags.src_dir = fullfile(root, 'matlab', 'src') ; flags.mex_dir = fullfile(root, 'matlab', 'mex') ; flags.bld_dir = fullfile(flags.mex_dir, '.build'); if ~exist(fullfile(flags.bld_dir,'bits','impl'), 'dir') mkdir(fullfile(flags.bld_dir,'bits','impl')) ; end % BASE: Base flags passed to `mex` and `nvcc` always. flags.base = {} ; if opts.enableGpu, flags.base{end+1} = '-DENABLE_GPU' ; end if opts.enableDouble, flags.base{end+1} = '-DENABLE_DOUBLE' ; end if opts.enableCudnn flags.base{end+1} = '-DENABLE_CUDNN' ; flags.base{end+1} = ['-I&quot;' opts.cudnnIncludeDir '&quot;'] ; end if opts.verbose &gt; 1, flags.base{end+1} = '-v' ; end if opts.debug flags.base{end+1} = '-g' ; flags.base{end+1} = '-DDEBUG' ; else flags.base{end+1} = '-O' ; flags.base{end+1} = '-DNDEBUG' ; end % MEX: Additional flags passed to `mex` for compiling C++ % code. CXX and CXXOPTIOM are passed directly to the encapsualted compiler. % &lt;### Changed by matconvnet_patch.py ###&gt; % flags.mex = {'-largeArrayDims'} ; flags.mex = opts.mexFlags; flags.cxx = {} ; flags.cxxoptim = {} ; if ~isempty(opts.mexConfig), flags.mex = horzcat(flags.mex, {'-f', opts.mexConfig}) ; end % MEX: Additional flags passed to `mex` for compiling CUDA % code. CXX and CXXOPTIOM are passed directly to the encapsualted compiler. % &lt;### Changed by matconvnet_patch.py ###&gt; % flags.mexcuda = {'-largeArrayDims'} ; flags.mexcuda = opts.mexFlags; flags.mexcuda_cxx = {} ; flags.mexcuda_cxxoptim = {} ; if ~isempty(opts.mexCudaConfig), flags.mexcuda = horzcat(flags.mexcuda, {'-f', opts.mexCudaConfig}) ; end % MEX_LINK: Additional flags passed to `mex` for linking. % &lt;### Changed by matconvnet_patch.py ###&gt; % flags.mexlink = {'-largeArrayDims','-lmwblas'} ; flags.mexlink = horzcat(opts.mexFlags, {'-lmwblas'}); flags.mexlink_ldflags = {} ; flags.mexlink_ldoptimflags = {} ; flags.mexlink_linklibs = {} ; % NVCC: Additional flags passed to `nvcc` for compiling CUDA code. flags.nvcc = {'-D_FORCE_INLINES', '--std=c++11', ... sprintf('-I&quot;%s&quot;',fullfile(matlabroot,'extern','include')), ... sprintf('-I&quot;%s&quot;',fullfile(toolboxdir('distcomp'),'gpu','extern','include')), ... opts.cudaArch} ; switch arch case {'maci64','glnxa64'} flags.cxx{end+1} = '--std=c++11' ; flags.nvcc{end+1} = '--compiler-options=-fPIC' ; if ~opts.debug flags.cxxoptim = horzcat(flags.cxxoptim,'-mssse3','-ffast-math') ; flags.mexcuda_cxxoptim{end+1} = '--compiler-options=-mssse3,-ffast-math' ; flags.nvcc{end+1} = '--compiler-options=-mssse3,-ffast-math' ; end case 'win64' % Visual Studio 2015 does C++11 without futher switches end if opts.enableGpu flags.mexlink = horzcat(flags.mexlink, ... {['-L&quot;' opts.cudaLibDir '&quot;'], '-lcudart', '-lcublas'}) ; switch arch case {'maci64', 'glnxa64'} flags.mexlink{end+1} = '-lmwgpu' ; case 'win64' flags.mexlink{end+1} = '-lgpu' ; end if opts.enableCudnn flags.mexlink{end+1} = ['-L&quot;' opts.cudnnLibDir '&quot;'] ; flags.mexlink{end+1} = '-lcudnn' ; end end switch arch case {'maci64'} flags.mex{end+1} = '-cxx' ; flags.nvcc{end+1} = '--compiler-options=-mmacosx-version-min=10.10' ; [s,r] = system('xcrun -f clang++') ; if s == 0 flags.nvcc{end+1} = sprintf('--compiler-bindir=&quot;%s&quot;',strtrim(r)) ; end if opts.enableGpu flags.mexlink_ldflags{end+1} = sprintf('-Wl,-rpath -Wl,&quot;%s&quot;', opts.cudaLibDir) ; end if opts.enableGpu &amp;&amp; opts.enableCudnn flags.mexlink_ldflags{end+1} = sprintf('-Wl,-rpath -Wl,&quot;%s&quot;', opts.cudnnLibDir) ; end case {'glnxa64'} flags.mex{end+1} = '-cxx' ; flags.mexlink{end+1} = '-lrt' ; if opts.enableGpu flags.mexlink_ldflags{end+1} = sprintf('-Wl,-rpath -Wl,&quot;%s&quot;', opts.cudaLibDir) ; end if opts.enableGpu &amp;&amp; opts.enableCudnn flags.mexlink_ldflags{end+1} = sprintf('-Wl,-rpath -Wl,&quot;%s&quot;', opts.cudnnLibDir) ; end case {'win64'} % VisualC does not pass this even if available in the CPU architecture flags.mex{end+1} = '-D__SSSE3__' ; cl_path = fileparts(check_clpath()); % check whether cl.exe in path flags.nvcc{end+1} = '--compiler-options=/MD' ; % flags.nvcc{end+1} = sprintf('--compiler-bindir=&quot;%s&quot;', cl_path) ; end if opts.enableImreadJpeg flags.mex = horzcat(flags.mex, opts.imageLibraryCompileFlags) ; flags.mexlink_linklibs = horzcat(flags.mexlink_linklibs, opts.imageLibraryLinkFlags) ; end if ~isempty(opts.preCompileFn) [opts, mex_src, lib_src, flags] = opts.preCompileFn(opts, mex_src, lib_src, flags) ; end % Compile intermediate object files srcs = horzcat(lib_src,mex_src) ; for i = 1:numel(horzcat(lib_src, mex_src)) [~,~,ext] = fileparts(srcs{i}) ; ext(1) = [] ; objfile = toobj(flags.bld_dir,srcs{i}); if strcmp(ext,'cu') if strcmp(opts.cudaMethod,'nvcc') nvcc_compile(opts, srcs{i}, objfile, flags) ; else mexcuda_compile(opts, srcs{i}, objfile, flags) ; end else mex_compile(opts, srcs{i}, objfile, flags) ; end assert(exist(objfile, 'file') ~= 0, 'Compilation of %s failed.', objfile); end % Link MEX files for i = 1:numel(mex_src) objs = toobj(flags.bld_dir, [mex_src(i), lib_src]) ; mex_link(opts, objs, flags.mex_dir, flags) ; end % Reset path adding the mex subdirectory just created vl_setupnn() ; if strcmp(arch, 'win64') &amp;&amp; opts.enableCudnn if opts.verbose(), fprintf('Copying CuDNN dll to mex folder.\n'); end copyfile(fullfile(opts.cudnnRoot, 'bin', '*.dll'), flags.mex_dir); end % Save the last compile flags to the build dir if isempty(opts.preCompileFn) save(fullfile(flags.bld_dir, 'last_compile_opts.mat'), '-struct', 'opts'); end % -------------------------------------------------------------------- % Utility functions % -------------------------------------------------------------------- % -------------------------------------------------------------------- function check_compability(arch) % -------------------------------------------------------------------- cc = mex.getCompilerConfigurations('C++'); if isempty(cc) error(['Mex is not configured.'... 'Run &quot;mex -setup&quot; to configure your compiler. See ',... 'http://www.mathworks.com/support/compilers ', ... 'for supported compilers for your platform.']); end switch arch case 'win64' clversion = str2double(cc.Version); if clversion &lt; 14 error('Unsupported VS C++ compiler, ver &gt;=14.0 required (VS 2015).'); end case 'maci64' case 'glnxa64' otherwise, error('Unsupported architecture ''%s''.', arch) ; end % -------------------------------------------------------------------- function done = check_deps(opts, tgt, src) % -------------------------------------------------------------------- done = false ; if ~iscell(src), src = {src} ; end if ~opts.continue, return ; end if ~exist(tgt,'file'), return ; end ttime = dir(tgt) ; ttime = ttime.datenum ; for i=1:numel(src) stime = dir(src{i}) ; stime = stime.datenum ; if stime &gt; ttime, return ; end end fprintf('%s: ''%s'' already there, skipping.\n', mfilename, tgt) ; done = true ; % -------------------------------------------------------------------- function objs = toobj(bld_dir, srcs) % -------------------------------------------------------------------- str = [filesep, 'src', filesep]; % NASTY. Do with regexp? multiple = iscell(srcs) ; if ~multiple, srcs = {srcs} ; end objs = cell(1, numel(srcs)); for t = 1:numel(srcs) i = strfind(srcs{t},str); i = i(end); % last occurence of '/src/' objs{t} = fullfile(bld_dir, srcs{t}(i+numel(str):end)) ; end if ~multiple, objs = objs{1} ; end objs = regexprep(objs,'.cpp$',['.' objext]) ; objs = regexprep(objs,'.cu$',['.' objext]) ; objs = regexprep(objs,'.c$',['.' objext]) ; % -------------------------------------------------------------------- function mex_compile(opts, src, tgt, flags) % -------------------------------------------------------------------- if check_deps(opts, tgt, src), return ; end args = horzcat({'-c', '-outdir', fileparts(tgt), src}, ... flags.base, flags.mex, ... {['CXXFLAGS=$CXXFLAGS ' strjoin(flags.cxx)]}, ... {['CXXOPTIMFLAGS=$CXXOPTIMFLAGS ' strjoin(flags.cxxoptim)]}) ; opts.verbose &amp;&amp; fprintf('%s: MEX CC: %s\n', mfilename, strjoin(args)) ; mex(args{:}) ; % -------------------------------------------------------------------- function mexcuda_compile(opts, src, tgt, flags) % -------------------------------------------------------------------- if check_deps(opts, tgt, src), return ; end % Hacky fix: In glnxa64 MATLAB includes the -ansi option by default, which % prevents -std=c++11 to work (an error?). This could be solved by editing the % mex configuration file; for convenience, we take it out here by % avoiding to append to the default flags. glue = '$CXXFLAGS' ; switch computer('arch') case {'glnxa64'} glue = '--compiler-options=-fexceptions,-fPIC,-fno-omit-frame-pointer,-pthread' ; end args = horzcat({'-c', '-outdir', fileparts(tgt), src}, ... flags.base, flags.mexcuda, ... {['CXXFLAGS=' glue ' ' strjoin(flags.mexcuda_cxx)]}, ... {['CXXOPTIMFLAGS=$CXXOPTIMFLAGS ' strjoin(flags.mexcuda_cxxoptim)]}) ; opts.verbose &amp;&amp; fprintf('%s: MEX CUDA: %s\n', mfilename, strjoin(args)) ; mexcuda(args{:}) ; % -------------------------------------------------------------------- function nvcc_compile(opts, src, tgt, flags) % -------------------------------------------------------------------- if check_deps(opts, tgt, src), return ; end nvcc_path = fullfile(opts.cudaRoot, 'bin', 'nvcc'); nvcc_cmd = sprintf('&quot;%s&quot; -c -o -v &quot;%s&quot; &quot;%s&quot; %s ', ... nvcc_path, tgt, src, ... strjoin(horzcat(flags.base,flags.nvcc))); opts.verbose &amp;&amp; fprintf('%s: NVCC CC: %s\n', mfilename, nvcc_cmd) ; status = system(nvcc_cmd); if status, error('Command %s failed.', nvcc_cmd); end; % -------------------------------------------------------------------- function mex_link(opts, objs, mex_dir, flags) % -------------------------------------------------------------------- args = horzcat({'-outdir', mex_dir}, ... flags.base, flags.mexlink, ... {['LDFLAGS=$LDFLAGS ' strjoin(flags.mexlink_ldflags)]}, ... {['LDOPTIMFLAGS=$LDOPTIMFLAGS ' strjoin(flags.mexlink_ldoptimflags)]}, ... {['LINKLIBS=' strjoin(flags.mexlink_linklibs) ' $LINKLIBS']}, ... objs) ; opts.verbose &amp;&amp; fprintf('%s: MEX LINK: %s\n', mfilename, strjoin(args)) ; mex(args{:}) ; % -------------------------------------------------------------------- function ext = objext() % -------------------------------------------------------------------- % Get the extension for an 'object' file for the current computer % architecture switch computer('arch') case 'win64', ext = 'obj'; case {'maci64', 'glnxa64'}, ext = 'o' ; otherwise, error('Unsupported architecture %s.', computer) ; end % -------------------------------------------------------------------- function cl_path = check_clpath() % -------------------------------------------------------------------- % Checks whether the cl.exe is in the path (needed for the nvcc). If % not, tries to guess the location out of mex configuration. cc = mex.getCompilerConfigurations('c++'); cl_path = fullfile(cc.Location, 'VC', 'bin', 'amd64'); [status, ~] = system('cl.exe -help'); if status == 1 % Add cl.exe to system path so that nvcc can find it. warning('CL.EXE not found in PATH. Trying to guess out of mex setup.'); prev_path = getenv('PATH'); setenv('PATH', [prev_path ';' cl_path]); status = system('cl.exe'); if status == 1 setenv('PATH', prev_path); error('Unable to find cl.exe'); else fprintf('Location of cl.exe (%s) successfully added to your PATH.\n', ... cl_path); end end % ------------------------------------------------------------------------- function paths = which_nvcc() % ------------------------------------------------------------------------- switch computer('arch') case 'win64' [~, paths] = system('where nvcc.exe'); paths = strtrim(paths); paths = paths(strfind(paths, '.exe')); case {'maci64', 'glnxa64'} [~, paths] = system('which nvcc'); paths = strtrim(paths) ; end % ------------------------------------------------------------------------- function cuda_root = search_cuda_devkit(opts) % ------------------------------------------------------------------------- % This function tries to to locate a working copy of the CUDA Devkit. opts.verbose &amp;&amp; fprintf(['%s:\tCUDA: searching for the CUDA Devkit' ... ' (use the option ''CudaRoot'' to override):\n'], mfilename); % Propose a number of candidate paths for NVCC paths = {getenv('MW_NVCC_PATH')} ; paths = [paths, which_nvcc()] ; for v = {'5.5', '6.0', '6.5', '7.0', '7.5', '8.0', '8.5', '9.0', '9.5', '10.0'} switch computer('arch') case 'glnxa64' paths{end+1} = sprintf('/usr/local/cuda-%s/bin/nvcc', char(v)) ; case 'maci64' paths{end+1} = sprintf('/Developer/NVIDIA/CUDA-%s/bin/nvcc', char(v)) ; case 'win64' paths{end+1} = sprintf('C:\\Program Files\\NVIDIA GPU Computing Toolkit\\CUDA\\v%s\\bin\\nvcc.exe', char(v)) ; end end paths{end+1} = sprintf('/usr/local/cuda/bin/nvcc') ; % Validate each candidate NVCC path for i=1:numel(paths) nvcc(i).path = paths{i} ; [nvcc(i).isvalid, nvcc(i).version] = validate_nvcc(paths{i}) ; end if opts.verbose fprintf('\t| %5s | %5s | %-70s |\n', 'valid', 'ver', 'NVCC path') ; for i=1:numel(paths) fprintf('\t| %5d | %5d | %-70s |\n', ... nvcc(i).isvalid, nvcc(i).version, nvcc(i).path) ; end end % Pick an entry index = find([nvcc.isvalid]) ; if isempty(index) error('Could not find a valid NVCC executable\n') ; end [~, newest] = max([nvcc(index).version]); nvcc = nvcc(index(newest)) ; cuda_root = fileparts(fileparts(nvcc.path)) ; if opts.verbose fprintf('%s:\tCUDA: choosing NVCC compiler ''%s'' (version %d)\n', ... mfilename, nvcc.path, nvcc.version) ; end % ------------------------------------------------------------------------- function [valid, cuver] = validate_nvcc(nvccPath) % ------------------------------------------------------------------------- [status, output] = system(sprintf('&quot;%s&quot; --version', nvccPath)) ; valid = (status == 0) ; if ~valid cuver = 0 ; return ; end match = regexp(output, 'V(\d+\.\d+\.\d+)', 'match') ; if isempty(match), valid = false ; return ; end cuver = [1e4 1e2 1] * sscanf(match{1}, 'V%d.%d.%d') ; % -------------------------------------------------------------------- function cuver = activate_nvcc(nvccPath) % -------------------------------------------------------------------- % Validate the NVCC compiler installation [valid, cuver] = validate_nvcc(nvccPath) ; if ~valid error('The NVCC compiler ''%s'' does not appear to be valid.', nvccPath) ; end % Make sure that NVCC is visible by MEX by setting the MW_NVCC_PATH % environment variable to the NVCC compiler path if ~strcmp(getenv('MW_NVCC_PATH'), nvccPath) warning('Setting the ''MW_NVCC_PATH'' environment variable to ''%s''', nvccPath) ; setenv('MW_NVCC_PATH', nvccPath) ; end % In some operating systems and MATLAB versions, NVCC must also be % available in the command line search path. Make sure that this is% % the case. [valid_, cuver_] = validate_nvcc('nvcc') ; if ~valid_ || cuver_ ~= cuver warning('NVCC not found in the command line path or the one found does not matches ''%s''.', nvccPath); nvccDir = fileparts(nvccPath) ; prevPath = getenv('PATH') ; switch computer case 'PCWIN64', separator = ';' ; case {'GLNXA64', 'MACI64'}, separator = ':' ; end setenv('PATH', [nvccDir separator prevPath]) ; [valid_, cuver_] = validate_nvcc('nvcc') ; if ~valid_ || cuver_ ~= cuver setenv('PATH', prevPath) ; error('Unable to set the command line path to point to ''%s'' correctly.', nvccPath) ; else fprintf('Location of NVCC (%s) added to your command search PATH.\n', nvccDir) ; end end % -------------------------------------------------------------------- function cudaArch = get_cuda_arch(opts) % -------------------------------------------------------------------- opts.verbose &amp;&amp; fprintf('%s:\tCUDA: determining GPU compute capability (use the ''CudaArch'' option to override)\n', mfilename); try gpu_device = gpuDevice(); arch = str2double(strrep(gpu_device.ComputeCapability, '.', '')); supparchs = get_nvcc_supported_archs(opts.nvccPath); [~, archi] = max(min(supparchs - arch, 0)); arch_code = num2str(supparchs(archi)); assert(~isempty(arch_code)); cudaArch = ... sprintf('-gencode=arch=compute_%s,code=\\\&quot;sm_%s,compute_%s\\\&quot; ', ... arch_code, arch_code, arch_code) ; catch opts.verbose &amp;&amp; fprintf(['%s:\tCUDA: cannot determine the capabilities of the installed GPU and/or CUDA; ' ... 'falling back to default\n'], mfilename); cudaArch = opts.defCudaArch; end % -------------------------------------------------------------------- function archs = get_nvcc_supported_archs(nvccPath) % -------------------------------------------------------------------- switch computer('arch') case {'win64'} [status, hstring] = system(sprintf('&quot;%s&quot; --help',nvccPath)); otherwise % fix possible output corruption (see manual) [status, hstring] = system(sprintf('&quot;%s&quot; --help &lt; /dev/null',nvccPath)) ; end archs = regexp(hstring, '''sm_(\d{2})''', 'tokens'); archs = cellfun(@(a) str2double(a{1}), archs); if status, error('NVCC command failed: %s', hstring); end; % Force the use of x64 architecture in the mex command if opts.enableGpu % If GPU is enabled, set the machine architecture explicitly to x64 flags.mexcc = [flags.mexcc ' -v COMPFLAGS=&quot;$COMPFLAGS /machine:x64&quot;']; else % For CPU-only, also ensure x64 architecture flags.mexcc = [flags.mexcc ' -v COMPFLAGS=&quot;$COMPFLAGS /machine:x64&quot;']; end % Ensure that no ARM64 libraries are used; focus on x64 libraries only. flags.link = [flags.link ' -L&quot;C:\path\to\x64\libs&quot;']; % Ensure correct x64 library paths flags.lib = [flags.lib ' -llibmx' ' -llibmex' ' -llibmat']; % Ensure correct x64 libraries are linked </code></pre>

https://stackoverflow.com/q/78868040 0 Reesh https://stackoverflow.com/users/26752151 2024-08-13T19:35:15Z 2024-08-18T14:27:08Z

<p>I am trying to follow this <a href="https://github.com/NIH-NEI/REShAPE" rel="nofollow noreferrer">github</a> and to open the code I am running into errors with the mex compiler:</p> <pre><code>% Apply pre-compilation modifier function to adjust the flags and % parameters. This can be used to add additional files to compile on the % fly. if ~isempty(opts.preCompileFn) [opts, mex_src, lib_src, flags] = opts.preCompileFn(opts, mex_src, lib_src, flags) ; end % Compile intermediate object files srcs = horzcat(lib_src,mex_src) ; for i = 1:numel(horzcat(lib_src, mex_src)) [~,~,ext] = fileparts(srcs{i}) ; ext(1) = [] ; objfile = toobj(flags.bld_dir,srcs{i}); if strcmp(ext,'cu') if strcmp(opts.cudaMethod,'nvcc') nvcc_compile(opts, srcs{i}, objfile, flags) ; else mexcuda_compile(opts, srcs{i}, objfile, flags) ; end else mex_compile(opts, srcs{i}, objfile, flags) ; end assert(exist(objfile, 'file') ~= 0, 'Compilation of %s failed.', objfile); end % Link MEX files for i = 1:numel(mex_src) objs = toobj(flags.bld_dir, [mex_src(i), lib_src]) ; mex_link(opts, objs, flags.mex_dir, flags) ; end % Reset path adding the mex subdirectory just created vl_setupnn() ; if strcmp(arch, 'win64') &amp;&amp; opts.enableCudnn if opts.verbose(), fprintf('Copying CuDNN dll to mex folder.\n'); end copyfile(fullfile(opts.cudnnRoot, 'bin', '*.dll'), flags.mex_dir); end % Save the last compile flags to the build dir if isempty(opts.preCompileFn) save(fullfile(flags.bld_dir, 'last_compile_opts.mat'), '-struct', 'opts'); end % -------------------------------------------------------------------- % Utility functions % -------------------------------------------------------------------- % -------------------------------------------------------------------- function check_compability(arch) % -------------------------------------------------------------------- cc = mex.getCompilerConfigurations('C++'); if isempty(cc) error(['Mex is not configured.'... 'Run &quot;mex -setup&quot; to configure your compiler. See ',... 'http://www.mathworks.com/support/compilers ', ... 'for supported compilers for your platform.']); end switch arch case 'win64' clversion = str2double(cc.Version); if clversion &lt; 14 error('Unsupported VS C++ compiler, ver &gt;=14.0 required (VS 2015).'); end case 'maci64' case 'glnxa64' otherwise, error('Unsupported architecture ''%s''.', arch) ; end % -------------------------------------------------------------------- function done = check_deps(opts, tgt, src) % -------------------------------------------------------------------- done = false ; if ~iscell(src), src = {src} ; end if ~opts.continue, return ; end if ~exist(tgt,'file'), return ; end ttime = dir(tgt) ; ttime = ttime.datenum ; for i=1:numel(src) stime = dir(src{i}) ; stime = stime.datenum ; if stime &gt; ttime, return ; end end fprintf('%s: ''%s'' already there, skipping.\n', mfilename, tgt) ; done = true ; % -------------------------------------------------------------------- function objs = toobj(bld_dir, srcs) % -------------------------------------------------------------------- str = [filesep, 'src', filesep]; % NASTY. Do with regexp? multiple = iscell(srcs) ; if ~multiple, srcs = {srcs} ; end objs = cell(1, numel(srcs)); for t = 1:numel(srcs) i = strfind(srcs{t},str); i = i(end); % last occurence of '/src/' objs{t} = fullfile(bld_dir, srcs{t}(i+numel(str):end)) ; end if ~multiple, objs = objs{1} ; end objs = regexprep(objs,'.cpp$',['.' objext]) ; objs = regexprep(objs,'.cu$',['.' objext]) ; objs = regexprep(objs,'.c$',['.' objext]) ; % -------------------------------------------------------------------- function mex_compile(opts, src, tgt, flags) % -------------------------------------------------------------------- if check_deps(opts, tgt, src), return ; end args = horzcat({'-c', '-outdir', fileparts(tgt), src}, ... flags.base, flags.mex, ... {['CXXFLAGS=$CXXFLAGS ' strjoin(flags.cxx)]}, ... {['CXXOPTIMFLAGS=$CXXOPTIMFLAGS ' strjoin(flags.cxxoptim)]}) ; opts.verbose &amp;&amp; fprintf('%s: MEX CC: %s\n', mfilename, strjoin(args)) ; mex(args{:}) ; </code></pre> <p>specifically this error in this block of code:</p> <pre><code>Error in vl_compilenn&gt;mex_compile (line 599) mex(args{:}) ; Error in vl_compilenn (line 508) mex_compile(opts, srcs{i}, objfile, flags) ; </code></pre> <p>How do I resolve this issue?</p> <p>We tried to use 'nvcc' and switched over to 'mex'. We haven't been able to figure it out for months. We also installed matlab compiler MinGW-w64 C/C++ compiler.</p> <p>Here is more of the error message: and that works here is more of the error message:</p> <pre><code>`:\Users\gelme\REShAPE\Segmentation\matconvnet\matlab\src\bits\impl\imread_gdiplus.cpp:14: C:/ProgramData/MATLAB/SupportPackages/R2019b/3P.instrset/mingw_w64.instrset/lib/gcc/x86_64-w64-mingw32/6.3.0/include/tmmintrin.h: In function '__m128i _mm_hadd_epi16(__m128i, __m128i)': C:/ProgramData/MATLAB/SupportPackages/R2019b/3P.instrset/mingw_w64.instrset/lib/gcc/x86_64-w64-mingw32/6.3.0/include/tmmintrin.h:42:70: error: '__builtin_ia32_phaddw128' was not declared in this scope return (__m128i) __builtin_ia32_phaddw128 ((__v8hi)__X, (__v8hi)__Y); ^ C:/ProgramData/MATLAB/SupportPackages/R2019b/3P.instrset/mingw_w64.instrset/lib/gcc/x86_64-w64-mingw32/6.3.0/include/tmmintrin.h: In function '__m128i _mm_hadd_epi32(__m128i, __m128i)': C:/ProgramData/MATLAB/SupportPackages/R2019b/3P.instrset/mingw_w64.instrset/lib/gcc/x86_64-w64-mingw32/6.3.0/include/tmmintrin.h:48:70: error: '__builtin_ia32_phaddd128' was not declared in this scope return (__m128i) __builtin_ia32_phaddd128 ((__v4si)__X, (__v4si)__Y); ^ C:/ProgramData/MATLAB/SupportPackages/R2019b/3P.instrset/mingw_w64.instrset/lib/gcc/x86_64-w64-mingw32/6.3.0/include/tmmintrin.h: In function '__m128i _mm_hadds_epi16(__m128i, __m128i)': C:/ProgramData/MATLAB/SupportPackages/R2019b/3P.instrset/mingw_w64.instrset/lib/gcc/x86_64-w64-mingw32/6.3.0/include/tmmintrin.h:54:71: error: '__builtin_ia32_phaddsw128' was not declared in this scope return (__m128i) __builtin_ia32_phaddsw128 ((__v8hi)__X, (__v8hi)__Y); ^ C:/ProgramData/MATLAB/SupportPackages/R2019b/3P.instrset/mingw_w64.instrset/lib/gcc/x86_64-w64-mingw32/6.3.0/include/tmmintrin.h: In function '__m64 _mm_hadd_pi16(__m64, __m64)': C:/ProgramData/MATLAB/SupportPackages/R2019b/3P.instrset/mingw_w64.instrset/lib/gcc/x86_64-w64-mingw32/6.3.0/include/tmmintrin.h:60:65: error: '__builtin_ia32_phaddw' was not declared in this scope return (__m64) __builtin_ia32_phaddw ((__v4hi)__X, (__v4hi)__Y);` ^ </code></pre>

https://stackoverflow.com/q/78439685 0 Jimbo https://stackoverflow.com/users/764365 2024-05-07T01:07:48Z 2024-08-14T11:01:49Z

<p>I have a matrix of bytes, or for this example we can assume I've converted to uint16. I am curious how to copy a subset of the data elements into a MATLABString (i.e., using the mex C++ interface). I asked ChatGPT but it is using a charArray and not a string. It also looks like it has an unnecessary data copy going on at the end (I could be wrong about that). It also is not clear why I need an iterator but that may be because I am used to C and working more with pointers. Is this the best approach? The original solution had a loop to copy/cast each element but I asked it to use <code>memcpy</code> as I'm under the impression that <code>memcpy</code> may be faster and also just seems like it should be the right way to copy a lot of data from one array to another.</p> <pre><code>#include &quot;mex.hpp&quot; #include &quot;mexAdapter.hpp&quot; class MexFunction : public matlab::mex::Function { public: void operator()(matlab::mex::ArgumentList outputs, matlab::mex::ArgumentList inputs) override { // Check input and output arguments checkArguments(outputs, inputs); // Get input arguments matlab::data::TypedArray&lt;uint16_t&gt; data = std::move(inputs[0]); double startIndex = inputs[1][0]; double endIndex = inputs[2][0]; size_t subsetLength = static_cast&lt;size_t&gt;(endIndex - startIndex + 1); // Check bounds size_t dataSize = data.getNumberOfElements(); if (startIndex &lt; 1 || endIndex &gt;= dataSize || startIndex &gt; endIndex) { matlab::engine::MATLABException::throwAsCaller(&quot;Subset indices are out of bounds.&quot;); } // Extract subset of uint16 array matlab::data::TypedArray&lt;uint16_t&gt;::iterator subsetData = data.begin() + static_cast&lt;size_t&gt;(startIndex) - 1; // Convert subset to character array using memcpy std::string charArray(subsetLength, '\0'); memcpy(&amp;charArray[0], subsetData, subsetLength * sizeof(uint16_t)); // Create MATLABString outputs[0] = factory.createCharArray(charArray); } private: matlab::data::ArrayFactory factory; void checkArguments(matlab::mex::ArgumentList outputs, matlab::mex::ArgumentList inputs) { if (inputs.size() != 3) { matlab::engine::MATLABException::throwAsCaller(&quot;Three input arguments required.&quot;); } if (outputs.size() != 1) { matlab::engine::MATLABException::throwAsCaller(&quot;One output argument required.&quot;); } if (!inputs[0].isUint16() || !inputs[1].isDouble() || !inputs[2].isDouble()) { matlab::engine::MATLABException::throwAsCaller(&quot;Inputs must be uint16 and double arrays.&quot;); } } }; </code></pre>

https://stackoverflow.com/q/77986687 1 Shinrei https://stackoverflow.com/users/7617574 2024-02-13T09:06:48Z 2024-02-13T09:06:48Z

<p>I want to call Fortran subroutines from Matlab and get the output from them. I'm using Windows 7 x64 and Matlab R2021b.</p> <p>I have tried compiling in Matlab using mex command. At</p> <pre><code>mex -setup FORTRAN </code></pre> <p>I have choosen Intel Parallel Studio XE 2019 for Fortran with Microsoft Visual Studio 2015.</p> <p>In mex settings: 'C:\Program Files\MATLAB\R2021b\bin\win64\mexopts\intel_fortran_19_vs2015.xml' I have removed \fixed option from COMPFLAGS. It appeared that it is not able to recognise f95 files:</p> <pre><code>Unknown file extension '.f95'. </code></pre> <p>After changing extension of fortran files from f95 to f90 and trying to compile I got error:</p> <pre><code>error #5078: Unrecognized token '@' skipped result_logic = FEXISTS@(filename,errorcode) </code></pre> <p>I'm not sure if this is the problem of compiler, which is probably not supporting Silver-frost FTN95, or I have to include FTN95 Library in some way.</p>

https://stackoverflow.com/q/77961724 1 Isabella Mariani https://stackoverflow.com/users/23365537 2024-02-08T12:38:41Z 2024-02-10T21:21:12Z

<p>I get this error. I've tried to fix it following the guidelines online but it seems impossible to solve. Can someone help me please?</p> <blockquote> <p>Error using spm_check_installation&gt;check_basic<br /> SPM uses a number of MEX files, which are compiled functions. These need to be compiled for the various platforms on which SPM is run. It seems that the compiled files for your computer platform are missing or not compatible. See <a href="https://en.wikibooks.org/wiki/SPM/Installation_on_64bit_Mac_OS_(Intel)" rel="nofollow noreferrer">https://en.wikibooks.org/wiki/SPM/Installation_on_64bit_Mac_OS_(Intel)</a> for information about how to compile MEX files for MACI64 in MATLAB 23.2.0.2485118 (R2023b) Update 6.</p> <pre><code>Error in spm_check_installation (line 28) check_basic; Error in spm (line 290) spm_check_installation('basic'); </code></pre> </blockquote> <p>I tried to follow <a href="https://en.wikibooks.org/wiki/SPM/Installation_on_64bit_Mac_OS_(Intel)and" rel="nofollow noreferrer">https://en.wikibooks.org/wiki/SPM/Installation_on_64bit_Mac_OS_(Intel)and</a> also this command <code>xattr -dr com.apple.quarantine *.mexmaci64</code> suggested by Chris Rorden.</p>

https://stackoverflow.com/q/77760595 1 Bryan https://stackoverflow.com/users/22869820 2024-01-04T18:58:43Z 2024-01-05T18:35:50Z

<p>My current MEX file is properly compiling, however, I would like to create let's say several MEX files:</p> <ul> <li>mex_a.cpp</li> <li>mex_b.cpp</li> <li>mex_c.cpp</li> </ul> <p>At this moment, mex_a.cpp contains the typical function definition:</p> <pre><code>class MexFunction : public matlab::mex::Function { private: // Get pointer to engine std::shared_ptr&lt;matlab::engine::MATLABEngine&gt; matlabPtr = getEngine(); // Get array factory std::shared_ptr&lt;matlab::data::ArrayFactory&gt; factoryPtr; private: void operator()(matlab::mex::ArgumentList outputs, matlab::mex::ArgumentList inputs) override { try { // Parse inputs.. mex_aux-&gt;mex_common_function(...) // Parse outputs.. } catch (int) { std::fprintf(stdout, &quot;ERROR: Something unexpected happened...&quot;); } } </code></pre> <p>In <code>mex_aux</code>, we can find the definition and implementation of <code>mex_common_function</code>. The purpose of this object is to be used by mex_a.cpp, mex_b.cpp and mex_c.cpp, basically what a common library do.</p> <p>My findings are:</p> <ul> <li>If mex_common_function is static, everything works.</li> <li>Otherwise, it doesn't work while using matlab objects.</li> </ul> <p>The .hpp file looks like:</p> <pre><code>#pragma once // Include system libraries # include &lt;string&gt; # include &lt;utility&gt; // MATLAB stuff #include &quot;mex.hpp&quot; class mex_aux { public: // Constructor mex_aux() = default; // Destructor ~mex_aux() = default; private: // Get pointer to engine std::shared_ptr&lt;matlab::engine::MATLABEngine&gt; matlabPtr; // Get array factory std::shared_ptr&lt;matlab::data::ArrayFactory&gt; factoryPtr; public: // Setters: void set_matlab_ptr(std::shared_ptr&lt;matlab::engine::MATLABEngine&gt;&amp; matlab_ptr) {this-&gt;matlabPtr = matlab_ptr; } void set_factory_ptr(std::shared_ptr&lt;matlab::data::ArrayFactory&gt;&amp; factory_ptr) {this-&gt;factoryPtr = factory_ptr; } public: // The used function, if this one is set as static (hence not using members of the class) it works. void mex_common_function(); </code></pre> <p>What I am doing now, and hence getting the following error in MATLAB (not during building time) is:</p> <ul> <li>Setting <code>matlabPtr</code> as an internal pointer in <code>mex_aux</code>.</li> <li>Same for <code>factoryPtr</code>.</li> </ul> <p>This both from the main mex function file (i.e., mex_a.cpp) after building <code>mex_aux</code>, the setters are used.</p> <p>Error in MATLAB:</p> <pre><code>Invalid MEX-file '/path/to/mex/mexfunction.mexa64': /path/to/external/libextern.so: undefined symbol: _ZN6matlab6engine12MATLABEngine5fevalERKNSt7__cxx1112basic_stringIDsSt11char_traitsIDsESaIDsEEEiRKSt6vectorINS_4data5ArrayESaISC_EERKSt10shared_ptrISt15basic_streambufIDsS5_EESM_ Error in mex_vsod (line 50) b(:,:,i) = mexfunction(...); </code></pre> <p>My CMakeLists.txt:</p> <pre><code># MEX AUX Library set(LIBRARY_MEX_AUX &quot;mex_aux&quot;) add_library(${LIBRARY_MEX_AUX} SHARED src/core/tools/mex_aux.cpp) target_link_libraries(${LIBRARY_MEX_AUX} ${Matlab_MEX_LIBRARY} ${Matlab_MX_LIBRARY}) set_target_properties(${LIBRARY_MEX_AUX} PROPERTIES COMPILE_FLAGS &quot;-fPIC&quot; LINK_FLAGS &quot;-Wl,-rpath,./&quot;) # To use relative paths in shared libs # MEX FILES if (BUILD_WITH_MATLAB) # MEX functions files matlab_add_mex(NAME mex_vsaod SRC src/main/mex_cpp/mex_a.cpp LINK_TO ${LIBRARY_target1} ... ${LIBRARY_targetN} ${LIBRARY_MEX_AUX}) endif () </code></pre> <p>I guess something wrong is happening within the memory due to a pointer from <code>mex_aux</code> pointing possibly to an un-accessible place in memory (only scoped to the mex file -mex_a.cpp?-)...</p> <p>Can someone give me a clue? I would really appreciate it.</p>

https://stackoverflow.com/q/77434788 0 Bryan https://stackoverflow.com/users/22869820 2023-11-06T23:37:49Z 2023-11-07T13:50:07Z

<p>I created a MEX executable and this one can be executed as a matlab function only if I <code>export LD_PRELOAD</code> all the needed libraries.</p> <p>I wanted to automate this procedure cause it is using a bunch of in house made libraries and it's a tedious task to write <code>export LD_PRELOAD</code> every time in the terminal.</p> <p>So I did a bash script:</p> <pre><code>#!/bin/bash echo &quot;PRELOAD MATLAB SCRIPT&quot; # Specify the folder containing the .so files so_folder=&quot;/path/to/lib&quot; # Check if the folder exists if [ ! -d &quot;$so_folder&quot; ]; then echo &quot;The folder '$so_folder' does not exist.&quot; exit 1 fi echo &quot;Searching in .so folder: $so_folder ...&quot; # STD library problem from LINUX users only export LD_PRELOAD=$LD_PRELOAD:/lib/x86_64-linux-gnu/libstdc++.so.6 matlab # Set DACE library export LD_PRELOAD=$LD_PRELOAD:/path/to/dace/lib/libdace.so matlab # Loop through .so files in the folder for so_file in &quot;$so_folder&quot;/*.so; do if [ -f &quot;$so_file&quot; ]; then # Set LD_PRELOAD and launch MATLAB for each .so file export LD_PRELOAD=&quot;$LD_PRELOAD:$so_file&quot; matlab echo &quot;Preloading library $so_file to MATLAB...&quot; fi done </code></pre> <p>Then, I call this script from the root of the project, where I also launch matlab afterwards:</p> <pre><code>$ bash scripts/preload_matlab_libs.sh </code></pre> <p>After it, I launch MATLAB and the .so files are not found as the error shows.</p> <pre><code>&gt;&gt; mex_test Invalid MEX-file '/path/to/project/mex/mex_vsod.mexa64': libdace.so.2: cannot open shared object file: No such file or directory Error in mex_test (line 13) b = mex_vsod(input1, input2) </code></pre> <p>No error is happening if I do the <code>export LD_PRELOAD</code>...</p> <p>I'm a bit new in MEX compiled files. Therefore my questions:</p> <ul> <li>Why is the bash script not working? How do I fix it?</li> <li>Is there any other better approach to solve this issue?</li> </ul> <p>Thanks!</p>

https://stackoverflow.com/q/77184611 1 kun zhao https://stackoverflow.com/users/12611025 2023-09-27T04:01:13Z 2023-09-27T20:20:45Z

<p>I am using Matlab 2023a to call Fortran source file using mex. I use the following compiler for Fortran language compilation:</p> <blockquote> <p>Intel oneAPI 2023 for Fortran with Microsoft Visual Studio 2022.</p> </blockquote> <p>The code I used to compile is</p> <blockquote> <p>mex ero_dep_fortran.f90 -output ero_dep_fortran</p> </blockquote> <p>However, I got the following erros:</p> <blockquote> <p>Error using mex C:\Users\40967\Desktop\Paper bank collapse and meander migration\Draft-version2\Code\ero_dep_fortran.lib �Ͷ��� C:\Users\40967\Desktop\Paper bank collapse and meander migration\Draft-version2\Code\ero_dep_fortran.exp ero_dep_fortran.obj : error LNK2019: �޷��������ⲿ���� NEW_EPOINTS_ASSIGNMENT���÷����ں��� MEXFUNCTION �б����� ero_dep_fortran.obj : error LNK2019: �޷��������ⲿ���� omp_get_wtime���÷����ں��� BANK_EROSION_CALCULATION_ENTRY �б����� C:\Users\40967\Desktop\Paper bank collapse and meander migration\�޸ĺ����\Draft-version2\Code\ero_dep_fortran.mexw64 : fatal error LNK1120: 2���޷��������ⲿ����</p> </blockquote> <p>A lot of messy code in the error. I found the main problem maybe</p> <blockquote> <p>ero_dep_fortran.obj : error LNK2019</p> </blockquote> <p>but I do not know how to slove it...This code can be run in other computers with lower version of Matlab and Intel Parallel Studio XE, Microsoft Visual Studio.</p> <p>Any suggestion is grateful....</p>

https://stackoverflow.com/q/77140911 0 MarryD https://stackoverflow.com/users/22598663 2023-09-20T09:06:44Z 2023-09-20T16:20:06Z

<p>I am trying to run example provided by field II on my computer. This example contains matlab functions as well as mex32 file. I have later version of matlab so i take initial .c file and converted it into mex64. However despite setting everything like in provided example my mex function does not work. I am not that familiar with c programming so i am not sure where exacly in the main function problem occurs. If anyone please could direct me what is going wrong in this .c program. Also, can i succesfully transform this c code to matlab function instead of using mex file or is there some nuances in the C code that can not be fully captured in the MATLAB translation?</p> <p><a href="https://i.sstatic.net/EJcDw.png" rel="nofollow noreferrer">Matlab error</a> <a href="https://field-ii.dk/examples/ftp_files/kidney/fast_int.c" rel="nofollow noreferrer">full original C code</a></p> <p>I was trying to write matlab function based on algorithm from C program but i am not sure if i am getting indexing correctly. I am aware that matlab indices starts from 1 not 0. What i don't understand is this env_pointer being set to be one variable's address and then is called like it contains N_samples+1 variables.</p> <pre><code>void make_interpolation (unsigned char *envelope_data, /* The envelope detected and log-compressed data */ int N_samples, /* Number of samples in one envelope line */ int *index_samp_line, /* Index for the data sample number */ int *image_index, /* Index for the image matrix */ float *weight_coef, /* The weight table */ int N_values, /* Number of values to calculate in the image */ unsigned char *image) /* The resulting image */ {int i; /* Integer loop counter */ int ij_index_coef; /* Index into coefficient array */ unsigned char *env_pointer; /* Pointer to the envelope data */ float *weight_pointer; /* Pointer to the weight coefficients */ ij_index_coef = 0; for (i=0; i&lt;N_values; i++) { weight_pointer = &amp;(weight_coef[ij_index_coef]); env_pointer = &amp;(envelope_data[index_samp_line[i]]); image[image_index[i]] = weight_pointer[0] * env_pointer[0] + weight_pointer[1] * env_pointer[1] + weight_pointer[2] * env_pointer[N_samples] + weight_pointer[3] * env_pointer[N_samples+1] + 0.5; ij_index_coef = ij_index_coef + 4; } } </code></pre>

https://stackoverflow.com/q/77129875 2 Vinicius Oliveira Fontes https://stackoverflow.com/users/22588044 2023-09-18T19:20:41Z 2023-09-19T01:06:03Z

<p>I would to ask how I can access a variable on Fortran from MATLAB without making a copy. Preferrably, I would like to access the variable a variable on Fortran (it may be 'read-only') that was declared and initialized on MATLAB and passed as an input to my Fortran-based MEX function as input (pretty much a pass-by-reference), and if possible, do it for a multidimensional array.</p> <p>The solutions that I have found so far require copying the variable, which is quite costly for large inputs.</p> <p>So far, I have tried the approach based on the solution given in the <a href="https://www.mathworks.com/help/matlab/matlab_external/create-fortran-source-mex-file.html" rel="nofollow noreferrer">MATLAB documentation</a>, where they implement a function in the form <code>y = times2(x)</code> which doubles my double scalar <code>x</code>. However, it requires copying the variable <code>x</code>, so there is the original <code>x</code> declared on MATLAB and another one on Fortran.</p> <p>For reference, MATLAB recommends first declaring a pointer and its associated variable</p> <p><code>mwPointer x_ptr real*8 x_input</code></p> <p>And then copying the content of the MATLAB variable <code>x</code> into the Fortran variable <code>x_input</code> via:</p> <p><code>call mxCopyPtrToReal8(x_ptr,x_input,size)</code></p> <p>Where <code>size</code> is 1, as <code>x</code> is a scalar, but in my application I will eventually use a multi-dimensional array.</p> <p>This solution is unappealing when one needs to operate on large arrays, as they would be copied, doubling my memory requirements. I have found <a href="https://www.mathworks.com/help/matlab/apiref/mexgetvariableptr.html" rel="nofollow noreferrer">this function</a> <code>mexGetVariablePtr</code> which returns a read-only MATLAB pointer that theoretically should give me access to the variable <code>x</code> without copying it.</p> <p>My question is: how do I access this variable, as the only function that I found in the API documentation, <code>mxCopyPtrToReal</code>, copies the variable.</p> <p>Additional information:</p> <ul> <li>I only want to read the contents of the variable 'x' and use as a normal Fortran variable without modifying it.</li> <li>The output of 'mexGetVariablePtr' is a 'mwPointer', which seems to not be easy to dereference like regular Fortran pointers.</li> </ul> <p>Thank you for your attention!</p>

https://stackoverflow.com/q/77105025 1 rinkert https://stackoverflow.com/users/7054628 2023-09-14T12:45:54Z 2023-09-15T14:03:55Z

<p>I ran into a weird case where Matlab closes without an error.</p> <p>I have the following mex file that uses OpenMP</p> <pre class="lang-cpp prettyprint-override"><code>// main.cpp #include &quot;mex.h&quot; #include &quot;omp.h&quot; void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { int num_threads; int kx; #pragma omp parallel for for (kx = 0; kx &lt; 100; kx++) { if (kx == 0) num_threads = omp_get_num_threads(); float tmp = (float)kx; } mexPrintf(&quot;Num threads: %i\n&quot;, num_threads); } </code></pre> <p>When I compile this with <code>mex COMPFLAGS=&quot;$COMPFLAGS /openmp&quot; main.cpp -output test_openmp</code> and run it, it works as expected</p> <pre><code>&gt;&gt; test_openmp() Num threads: 8 </code></pre> <p>However, if I immediately add a <code>clear mex</code> after the call, Matlab crashes and immediately closes without error report.</p> <pre><code>&gt;&gt; test_openmp(); clear mex </code></pre> <p>If I add a <code>pause(1)</code>, there is no crash</p> <pre><code>&gt;&gt; test_openmp(); pause(1); clear mex </code></pre> <p>This is a problem in another program where call a mex program can case an error and I need to do:</p> <pre><code>try mexprogram() catch ME clear mexprogram end </code></pre> <p>Does anybody know what is going on? When I run Matlab from the command line, I see the following return code: -1073741819. Attaching the debugger does not give any info.</p> <p>I am using MSVC 2022 and Matlab R2021b.</p> <hr /> <p><strong>Update</strong></p> <p>Linking against the intel openMP implementation <code>libiomp5md.lib</code> fixes the issue! Compiling from MATLAB with</p> <pre><code>mex -v -g COMPFLAGS=&quot;$COMPFLAGS /openmp&quot; LINKFLAGS=&quot;$LINKFLAGS /nodefaultlib:vcomp&quot; '-LC:\Program Files\MATLAB\R2021b\bin\win64\' -llibiomp5md.lib main.cpp -output test_openmp_intel </code></pre>

https://stackoverflow.com/q/76672265 2 zhaokg https://stackoverflow.com/users/17431200 2023-07-12T15:34:19Z 2023-07-12T20:18:25Z

<p>In Matlab, we can use <code>mex</code> to compile a library from multiple files (e.g., 100+ .c files) using something like</p> <pre><code>mex *.c -o mylib.mex </code></pre> <p>But in Octave, it seems like that only one input source file accepted, it always gives an error</p> <pre><code>cc1.exe: fatal error: *.c: Invalid argument </code></pre> <p>Tried googing a lot and find no solution. Any pointer to getting this problem solved will be appreciated.</p>

https://stackoverflow.com/q/76474884 0 drakon101 https://stackoverflow.com/users/17059713 2023-06-14T14:54:17Z 2023-06-14T14:54:17Z

<p>I am using a <a href="https://marketplace.visualstudio.com/items?itemName=gharveymn.MEXFunctionTemplates" rel="nofollow noreferrer">visual studio template</a> to produce Matlab MEX functions. Using this template alone, I can compile functions from Visual Studio using either Visual Studio's compiler or one of the Intel C++/DPC++ compilers. However, when I include mkl libraries, I encounter an error at the include line saying &quot;cannot open source file mkl.h&quot;</p> <p>I have included and linked all the relevant directories and libraries. Here are the command lines from the configuration properties in visual studio. Why could this be happening?</p> <p>C/C++ Properties:</p> <pre><code>/GS /TP /W3 /I&quot;C:\Program Files (x86)\Intel\oneAPI\mkl\latest\include&quot; /I&quot;C:\Program Files\MATLAB\R2023a\extern\include\&quot; /I&quot;C:\Program Files\MATLAB\R2023a\simulink\include\&quot; /Z7 /O2 /Zp8 /D &quot;USE_MEX_CMD&quot; /D &quot;_CRT_SECURE_NO_DEPRECATE&quot; /D &quot;_SCL_SECURE_NO_DEPRECATE&quot; /D &quot;_SECURE_SCL=0&quot; /D &quot;MATLAB_MEX_FILE&quot; /D &quot;__INTEL_LLVM_COMPILER=20230100&quot; /D &quot;_WINDLL&quot; /Zc:forScope /GR /MD /Fa&quot;C:\Users\nikun\source\repos\JCFMKL\JCFMKL\build\x64\Release\&quot; /EHs /nologo /Fo&quot;C:\Users\nikun\source\repos\JCFMKL\JCFMKL\build\x64\Release\&quot; //fprofile-instr-use &quot;C:\Users\nikun\source\repos\JCFMKL\JCFMKL\build\x64\Release\&quot; /Fp&quot;C:\Users\nikun\source\repos\JCFMKL\JCFMKL\build\x64\Release\JCFMKL.pch&quot; </code></pre> <p>Linker Properties</p> <pre><code>/OUT:&quot;C:\Users\nikun\source\repos\JCFMKL\bin\x64\Release\JCFMKL.mexw64&quot; /MANIFEST /NXCOMPAT /PDB:&quot;C:\Users\nikun\source\repos\JCFMKL\bin\x64\Release\JCFMKL.mexw64.pdb&quot; /DYNAMICBASE &quot;mkl_rt.lib&quot; &quot;mkl_intel_c.lib&quot; &quot;mkl_intel_thread.lib&quot; &quot;mkl_intel_lp64.lib&quot; &quot;mkl_core.lib&quot; &quot;libiomp5md.lib&quot; &quot;libmx.lib&quot; &quot;libmex.lib&quot; &quot;libmat.lib&quot; &quot;libMatlabDataArray.lib&quot; &quot;libMatlabEngine.lib&quot; &quot;kernel32.lib&quot; &quot;user32.lib&quot; &quot;gdi32.lib&quot; &quot;winspool.lib&quot; &quot;comdlg32.lib&quot; &quot;advapi32.lib&quot; &quot;shell32.lib&quot; &quot;ole32.lib&quot; &quot;oleaut32.lib&quot; &quot;uuid.lib&quot; &quot;odbc32.lib&quot; &quot;odbccp32.lib&quot; /IMPLIB:&quot;C:\Users\nikun\source\repos\JCFMKL\bin\x64\Release\JCFMKL.lib&quot; /DLL /MACHINE:X64 /INCREMENTAL:NO /MANIFESTUAC:&quot;level='asInvoker' uiAccess='false'&quot; /ManifestFile:&quot;C:\Users\nikun\source\repos\JCFMKL\JCFMKL\build\x64\Release\JCFMKL.mexw64.int ermediate.manifest&quot; /NOLOGO /LIBPATH:&quot;C:\Program Files (x86)\Intel\oneAPI\mkl\latest\lib\intel64&quot; /LIBPATH:&quot;C:\Program Files (x86)\Intel\oneAPI\mkl\latest\lib&quot; /LIBPATH:&quot;C:\Program Files\MATLAB\R2023a\extern\lib\win64\microsoft\&quot; /TLBID:1 </code></pre>

https://stackoverflow.com/q/76393248 2 drakon101 https://stackoverflow.com/users/17059713 2023-06-02T20:18:49Z 2023-06-03T17:55:31Z

<p>This question is a follow up to a previous <a href="https://stackoverflow.com/questions/76353633/efficient-method-to-pass-large-arrays-from-matlab-to-c-eigen-using-map-class/76353749#76353749">question</a> about building an Eigen::Map to an input matlab::data::Array object in the Matlab Data API for C++. Now I want to know how I can implement something similar but to an output array which I will populate with data as I loop through it.</p> <p>In the previous question, I implemented the linked solution from <a href="https://www.mathworks.com/matlabcentral/answers/1573713-c-mex-data-api-how-can-i-get-the-raw-pointer-of-an-array-input-without-copying-it" rel="nofollow noreferrer">Matlab's forums</a>. Here, I tried using the following to get a non-const pointer:</p> <pre class="lang-cpp prettyprint-override"><code>template &lt;typename T&gt; T* getOutDataPtr(matlab::data::Array arr) { matlab::data::TypedArray&lt;T&gt; arr_t = arr; matlab::data::TypedIterator&lt;T&gt; it(arr_t.begin()); return it.operator-&gt;(); } </code></pre> <p>Which is used like so:</p> <pre class="lang-cpp prettyprint-override"><code>// Initialize input parameters long int numEl = inputs[0][0]; long int na = inputs[0][1]; long int nPoints = inputs[0][2]; // Preallocate output and define map to output array outputs[0] = factory.createArray&lt;double&gt; ({static_cast&lt;size_t&gt;(nPoints),static_cast&lt;size_t&gt;(numEl*na)}); auto ptrRecon = getOutDataPtr&lt;double&gt;(outputs[0]); Eigen::Map&lt;Eigen::MatrixXd&gt; Recon(ptrRecon,nPoints,numEl*na); for(int i = 0; i &lt; nPoints; i++) { for(int n = 0; n &lt; na; n++) { for (int j = 0; j &lt; numEl; j++) { Recon(i,n*numEl + j) = 5.0; } } } </code></pre> <p>Even with a simple assignment like that however, I encounter runtime errors after compiling and Matlab crashes. What exactly is going wrong?</p> <p><strong>Additional Trials</strong> I also tried the following for the inner loop which does seem to work:</p> <pre class="lang-cpp prettyprint-override"><code>for(int i = 0; i &lt; nPoints; i++) { for(int n = 0; n &lt; na; n++) { for (int j = 0; j &lt; numEl; j++) { outputs[0][i][n*numEl + j] = 5.0; } } } </code></pre> <p>I suspect this means there's some issue with the way I am getting the pointer to the output matlab data array?</p>

https://stackoverflow.com/q/76316325 0 andro1666 https://stackoverflow.com/users/21947713 2023-05-23T15:19:31Z 2024-03-18T15:54:43Z

<p>I am running a Monte-Carlo experiment and evaluate the absolute loss function described below. Since it is very computationally intensive I would like to optimise my code and improve the speed further. My main code is in MATLAB, but I am evaluating the function in C using the MEX functionality of MATLAB.</p> <p>The mathematical problem is as follows: I have a matrix D with dimensions (M times N). Usually M is around 20,000 and N takes values around {10, 30, 144}.</p> <p><a href="https://i.sstatic.net/0J8wh.png" rel="nofollow noreferrer">Definition of matrix D</a></p> <p>Effectively, I need to obtain L column vector with dimensions (M times 1) defined as</p> <p><a href="https://i.sstatic.net/5N60v.png" rel="nofollow noreferrer">Definition of matrix L</a></p> <p>My C function looks like this:</p> <pre><code>void absolute_loss(double *D, double *L, mwSize cols, mwSize rows) { double aux; int i; int j; int k; for (i = 0; i &lt; rows; i++) { for (j = 0; j &lt; rows; j++){ aux = 0; for (k = 0; k &lt; cols; k++) { aux = aux + fabs(D[j + rows * k] - D[i + rows * k]); } L[i] = L[i] + aux; } } for (i = 0; i &lt; rows; i++) { L[i] /= rows; } } </code></pre>

https://stackoverflow.com/q/75793203 0 drakon101 https://stackoverflow.com/users/17059713 2023-03-20T17:00:20Z 2023-03-20T19:00:20Z

<p>I have the following code snippet:</p> <pre class="lang-cpp prettyprint-override"><code>Eigen::MatrixXcd Recon = Eigen::MatrixXcd::Zero(p.nPoints*p.na, p.numEl); // Do some functions to add data to Recon outputs[0] = factory.createArray&lt;std::complex&lt;double&gt;&gt;({static_cast&lt;size_t&gt;(p.numEl),static_cast&lt;size_t&gt;(p.na),static_cast&lt;size_t&gt;(p.nPoints)}); for(int i = 0; i &lt; p.na; i++) { for (int j = 0; j &lt; p.nPoints; j++) { for (int n = 0; n &lt; p.numEl; n++) { outputs[n][i][j] = Recon( i*p.nPoints + j ,n); } } } </code></pre> <p>How can I adjust my code above so that I am mapping the data in outputs to the data in Recon instead of copying it element by element?</p>

https://stackoverflow.com/q/75777807 1 Nukesub https://stackoverflow.com/users/5154085 2023-03-18T18:15:51Z 2023-03-19T07:19:21Z

<p>For some context, I am familiar with C programming (not an expert but a daily user) and also MATLAB. I have used C++ minimally and not recently.</p> <p>I had recently started looking into the MEX functionality of MATLAB as I hadn't ever looked at it before. While I was working through a <a href="https://www.mathworks.com/help/matlab/matlab_external/standalone-example.html" rel="nofollow noreferrer">basic example</a> I noticed something I had not seen before in the entry point function <code>mexFunction</code>, copy/pasted from the docs below.</p> <p>The output and input arguments are returned as pointers of type <code>mxArray</code> but with an empty bracket (e.g., <code>mxArray *plhs[]</code>). Coming from C I would have just expected <code>mxArray *plhs</code>. I can also understand you might want to return multiple types like scalars, matrices, etc. so maybe you want a pointer to pointers so wouldn't you do <code>mxArray **plhs</code>?</p> <p>From this post <a href="https://stackoverflow.com/questions/10350461/c-asterisk-and-bracket-operators-used-together">here</a> I understand the syntax <code>mxArray *plhs[]</code> means an array of pointers. Is that different than doing <code>mxArray **plhs</code>?</p> <pre class="lang-matlab prettyprint-override"><code>/* The gateway function */ void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { /* variable declarations here */ /* code here */ } </code></pre>

https://stackoverflow.com/q/75693305 0 G. B. https://stackoverflow.com/users/9963147 2023-03-10T07:31:18Z 2023-03-10T09:05:24Z

<p>We are building a matlab executable (.mexw64) from C++ using MSVC compiler (we tested compiler versions 15 and 16) and a cmake toolchain (not with the mex command). The mexw64 contains functionality, that we can call from Matlab, and these functions need to log. My problem is, that when building in debug, the std::cout is not routed to the matlab console (we display log messages this way).</p> <p>Details:</p> <ul> <li>if built in release mode, the log messages display correctly in R2017 and R2021</li> <li>if build in debug mode, the log messages are only displayed in R2017 and not displayed R2021</li> <li>MSVC compiler version 15 and 16 work the same way</li> </ul> <p>What was changed in this regard between R2017 and R2021? How can I solve this issue? Changing the logging from std::cout to something else (e.g. mex logging function or std::cerr) is unfortunately not an option.</p> <p>Note: I noticed in the build log, that the option <code>-DMATLAB_DEFAULT_RELEASE=R2017b</code>is set. What does this do? Is it ok, to have this value, even if my Matlab version is R2021b? The option is set in a cmake function:</p> <pre><code> # For 9.4 (R2018a) and newer, add API macro. # Add it for unknown versions too, just in case. if(NOT Matlab_VERSION_STRING VERSION_LESS &quot;9.4&quot; OR Matlab_VERSION_STRING STREQUAL &quot;unknown&quot;) if(${${prefix}_R2018a}) set(MEX_API_MACRO &quot;MATLAB_DEFAULT_RELEASE=R2018a&quot;) else() set(MEX_API_MACRO &quot;MATLAB_DEFAULT_RELEASE=R2017b&quot;) endif() endif() </code></pre> <p>Edit: I did not mention that we correctly <a href="https://stackoverflow.com/questions/243696/correctly-over-loading-a-stringbuf-to-replace-cout-in-a-matlab-mex-file">overload the std::streambuf</a> as @Cris Luengo pointed out. I do not copy the code here, because it's long and proprietary, but it works with R2017b and also with R2021b in release mode.</p>

https://stackoverflow.com/q/74913843 0 SheppLogan https://stackoverflow.com/users/7514181 2022-12-25T14:22:26Z 2022-12-26T04:30:46Z

<p>I have a Windows Matlab R2020a and was able to generate C source code for my Matlab function. However, I need to run a MEX file (for speed) on a Linux platform which has a Matlab installed (without Coder) and no way to install it (even trial version).</p> <p>Please, provide the steps to compile this Matlab-generated C code on the linux (Ubuntu) system. I guess it can be done with gcc. I am unable to find a clear description on the web, as this process seems extremely complicated due to the obscure command line args needed to give to gcc (and i only rarely compile C code).</p>

https://stackoverflow.com/q/74722979 1 mfeuling https://stackoverflow.com/users/5763459 2022-12-07T21:25:43Z 2022-12-08T19:24:48Z

<p>I am putting together a minimal example leveraging <a href="https://en.cppreference.com/w/cpp/experimental/parallelism" rel="nofollow noreferrer">parallelism features in C++17/20</a> within Matlab MEX functions. I am able to compile and run the mex function from Matlab, but when I set the <a href="https://en.cppreference.com/w/cpp/experimental/execution_policy_tag_t" rel="nofollow noreferrer">execution policy</a> of my C++ STL function to &quot;par&quot; instead of &quot;seq&quot;, Matlab gives a runtime linkage complaint. Code and error message follows:</p> <p>test.m (Matlab top-level script):</p> <pre><code>vec_in = zeros(5); coeff = 0.05; vec_out = test_mex_gateway(vec_in, coeff); </code></pre> <p>test_mex_gateway.cpp (C++ interface to Matlab):</p> <pre><code>#include &quot;mex.h&quot; extern void test_execute(float *array_in, float *array_out, const size_t vec_size, const float coeff); void mexFunction( int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[] ) { // Check for proper number of input and output arguments if( nrhs != 2 ) { mexErrMsgTxt( &quot;3 input arguments required: input_data, coeff&quot; ); } if( nlhs &gt; 2 ) { mexErrMsgTxt( &quot;Too many output arguments.&quot; ); } const mwSize *matlab_data_dims_in; mwSize matlab_data_dims_out[1]; // Input Parameters float *input_data = (float *) mxGetData(prhs[0]); float coeff = mxGetScalar(prhs[1]); // Get dimensions matlab_data_dims_in = mxGetDimensions(prhs[0]); const int vec_len = matlab_data_dims_in[1]; // Set output data dimension matlab_data_dims_out[0] = vec_len; // Output data plhs[0] = mxCreateNumericArray(1, matlab_data_dims_out, mxSINGLE_CLASS, mxREAL); float *output_data = (float *) mxGetData(plhs[0]); test_execute(input_data, output_data, vec_len, coeff); } </code></pre> <p>test_execute.cpp (This is where the actual C++ STL call is made):</p> <pre><code>#include &lt;execution&gt; // std::execution::* #include &lt;numeric&gt; // std::exclusive_scan() void test_execute(float *array_in, float *array_out, const size_t vec_size, const float coeff) { std::exclusive_scan ( std::execution::par, // std::execution::seq works here for Mex call, par does not array_in, array_in + vec_size, array_out, 0.0f, [coeff](float a, float b) { float ret = a + b + coeff; return ret; } ); } </code></pre> <p>I also have a stand-alone main function to replace the Mex wrapper to do a pure C++ test, test_standalone.cpp:</p> <pre><code>#include &lt;vector&gt; #include &lt;iostream&gt; size_t VEC_NUM_ELEM = 10; extern void test_execute(float *array_in, float *array_out, const size_t vec_size, const float coeff); int main(int argc, char **argv) { if (argc != 2) { std::cout &lt;&lt; &quot;Try: &quot; &lt;&lt; argv[0] &lt;&lt; &quot;&lt;coeff&gt;&quot; &lt;&lt; std::endl; return -1; } const float coeff = std::stof(argv[1]); std::cout &lt;&lt; &quot;Coeff: &quot; &lt;&lt; coeff &lt;&lt; std::endl; float __attribute__ ((aligned (64))) *vec1_array = (float *)malloc(VEC_NUM_ELEM * sizeof(float)); float __attribute__ ((aligned (64))) *vec2_array = (float *)malloc(VEC_NUM_ELEM * sizeof(float)); for (unsigned i = 0; i &lt; VEC_NUM_ELEM; i++) { vec1_array[i] = static_cast&lt;float&gt;(i); } test_execute(vec1_array, vec2_array, VEC_NUM_ELEM, coeff); return 0; } </code></pre> <p>Here is how I am building and linking, build.sh:</p> <pre><code>#!/bin/bash rm *.o rm *.exe rm *.mexa64 cstd=c++17 gpp910=/home/m/compilers/bin/g++ tbblib=/home/m/reqs/tbb/lib/intel64/gcc4.8 echo &quot;Building test_execute.cpp...&quot; $gpp910 -std=$cstd -I/home/m/reqs/tbb/include -L$tbblib -ltbb -Wl,rpath=$tbblib -c test_execute.cpp -fPIC echo &quot;Building test_standalone.cpp...&quot; $gpp910 -std=$cstd -L$tbblib test_execute.o test_standalone.cpp -o test_standalone.exe -ltbb echo &quot;Building test_mex_gateway.cpp...&quot; mex test_execute.o test_mex_gateway.cpp -L$tbblib -ltbb </code></pre> <p>The parallel STL calls has a requirement to link against the Intel TBB (Threading Building Blocks), so before I run Matlab to call test.m OR before I run my test_standalone.exe, I run:</p> <p>export LD_LIBRARY_PATH=/home/m/reqs/tbb/lib/intel64/gcc4.8:$LD_LIBRARY_PATH</p> <p>I also make sure to make the the C++ library associated with the version of GCC we built with available at runtime:</p> <p>export LD_LIBRARY_PATH=/home/m/compilers/lib64:$LD_LIBRARY_PATH</p> <p>When I run test_standalone.exe, everything behaves normally whether I have the execution policy set to &quot;par&quot; or &quot;seq&quot; on std::exclusive_scan. When run test.m, if &quot;seq&quot; was compiled, I can run with no errors. If &quot;par&quot; was compiled, Matlab complains at runtime about a linkage issue:</p> <p><em>Invalid MEX-file 'test_mex_gateway.mexa64': test_mex_gateway.mexa64: undefined symbol: _ZN3tbb10interface78internal20isolate_within_arenaERNS1_13delegate_baseEl</em></p> <p>I suspect this was a function that was supposed to be linked from TBB, which I confirmed:</p> <p>$ nm /home/m/reqs/tbb/lib/intel64/gcc4.8/libtbb.so.2 | grep baseEl</p> <p>0000000000028a30 T _ZN3tbb10interface78internal20isolate_within_arenaERNS1_13delegate_baseEl</p> <p>000000000005ed70 r _ZN3tbb10interface78internal20isolate_within_arenaERNS1_13delegate_baseEl$$LSDA</p> <p>I confirmed Matlab's LD_LIBRARY_PATH has the path I supplied in the above &quot;export ..&quot; to this library.</p> <p>I tried making sure my libraries came before the many Matlab-centric paths Matlab adds to LD_LIBRARY_PATH after it launches from the terminal.</p> <p>I tried baking the path to the linked libraries via a -Wl,rpath=&lt;path_to_tbb.so&gt; passage to the linker.</p> <p>After almost two days, I can't figure out why Matlab is having this very specific runtime issue, especially when the pure C++ version is not. Any help would be appreciated.</p> <p><strong>RHEL 7.9</strong></p> <p><strong>Matlab R2020a</strong></p> <p><strong>GCC 9.1.0</strong></p> <p><strong>TBB (Intel Thread Building Blocks) 2020.3</strong></p>

https://stackoverflow.com/q/74041072 0 Gagandeep Singh https://stackoverflow.com/users/17960694 2022-10-12T11:29:28Z 2022-10-12T15:34:38Z

<p><strong>mex_example0.cpp</strong></p> <pre><code>#include &quot;mex.hpp&quot; #include &quot;mexAdapter.hpp&quot; #include &lt;iostream&gt; using namespace matlab::data; using matlab::mex::ArgumentList; class MexFunction : public matlab::mex::Function { public: void operator()(ArgumentList outputs, ArgumentList inputs) { ArrayFactory factory; TypedArray&lt;double&gt; x = std:: move(inputs[0]); for (auto&amp; i: x){ i += 5; } outputs[0] = x; std::cout&lt;&lt;&quot;Hello MEX&quot;&lt;&lt;std::endl; } }; </code></pre> <p><strong>MATLAB code</strong></p> <pre><code>mex mex_example0.cpp </code></pre> <p><strong>Error</strong></p> <pre><code>Building with 'Microsoft Visual C++ 2012'. Error using mex mex_example0.cpp C:\Users\iisermbionap\Downloads\mex_example0.cpp(1) : fatal error C1083: Cannot open include file: 'mex.hpp': No such file or directory </code></pre> <p>I am using MATLAB 2015a and MVC 2012 professional.<br /> Can somebody please help me how to rectify this error.</p>

https://stackoverflow.com/q/73947526 0 Hadi_IronSide https://stackoverflow.com/users/20158033 2022-10-04T11:57:34Z 2022-10-06T15:38:24Z

<p>I have a question about display C++ output in MATLAB. I wrote a mex function to include my C++ algorithm in MATLAB but I have a problem with displaying my output in MATLAB. The c++ mex function calls other classes where my output is generated. To keep the matter simple and clear, here is an example.</p> <p>Mexfunction.cpp:</p> <pre class="lang-cpp prettyprint-override"><code>#include &quot;mex.hpp&quot; #include &quot;mexAdapter.hpp&quot; #include &quot;Helper.h&quot; class MexFunction : public matlab::mex::Function { public: void operator()(matlab::mex::ArgumentList outputs, matlab::mex::ArgumentList inputs) { // Function implementation ... Helper(); } }; </code></pre> <p>(Helper Class, that is called in Mexfunction.cpp)</p> <pre class="lang-cpp prettyprint-override"><code>#ifndef HELPER_H #define HELPER_H #include &lt;iostream&gt; class Helper { public: Helper() { std::cout &lt;&lt; &quot;Hello World&quot; &lt;&lt; std::endl; } }; #endif </code></pre> <p>How I can display the string <code>&quot;Hello world&quot;</code> in the MATLAB command window?</p>

https://stackoverflow.com/q/73928582 1 Paul Aner https://stackoverflow.com/users/16975921 2022-10-02T19:06:27Z 2022-10-02T21:34:53Z

<p>I want to do some simulations that need C++-Code in Matlab. This works like a charm, but as soon as I put in some additional <code>#include</code>'s in the C++/MEX-Code I get a ton of error messages like this:</p> <pre><code> C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(19) : error C2054: expected '(' to follow 'using' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(19) : error C2061: syntax error : identifier 'using' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(21) : error C2061: syntax error : identifier 'using' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(21) : error C2061: syntax error : identifier 'abs' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(21) : error C2059: syntax error : ';' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(21) : error C2061: syntax error : identifier 'atexit' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(22) : error C2061: syntax error : identifier 'atof' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(22) : error C2059: syntax error : ';' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(22) : error C2061: syntax error : identifier 'atoi' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(22) : error C2061: syntax error : identifier 'atol' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(23) : error C2061: syntax error : identifier 'bsearch' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(23) : error C2059: syntax error : ';' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(23) : error C2061: syntax error : identifier 'calloc' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(23) : error C2061: syntax error : identifier 'div' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(24) : error C2061: syntax error : identifier 'exit' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(24) : error C2059: syntax error : ';' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(24) : error C2061: syntax error : identifier 'free' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(25) : error C2061: syntax error : identifier 'labs' C:\Program Files (x86)\Microsoft Visual Studio 11.0\VC\INCLUDE\cstdlib(25) : error C2059: syntax error : ';' ... </code></pre> <p>I tried different versions of Matlab, Visual Studio and different C++-Compilers. I have code that works perfectly fine and also tried this with example-code from Mathworks that also works fine.</p> <p>For testing-purposes I simply put <code>#include &lt;random&gt;</code> or <code>include &lt;algorithm&gt;</code> on top of working examples and then those error-messages occur. No &quot;file not found&quot;-errors.</p> <p>Can anybody give me a hint what I'm doing wrong here? I didn't change anything in <strong>mexopts.bat</strong> and I believe the right include-directories are in there. Programming in Visual Studio itself works without any problems with those <code>#include</code>'s...</p>

https://stackoverflow.com/q/73885145 0 Naveen https://stackoverflow.com/users/3565667 2022-09-28T16:49:20Z 2022-09-28T17:08:59Z

<p>I have an octave DLD function say main.cpp which calls several functions declared and defined(in .h and .cpp files resp) in the directory mathlib/ which is in the same location as main.cpp.</p> <p>I'm referring to <a href="https://docs.octave.org/v4.2.0/Getting-Started-with-Oct_002dFiles.html" rel="nofollow noreferrer">this</a> octave documentation page for compile options for including and linking libraries</p> <p>I'm getting the following errors:</p> <pre><code>&gt;&gt; mkoctfile -Imathlib \ main.cpp g++: warning: main.cpp -o &lt;PATH&gt;\AppData\Local\Temp/oct-uoZEQc.o: linker input file unused because linking not done g++: error: main.cpp -o &lt;PATH&gt;\AppData\Local\Temp/oct-uoZEQc.o: linker input file not found: Invalid argument error: mkoctfile: building exited with failure status </code></pre> <p>What are the other compile options and steps needed to get the linking to work and compile main.cpp successfully?</p>

https://stackoverflow.com/q/73827088 1 Naveen https://stackoverflow.com/users/3565667 2022-09-23T11:35:47Z 2022-09-23T14:32:58Z

<p>Consider the following simplified example:</p> <p>File 1 (The Octave .m file)</p> <pre><code>% test_multiplyCxArrays.m clc; clear all; close all; pkg load signal mkoctfile multiplyCxArrays.cpp disp('Compiled Octave-C++ function multiplyCxArrays.cpp'); cx_one = 1+1i; a = [cx_one; cx_one*2; cx_one*3; cx_one*4; cx_one*5]; b = [cx_one; cx_one/2; cx_one/3; cx_one/4; cx_one/5]; disp('First input vector: '); disp(a); disp('First input vector inside DLD: '); cx_out = multiplyCxArrays(a,b); disp('Output vector returned DLD: '); disp(cx_out); </code></pre> <p>File 2 (The C++ DLD file)</p> <pre><code>// multiplyCxArrays.cpp #include &lt;iostream&gt; #include &lt;octave/oct.h&gt; #include &lt;iostream&gt; #include &lt;complex&gt; #include &lt;vector&gt; #include &lt;fstream&gt; #include &lt;cmath&gt; #include &lt;math.h&gt; #include &quot;CNDArray.h&quot; using namespace std; DEFUN_DLD (multiplyCxArrays, args, nargout, &quot;multiplyCxArrays is a c++ compiled example to illustrate dynamically loaded functions in octave &quot;) { ComplexNDArray cxIn1 = args(0).array_value(); ComplexNDArray cxIn2 = args(1).array_value(); const int cxIn1_length = cxIn1.rows(); const int cxIn2_length = cxIn2.rows(); std::vector&lt;std::complex&lt;double&gt;&gt; cxVec; std::complex&lt;double&gt; cx_elem1, cx_elem2; for(octave_idx_type i = 0; i &lt; cxIn1_length; ++i) { cout &lt;&lt; &quot;i = &quot; &lt;&lt; i &lt;&lt; &quot;, cxIn1(i) = &quot; &lt;&lt; cxIn1(i) &lt;&lt; endl; cx_elem1 = cxIn1(i); cx_elem2 = cxIn2(i); cxVec.push_back(cx_elem1*cx_elem2); } const int cxOut_length = (int)cxVec.size(); dim_vector dv (cxOut_length, 1); ComplexNDArray cxOut(dv); dv = cxOut.dims(); for (octave_idx_type i = 0; i &lt; cxOut_length; ++i) { cxOut(i) = cxVec[(int)i]; } octave_value_list retval(nargout); retval(0) = octave_value(cxOut); return retval; } </code></pre> <p>OUTPUT printout in Octave command window:</p> <pre><code>First input vector: 1 + 1i 2 + 2i 3 + 3i 4 + 4i 5 + 5i First input vector inside DLD: i = 0, cxIn1(i) = (1,0) i = 1, cxIn1(i) = (2,0) i = 2, cxIn1(i) = (3,0) i = 3, cxIn1(i) = (4,0) i = 4, cxIn1(i) = (5,0) Output vector returned DLD: 1 1 1 1 1 </code></pre> <p>For some reason the first argument of the DLD <strong>args(0).array_value()</strong> is only returning the real part of the input complex array and zeroing out the imaginary part as you can see from the output printout for <strong>cxIn1(i)</strong>. Is there any way to retrieve the full complex array from the input argument list?</p>