Compare commits

3 Commits

Author SHA1 Message Date
b6f452e74c feat(libconfig): new version of libconfig 2025-12-06 11:41:57 -05:00
7242c765f3 build(wasm): major progress on gridfire compiling to wasm 2025-12-03 11:38:08 -05:00
d852ee43fe perf(precomputation): cleaned up allocations
recovered about 5% execution time
2025-12-02 13:09:19 -05:00
35 changed files with 437 additions and 280 deletions

3
.gitignore vendored
View File

@@ -82,6 +82,7 @@ subprojects/cvode-*/
subprojects/kinsol-*/ subprojects/kinsol-*/
subprojects/CLI11-*/ subprojects/CLI11-*/
subprojects/openssl-*/ subprojects/openssl-*/
subprojects/tomlplusplus-*/
*.fbundle *.fbundle
*.wraplock *.wraplock
@@ -98,6 +99,8 @@ liblogging.wrap
libplugin.wrap libplugin.wrap
minizip-ng.wrap minizip-ng.wrap
openssl.wrap openssl.wrap
glaze.wrap
tomlplusplus.wrap
.vscode/ .vscode/

View File

View File

@@ -0,0 +1,33 @@
cppc = meson.get_compiler('cpp')
if cppc.get_id() == 'clang'
message('disabling bitwise-instead-of-logical warnings for clang')
add_project_arguments('-Wno-bitwise-instead-of-logical', language: 'cpp')
endif
if cppc.get_id() == 'gcc'
message('disabling psabi warnings for gcc')
add_project_arguments('-Wno-psabi', language: 'cpp')
if (cppc.version().version_compare('<14.0'))
error('g++ version must be at least 14.0, found ' + cppc.version())
endif
endif
if not cppc.has_header('print')
error('C++ standard library header <print> not found. Please ensure your compiler and standard library supports C++23. We have already validated your compiler version so this is likely an issue with your standard library installation.')
endif
if not cppc.has_header('format')
error('C++ standard library header <format> not found. Please ensure your compiler and standard library supports C++23. We have already validated your compiler version so this is likely an issue with your standard library installation.')
endif
# For Eigen
add_project_arguments('-Wno-deprecated-declarations', language: 'cpp')
if get_option('build_python')
message('enabling hidden visibility for C++ symbols when building Python extension. This reduces the size of the resulting shared library.')
add_project_arguments('-fvisibility=hidden', language: 'cpp')
else
message('enabling default visibility for C++ symbols')
add_project_arguments('-fvisibility=default', language: 'cpp')
endif

View File

@@ -0,0 +1,15 @@
if get_option('build_fortran')
add_languages('fortran', native: true)
message('Found FORTRAN compiler: ' + meson.get_compiler('fortran').get_id())
message('Fortran standard set to: ' + get_option('fortran_std'))
message('Building fortran module (gridfire_mod.mod)')
fc = meson.get_compiler('fortran')
if not get_option('unsafe_fortran')
if fc.get_id() != 'gcc'
error('The only supported fortran compiler for GridFire is gfortran (version >= 14.0), found ' + fc + '. GridFire has not been tested with any other compilers. You can disable this check with the -Dunsafe-fortran=true flag to try other compilers')
endif
endif
if (fc.version().version_compare('<14.0'))
error('gfortran version must be at least 14.0, found ' + fc.version())
endif
endif

15
build-check/meson.build Normal file
View File

@@ -0,0 +1,15 @@
message('Found CXX compiler: ' + meson.get_compiler('cpp').get_id())
message('C++ standard set to: ' + get_option('cpp_std'))
cc = meson.get_compiler('c')
ignore_unused_args = '-Wno-unused-command-line-argument'
add_global_arguments(ignore_unused_args, language: 'cpp')
add_global_arguments(ignore_unused_args, language: 'c')
subdir('CPPC')
subdir('FC')

View File

@@ -1,17 +1,21 @@
# bring in all of the fourdst utility repositories # bring in all of the fourdst utility repositories
fourdst_build_lib_all = true fourdst_build_lib_all = true
if get_option('unity-safe') if not get_option('plugin_support')
fourdst_build_lib_all=false fourdst_build_lib_all=false
message('Disabling fourdst plugin support as per user request.')
endif endif
fourdst_sp = subproject('fourdst', fourdst_sp = subproject('fourdst',
default_options: default_options:
['build-tests=' + get_option('build-tests').to_string(), ['build_tests=' + get_option('build_tests').to_string(),
'build-python=' + get_option('build-python').to_string(), 'build_python=' + get_option('build_python').to_string(),
'build-lib-all=' + fourdst_build_lib_all.to_string(), 'build_lib_all=' + fourdst_build_lib_all.to_string(),
'pkg-config=' + get_option('pkg-config').to_string(), 'build_lib_comp=true',
'build-lib-log=true' 'build_lib_config=true',
'build_lib_log=true',
'build_lib_const=true',
'pkg_config=' + get_option('pkg_config').to_string(),
] ]
) )
@@ -19,15 +23,16 @@ composition_dep = fourdst_sp.get_variable('composition_dep')
log_dep = fourdst_sp.get_variable('log_dep') log_dep = fourdst_sp.get_variable('log_dep')
const_dep = fourdst_sp.get_variable('const_dep') const_dep = fourdst_sp.get_variable('const_dep')
config_dep = fourdst_sp.get_variable('config_dep') config_dep = fourdst_sp.get_variable('config_dep')
if not get_option('unity-safe') if get_option('plugin_support')
warning('Including plugin library from fourdst. Note this will bring in minizip-ng and openssl, which can cause build issues with cross compilation due to their complexity.')
plugin_dep = fourdst_sp.get_variable('plugin_dep') plugin_dep = fourdst_sp.get_variable('plugin_dep')
endif endif
libcomposition = fourdst_sp.get_variable('libcomposition') libcomposition = fourdst_sp.get_variable('libcomposition')
libconst = fourdst_sp.get_variable('libconst') libconst = fourdst_sp.get_variable('libconst')
libconfig = fourdst_sp.get_variable('libconfig')
liblogging = fourdst_sp.get_variable('liblogging') liblogging = fourdst_sp.get_variable('liblogging')
if not get_option('unity-safe') if get_option('plugin_support')
warning('Including plugin library from fourdst. Note this will bring in minizip-ng and openssl, which can cause build issues with cross compilation due to their complexity.')
libplugin = fourdst_sp.get_variable('libplugin') libplugin = fourdst_sp.get_variable('libplugin')
endif endif

View File

@@ -1,7 +1,9 @@
cmake = import('cmake') cmake = import('cmake')
if get_option('build_python')
subdir('python') subdir('python')
subdir('pybind')
endif
subdir('fourdst') subdir('fourdst')
subdir('sundials') subdir('sundials')
@@ -11,6 +13,5 @@ subdir('eigen')
subdir('json') subdir('json')
subdir('pybind')
subdir('CLI11') subdir('CLI11')

View File

@@ -6,7 +6,7 @@ cvode_cmake_options.add_cmake_defines({
'CMAKE_C_FLAGS' : '-Wno-deprecated-declarations', 'CMAKE_C_FLAGS' : '-Wno-deprecated-declarations',
'BUILD_SHARED_LIBS' : 'OFF', 'BUILD_SHARED_LIBS' : 'OFF',
'BUILD_STATIC_LIBS' : 'ON', 'BUILD_STATIC_LIBS' : 'ON',
'EXAMPLES_ENABLE_C': 'OFF', 'EXAMPLES_ENABLE_C' : 'OFF',
'CMAKE_POSITION_INDEPENDENT_CODE': true 'CMAKE_POSITION_INDEPENDENT_CODE': true
}) })
@@ -16,6 +16,15 @@ cvode_cmake_options.add_cmake_defines({
'CMAKE_INSTALL_INCLUDEDIR': get_option('includedir') 'CMAKE_INSTALL_INCLUDEDIR': get_option('includedir')
}) })
if meson.is_cross_build() and host_machine.system() == 'emscripten'
cvode_cmake_options.add_cmake_defines({
'CMAKE_C_FLAGS': '-s MEMORY64=1 -s ALLOW_MEMORY_GROWTH=1',
'CMAKE_CXX_FLAGS': '-s MEMORY64=1 -s ALLOW_MEMORY_GROWTH=1',
'CMAKE_SHARED_LINKER_FLAGS': '-s MEMORY64=1 -s ALLOW_MEMORY_GROWTH=1',
'CMAKE_EXE_LINKER_FLAGS': '-s MEMORY64=1 -s ALLOW_MEMORY_GROWTH=1'
})
endif
cvode_sp = cmake.subproject( cvode_sp = cmake.subproject(
'cvode', 'cvode',
options: cvode_cmake_options, options: cvode_cmake_options,

View File

@@ -0,0 +1,32 @@
llevel = get_option('log_level')
logbase='QUILL_COMPILE_ACTIVE_LOG_LEVEL_'
if (llevel == 'traceL3')
message('Setting log level to TRACE_L3')
log_argument = logbase + 'TRACE_L3'
elif (llevel == 'traceL2')
message('Setting log level to TRACE_L2')
log_argument = logbase + 'TRACE_L2'
elif (llevel == 'traceL1')
message('Setting log level to TRACE_L1')
log_argument = logbase + 'TRACE_L1'
elif (llevel == 'debug')
message('Setting log level to DEBUG')
log_argument = logbase + 'DEBUG'
elif (llevel == 'info')
message('Setting log level to INFO')
log_argument = logbase + 'INFO'
elif (llevel == 'warning')
message('Setting log level to WARNING')
log_argument = logbase + 'WARNING'
elif (llevel == 'error')
message('Setting log level to ERROR')
log_argument = logbase + 'ERROR'
elif (llevel == 'critical')
message('Setting log level to CRITICAL')
log_argument = logbase + 'CRITICAL'
endif
log_argument = '-DQUILL_COMPILE_ACTIVE_LOG_LEVEL=' + log_argument
add_project_arguments(log_argument, language: 'cpp')

View File

@@ -0,0 +1,18 @@
if get_option('pkg_config')
message('Generating pkg-config file for GridFire...')
pkg = import('pkgconfig')
pkg.generate(
name: 'gridfire',
description: 'GridFire nuclear reaction network solver',
version: meson.project_version(),
libraries: [
libgridfire,
libcomposition,
libconst,
liblogging
],
subdirs: ['gridfire'],
filebase: 'gridfire',
install_dir: join_paths(get_option('libdir'), 'pkgconfig')
)
endif

View File

@@ -1,90 +1,95 @@
if get_option('build_python')
message('Building Python bindings...')
gridfire_py_deps = [ gridfire_py_deps = [
py_dep, py_dep,
pybind11_dep, pybind11_dep,
const_dep, const_dep,
config_dep, config_dep,
composition_dep, composition_dep,
gridfire_dep gridfire_dep
] ]
py_sources = [ py_sources = [
meson.project_source_root() + '/src/python/bindings.cpp', meson.project_source_root() + '/src/python/bindings.cpp',
meson.project_source_root() + '/src/python/types/bindings.cpp', meson.project_source_root() + '/src/python/types/bindings.cpp',
meson.project_source_root() + '/src/python/partition/bindings.cpp', meson.project_source_root() + '/src/python/partition/bindings.cpp',
meson.project_source_root() + '/src/python/partition/trampoline/py_partition.cpp', meson.project_source_root() + '/src/python/partition/trampoline/py_partition.cpp',
meson.project_source_root() + '/src/python/reaction/bindings.cpp', meson.project_source_root() + '/src/python/reaction/bindings.cpp',
meson.project_source_root() + '/src/python/screening/bindings.cpp', meson.project_source_root() + '/src/python/screening/bindings.cpp',
meson.project_source_root() + '/src/python/screening/trampoline/py_screening.cpp', meson.project_source_root() + '/src/python/screening/trampoline/py_screening.cpp',
meson.project_source_root() + '/src/python/io/bindings.cpp', meson.project_source_root() + '/src/python/io/bindings.cpp',
meson.project_source_root() + '/src/python/io/trampoline/py_io.cpp', meson.project_source_root() + '/src/python/io/trampoline/py_io.cpp',
meson.project_source_root() + '/src/python/exceptions/bindings.cpp', meson.project_source_root() + '/src/python/exceptions/bindings.cpp',
meson.project_source_root() + '/src/python/engine/bindings.cpp', meson.project_source_root() + '/src/python/engine/bindings.cpp',
meson.project_source_root() + '/src/python/engine/trampoline/py_engine.cpp', meson.project_source_root() + '/src/python/engine/trampoline/py_engine.cpp',
meson.project_source_root() + '/src/python/solver/bindings.cpp', meson.project_source_root() + '/src/python/solver/bindings.cpp',
meson.project_source_root() + '/src/python/solver/trampoline/py_solver.cpp', meson.project_source_root() + '/src/python/solver/trampoline/py_solver.cpp',
meson.project_source_root() + '/src/python/policy/bindings.cpp', meson.project_source_root() + '/src/python/policy/bindings.cpp',
meson.project_source_root() + '/src/python/policy/trampoline/py_policy.cpp', meson.project_source_root() + '/src/python/policy/trampoline/py_policy.cpp',
meson.project_source_root() + '/src/python/utils/bindings.cpp', meson.project_source_root() + '/src/python/utils/bindings.cpp',
] ]
if meson.is_cross_build() and host_machine.system() == 'darwin' if meson.is_cross_build() and host_machine.system() == 'darwin'
py_mod = shared_module( py_mod = shared_module(
'_gridfire', '_gridfire',
sources: py_sources,
dependencies: gridfire_py_deps,
name_prefix: '',
name_suffix: 'so',
install: true,
install_dir: py_installation.get_install_dir() + '/gridfire'
)
else
py_mod = py_installation.extension_module(
'_gridfire', # Name of the generated .so/.pyd file (without extension)
sources: py_sources, sources: py_sources,
dependencies : gridfire_py_deps, dependencies: gridfire_py_deps,
install : true, name_prefix: '',
name_suffix: 'so',
install: true,
install_dir: py_installation.get_install_dir() + '/gridfire'
)
else
py_mod = py_installation.extension_module(
'_gridfire', # Name of the generated .so/.pyd file (without extension)
sources: py_sources,
dependencies : gridfire_py_deps,
install : true,
subdir: 'gridfire',
)
endif
py_installation.install_sources(
files(
meson.project_source_root() + '/src/python/gridfire/__init__.py',
meson.project_source_root() + '/stubs/gridfire/_gridfire/__init__.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/exceptions.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/partition.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/reaction.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/screening.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/io.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/solver.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/policy.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/type.pyi'
),
subdir: 'gridfire', subdir: 'gridfire',
) )
py_installation.install_sources(
files(
meson.project_source_root() + '/stubs/gridfire/_gridfire/engine/__init__.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/engine/diagnostics.pyi',
),
subdir: 'gridfire/engine',
)
py_installation.install_sources(
files(
meson.project_source_root() + '/stubs/gridfire/_gridfire/utils/__init__.pyi',
),
subdir: 'gridfire/utils',
)
py_installation.install_sources(
files(
meson.project_source_root() + '/stubs/gridfire/_gridfire/utils/hashing/__init__.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/utils/hashing/reaction.pyi',
),
subdir: 'gridfire/utils/hashing',
)
else
message('Python bindings disabled')
endif endif
py_installation.install_sources(
files(
meson.project_source_root() + '/src/python/gridfire/__init__.py',
meson.project_source_root() + '/stubs/gridfire/_gridfire/__init__.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/exceptions.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/partition.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/reaction.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/screening.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/io.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/solver.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/policy.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/type.pyi'
),
subdir: 'gridfire',
)
py_installation.install_sources(
files(
meson.project_source_root() + '/stubs/gridfire/_gridfire/engine/__init__.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/engine/diagnostics.pyi',
),
subdir: 'gridfire/engine',
)
py_installation.install_sources(
files(
meson.project_source_root() + '/stubs/gridfire/_gridfire/utils/__init__.pyi',
),
subdir: 'gridfire/utils',
)
py_installation.install_sources(
files(
meson.project_source_root() + '/stubs/gridfire/_gridfire/utils/hashing/__init__.pyi',
meson.project_source_root() + '/stubs/gridfire/_gridfire/utils/hashing/reaction.pyi',
),
subdir: 'gridfire/utils/hashing',
)

23
cross/wasm.ini Normal file
View File

@@ -0,0 +1,23 @@
[binaries]
c = 'emcc'
cpp = 'em++'
ar = 'emar'
strip = 'emstrip'
exec_wrapper = 'node'
[built-in options]
c_args = ['-Dpkg_config=false', '-Dbuild_tests=false', '-Dbuild_examples=true', '-Dbuild_fortran=falase', '-Dplugin_support=false', '-s', 'MEMORY64=1', '-pthread', '-DQUILL_NO_THREAD_NAME_SUPPORT', '-DQUILL_IMMEDIATE_FLUSH']
cpp_args = ['-Dpkg_config=false', '-Dbuild_tests=false', '-Dbuild_examples=true', '-Dbuild_fortran=falase', '-Dplugin_support=false', '-s', 'MEMORY64=1', '-pthread', '-DQUILL_NO_THREAD_NAME_SUPPORT', '-DQUILL_IMMEDIATE_FLUSH']
c_link_args = ['-s', 'WASM=1', '-s', 'ALLOW_MEMORY_GROWTH=1', '-s', 'MEMORY64=1', '-fwasm-exceptions', '-pthread', '-s', 'EXPORTED_RUNTIME_METHODS=["FS", "callMain"]', '-s', 'STACK_SIZE=10485760']
cpp_link_args = ['-s', 'WASM=1', '-s', 'ALLOW_MEMORY_GROWTH=1', '-s', 'MEMORY64=1', '-fwasm-exceptions', '-pthread', '-s', 'EXPORTED_RUNTIME_METHODS=["FS", "callMain"]', '-s', 'STACK_SIZE=10485760']
[host_machine]
system = 'emscripten'
cpu_family = 'wasm64'
cpu = 'wasm64'
endian = 'little'
[properties]
cmake_toolchain_file = '/home/tboudreaux/Programming/emsdk/upstream/emscripten/cmake/Modules/Platform/Emscripten.cmake'

View File

@@ -20,141 +20,27 @@
# *********************************************************************** # # *********************************************************************** #
project('GridFire', ['c', 'cpp'], version: 'v0.7.4_rc2', default_options: ['cpp_std=c++23'], meson_version: '>=1.5.0') project('GridFire', ['c', 'cpp'], version: 'v0.7.4_rc2', default_options: ['cpp_std=c++23'], meson_version: '>=1.5.0')
if get_option('build-python') # Start by running the code which validates the build environment
add_project_arguments('-fvisibility=hidden', language: 'cpp') subdir('build-check')
else
add_project_arguments('-fvisibility=default', language: 'cpp')
endif
message('Found CXX compiler: ' + meson.get_compiler('cpp').get_id()) # Configure the logging level
message('C++ standard set to: ' + get_option('cpp_std')) subdir('build-extra/log-level')
cppc = meson.get_compiler('cpp')
cc = meson.get_compiler('c')
if cppc.get_id() == 'clang'
message('disabling bitwise-instead-of-logical warnings for clang')
add_project_arguments('-Wno-bitwise-instead-of-logical', language: 'cpp')
endif
if cppc.get_id() == 'gcc'
message('disabling psabi warnings for gcc')
add_project_arguments('-Wno-psabi', language: 'cpp')
if (cppc.version().version_compare('<14.0'))
error('g++ version must be at least 14.0, found ' + cppc.version())
endif
endif
build_fortran = get_option('build-fortran')
if (build_fortran)
add_languages('fortran', native: true)
message('Found FORTRAN compiler: ' + meson.get_compiler('fortran').get_id())
message('Fortran standard set to: ' + get_option('fortran_std'))
message('Building fortran module (gridfire_mod.mod)')
fc = meson.get_compiler('fortran')
if not get_option('unsafe-fortran')
if fc.get_id() != 'gcc'
error('The only supported fortran compiler for GridFire is gfortran (version >= 14.0), found ' + fc + '. GridFire has not been tested with any other compilers. You can disable this check with the -Dunsafe-fortran=true flag to try other compilers')
endif
endif
if (fc.version().version_compare('<14.0'))
error('gfortran version must be at least 14.0, found ' + fc.version())
endif
endif
if not cppc.has_header('print')
error('C++ standard library header <print> not found. Please ensure your compiler and standard library supports C++23. We have already validated your compiler version so this is likely an issue with your standard library installation.')
endif
if not cppc.has_header('format')
error('C++ standard library header <format> not found. Please ensure your compiler and standard library supports C++23. We have already validated your compiler version so this is likely an issue with your standard library installation.')
endif
ignore_unused_args = '-Wno-unused-command-line-argument'
add_global_arguments(ignore_unused_args, language: 'cpp')
add_global_arguments(ignore_unused_args, language: 'c')
# For Eigen
add_project_arguments('-Wno-deprecated-declarations', language: 'cpp')
llevel = get_option('log-level')
logbase='QUILL_COMPILE_ACTIVE_LOG_LEVEL_'
if (llevel == 'traceL3')
message('Setting log level to TRACE_L3')
log_argument = logbase + 'TRACE_L3'
elif (llevel == 'traceL2')
message('Setting log level to TRACE_L2')
log_argument = logbase + 'TRACE_L2'
elif (llevel == 'traceL1')
message('Setting log level to TRACE_L1')
log_argument = logbase + 'TRACE_L1'
elif (llevel == 'debug')
message('Setting log level to DEBUG')
log_argument = logbase + 'DEBUG'
elif (llevel == 'info')
message('Setting log level to INFO')
log_argument = logbase + 'INFO'
elif (llevel == 'warning')
message('Setting log level to WARNING')
log_argument = logbase + 'WARNING'
elif (llevel == 'error')
message('Setting log level to ERROR')
log_argument = logbase + 'ERROR'
elif (llevel == 'critical')
message('Setting log level to CRITICAL')
log_argument = logbase + 'CRITICAL'
endif
log_argument = '-DQUILL_COMPILE_ACTIVE_LOG_LEVEL=' + log_argument
add_project_arguments(log_argument, language: 'cpp')
cpp = meson.get_compiler('cpp')
# Then build the external dependencies
subdir('build-config') subdir('build-config')
# Build the main source code
subdir('src') subdir('src')
if get_option('build-python') # Build the Python bindings
message('Configuring Python bindings...') subdir('build-python')
subdir('build-python')
else
message('Skipping Python bindings...')
endif
if get_option('build-tests') # Buil the test suite
message('Setting up tests for GridFire...') subdir('tests')
subdir('tests')
else subdir('tools')
message('Skipping tests for GridFire...')
endif # Build the pkg-config file
subdir('build-extra/pkg-config')
if get_option('pkg-config')
message('Generating pkg-config file for GridFire...')
pkg = import('pkgconfig')
pkg.generate(
name: 'gridfire',
description: 'GridFire nuclear reaction network solver',
version: meson.project_version(),
libraries: [
libgridfire,
libcomposition,
libconfig,
libconst,
liblogging
],
subdirs: ['gridfire'],
filebase: 'gridfire',
install_dir: join_paths(get_option('libdir'), 'pkgconfig')
)
endif

View File

@@ -1,8 +1,11 @@
option('log-level', type: 'combo', choices: ['traceL3', 'traceL2', 'traceL1', 'debug', 'info', 'warning', 'error', 'critial'], value: 'info', description: 'Set the log level for the GridFire library') option('log_level', type: 'combo', choices: ['traceL3', 'traceL2', 'traceL1', 'debug', 'info', 'warning', 'error', 'critial'], value: 'info', description: 'Set the log level for the GridFire library')
option('pkg-config', type: 'boolean', value: true, description: 'generate pkg-config file for GridFire (gridfire.pc)') option('pkg_config', type: 'boolean', value: true, description: 'generate pkg-config file for GridFire (gridfire.pc)')
option('build-python', type: 'boolean', value: false, description: 'build the python bindings so you can use GridFire from python') option('build_python', type: 'boolean', value: false, description: 'build the python bindings so you can use GridFire from python')
option('build-tests', type: 'boolean', value: true, description: 'build the test suite') option('build_tests', type: 'boolean', value: true, description: 'build the test suite')
option('build-fortran', type: 'boolean', value: false, description: 'build fortran module support') option('build_examples', type: 'boolean', value: true, description: 'build example code')
option('unsafe-fortran', type: 'boolean', value: false, description: 'Allow untested fortran compilers (compilers other than gfortran)') option('build_fortran', type: 'boolean', value: false, description: 'build fortran module support')
option('unity-safe', type: 'boolean', value: false, description: 'Enable safe unity builds for better compatibility across different compilers and platforms') option('unsafe_fortran', type: 'boolean', value: false, description: 'Allow untested fortran compilers (compilers other than gfortran)')
option('python-target-version', type: 'string', value: '3.13', description: 'Target version for python compilation, only used for cross compilation') option('plugin_support', type: 'boolean', value: false, description: 'Enable support for libplugin plugins')
option('python_target_version', type: 'string', value: '3.13', description: 'Target version for python compilation, only used for cross compilation')
option('build_c_api', type: 'boolean', value: true, description: 'compile the C API')
option('build_tools', type: 'boolean', value: true, description: 'build the GridFire command line tools')

View File

@@ -23,7 +23,7 @@ gridfire_extern_dep = declare_dependency(
install_subdir('include/gridfire', install_dir: get_option('includedir')) install_subdir('include/gridfire', install_dir: get_option('includedir'))
if get_option('build-fortran') if get_option('build_fortran')
message('Configuring Fortran bindings...') message('Configuring Fortran bindings...')
subdir('fortran') subdir('fortran')
endif endif

View File

@@ -0,0 +1,35 @@
#pragma once
#include "fourdst/config/config.h"
namespace gridfire::config {
struct CVODESolverConfig {
double absTol = 1.0e-8;
double relTol = 1.0e-5;
};
struct SolverConfig {
CVODESolverConfig cvode;
};
struct AdaptiveEngineViewConfig {
double relativeCullingThreshold = 1.0e-75;
};
struct EngineViewConfig {
AdaptiveEngineViewConfig adaptiveEngineView;
};
struct EngineConfig {
EngineViewConfig views;
};
struct GridFireConfig {
SolverConfig solver;
EngineConfig engine;
};
}

View File

@@ -53,7 +53,7 @@ namespace gridfire::engine {
struct StepDerivatives { struct StepDerivatives {
std::map<fourdst::atomic::Species, T> dydt{}; ///< Derivatives of abundances (dY/dt for each species). std::map<fourdst::atomic::Species, T> dydt{}; ///< Derivatives of abundances (dY/dt for each species).
T nuclearEnergyGenerationRate = T(0.0); ///< Specific energy generation rate (e.g., erg/g/s). T nuclearEnergyGenerationRate = T(0.0); ///< Specific energy generation rate (e.g., erg/g/s).
std::map<fourdst::atomic::Species, std::unordered_map<std::string, T>> reactionContributions{}; std::optional<std::map<fourdst::atomic::Species, std::unordered_map<std::string, T>>> reactionContributions = std::nullopt;
T neutrinoEnergyLossRate = T(0.0); // (erg/g/s) T neutrinoEnergyLossRate = T(0.0); // (erg/g/s)
T totalNeutrinoFlux = T(0.0); // (neutrinos/g/s) T totalNeutrinoFlux = T(0.0); // (neutrinos/g/s)

View File

@@ -12,6 +12,7 @@
#include "gridfire/screening/screening_types.h" #include "gridfire/screening/screening_types.h"
#include "gridfire/partition/partition_abstract.h" #include "gridfire/partition/partition_abstract.h"
#include "gridfire/engine/procedures/construction.h" #include "gridfire/engine/procedures/construction.h"
#include "gridfire/config/config.h"
#include <string> #include <string>
#include <unordered_map> #include <unordered_map>
@@ -753,6 +754,14 @@ namespace gridfire::engine {
[[nodiscard]] [[nodiscard]]
SpeciesStatus getSpeciesStatus(const fourdst::atomic::Species &species) const override; SpeciesStatus getSpeciesStatus(const fourdst::atomic::Species &species) const override;
[[nodiscard]] bool get_store_intermediate_reaction_contributions() const {
return m_store_intermediate_reaction_contributions;
}
void set_store_intermediate_reaction_contributions(const bool value) {
m_store_intermediate_reaction_contributions = value;
}
private: private:
struct PrecomputedReaction { struct PrecomputedReaction {
@@ -847,7 +856,7 @@ namespace gridfire::engine {
const GraphEngine& m_engine; const GraphEngine& m_engine;
}; };
private: private:
Config& m_config = Config::getInstance(); Config<config::GridFireConfig> m_config;
quill::Logger* m_logger = LogManager::getInstance().getLogger("log"); quill::Logger* m_logger = LogManager::getInstance().getLogger("log");
constants m_constants; constants m_constants;
@@ -879,6 +888,7 @@ namespace gridfire::engine {
bool m_usePrecomputation = true; ///< Flag to enable or disable using precomputed reactions for efficiency. Mathematically, this should not change the results. Generally end users should not need to change this. bool m_usePrecomputation = true; ///< Flag to enable or disable using precomputed reactions for efficiency. Mathematically, this should not change the results. Generally end users should not need to change this.
bool m_useReverseReactions = true; ///< Flag to enable or disable reverse reactions. If false, only forward reactions are considered. bool m_useReverseReactions = true; ///< Flag to enable or disable reverse reactions. If false, only forward reactions are considered.
bool m_store_intermediate_reaction_contributions = false; ///< Flag to enable or disable storing intermediate reaction contributions for debugging.
BuildDepthType m_depth; BuildDepthType m_depth;
@@ -1207,7 +1217,10 @@ namespace gridfire::engine {
const T nu_ij = static_cast<T>(reaction.stoichiometry(species)); const T nu_ij = static_cast<T>(reaction.stoichiometry(species));
const T dydt_increment = threshold_flag * molarReactionFlow * nu_ij; const T dydt_increment = threshold_flag * molarReactionFlow * nu_ij;
dydt_vec[speciesIdx] += dydt_increment; dydt_vec[speciesIdx] += dydt_increment;
result.reactionContributions[species][std::string(reaction.id())] = dydt_increment;
if (m_store_intermediate_reaction_contributions) {
result.reactionContributions.value()[species][std::string(reaction.id())] = dydt_increment;
}
} }
} }

View File

@@ -4,6 +4,7 @@
#include "gridfire/screening/screening_abstract.h" #include "gridfire/screening/screening_abstract.h"
#include "gridfire/screening/screening_types.h" #include "gridfire/screening/screening_types.h"
#include "gridfire/types/types.h" #include "gridfire/types/types.h"
#include "gridfire/config/config.h"
#include "fourdst/atomic/atomicSpecies.h" #include "fourdst/atomic/atomicSpecies.h"
#include "fourdst/config/config.h" #include "fourdst/config/config.h"
@@ -386,10 +387,10 @@ namespace gridfire::engine {
*/ */
[[nodiscard]] SpeciesStatus getSpeciesStatus(const fourdst::atomic::Species &species) const override; [[nodiscard]] SpeciesStatus getSpeciesStatus(const fourdst::atomic::Species &species) const override;
private: private:
using Config = fourdst::config::Config;
using LogManager = fourdst::logging::LogManager; using LogManager = fourdst::logging::LogManager;
/** @brief A reference to the singleton Config instance, used for retrieving configuration parameters. */
Config& m_config = Config::getInstance(); fourdst::config::Config<config::GridFireConfig> m_config;
/** @brief A pointer to the logger instance, used for logging messages. */ /** @brief A pointer to the logger instance, used for logging messages. */
quill::Logger* m_logger = LogManager::getInstance().getLogger("log"); quill::Logger* m_logger = LogManager::getInstance().getLogger("log");

View File

@@ -6,6 +6,8 @@
#include "gridfire/io/network_file.h" #include "gridfire/io/network_file.h"
#include "gridfire/types/types.h" #include "gridfire/types/types.h"
#include "gridfire/config/config.h"
#include "fourdst/config/config.h" #include "fourdst/config/config.h"
#include "fourdst/logging/logging.h" #include "fourdst/logging/logging.h"
@@ -365,9 +367,9 @@ namespace gridfire::engine {
[[nodiscard]] std::string getNetworkFile() const { return m_fileName; } [[nodiscard]] std::string getNetworkFile() const { return m_fileName; }
[[nodiscard]] const io::NetworkFileParser& getParser() const { return m_parser; } [[nodiscard]] const io::NetworkFileParser& getParser() const { return m_parser; }
private: private:
using Config = fourdst::config::Config; using LogManager = LogManager;
using LogManager = fourdst::logging::LogManager; Config<config::GridFireConfig> m_config;
Config& m_config = Config::getInstance();
quill::Logger* m_logger = LogManager::getInstance().getLogger("log"); quill::Logger* m_logger = LogManager::getInstance().getLogger("log");
std::string m_fileName; std::string m_fileName;
///< Parser for the network file. ///< Parser for the network file.

View File

@@ -2,6 +2,7 @@
#include "fourdst/config/config.h" #include "fourdst/config/config.h"
#include "fourdst/logging/logging.h" #include "fourdst/logging/logging.h"
#include "gridfire/config/config.h"
#include "quill/Logger.h" #include "quill/Logger.h"
@@ -101,9 +102,8 @@ namespace gridfire::io {
*/ */
[[nodiscard]] ParsedNetworkData parse(const std::string& filename) const override; [[nodiscard]] ParsedNetworkData parse(const std::string& filename) const override;
private: private:
using Config = fourdst::config::Config;
using LogManager = fourdst::logging::LogManager; using LogManager = fourdst::logging::LogManager;
Config& m_config = Config::getInstance(); fourdst::config::Config<config::GridFireConfig> m_config;
quill::Logger* m_logger = LogManager::getInstance().getLogger("log"); quill::Logger* m_logger = LogManager::getInstance().getLogger("log");
}; };
@@ -141,9 +141,8 @@ namespace gridfire::io {
*/ */
[[nodiscard]] ParsedNetworkData parse(const std::string& filename) const override; [[nodiscard]] ParsedNetworkData parse(const std::string& filename) const override;
private: private:
using Config = fourdst::config::Config;
using LogManager = fourdst::logging::LogManager; using LogManager = fourdst::logging::LogManager;
Config& m_config = Config::getInstance(); fourdst::config::Config<config::GridFireConfig> m_config;
quill::Logger* m_logger = LogManager::getInstance().getLogger("log"); quill::Logger* m_logger = LogManager::getInstance().getLogger("log");
std::string m_filename; std::string m_filename;

View File

@@ -4,6 +4,7 @@
#include "gridfire/engine/engine_abstract.h" #include "gridfire/engine/engine_abstract.h"
#include "gridfire/types/types.h" #include "gridfire/types/types.h"
#include "gridfire/exceptions/exceptions.h" #include "gridfire/exceptions/exceptions.h"
#include "gridfire/config/config.h"
#include "fourdst/atomic/atomicSpecies.h" #include "fourdst/atomic/atomicSpecies.h"
#include "fourdst/config/config.h" #include "fourdst/config/config.h"
@@ -237,13 +238,13 @@ namespace gridfire::solver {
}; };
struct CVODERHSOutputData { struct CVODERHSOutputData {
std::map<fourdst::atomic::Species, std::unordered_map<std::string, double>> reaction_contribution_map; std::optional<std::map<fourdst::atomic::Species, std::unordered_map<std::string, double>>> reaction_contribution_map;
double neutrino_energy_loss_rate; double neutrino_energy_loss_rate;
double total_neutrino_flux; double total_neutrino_flux;
}; };
private: private:
fourdst::config::Config& m_config = fourdst::config::Config::getInstance(); fourdst::config::Config<config::GridFireConfig> m_config;
quill::Logger* m_logger = fourdst::logging::LogManager::getInstance().getLogger("log"); quill::Logger* m_logger = fourdst::logging::LogManager::getInstance().getLogger("log");
/** /**
* @brief CVODE RHS C-wrapper that delegates to calculate_rhs and captures exceptions. * @brief CVODE RHS C-wrapper that delegates to calculate_rhs and captures exceptions.

View File

@@ -684,7 +684,7 @@ namespace gridfire::engine {
// --- Efficient lookup of only the active reactions --- // --- Efficient lookup of only the active reactions ---
uint64_t reactionHash = utils::hash_reaction(*reaction); uint64_t reactionHash = utils::hash_reaction(*reaction);
const size_t reactionIndex = m_precomputedReactionIndexMap.at(reactionHash); const size_t reactionIndex = m_precomputedReactionIndexMap.at(reactionHash);
PrecomputedReaction precomputedReaction = m_precomputedReactions[reactionIndex]; const PrecomputedReaction& precomputedReaction = m_precomputedReactions[reactionIndex];
// --- Forward abundance product --- // --- Forward abundance product ---
double forwardAbundanceProduct = 1.0; double forwardAbundanceProduct = 1.0;
@@ -697,12 +697,12 @@ namespace gridfire::engine {
forwardAbundanceProduct = 0.0; forwardAbundanceProduct = 0.0;
break; // No need to continue if one of the reactants has zero abundance break; // No need to continue if one of the reactants has zero abundance
} }
double factor = std::pow(comp.getMolarAbundance(reactant), power); const double factor = std::pow(comp.getMolarAbundance(reactant), power);
if (!std::isfinite(factor)) { if (!std::isfinite(factor)) {
LOG_CRITICAL(m_logger, "Non-finite factor encountered in forward abundance product for reaction '{}'. Check input abundances for validity.", reaction->id()); LOG_CRITICAL(m_logger, "Non-finite factor encountered in forward abundance product for reaction '{}'. Check input abundances for validity.", reaction->id());
throw exceptions::BadRHSEngineError("Non-finite factor encountered in forward abundance product."); throw exceptions::BadRHSEngineError("Non-finite factor encountered in forward abundance product.");
} }
forwardAbundanceProduct *= std::pow(comp.getMolarAbundance(reactant), power); forwardAbundanceProduct *= factor;
} }
const double bare_rate = bare_rates.at(reactionCounter); const double bare_rate = bare_rates.at(reactionCounter);
@@ -764,8 +764,8 @@ namespace gridfire::engine {
default: ; default: ;
} }
double local_neutrino_loss = molarReactionFlows.back() * q_abs * neutrino_loss_fraction * m_constants.Na * m_constants.MeV_to_erg; const double local_neutrino_loss = molarReactionFlows.back() * q_abs * neutrino_loss_fraction * m_constants.Na * m_constants.MeV_to_erg;
double local_neutrino_flux = molarReactionFlows.back() * m_constants.Na; const double local_neutrino_flux = molarReactionFlows.back() * m_constants.Na;
result.totalNeutrinoFlux += local_neutrino_flux; result.totalNeutrinoFlux += local_neutrino_flux;
result.neutrinoEnergyLossRate += local_neutrino_loss; result.neutrinoEnergyLossRate += local_neutrino_loss;
@@ -782,7 +782,7 @@ namespace gridfire::engine {
reactionCounter = 0; reactionCounter = 0;
for (const auto& reaction: activeReactions) { for (const auto& reaction: activeReactions) {
size_t j = m_precomputedReactionIndexMap.at(utils::hash_reaction(*reaction)); const size_t j = m_precomputedReactionIndexMap.at(utils::hash_reaction(*reaction));
const auto& precomp = m_precomputedReactions[j]; const auto& precomp = m_precomputedReactions[j];
const double R_j = molarReactionFlows[reactionCounter]; const double R_j = molarReactionFlows[reactionCounter];
@@ -793,9 +793,12 @@ namespace gridfire::engine {
const int stoichiometricCoefficient = precomp.stoichiometric_coefficients[i]; const int stoichiometricCoefficient = precomp.stoichiometric_coefficients[i];
// Update the derivative for this species // Update the derivative for this species
double dydt_increment = static_cast<double>(stoichiometricCoefficient) * R_j; const double dydt_increment = static_cast<double>(stoichiometricCoefficient) * R_j;
result.dydt.at(species) += dydt_increment; result.dydt.at(species) += dydt_increment;
result.reactionContributions[species][std::string(reaction->id())] = dydt_increment;
if (m_store_intermediate_reaction_contributions) {
result.reactionContributions.value()[species][std::string(reaction->id())] = dydt_increment;
}
} }
reactionCounter++; reactionCounter++;
} }

View File

@@ -394,7 +394,9 @@ namespace gridfire::engine {
const double maxFlow const double maxFlow
) const { ) const {
LOG_TRACE_L1(m_logger, "Culling reactions based on flow rates..."); LOG_TRACE_L1(m_logger, "Culling reactions based on flow rates...");
const auto relative_culling_threshold = m_config.get<double>("gridfire:AdaptiveEngineView:RelativeCullingThreshold", 1e-75);
const auto relative_culling_threshold = m_config->engine.views.adaptiveEngineView.relativeCullingThreshold;
double absoluteCullingThreshold = relative_culling_threshold * maxFlow; double absoluteCullingThreshold = relative_culling_threshold * maxFlow;
LOG_DEBUG(m_logger, "Relative culling threshold: {:7.3E} ({:7.3E})", relative_culling_threshold, absoluteCullingThreshold); LOG_DEBUG(m_logger, "Relative culling threshold: {:7.3E} ({:7.3E})", relative_culling_threshold, absoluteCullingThreshold);
std::vector<const reaction::Reaction*> culledReactions; std::vector<const reaction::Reaction*> culledReactions;

View File

@@ -112,8 +112,8 @@ namespace gridfire::solver {
// 2. If the user has set tolerances in code, those override the config // 2. If the user has set tolerances in code, those override the config
// 3. If the user has not set tolerances in code and the config does not have them, use hardcoded defaults // 3. If the user has not set tolerances in code and the config does not have them, use hardcoded defaults
auto absTol = m_config.get<double>("gridfire:solver:CVODESolverStrategy:absTol", 1.0e-8); auto absTol = m_config->solver.cvode.absTol;
auto relTol = m_config.get<double>("gridfire:solver:CVODESolverStrategy:relTol", 1.0e-5); auto relTol = m_config->solver.cvode.relTol;
if (m_absTol) { if (m_absTol) {
absTol = *m_absTol; absTol = *m_absTol;
@@ -935,8 +935,8 @@ namespace gridfire::solver {
sunrealtype *y_data = N_VGetArrayPointer(m_Y); sunrealtype *y_data = N_VGetArrayPointer(m_Y);
sunrealtype *y_err_data = N_VGetArrayPointer(m_YErr); sunrealtype *y_err_data = N_VGetArrayPointer(m_YErr);
const auto absTol = m_config.get<double>("gridfire:solver:CVODESolverStrategy:absTol", 1.0e-8); const auto absTol = m_config->solver.cvode.absTol;
const auto relTol = m_config.get<double>("gridfire:solver:CVODESolverStrategy:relTol", 1.0e-8); const auto relTol = m_config->solver.cvode.relTol;
std::vector<double> err_ratios; std::vector<double> err_ratios;
const size_t num_components = N_VGetLength(m_Y); const size_t num_components = N_VGetLength(m_Y);

View File

@@ -42,7 +42,7 @@ gridfire_build_dependencies = [
json_dep, json_dep,
] ]
if not get_option('unity-safe') if get_option('plugin_support')
gridfire_build_dependencies += [plugin_dep] gridfire_build_dependencies += [plugin_dep]
endif endif
@@ -63,12 +63,11 @@ gridfire_dep = declare_dependency(
install_subdir('include/gridfire', install_dir: get_option('includedir')) install_subdir('include/gridfire', install_dir: get_option('includedir'))
message('Configuring C API...')
subdir('extern') if not get_option('build_c_api') and get_option('build_fortran')
# error('Cannot build fortran without C API. Set -Dbuild-c-api=true and -Dbuild-fortran=true')
#if get_option('build-python') endif
# message('Configuring Python bindings...') if get_option('build_c_api')
# subdir('python') message('Configuring C API...')
#else subdir('extern')
# message('Skipping Python bindings...') endif
#endif

View File

@@ -1,4 +1,4 @@
[wrap-git] [wrap-git]
url = https://github.com/4D-STAR/fourdst url = https://github.com/4D-STAR/fourdst
revision = v0.9.11 revision = v0.9.14
depth = 1 depth = 1

View File

@@ -1,5 +1,7 @@
subdir('C') if get_option('build_c_api')
subdir('C')
endif
if get_option('build-fortran') if get_option('build_fortran')
subdir('fortran') subdir('fortran')
endif endif

View File

@@ -1,5 +1,5 @@
executable( executable(
'graphnet_sandbox', 'graphnet_sandbox',
'main.cpp', 'main.cpp',
dependencies: [gridfire_dep, composition_dep, cli11_dep], dependencies: [gridfire_dep, cli11_dep],
) )

View File

@@ -1,7 +1,7 @@
# Google Test dependency # Google Test dependency
gtest_dep = dependency('gtest', main: true, required : true) #gtest_dep = dependency('gtest', main: true, required : true)
gtest_main = dependency('gtest_main', required: true) #gtest_main = dependency('gtest_main', required: true)
gtest_nomain_dep = dependency('gtest', main: false, required : true) #gtest_nomain_dep = dependency('gtest', main: false, required : true)
# Subdirectories for unit and integration tests # Subdirectories for unit and integration tests
subdir('graphnet_sandbox') subdir('graphnet_sandbox')

View File

@@ -0,0 +1,48 @@
#include "fourdst/config/config.h"
#include "gridfire/config/config.h"
#include <source_location>
#include <filesystem>
#include "CLI/CLI.hpp"
consteval std::string_view strip_namespaces(const std::string_view fullName) {
const size_t pos = fullName.rfind("::");
if (pos == std::string_view::npos) {
return fullName;
}
return fullName.substr(pos + 2);
}
template <typename T>
consteval std::string_view get_type_name() {
constexpr std::string_view name = std::source_location::current().function_name();
const auto pos = name.find("T = ");
if (pos == std::string_view::npos) return name;
const auto start = pos + 4;
const auto end = name.rfind(']');
return name.substr(start, end - start);
}
int main(int argc, char** argv) {
CLI::App app{"GridFire Sandbox Application."};
std::string outputPath = ".";
app.add_option("-p,--path", outputPath, "path to save generated config files (default: current directory)");
CLI11_PARSE(app, argc, argv);
const std::filesystem::path outPath(outputPath);
if (!std::filesystem::exists(outPath)) {
std::cerr << "Error: The specified path does not exist: " << outputPath << std::endl;
return 1;
}
fourdst::config::Config<gridfire::config::GridFireConfig> configConfig;
const std::string_view name = strip_namespaces(get_type_name<gridfire::config::GridFireConfig>());
const std::string defaultConfigFilePath = (outPath / (std::string(name) + ".toml")).string();
const std::string schemaFilePath = (outPath / (std::string(name) + ".schema.json")).string();
configConfig.save(defaultConfigFilePath);
configConfig.save_schema(schemaFilePath);
}

1
tools/config/meson.build Normal file
View File

@@ -0,0 +1 @@
executable('gf_generate_config_file', 'generate_config_files.cpp', dependencies: [gridfire_dep, cli11_dep], install: true)

3
tools/meson.build Normal file
View File

@@ -0,0 +1,3 @@
if get_option('build_tools')
subdir('config')
endif