Initial commit

.clang-format

+---
+# https://clang.llvm.org/docs/ClangFormatStyleOptions.html#
+# the "Visual Studio" style is similar to:
+BasedOnStyle: LLVM
+UseTab: Never
+IndentWidth: 4
+BreakBeforeBraces: Allman
+AccessModifierOffset: -4
+SortIncludes: false
+ColumnLimit: 0
+...

.gitignore

+# Compiled Object files
+*.slo
+*.lo
+*.o
+
+# Compiled Dynamic libraries
+*.so
+*.dylib
+
+# Compiled Static libraries
+*.lai
+*.la
+*.a
+
+# Misc
+core
+*.pyc
+*.*~
+
+# OS
+.DS_Store
+*.swp
+*.bak
+
+# Workspace
+workspace
+*.tar.gz
+*.tgz
+
+# build dir
+build
+
+# gitlab-ci / clang-format
+patches
+
+# Gmsh
+*.pos
+*.msh # do not commit the mesh...
+!*_lfs.msh # ... except in the lfs
+
+# IDE
+.project
+.pydevproject
+.settings
+.vscode
+.vscode/*

.gitlab-ci.yml

+default:
+    image: rboman/waves-py3:2020.3
+    before_script:
+        - source /opt/intel/mkl/bin/mklvars.sh intel64
+        - source /opt/intel/tbb/bin/tbbvars.sh intel64
+        - echo $(nproc)
+        - printenv | sort
+
+.global_tag: &global_tag_def
+    tags:
+        - mn2l # you can choose a set of runners here
+
+variables:
+    GIT_STRATEGY: clone # workaround full clone for each pipeline (https://gitlab.com/gitlab-org/gitlab-runner/-/issues/26993)
+    GIT_LFS_SKIP_SMUDGE: 1 # do not pull LFS
+
+stages:
+    - build
+    - test
+
+format:
+    <<: *global_tag_def
+    stage: build
+    script:
+        - clang-format --version # we use clang-format-10 exclusively
+        - ./scripts/format_code.py
+        - mkdir -p patches
+        - if git diff --patch --exit-code > patches/clang-format.patch; then echo "Clang format changed nothing"; else echo "Clang format found changes to make!"; false; fi
+    artifacts:
+        paths:
+            - patches/
+        expire_in: 1 day
+        when: on_failure
+    allow_failure: true
+
+build:
+    <<: *global_tag_def
+    stage: build
+    script:
+        - rm -rf build workspace
+        - mkdir build
+        - cd build
+        - cmake -Wno-dev ..
+        - make -j 8
+    artifacts:
+        paths:
+            - build/
+        expire_in: 1 day
+
+doxygen:
+    <<: *global_tag_def
+    stage: test
+    script:
+        - cd build
+        - make dox
+    artifacts:
+        paths:
+            - build/doxygen/
+        expire_in: 1 week
+    dependencies:
+        - build
+
+ctest:
+    <<: *global_tag_def
+    stage: test
+    script:
+        - cd build
+        - ctest --output-on-failure -j 8
+    #timeout: 10 hours  # will be available in 12.3
+    dependencies:
+        - build

CMake/FindEIGEN.cmake

+# Copyright 2020 University of Liège
+# 
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+# 
+#     http://www.apache.org/licenses/LICENSE-2.0
+# 
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+# FindEIGEN.cmake - try to find Eigen 3 headers
+# The logic that checks the version comes from FindEigen3.cmake (source dir)
+# ----------------------------------------------------------------------------   
+# output:
+#    EIGEN_FOUND             : TRUE/FALSE
+#    EIGEN_INCLUDE_DIRS      : where the Eigen/*.h are              [cached]
+#    EIGEN_VERSION           : version number of Eigen
+# ----------------------------------------------------------------------------
+# autodetection:
+#    utiliser "CMAKE_PREFIX_PATH=c:\local"
+#          ou "CMAKE_INCLUDE_PATH=c:\local\include"
+#          ou "INCLUDE=c:\local\include"
+# ----------------------------------------------------------------------------   
+
+# Check version (only for Eigen 3)
+macro(_eigen3_check_version)
+    file(READ "${EIGEN_INCLUDE_DIRS}/Eigen/src/Core/util/Macros.h" _eigen3_version_header)
+
+    string(REGEX MATCH "define[ \t]+EIGEN_WORLD_VERSION[ \t]+([0-9]+)" _eigen3_world_version_match "${_eigen3_version_header}")
+    set(EIGEN_WORLD_VERSION "${CMAKE_MATCH_1}")
+    string(REGEX MATCH "define[ \t]+EIGEN_MAJOR_VERSION[ \t]+([0-9]+)" _eigen3_major_version_match "${_eigen3_version_header}")
+    set(EIGEN_MAJOR_VERSION "${CMAKE_MATCH_1}")
+    string(REGEX MATCH "define[ \t]+EIGEN_MINOR_VERSION[ \t]+([0-9]+)" _eigen3_minor_version_match "${_eigen3_version_header}")
+    set(EIGEN_MINOR_VERSION "${CMAKE_MATCH_1}")
+
+    set(EIGEN_VERSION ${EIGEN_WORLD_VERSION}.${EIGEN_MAJOR_VERSION}.${EIGEN_MINOR_VERSION})
+
+    if(${EIGEN_VERSION} VERSION_LESS ${EIGEN_FIND_VERSION})
+        message(FATAL_ERROR "Eigen version ${EIGEN_VERSION} found in ${EIGEN_INCLUDE_DIRS}, but at least version ${EIGEN_FIND_VERSION} is required!")
+    endif()
+endmacro()
+
+# Set a dummy version if not provided
+if(NOT EIGEN_FIND_VERSION)
+    if(NOT EIGEN_FIND_VERSION_MAJOR)
+        set(EIGEN_FIND_VERSION_MAJOR 3)
+    endif()
+    if(NOT EIGEN_FIND_VERSION_MINOR)
+        set(EIGEN_FIND_VERSION_MINOR 3)
+    endif()
+    if(NOT EIGEN_FIND_VERSION_PATCH)
+        set(EIGEN_FIND_VERSION_PATCH 4)
+    endif()
+
+    set(EIGEN_FIND_VERSION "${EIGEN_FIND_VERSION_MAJOR}.${EIGEN_FIND_VERSION_MINOR}.${EIGEN_FIND_VERSION_PATCH}")
+endif()
+
+# Find the header and check the version
+find_path(EIGEN_INCLUDE_DIRS "Eigen/Dense" PATHS "/usr/include/eigen3")
+if (EIGEN_INCLUDE_DIRS)
+    _eigen3_check_version()
+else()
+    message(FATAL_ERROR "Eigen3 headers not found! Define the path in the INCLUDE environment variable.")
+endif()
+
+# handle the QUIETLY and REQUIRED arguments and set EIGEN_FOUND to TRUE
+# if all listed variables are TRUE
+include(FindPackageHandleStandardArgs)
+find_package_handle_standard_args(EIGEN
+                                  FOUND_VAR EIGEN_FOUND
+                                  REQUIRED_VARS EIGEN_INCLUDE_DIRS
+                                  VERSION_VAR EIGEN_VERSION)
+
+#MESSAGE(STATUS "EIGEN_FOUND = ${EIGEN_FOUND}")

CMake/sdpmMacros.cmake

+# Copyright 2023 University of Liège
+# 
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+# 
+#     http://www.apache.org/licenses/LICENSE-2.0
+# 
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+# useful macros/fcts
+
+MACRO(MACRO_AddTest srcDir)
+    file(GLOB tfiles RELATIVE ${srcDir} ${srcDir}/*)
+    foreach(tfile ${tfiles})
+        set(spath ${srcDir}/${tfile})
+        if(NOT IS_DIRECTORY ${spath} AND 
+           ${spath} MATCHES ".+\\.py$" AND 
+           NOT ${tfile} STREQUAL "__init__.py")
+            string(REPLACE "${PROJECT_SOURCE_DIR}/" "" strip ${spath}) 
+            message(STATUS "Adding test ${strip}")
+            add_test(NAME ${strip} 
+                    WORKING_DIRECTORY ${PROJECT_SOURCE_DIR} 
+                    COMMAND ${Python3_EXECUTABLE} "${PROJECT_SOURCE_DIR}/run.py" --nogui --clean ${strip})            
+        else()
+            MACRO_AddTest(${srcDir}/${tfile})
+        endif()
+    endforeach()
+ENDMACRO()
+
+MACRO(MACRO_DebugPostfix libname)
+    IF(MSVC)
+        SET_TARGET_PROPERTIES(${libname} PROPERTIES DEBUG_POSTFIX "_d")
+    ENDIF(MSVC)
+ENDMACRO()

CMakeLists.txt

+# Copyright 2020 University of Liège
+# 
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+# 
+#     http://www.apache.org/licenses/LICENSE-2.0
+# 
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+# ----------------------------------------------------------------------------
+PROJECT(SDPM)
+# ----------------------------------------------------------------------------
+CMAKE_MINIMUM_REQUIRED(VERSION 3.14)
+
+# -- I/O
+# Lib/Exe dir
+SET(LIBRARY_OUTPUT_PATH ${PROJECT_BINARY_DIR}/bin CACHE PATH
+                        "Single output directory for building all libraries.")
+SET(EXECUTABLE_OUTPUT_PATH ${PROJECT_BINARY_DIR}/bin CACHE PATH
+                        "Single output directory for building all executables.")
+MARK_AS_ADVANCED(LIBRARY_OUTPUT_PATH EXECUTABLE_OUTPUT_PATH)
+
+# Build type
+IF(NOT CMAKE_BUILD_TYPE)
+    SET( CMAKE_BUILD_TYPE "Release" CACHE STRING 
+         "Choose the type of build, options are: None Debug Release RelWithDebInfo MinSizeRel."
+         FORCE)
+ENDIF(NOT CMAKE_BUILD_TYPE)
+
+# Additional modules and macros
+LIST(APPEND CMAKE_MODULE_PATH "${PROJECT_SOURCE_DIR}/CMake")
+INCLUDE(sdpmMacros)
+
+# -- C/C++
+# Set specific languages flags
+SET(CMAKE_CXX_STANDARD 11) # newer way to set C++11 (requires cmake>=3.1)
+SET(CMAKE_CXX_STANDARD_REQUIRED ON)
+IF((CMAKE_CXX_COMPILER_ID MATCHES "GNU") OR (CMAKE_CXX_COMPILER_ID MATCHES "Intel"))
+    IF(NOT APPLE)
+        SET(CMAKE_SHARED_LINKER_FLAGS "-Wl,--no-as-needed")
+    ENDIF()
+    SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Wall") # add verbosity
+ELSEIF(CMAKE_CXX_COMPILER_ID MATCHES "MSVC")
+    ADD_DEFINITIONS(-D_CRT_SECURE_NO_WARNINGS -D_SCL_SECURE_NO_DEPRECATE)
+    ADD_DEFINITIONS(-D_USE_MATH_DEFINES) # otherwise M_PI is undefined
+    SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} /MP")  # parallel build with MSVC
+    #SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} /W4")  # add verbosity
+ELSEIF(CMAKE_CXX_COMPILER_ID MATCHES "Clang")
+    SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Wno-deprecated-register")
+    #SET(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -Weverything") # add verbosity
+ENDIF()
+
+# -- DEPENDENCIES
+# Python
+# use Python3_ROOT_DIR if wrong python found (e.g. anaconda)
+FIND_PACKAGE(Python3 COMPONENTS Interpreter Development REQUIRED)   
+
+# SWIG
+FIND_PACKAGE(SWIG REQUIRED)
+IF(CMAKE_GENERATOR MATCHES "Visual Studio") # not MSVC because of nmake & jom
+    SET(CMAKE_SWIG_OUTDIR "${EXECUTABLE_OUTPUT_PATH}/$(Configuration)/")
+ELSE()
+    SET(CMAKE_SWIG_OUTDIR "${EXECUTABLE_OUTPUT_PATH}")
+ENDIF()
+
+# Eigen
+FIND_PACKAGE(EIGEN 3.3.4 REQUIRED)
+
+# Doxygen (https://vicrucann.github.io/tutorials/quick-cmake-doxygen/)
+FIND_PACKAGE(Doxygen) # check if Doxygen is installed
+IF(DOXYGEN_FOUND)
+    # set input and output files
+    SET(DOXYGEN_IN ${PROJECT_SOURCE_DIR}/Doxyfile.in)
+    SET(DOXYGEN_OUT ${PROJECT_BINARY_DIR}/Doxyfile)
+    # request to configure the file
+    CONFIGURE_FILE(${DOXYGEN_IN} ${DOXYGEN_OUT} @ONLY)
+    # note the option ALL which allows to build the docs together with the application
+    # "make dox" rebuilds the doc
+    ADD_CUSTOM_TARGET( dox #ALL
+        COMMAND ${DOXYGEN_EXECUTABLE} ${DOXYGEN_OUT}
+        WORKING_DIRECTORY ${CMAKE_CURRENT_BINARY_DIR}
+        COMMENT "Generating API documentation with Doxygen"
+        VERBATIM )
+ELSE()
+    MESSAGE("Doxygen needs to be installed to generate the doxygen documentation")
+ENDIF()
+
+# -- CTest
+ENABLE_TESTING()
+
+# -- INSTALL
+IF(CMAKE_INSTALL_PREFIX_INITIALIZED_TO_DEFAULT)
+    EXECUTE_PROCESS(COMMAND ${Python3_EXECUTABLE} -m site --user-site OUTPUT_VARIABLE PY_SITE OUTPUT_STRIP_TRAILING_WHITESPACE)
+    STRING(REGEX REPLACE "\\\\" "/" PY_SITE ${PY_SITE})
+    SET(CMAKE_INSTALL_PREFIX "${PY_SITE}/sdpm" CACHE STRING "Install location" FORCE)
+    SET(CMAKE_INSTALL_PREFIX_INITIALIZED_TO_DEFAULT FALSE)
+ELSE()
+    IF(NOT(CMAKE_INSTALL_PREFIX MATCHES "sdpm$"))
+        SET(CMAKE_INSTALL_PREFIX "${CMAKE_INSTALL_PREFIX}/sdpm" CACHE STRING "Install location" FORCE)
+    ENDIF()
+ENDIF()
+# RPATH for Linux
+IF(UNIX AND NOT APPLE)
+    SET(CMAKE_INSTALL_RPATH "${CMAKE_INSTALL_PREFIX}")
+ENDIF()
+# Copy init.py as __init__.py to install location
+INSTALL(FILES ${CMAKE_CURRENT_SOURCE_DIR}/init.py DESTINATION ${CMAKE_INSTALL_PREFIX} RENAME __init__.py)
+
+# -- Sub directories
+ADD_SUBDIRECTORY( sdpm )
+
+# -- FINAL
+MESSAGE(STATUS "PROJECT: ${CMAKE_PROJECT_NAME}")
+MESSAGE(STATUS "* SYSTEM NAME=\"${CMAKE_SYSTEM_NAME}\"")
+MESSAGE(STATUS "* CXX COMPILER: ${CMAKE_CXX_COMPILER_ID}")
+MESSAGE(STATUS "* CXX STANDARD: ${CMAKE_CXX_STANDARD}")
+MESSAGE(STATUS "* INSTALL DIR: ${CMAKE_INSTALL_PREFIX}")
+MESSAGE(STATUS "* BUILD TYPE: ${CMAKE_BUILD_TYPE}")

CODEOWNERS

+# This is a comment.
+# Each line is a file pattern followed by one or more owners.
+
+# These owners will be the default owners for everything in
+# the repo. Unless a later match takes precedence, they
+# will be requested for review when someone opens a pull request.
+*       @rboman
+*       @acrovato
+
+# These owners own a directory and nested subdirectories
+/sdpm/ @acrovato
+
+# You can also use email addresses if you prefer. They'll be
+# used to look up users just like we do for commit author
+# emails.
+#*.go docs@example.com
+

CREDITS.md

+SDPM is Copyright (C) of University of Liège.
+
+The main author of SDPM is Adrien Crovato.
+
+Direct and indirect contributions have been made to SDPM by the following people:
+- Grigorios Dimitriadis
+- Romain Boman
+- Luc Papeleux

LICENSE

+                                 Apache License
+                           Version 2.0, January 2004
+                        http://www.apache.org/licenses/
+
+   TERMS AND CONDITIONS FOR USE, REPRODUCTION, AND DISTRIBUTION
+
+   1. Definitions.
+
+      "License" shall mean the terms and conditions for use, reproduction,
+      and distribution as defined by Sections 1 through 9 of this document.
+
+      "Licensor" shall mean the copyright owner or entity authorized by
+      the copyright owner that is granting the License.
+
+      "Legal Entity" shall mean the union of the acting entity and all
+      other entities that control, are controlled by, or are under common
+      control with that entity. For the purposes of this definition,
+      "control" means (i) the power, direct or indirect, to cause the
+      direction or management of such entity, whether by contract or
+      otherwise, or (ii) ownership of fifty percent (50%) or more of the
+      outstanding shares, or (iii) beneficial ownership of such entity.
+
+      "You" (or "Your") shall mean an individual or Legal Entity
+      exercising permissions granted by this License.
+
+      "Source" form shall mean the preferred form for making modifications,
+      including but not limited to software source code, documentation
+      source, and configuration files.
+
+      "Object" form shall mean any form resulting from mechanical
+      transformation or translation of a Source form, including but
+      not limited to compiled object code, generated documentation,
+      and conversions to other media types.
+
+      "Work" shall mean the work of authorship, whether in Source or
+      Object form, made available under the License, as indicated by a
+      copyright notice that is included in or attached to the work
+      (an example is provided in the Appendix below).
+
+      "Derivative Works" shall mean any work, whether in Source or Object
+      form, that is based on (or derived from) the Work and for which the
+      editorial revisions, annotations, elaborations, or other modifications
+      represent, as a whole, an original work of authorship. For the purposes
+      of this License, Derivative Works shall not include works that remain
+      separable from, or merely link (or bind by name) to the interfaces of,
+      the Work and Derivative Works thereof.
+
+      "Contribution" shall mean any work of authorship, including
+      the original version of the Work and any modifications or additions
+      to that Work or Derivative Works thereof, that is intentionally
+      submitted to Licensor for inclusion in the Work by the copyright owner
+      or by an individual or Legal Entity authorized to submit on behalf of
+      the copyright owner. For the purposes of this definition, "submitted"
+      means any form of electronic, verbal, or written communication sent
+      to the Licensor or its representatives, including but not limited to
+      communication on electronic mailing lists, source code control systems,
+      and issue tracking systems that are managed by, or on behalf of, the
+      Licensor for the purpose of discussing and improving the Work, but
+      excluding communication that is conspicuously marked or otherwise
+      designated in writing by the copyright owner as "Not a Contribution."
+
+      "Contributor" shall mean Licensor and any individual or Legal Entity
+      on behalf of whom a Contribution has been received by Licensor and
+      subsequently incorporated within the Work.
+
+   2. Grant of Copyright License. Subject to the terms and conditions of
+      this License, each Contributor hereby grants to You a perpetual,
+      worldwide, non-exclusive, no-charge, royalty-free, irrevocable
+      copyright license to reproduce, prepare Derivative Works of,
+      publicly display, publicly perform, sublicense, and distribute the
+      Work and such Derivative Works in Source or Object form.
+
+   3. Grant of Patent License. Subject to the terms and conditions of
+      this License, each Contributor hereby grants to You a perpetual,
+      worldwide, non-exclusive, no-charge, royalty-free, irrevocable
+      (except as stated in this section) patent license to make, have made,
+      use, offer to sell, sell, import, and otherwise transfer the Work,
+      where such license applies only to those patent claims licensable
+      by such Contributor that are necessarily infringed by their
+      Contribution(s) alone or by combination of their Contribution(s)
+      with the Work to which such Contribution(s) was submitted. If You
+      institute patent litigation against any entity (including a
+      cross-claim or counterclaim in a lawsuit) alleging that the Work
+      or a Contribution incorporated within the Work constitutes direct
+      or contributory patent infringement, then any patent licenses
+      granted to You under this License for that Work shall terminate
+      as of the date such litigation is filed.
+
+   4. Redistribution. You may reproduce and distribute copies of the
+      Work or Derivative Works thereof in any medium, with or without
+      modifications, and in Source or Object form, provided that You
+      meet the following conditions:
+
+      (a) You must give any other recipients of the Work or
+          Derivative Works a copy of this License; and
+
+      (b) You must cause any modified files to carry prominent notices
+          stating that You changed the files; and
+
+      (c) You must retain, in the Source form of any Derivative Works
+          that You distribute, all copyright, patent, trademark, and
+          attribution notices from the Source form of the Work,
+          excluding those notices that do not pertain to any part of
+          the Derivative Works; and
+
+      (d) If the Work includes a "NOTICE" text file as part of its
+          distribution, then any Derivative Works that You distribute must
+          include a readable copy of the attribution notices contained
+          within such NOTICE file, excluding those notices that do not
+          pertain to any part of the Derivative Works, in at least one
+          of the following places: within a NOTICE text file distributed
+          as part of the Derivative Works; within the Source form or
+          documentation, if provided along with the Derivative Works; or,
+          within a display generated by the Derivative Works, if and
+          wherever such third-party notices normally appear. The contents
+          of the NOTICE file are for informational purposes only and
+          do not modify the License. You may add Your own attribution
+          notices within Derivative Works that You distribute, alongside
+          or as an addendum to the NOTICE text from the Work, provided
+          that such additional attribution notices cannot be construed
+          as modifying the License.
+
+      You may add Your own copyright statement to Your modifications and
+      may provide additional or different license terms and conditions
+      for use, reproduction, or distribution of Your modifications, or
+      for any such Derivative Works as a whole, provided Your use,
+      reproduction, and distribution of the Work otherwise complies with
+      the conditions stated in this License.
+
+   5. Submission of Contributions. Unless You explicitly state otherwise,
+      any Contribution intentionally submitted for inclusion in the Work
+      by You to the Licensor shall be under the terms and conditions of
+      this License, without any additional terms or conditions.
+      Notwithstanding the above, nothing herein shall supersede or modify
+      the terms of any separate license agreement you may have executed
+      with Licensor regarding such Contributions.
+
+   6. Trademarks. This License does not grant permission to use the trade
+      names, trademarks, service marks, or product names of the Licensor,
+      except as required for reasonable and customary use in describing the
+      origin of the Work and reproducing the content of the NOTICE file.
+
+   7. Disclaimer of Warranty. Unless required by applicable law or
+      agreed to in writing, Licensor provides the Work (and each
+      Contributor provides its Contributions) on an "AS IS" BASIS,
+      WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or
+      implied, including, without limitation, any warranties or conditions
+      of TITLE, NON-INFRINGEMENT, MERCHANTABILITY, or FITNESS FOR A
+      PARTICULAR PURPOSE. You are solely responsible for determining the
+      appropriateness of using or redistributing the Work and assume any
+      risks associated with Your exercise of permissions under this License.
+
+   8. Limitation of Liability. In no event and under no legal theory,
+      whether in tort (including negligence), contract, or otherwise,
+      unless required by applicable law (such as deliberate and grossly
+      negligent acts) or agreed to in writing, shall any Contributor be
+      liable to You for damages, including any direct, indirect, special,
+      incidental, or consequential damages of any character arising as a
+      result of this License or out of the use or inability to use the
+      Work (including but not limited to damages for loss of goodwill,
+      work stoppage, computer failure or malfunction, or any and all
+      other commercial damages or losses), even if such Contributor
+      has been advised of the possibility of such damages.
+
+   9. Accepting Warranty or Additional Liability. While redistributing
+      the Work or Derivative Works thereof, You may choose to offer,
+      and charge a fee for, acceptance of support, warranty, indemnity,
+      or other liability obligations and/or rights consistent with this
+      License. However, in accepting such obligations, You may act only
+      on Your own behalf and on Your sole responsibility, not on behalf
+      of any other Contributor, and only if You agree to indemnify,
+      defend, and hold each Contributor harmless for any liability
+      incurred by, or claims asserted against, such Contributor by reason
+      of your accepting any such warranty or additional liability.
+
+   END OF TERMS AND CONDITIONS
+
+   APPENDIX: How to apply the Apache License to your work.
+
+      To apply the Apache License to your work, attach the following
+      boilerplate notice, with the fields enclosed by brackets "{}"
+      replaced with your own identifying information. (Don't include
+      the brackets!)  The text should be enclosed in the appropriate
+      comment syntax for the file format. We also recommend that a
+      file or class name and description of purpose be included on the
+      same "printed page" as the copyright notice for easier
+      identification within third-party archives.
+
+   Copyright {yyyy} {name of copyright owner}
+
+   Licensed under the Apache License, Version 2.0 (the "License");
+   you may not use this file except in compliance with the License.
+   You may obtain a copy of the License at
+
+       http://www.apache.org/licenses/LICENSE-2.0
+
+   Unless required by applicable law or agreed to in writing, software
+   distributed under the License is distributed on an "AS IS" BASIS,
+   WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+   See the License for the specific language governing permissions and
+   limitations under the License.

README.md

+# SDPM
+SDPM is an open-source C++/python, Source-Doublet Panel Method, developed at the University of Liège by Adrien Crovato and Grigorios Dimitriadis, during the academic years 2022-2023.  
+
+## Main features
+* Cross platform (Windows and Unix) C++/python code
+
+## Documentation
+Detailed build and use instructions can be found in the [wiki](https://gitlab.uliege.be/am-dept/sdpm/wikis/home).
+
+## References
+Sanchez Martinez Mariano and Dimitriadis Grigorios, [Subsonic source and doublet panel methods](https://hdl.handle.net/2268/291174), Journal of Fluids and Structures, 2022.

init.py

+# THIS FILE IS INTENDED FOR INSTALLATION ONLY (copied as __init__.py), NOT FOR DEVELOPMENT
+from .sdpmw import *

run.py

+#!/usr/bin/env python3
+# -*- coding: utf8 -*-
+# test encoding: à-é-è-ô-ï-€
+
+# Copyright 2023 University of Liège
+# 
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+# 
+#     http://www.apache.org/licenses/LICENSE-2.0
+# 
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+# Run a script as if sdpm was installed
+
+def setup_workdir(testname, clean, verb=True):
+    """Create a single directory for the given test
+    """
+    import os, os.path
+
+    # build the name of the workspace folder
+    dir1 = os.path.abspath(os.path.dirname(__file__) + os.sep + '..') + os.sep
+    if verb: print(f'Setting workspace for "{testname}"')
+    common = os.path.commonprefix((testname, dir1))
+    resdir = testname[len(common):].replace(os.sep, "_")
+    resdir = os.path.splitext(resdir)[0] # remove ".py"
+    wdir = os.path.join('workspace', resdir)
+
+    # create the directory
+    if not os.path.isdir(wdir):
+        if verb: print('- creating', wdir)
+        os.makedirs(wdir)
+    elif os.path.isdir(wdir) and clean:
+        if verb: print('- cleaning', wdir)
+        import shutil
+        for f in os.listdir(wdir):
+            fpth = os.path.join(wdir, f)
+            if os.path.isfile(fpth):
+                os.remove(fpth)
+            elif os.path.isdir(fpth):
+                shutil.rmtree(fpth)
+    # change dir
+    if verb: print('- changing to', wdir)
+    os.chdir(wdir)
+
+def main():
+    # adds '.' and libs to the pythonpath
+    import os, sys, time, socket
+    thisdir = os.path.dirname(os.path.realpath(__file__))
+    sys.path.append(os.path.join(thisdir, 'build', 'bin'))
+    sys.path.append(os.path.join(thisdir, 'build', 'bin', 'Release'))
+    sys.path.append(os.path.join(thisdir, 'build', 'bin', 'Debug'))
+
+    # import module
+    import sdpm
+    from sdpm.utils import parse_args
+    from sdpm.streams import Tee
+    # redirect C++ streams to Python
+    redirect = sdpm.StdOutErr2Py()
+
+    # prepare run
+    args = parse_args()
+    testname = os.path.abspath(args.file)
+    if not os.path.isfile(testname):
+        raise Exception(f'File not found: {testname}')
+    setup_workdir(testname, args.clean)
+
+    # create log and run
+    tee = Tee('log') # split streams
+    print('Starting test', testname)
+    print('Time:', time.strftime('%c'))
+    print('Hostname:', socket.gethostname())
+    script = open(testname, 'r', encoding='utf-8').read()
+    exec(compile(script, testname, 'exec'), {'__file__': testname, '__name__':'__main__'}) # compile to get clean backtrace
+
+if __name__ == '__main__':
+    main()

scripts/format_code.py

+#!/usr/bin/env python3
+# -*- coding: utf-8 -*-
+#
+# This script runs "clang-format" recursively on all C/C++ files
+
+import sys
+import os
+import fnmatch
+import re
+import subprocess
+
+def all_files(root,
+              patterns='*',
+              skips='*.git*;*build*',
+              single_level=False,
+              yield_folders=False):
+    # self.checkPath(root)
+    patterns = patterns.split(';')
+    skips = skips.split(';')
+    for path, subdirs, files in os.walk(root):
+        if yield_folders:
+            files.extend(subdirs)
+        files.sort()
+        for name in files:
+            for pattern in patterns:
+                if fnmatch.fnmatch(name, pattern):
+                    fullname = os.path.join(path, name)
+                    ok = True
+                    for skip in skips:
+                        if fnmatch.fnmatch(os.path.relpath(fullname, start=root), skip):
+                            ok = False
+                    if ok:
+                        yield fullname
+                        break
+        if single_level:
+            break
+
+def main():
+    # loop over all files and format them
+    encs = {}
+    for f in all_files(os.getcwd(), patterns='*.cpp;*.c;*.h;*.hpp'):
+        # print(f)
+        cmd = ['clang-format-10', "-style=file", "-i", f]
+        retcode = subprocess.call(cmd)
+        if retcode != 0:
+            print(f'ERROR: retcode = {retcode}')
+            break
+
+
+if __name__ == "__main__":
+    # print('running format_code.py...')
+    main()

sdpm/CMakeLists.txt

+# Copyright 2023 University of Liège
+# 
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+# 
+#     http://www.apache.org/licenses/LICENSE-2.0
+# 
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+# Add source dir
+ADD_SUBDIRECTORY( src )
+ADD_SUBDIRECTORY( _src )
+
+# Add test dir
+MACRO_AddTest(${CMAKE_CURRENT_SOURCE_DIR}/tests)
+
+# Add to install
+#FILE(GLOB files "${CMAKE_CURRENT_SOURCE_DIR}/utils.py")
+#INSTALL(FILES ${files} DESTINATION ${CMAKE_INSTALL_PREFIX})
+INSTALL(FILES "${CMAKE_CURRENT_SOURCE_DIR}/utils.py" DESTINATION ${CMAKE_INSTALL_PREFIX})
+INSTALL(FILES "${CMAKE_CURRENT_SOURCE_DIR}/gmsh.py" DESTINATION ${CMAKE_INSTALL_PREFIX})

sdpm/__init__.py

+# -*- coding: utf-8 -*-
+# fwk MODULE initialization file
+
+# Copyright 2023 University of Liège
+# 
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+# 
+#     http://www.apache.org/licenses/LICENSE-2.0
+# 
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+from sdpmw import *

sdpm/_src/CMakeLists.txt

+# Copyright 2023 University of Liège
+# 
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+# 
+#     http://www.apache.org/licenses/LICENSE-2.0
+# 
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+INCLUDE(${SWIG_USE_FILE})
+
+FILE(GLOB SRCS *.h *.cpp *.inl *.swg)
+FILE(GLOB ISRCS *.i)
+
+SET_SOURCE_FILES_PROPERTIES(${ISRCS} PROPERTIES CPLUSPLUS ON)
+SET(CMAKE_SWIG_FLAGS ${CMAKE_SWIG_FLAGS} "-interface" "_sdpmw") # avoids "import _module_d" with MSVC/Debug
+SWIG_ADD_LIBRARY(sdpmw LANGUAGE python SOURCES ${ISRCS} ${SRCS})
+SET_PROPERTY(TARGET sdpmw PROPERTY SWIG_USE_TARGET_INCLUDE_DIRECTORIES ON)
+MACRO_DebugPostfix(sdpmw)
+
+TARGET_INCLUDE_DIRECTORIES(sdpmw PRIVATE ${PROJECT_SOURCE_DIR}/sdpm/_src
+                                        ${Python3_INCLUDE_DIRS})
+TARGET_LINK_LIBRARIES(sdpmw PRIVATE sdpm ${Python3_LIBRARIES})
+
+INSTALL(FILES "${EXECUTABLE_OUTPUT_PATH}/\${BUILD_TYPE}/sdpmw.py"
+        DESTINATION ${CMAKE_INSTALL_PREFIX})
+INSTALL(TARGETS sdpmw DESTINATION ${CMAKE_INSTALL_PREFIX})

sdpm/_src/sdpmCppBuf2Py.cpp

+/*
+ * Copyright 2020 University of Liège
+ *
+ * Licensed under the Apache License, Version 2.0 (the "License");
+ * you may not use this file except in compliance with the License.
+ * You may obtain a copy of the License at
+ *
+ *     http://www.apache.org/licenses/LICENSE-2.0
+ *
+ * Unless required by applicable law or agreed to in writing, software
+ * distributed under the License is distributed on an "AS IS" BASIS,
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ * See the License for the specific language governing permissions and
+ * limitations under the License.
+ */
+
+#include <Python.h>
+#include "sdpmCppBuf2Py.h"
+#include <stdio.h>
+#include <algorithm>
+
+CppBuf2Py::CppBuf2Py(std::ostream &o1, bool err) : stream(o1), errstream(err)
+{
+    oldbuf = stream.rdbuf();
+    stream.rdbuf(this);
+    verb = false;
+}
+
+CppBuf2Py::~CppBuf2Py()
+{
+    stream.rdbuf(oldbuf); // restore the ostream
+}
+
+std::streamsize
+CppBuf2Py::xsputn(const char *s, std::streamsize n)
+{
+    if (verb)
+        printf("CppBuf2Py::xsputn(%ld)\n", static_cast<long>(n));
+    // PySys_WriteStdout truncates to 1000 chars
+    static const std::streamsize MAXSIZE = 1000;
+    std::streamsize written = std::min(n, MAXSIZE);
+    if (verb)
+        printf("written=%ld\n", static_cast<long>(written));
+    std::string str(s, n);
+    if (verb)
+        printf("string=%s\n", str.c_str());
+
+    PyGILState_STATE gstate = PyGILState_Ensure();
+    if (errstream)
+        PySys_WriteStderr("%s", str.c_str());
+    else
+        PySys_WriteStdout("%s", str.c_str());
+    PyGILState_Release(gstate);
+
+    return written;
+}
+
+int CppBuf2Py::overflow(int c)
+{
+    if (verb)
+        printf("CppBuf2Py::overflow(%d)\n", c);
+
+    if (c != EOF)
+    {
+        PyGILState_STATE gstate = PyGILState_Ensure();
+        if (errstream)
+            PySys_WriteStderr("%c", c);
+        else
+            PySys_WriteStdout("%c", c);
+        PyGILState_Release(gstate);
+    }
+    return c;
+}
+
+void CppBuf2Py::test()
+{
+    CppBuf2Py buf(std::cout);
+    for (int i = 0; i < 1000; ++i)
+        std::cout << "(" << i << ")"
+                  << " this is a test - ";
+}
+
+// =============================================================================
+
+StdOutErr2Py::StdOutErr2Py() : out(std::cout), err(std::cerr, true)
+{
+}
+
+void StdOutErr2Py::test()
+{
+    StdOutErr2Py redirector;
+
+    for (int i = 0; i < 50; ++i)
+    {
+        std::cout << "(" << i << ")"
+                  << " this is a test - ";
+        std::cerr << "(" << i << ")"
+                  << " error - ";
+    }
+    std::cout << '\n';
+    std::cerr << '\n';
+}

sdpm/_src/sdpmCppBuf2Py.h

+/*
+ * Copyright 2020 University of Liège
+ *
+ * Licensed under the Apache License, Version 2.0 (the "License");
+ * you may not use this file except in compliance with the License.
+ * You may obtain a copy of the License at
+ *
+ *     http://www.apache.org/licenses/LICENSE-2.0
+ *
+ * Unless required by applicable law or agreed to in writing, software
+ * distributed under the License is distributed on an "AS IS" BASIS,
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ * See the License for the specific language governing permissions and
+ * limitations under the License.
+ */
+
+#ifndef SDPMCPPBUF2PY_H
+#define SDPMCPPBUF2PY_H
+
+#include "sdpm.h"
+#include <iostream>
+
+// http://bo-peng.blogspot.be/2004/10/how-to-re-direct-cout-to-python_05.html
+// https://docs.python.org/2/c-api/sys.html#system-functions
+
+/**
+ * @brief redirect a std::stream to python
+ */
+#ifndef SWIG
+class CppBuf2Py : public std::basic_streambuf<char>
+{
+public:
+    typedef std::char_traits<char> traits_type;
+    typedef traits_type::int_type int_type;
+
+private:
+    std::ostream &stream;
+    std::streambuf *oldbuf;
+    bool verb;
+    bool errstream;
+
+public:
+    CppBuf2Py(std::ostream &o1, bool err = false);
+    ~CppBuf2Py();
+
+    static void test();
+
+protected:
+    virtual std::streamsize xsputn(const char *s, std::streamsize n);
+    virtual int_type overflow(int c);
+};
+#endif // SWIG
+
+/**
+ * @brief redirect std::cout/std::cerr to python sys.stdout/sys.stderr
+ */
+
+class StdOutErr2Py
+{
+    CppBuf2Py out;
+    CppBuf2Py err;
+
+public:
+    StdOutErr2Py();
+
+    static void test();
+};
+
+#endif // SDPMCPPBUF2PY_H

sdpm/_src/sdpmw.i

+/*
+ * Copyright 2023 University of Liège
+ *
+ * Licensed under the Apache License, Version 2.0 (the "License");
+ * you may not use this file except in compliance with the License.
+ * You may obtain a copy of the License at
+ *
+ *     http://www.apache.org/licenses/LICENSE-2.0
+ *
+ * Unless required by applicable law or agreed to in writing, software
+ * distributed under the License is distributed on an "AS IS" BASIS,
+ * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+ * See the License for the specific language governing permissions and
+ * limitations under the License.
+ */
+
+// fwk.i: SWIG input file for python interface
+
+%feature("autodoc","1");
+
+%module(docstring=
+"'sdpmw' module: Source-Doublet Panel Method
+(c) ULg - A&M",
+threads="1"
+) sdpmw
+%{
+#include "sdpm.h"
+#include "sdpmObject.h"
+#include "sdpmTimers.h"
+#include "sdpmMesh.h"
+#include "sdpmNode.h"
+#include "sdpmElement.h"
+#include "sdpmLine2.h"
+#include "sdpmQuad4.h"
+#include "sdpmTag.h"
+#include "sdpmGmshImport.h"
+#include "sdpmGmshExport.h"
+#include "sdpmUnitTest.h"
+
+#include "sdpmCppBuf2Py.h"
+
+#include <sstream>
+%}
+
+// COMMON
+%include "std_string.i"
+%include "std_vector.i"
+%include "std_map.i"
+
+// EXCEPTION
+// from: http://swig.10945.n7.nabble.com/Trapping-Swig-DirectorException-td6013.html
+// le code suivant permet de voir la call stack dans les appels C++ => python
+%include "exception.i"
+%{ 
+   static void handle_exception(void) { 
+     try { 
+       throw; 
+     } catch (std::exception &e) { 
+        std::stringstream txt; 
+        txt << e.what(); // << ' ' << typeid(e).name();
+        PyErr_SetString(PyExc_RuntimeError, e.what()); 
+     } 
+     catch(...) 
+     {
+        PyErr_SetString(PyExc_RuntimeError, "Unknown C++ Runtime Error");
+     } 
+   } 
+%} 
+%exception { 
+   try { 
+     $action 
+   } catch (...) { 
+     // Note that if a director method failed, the Python error indicator 
+     // already contains full details of the exception, and it will be 
+     // reraised when we go to SWIG_fail; so no need to convert the C++ 
+     // exception back to a Python one 
+     if (!PyErr_Occurred()) { 
+       handle_exception(); 
+     } 
+     SWIG_fail; 
+   } 
+}
+
+// SDPM
+%include "sdpm.h"
+%include "sdpmObject.h"
+%include "sdpmTimers.h"
+%include "sdpmMesh.h"
+%include "sdpmNode.h"
+%nodefault sdpm::Element; // Element is pure virtual
+%include "sdpmElement.h"
+%include "sdpmLine2.h"
+%include "sdpmQuad4.h"
+%include "sdpmTag.h"
+%include "sdpmGmshImport.h"
+%include "sdpmGmshExport.h"
+%include "sdpmUnitTest.h"
+
+%warnfilter(401); //Nothing known about base class 'std::basic_streambuf< char >'
+%include "sdpmCppBuf2Py.h"
+%warnfilter(+401);
+
+// SDPM template
+%template(std_vector_Nodep) std::vector<sdpm::Node *>;
+%template(std_vector_Elementp) std::vector<sdpm::Element *>;
+%template(std_vector_Tagp) std::vector<sdpm::Tag *>;
+%template(std_map_string_Tagp) std::map<std::string, sdpm::Tag *>;
+
+namespace sdpm {
+
+// So that print(obj) use << *obj
+%extend sdpmObject {
+    std::string __str__() {
+        std::stringstream str;
+        str << *self;
+        return str.str();
+    }
+}
+
+// To get a Timer from Timers[]
+%extend Timers {
+    Timer &getitem(std::string const &name)
+    {
+        return (*self)[name];
+    }
+
+    %pythoncode {
+def __getitem__(self, name):
+    return self.getitem(name)
+    }
+}
+};