Annotation raytracing as overlay

This was the initial attempt at non-scene-geometry annotation raytracing.
Squashed down to one commit for simpler diff review.

Clear Werror=subobject-linkage issue
use bu_strcmp
Make annotation intersections scale independent
Use relative parallel and determinant tests so valid annotations are not rejected after small model transforms.  Carry the caller's model tolerance into visibility decisions and cover both behaviors in the focused annotation semantics test.
Use a strict annotation hit ordering
Compare hit distances exactly before using the path as a tie breaker.  Approximate equality in a std::sort predicate is not transitive and can violate the strict weak ordering required by the C++ standard library.
Preserve hit distance in simple rt lightmodels
Record the first partition distance in viewit, matching colorview.  Annotation compositing uses this value to keep model-space marks behind foreground geometry for lightmodels 1, 2, 4, 5, and 6.
Defer NIRT annotation reporting
Remove the provisional useannot mode instead of emitting records outside NIRT's fmt system.  This preserves existing structured output contracts and avoids extending the public nirt_opt_vals layout until annotation records have a complete formatter design.
Use annotation line styles for render width
Remove the redundant object-wide render_width field, ANR2 payload, and --stroke-width option.  Ray renderers now rely on the established per-segment line_width style over an automatic base width, preserving the public primitive layout and avoiding another database extension.  Document and test the ray coverage behavior.
Name the ADRT model coordinate scale
Replace repeated 1000-to-1 load conversions with ADRT_MODEL_UNITS_PER_TIE_UNIT so geometry and annotation compositing share one documented coordinate convention.
Keep annotations out of the librender camera ABI
Replace camera-resident annotation fields with an additive render_camera_render_annotations entry point.  Per-render private state now carries the scene, coordinate scale, and derived view while the public render_camera_t and thread-data layouts remain unchanged.  Preserve primary rays before shaders mutate them so annotation compositing evaluates the intended sample.
Define ISST annotation scene ownership
Make GFile release its annotation scene and TIE data, destroy the render window before its GFile owner, and clear the renderer under its mutex before reloading.  Queued renders now reject cleared scenes after acquiring the context, preventing use-after-free during reload and shutdown.
Cache rt annotation preparation state
Keep annotation scene and view state private to do.c, expose read-only accessors to workers, and reuse prepared coverage when the database, selected paths, and tolerances are unchanged.  Explicitly invalidate on animation changes and centralize cleanup.
Add a NIRT annotation output format
Add opt-in NIRT annotation reporting
Document NIRT annotation diagnostics
Document annotation and view-object roles
Remove the unused ISST application meta-object
ISSTApp declares no signals, slots, properties, or invokable methods.  Avoid generating a redundant derived meta-object, which also prevents Qt 6 from rejecting the member main-window construction during application startup.
Use Qt 6 OpenGL dependencies for ISST
Qt6::OpenGL and Qt6::OpenGLWidgets already carry the platform OpenGL implementation.  Do not also link the independently discovered host OpenGL target, which can select a different GL ABI and bypass Qt's configured backend.
Keep ISST annotations out of NMG tessellation
ISST renders annotations through its dedicated annotation scene.  Exclude annotation-only roots from the geometry walk and replace annotation leaves in mixed trees with empty nodes so the independent NMG path does not report unsupported tessellation.
Share ADRT annotation tessellation filtering
Frame ISST views using annotation bounds
Gate ART annotation frame processing
Honor inherited colors in annotation renders
Make RTree ray queries const-safe
Composite ordered annotation coverage with RTree indices
Blend all annotation layers in image renderers
Test ordered annotation coverage and alpha blending
Test ADRT annotation filtering and scene framing
Document annotation layering and lookup design
Tighten annotation scene preparation
Test annotation combination color inheritance

.github/workflows/push.yml

@@ -384,7 +384,7 @@ jobs:
   osx:
     name: macOS Clang
     needs: [linux]
-    runs-on: macos-15
+    runs-on: macos-14
     strategy:
       fail-fast: true
     steps:
@@ -403,7 +403,7 @@ jobs:
       - name: Prepare bext dependencies
         uses: ./brlcad/.github/actions/prepare-bext
         with:
-          platform-channel: macos-15-${{ runner.arch }}-release
+          platform-channel: macos-14-${{ runner.arch }}-release
           source-dir: ${{ github.workspace }}/brlcad
           output-dir: ${{ runner.temp }}/bext_output
           build-type: Release

.github/workflows/scheduled.yml

@@ -177,7 +177,7 @@ jobs:
   osx_bext:
     name: macOS - BRL-CAD bext build
     needs: [linux_bext]
-    runs-on: macos-15
+    runs-on: macos-14
     strategy:
       fail-fast: true
     if: ${{ !cancelled() }}

CMakeLists.txt

@@ -231,7 +231,7 @@ set(CMAKE_FIND_USE_PACKAGE_REGISTRY FALSE)
 # Prepare to manage and build BRL-CAD's external dependencies. Specifies a
 # specific commit hash to use in the 'bext' repository.
 include(BRLCAD_ExternalDeps)
-brlcad_bext_init(42d8d7ca57f9d2f61192b41f0b41a731f1ba6c47)
+brlcad_bext_init(12115e6c5922e2db27a9b88e00dd714af5315f69)
 
 #---------------------------------------------------------------------
 # Define various utilities.
@@ -352,10 +352,9 @@ set(CMAKE_SKIP_BUILD_RPATH FALSE)
 # it should only be set to the installation value when actually installed.
 set(CMAKE_BUILD_WITH_INSTALL_RPATH FALSE)
 
-# Installed runtime paths are defined explicitly below.  Do not append
-# directories discovered on the build host, since those paths may not exist
-# or may resolve to unrelated libraries on the installed system.
-set(CMAKE_INSTALL_RPATH_USE_LINK_PATH FALSE)
+# Add the automatically determined parts of the RPATH which point to
+# directories outside the build tree to the install RPATH
+set(CMAKE_INSTALL_RPATH_USE_LINK_PATH TRUE)
 
 # The following logic is what allows binaries to run successfully in the build
 # directory AND install directory.
@@ -679,13 +678,8 @@ endif(MSVC)
 # before use so platforms without GNU/LLVM-style linker options are skipped.
 brlcad_check_linker_flag(C "LINKER:--no-undefined" BRLCAD_LINKER_NO_UNDEFINED_C_FLAG_FOUND)
 brlcad_check_linker_flag(CXX "LINKER:--no-undefined" BRLCAD_LINKER_NO_UNDEFINED_CXX_FLAG_FOUND)
-if(BRLCAD_LINKER_NO_UNDEFINED_C_FLAG_FOUND AND BRLCAD_LINKER_NO_UNDEFINED_CXX_FLAG_FOUND
-   AND NOT BRLCAD_ENABLE_ADDRESS_SANITIZER)
+if(BRLCAD_LINKER_NO_UNDEFINED_C_FLAG_FOUND AND BRLCAD_LINKER_NO_UNDEFINED_CXX_FLAG_FOUND)
   set(BRLCAD_LINKER_NO_UNDEFINED_FLAG "LINKER:--no-undefined" CACHE STRING "Linker flag requiring resolved shared library symbols" FORCE)
-else()
-  # Clang links the AddressSanitizer runtime into executables, leaving its
-  # shared-library instrumentation symbols unresolved until then.
-  set(BRLCAD_LINKER_NO_UNDEFINED_FLAG "" CACHE STRING "Linker flag requiring resolved shared library symbols" FORCE)
 endif()
 brlcad_check_linker_flags(C BRLCAD_LINKER_WHOLE_ARCHIVE_GROUP_C_FLAG_FOUND "LINKER:--whole-archive" "LINKER:--no-whole-archive")
 brlcad_check_linker_flags(CXX BRLCAD_LINKER_WHOLE_ARCHIVE_GROUP_CXX_FLAG_FOUND "LINKER:--whole-archive" "LINKER:--no-whole-archive")
@@ -1648,31 +1642,21 @@ if(NOT MSVC)
   endif(BRLCAD_WARNINGS OR BRLCAD_ENABLE_STRICT)
 
   if(BRLCAD_ENABLE_COVERAGE)
-    if(CMAKE_C_COMPILER_ID MATCHES "Clang" AND CMAKE_CXX_COMPILER_ID MATCHES "Clang")
-      check_c_flag(fprofile-instr-generate)
-      check_cxx_flag(fprofile-instr-generate)
-      check_c_flag(fcoverage-mapping)
-      check_cxx_flag(fcoverage-mapping)
-      if(NOT FPROFILE_INSTR_GENERATE_C_FLAG_FOUND OR
-         NOT FPROFILE_INSTR_GENERATE_CXX_FLAG_FOUND OR
-         NOT FCOVERAGE_MAPPING_C_FLAG_FOUND OR
-         NOT FCOVERAGE_MAPPING_CXX_FLAG_FOUND)
-        message(FATAL_ERROR "Clang source-based coverage is not supported by this toolchain.")
-      endif()
-      add_compile_options(-fprofile-instr-generate -fcoverage-mapping)
-      add_link_options(-fprofile-instr-generate)
-    else()
-      check_c_flag(ftest-coverage)
-      check_cxx_flag(ftest-coverage)
-      if(FTEST_COVERAGE_C_FLAG_FOUND AND FTEST_COVERAGE_CXX_FLAG_FOUND)
-        set(CMAKE_C_FLAGS "${CMAKE_C_FLAGS} -ftest-coverage -fprofile-arcs -fprofile-update=atomic")
-        set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -ftest-coverage -fprofile-arcs -fprofile-update=atomic")
-        set(CMAKE_SHARED_LINKER_FLAGS "${CMAKE_SHARED_LINKER_FLAGS} -lgcov -ftest-coverage -fprofile-arcs -fprofile-update=atomic")
-        set(CMAKE_EXE_LINKER_FLAGS "${CMAKE_EXE_LINKER_FLAGS} -lgcov -ftest-coverage -fprofile-arcs -fprofile-update=atomic")
-      else()
-        message(FATAL_ERROR "GCC coverage is not supported by this toolchain.")
-      endif()
-    endif()
+    # TODO: These seem to GCC specific flags - should probably set up for clang
+    # as well if that's our compiler...
+    # https://clang.llvm.org/docs/SourceBasedCodeCoverage.html
+    # If we need different settings for different tools, should
+    # encapsulate each tool's setup in a function
+    check_c_flag(ftest-coverage)
+    check_cxx_flag(ftest-coverage)
+    if(FTEST_COVERAGE_C_FLAG_FOUND)
+      set(CMAKE_C_FLAGS "${CMAKE_C_FLAGS} -ftest-coverage -fprofile-arcs -fprofile-update=atomic")
+      set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -ftest-coverage -fprofile-arcs -fprofile-update=atomic")
+      set(CMAKE_SHARED_LINKER_FLAGS "${CMAKE_SHARED_LINKER_FLAGS} -lgcov -ftest-coverage -fprofile-arcs -fprofile-update=atomic")
+      set(CMAKE_EXE_LINKER_FLAGS "${CMAKE_EXE_LINKER_FLAGS} -lgcov -ftest-coverage -fprofile-arcs -fprofile-update=atomic")
+    else(FTEST_COVERAGE_C_FLAG_FOUND)
+      message(SEND_ERROR "Building with coverage is not supported by BRL-CAD on this platform.")
+    endif(FTEST_COVERAGE_C_FLAG_FOUND)
   endif(BRLCAD_ENABLE_COVERAGE)
 
   if(BRLCAD_ENABLE_STRICT)

INSTALL

@@ -151,13 +151,10 @@ recommended for serious developers.
 
 Various useful configurations for working with bext:
 
-* Acquire bext and its dependencies from another repository host
+* Get a git clone that has everything which might be needed (large,
+  but useful for offline building)
 
-cmake .. -DBRLCAD_REPOSITORY_ROOT=https://git.example.org/brlcad
-
-* Acquire bext and its dependencies from local bare Git mirrors
-
-cmake .. "-DBRLCAD_REPOSITORY_ROOT=/path/to/mirror/BRL-CAD"
+git clone --recurse-submodules https://github.com/BRL-CAD/brlcad
 
 * Have prebuilt bext outputs to reuse - BRL-CAD will simply copy and prepare
   the pre-existing files
@@ -401,40 +398,10 @@ Aliases:  ENABLE_RUNTIME_DEBUG, ENABLE_RUN_TIME_DEBUG, ENABLE_RUNTIME_DEBUGGING
           ENABLE_RUN_TIME_DEBUGGING
 
 
---- BRLCAD_REPOSITORY_ROOT ---
-
-Repository URL or local directory containing bext and its dependency
-repositories. Defaults to https://github.com/BRL-CAD. The root directly
-contains repositories such as bext.git and astyle.git; BRL-CAD does not
-append BRL-CAD to the root.
-
-Remote Git URLs (including SSH), native paths, and file:// URLs are supported.
-Relative native paths are resolved against the BRL-CAD source directory.
-Quote the entire -D argument for paths containing spaces. The setting is
-cached and is also passed as BEXT_REPOSITORY_ROOT when BRL-CAD configures
-bext, including when BRLCAD_EXT_SOURCE_DIR supplies an existing checkout.
-Alternate roots require bext sources with BEXT_REPOSITORY_ROOT support.
-
-Local roots restrict CMake-driven Git acquisition and queries to file
-transport, with no online fallback. Mirrors must contain all required pinned
-commits and recursive submodules. Dependencies outside bext's canonical URL
-prefix require their own Git URL mappings. This setting does not restrict
-network access by dependency build scripts or change persistent Git settings.
-BRLCAD_EXT_DIR still selects prebuilt outputs, and changing the repository
-root alone does not force those outputs to be rebuilt.
-
-To return to the default source:
-
-cmake .. -DBRLCAD_REPOSITORY_ROOT=https://github.com/BRL-CAD
-
-
 --- BRLCAD_EXT_DIR ---
 
 If specified, BRL-CAD will look in this directory for bext compiled outputs to
-bundle and use during the build. Relative paths are resolved against the
-BRL-CAD source directory. Distcheck subconfigures reuse this resolved output
-directory and preserve BRLCAD_REPOSITORY_ROOT, BRLCAD_EXT_SOURCE_DIR, and
-BRLCAD_EXT_PARALLEL when specified.
+bundle and use during the build.
 
 
 --- BRLCAD_EXT_SOURCE_DIR ---
@@ -442,8 +409,7 @@ BRLCAD_EXT_PARALLEL when specified.
 If specified, BRL-CAD will look in this directory for bext source files. Use
 this when you have a pre-existing bext clone and don't want BRL-CAD's configure
 process to pull down another copy, but are still looking to have the BRL-CAD
-configure process manage producing bext build outputs. Relative paths are
-resolved against the BRL-CAD source directory.
+configure process manage producing bext build outputs.
 
 
 --- BRLCAD_EXT_PARALLEL ---

NEWS

@@ -10,112 +10,78 @@ descriptive text that emphasizes or further describes significant
 changes made.  See document footer for additional details.
 
 ----------------------------------------------------------------------
---- 202X-XX-XX  Release 7.48.0                                     ---
+--- 20XX-XX-XX  Release 7.46.0                                     ---
 ----------------------------------------------------------------------
 
-
-----------------------------------------------------------------------
---- 2026-10-12  Release 7.46.0                                     ---
-----------------------------------------------------------------------
-
-This release expands geometry conversion and improves modeling
-reliability.  STEP support adds AP203e2, AP214, and AP242
-import/export, native CSG round trips, and AP242 product and
-manufacturing information (PMI) import.  IGES and JT conversion is
-expanded, and new OpenSCAD CSG import/export tools are available.
-
-Facetization is more reliable, NMG solid ray tracing is 10-20 times
-faster, CSG ray tracing is faster, and BoT, DSP, and point-cloud ray
-intersections are corrected.  MGED gains numerous fixes to primitive,
-sketch, and combination editing, drawing, dialogs, and preferences.
-
-New tools support object nesting with arrange, building generation
-with mkbuilding, annotation editing, and primitive repair.  Geometry
-analysis, image utilities, Windows support, documentation, and bundled
-dependencies also receive substantial updates.
-
-* corrected Archer logic for saving window size - Cliff Yapp
-* fixed inside BoT rays, DSP surface normals, PNTS bounds - Cliff Yapp
-* improved mged fallback HTML manual viewer - Cliff Yapp
-* fixed gqa undefined densities and exit status - Chris McGregor
-* added ARBN support to the mged inside command - Cliff Yapp
-* improved HYP neck-ratio labels and BoT face selection - Cliff Yapp
-* fixed mged view panning, edit rotation, & ray queries - Cliff Yapp
-* preserved ARB, CLINE, TGC, & extrusion geometry edits - Cliff Yapp
-* restored mged command prompt focus at GUI startup - Cliff Yapp
-* improved pixscale & bwscale scaling and image piping - Sean Morrison
-* restored pixdiff behavior on .pix, added -b for .bw - Sean Morrison
-* improved facetize fallback methods and reliability - Cliff Yapp
-* updated mged and Archer for Tcl 8.6.18 and Itcl/Itk 4 - Cliff Yapp
+* updated MGED and Archer for Tcl 8.6.18 and Itcl/Itk 4 - Cliff Yapp
 * fixed writes to fileless geometry databases - Cliff Yapp
 * fixed nastran-g handling of large property IDs - Rupesh, Cliff Yapp
-* fixed dxf-g invalid polyface idx crash - Meghna45-code, Cliff Yapp
+* fixed dxf-g invalid polyface index crashes - Meghna45-code, Cliff Yapp
+* fixed Archer primitive editors under Itcl 4 - Cliff Yapp
+* fixed Archer shutdown menu cleanup errors - Cliff Yapp
+* fixed MGED Build Pattern dialog under Itcl 4 - Cliff Yapp
 * fixed TclCAD framebuffer listing and deletion - Cliff Yapp
-* fixed mged Build Pattern progress updates - Cliff Yapp
-* fixed mged search -exec reentry and Ctrl-C crashes - Cliff Yapp
-* fixed mged asynchronous commands and clone progress - Cliff Yapp
-* accelerated large and repeated mged path drawing - Cliff Yapp
+* fixed MGED Build Pattern progress updates - Cliff Yapp
+* fixed MGED search -exec reentry and Ctrl-C crashes - Cliff Yapp
+* fixed MGED asynchronous commands and clone progress - Cliff Yapp
+* accelerated large and repeated MGED path drawing - Cliff Yapp
 * fixed ARB7 inside solid construction - Cliff Yapp
 * fixed direct ARB face and edge-direction editing - Cliff Yapp
 * fixed rtwizard startup with framebuffer contention - Cliff Yapp
 * fixed Archer annotation line display updates - Cliff Yapp
 * fixed Archer sketch scale drift in non-mm databases - Cliff Yapp
 * fixed toggle arrow callbacks during setup - Cliff Yapp
-* fixed Tcl dialogs with multi-word titles and msgs - Cliff Yapp
+* fixed Tcl dialogs with multi-word titles and messages - Cliff Yapp
 * fixed Archer Tank Wizard gun generation - Cliff Yapp
-* fixed mged and Archer sketch arc tangency editing - Cliff Yapp
+* fixed MGED and Archer sketch arc tangency editing - Cliff Yapp
 * fixed Archer ARB4 point 3 editing - Cliff Yapp
 * fixed ARB4 storage from make and mk_arb4 - Cliff Yapp
 * fixed ARB4 and ARB6 center point calculations - Cliff Yapp
-* fixed concurrent mged and rtwizard framebuffer use - Cliff Yapp
-* fixed mged ShotVis editing and reload failures - Cliff Yapp
+* fixed concurrent MGED and rtwizard framebuffer use - Cliff Yapp
+* fixed MGED ShotVis editing and reload failures - Cliff Yapp
 * fixed Tcl widget geometry and selection defects - Cliff Yapp
-* fixed mged BoT editor property handling - Cliff Yapp
-* fixed mged LOD dialog and adaptive BoT wireframes - Cliff Yapp
-* fixed mged Plot, PostScript, and manual launch tools - Cliff Yapp
-* fixed mged Overlap Tool multi-object result handling - Cliff Yapp
+* fixed MGED BoT editor property handling - Cliff Yapp
+* fixed MGED LOD dialog and adaptive BoT wireframes - Cliff Yapp
+* fixed MGED Plot, PostScript, and manual launch tools - Cliff Yapp
+* fixed MGED Overlap Tool multi-object result handling - Cliff Yapp
 * eliminated 16-fold redundant gchecker grid analysis - Cliff Yapp
-* fixed submodel transforms & relative database paths - Sean Morrison
-* fixed Archer tree handling for names with spaces - Cliff Yapp
+* fixed Archer tree handling for object names with spaces - Cliff Yapp
 * fixed RtWizard multi-object edge occlusion handling - Cliff Yapp
 * fixed RtWizard title, About, and preference handling - Cliff Yapp
-* fixed mged startup file false legacy-syntax warnings - Cliff Yapp
-* fixed mged saved-view deletion and collaboration sync - Cliff Yapp
-* fixed mged font preferences apply and persistence - Cliff Yapp
+* fixed MGED startup file false legacy-syntax warnings - Cliff Yapp
+* fixed MGED saved-view deletion and collaboration sync - Cliff Yapp
+* fixed MGED font preferences apply and persistence - Cliff Yapp
 * fixed Archer combination edits changing Boolean trees - Cliff Yapp
 * fixed Archer command-line help startup crash - Cliff Yapp
-* fixed mged Geometry Browser menu and close behavior - Cliff Yapp
-* fixed mged collaboration view unshare crash - Cliff Yapp
-* added override for Archer's preference file location - Cliff Yapp
+* fixed MGED Geometry Browser menu and close behavior - Cliff Yapp
+* fixed MGED collaboration view unshare crash - Cliff Yapp
+* added an override for Archer's preference file location - Cliff Yapp
 * fixed embedded Archer display and mode menus - Cliff Yapp
 * added Archer -h and --help command-line options - Cliff Yapp
 * fixed Archer Help opening the documentation menu - Sean Morrison
-* added arrange nest command for packing db objects - Sean Morrison
-* added mkbuilding for parametric/OSM/CityJSON models - Sean Morrison
+* added arrange nest command for packing database objects - Sean Morrison
+* new mkbuilding tool for parametric/OSM/CityJSON models - Sean Morrison
 * cross-platform BINUNIF and DSP byte-order handling - Sean Morrison
-* improved several cv default behaviors (s,32,64,n) - Sean Morrison
+* improved several cv default behaviors (s, 32, 64, n) - Sean Morrison
 * fixed invalid bot handling bug applying matrices - Sean Morrison
 * improved glTF import robustness via gcv - Sean Morrison
 * fixed bot dump glTF/GLB export; -F is now rejected - Cliff Yapp
 * faster NMG solid ray tracing (10-20x faster) - Sean Morrison
-* preserved prim colors in shaded BoT and B-Rep drawing - Cliff Yapp
-* preserved B-Rep attrs through in-place geometry edits - Cliff Yapp
-* fixed gcv filter arguments w/ one-character values - Cliff Yapp
+* preserved primitive colors in shaded BoT and B-Rep drawing - Cliff Yapp
+* preserved B-Rep attributes through in-place geometry edits - Cliff Yapp
+* fixed gcv filter arguments containing one-character values - Cliff Yapp
 * improved iges-g importer robustness - Cliff Yapp
 * expanded g-iges native CSG and B-Rep export support - Cliff Yapp
 * added IGES import and export support to gcv - Cliff Yapp
-* rt -b query option now supports a single value idx - Sean Morrison
+* rt -b query option now supports a single value index - Sean Morrison
 * fixed .dpix file format portability and icv support - Sean Morrison
 * expanded pixscale support for piped image data - Sean Morrison
 * fixed HYP volume and ARBN totals reported by analyze - Cliff Yapp
-* fixed VOL data-file lookup and missing-data errors - Chris McGregor
-* corrected rtarea sampling/exposed-area calculations - Cliff Yapp
+* corrected rtarea sampling and exposed-area calculations - Cliff Yapp
 * bot split now groups results and preserves mesh data - Cliff Yapp
-* bot decimate now preserves nrmls, UVs, and plate data - Cliff Yapp
+* bot decimate now preserves normals, UVs, and plate data - Cliff Yapp
 * added raytracing support for annotation objects - Cliff Yapp
-* fixed Windows subprocess quoting and empty arguments - Cliff Yapp
-* fixed mged X embedded framebuffer display manager crash - Cliff Yapp
-* fixed Windows drive-letter path handling - Cliff Yapp
+* fixed MGED X embedded framebuffer display manager crash - Cliff Yapp
 * added annotate command for managing annotation objects - Cliff Yapp
 * added repair command for supported invalid primitives - Cliff Yapp
 * changed default new-region material code from 1 to 0 - Keith Bowman
@@ -137,34 +103,26 @@ dependencies also receive substantial updates.
 * added option validation to image & framebuffer tools - Sean Morrison
 * added contributor, API, command, & raytracing guides - Sean Morrison
 * added and updated command and library manual pages - Sean Morrison
-* added pnts write -f xyz export option - Chris McGregor
-* fixed mged ARB5 point editing - Cliff Yapp
-* fixed mged gqa and rtcheck stale overlap plots - Cliff Yapp
-* fixed procedural tree generator crashes - Sean Morrison, Cliff Yapp
-* added gaia Earth models with day/dusk/night textures - Sean Morrison
-* fixed fbcolor ramps & Windows input - Cliff Yapp, Sean Morrison
 * added parashot gravity trajectory path generator - Sean Morrison
-* added metaball pnt scaling (via metaball_move_pnt) - Sean Morrison
+* added metaball point scaling (via metaball_move_pnt) - Sean Morrison
 * all raytracing apps now report object prep failures - Sean Morrison
 * improved fast4-g robustness and error reporting - Cliff Yapp
 * improved robustness of g-vrml plate mode handling - Cliff Yapp
-* added support for tick drawing on mged model axis - Cliff Yapp
+* fixed MGED ARB5 point, pipe, and sketch editing - Cliff Yapp
+* added pnts write -f xyz export option - Chris McGregor
+* added support for tick drawing on MGED model axis - Cliff Yapp
 * added naca456 left and paired wing generation - Cliff Yapp
 * added patch-g -j option, specify object thicknesses - Sean Morrison
-* corrected bwstat median and input error handling - Sean Morrison
+* fixed fbcolor color ramp bug drawing the last color - Sean Morrison
 * fixed gist -Q quick preview option - Sean Morrison
 * rt preserves eye_pt setting in command scripts - Sean Morrison
-* removed trailing whitespace from nirt output - Sean Morrison
 * rt -l8 heat graphs now measure CPU time - Sean Morrison
 * documented gqa centroid & moments of inertia options - Sean Morrison
 * pixdiff now calculates pixel differences on bw files - Sean Morrison
-* fixed bwcrop single-axis crops & pixembed edge fill - Sean Morrison
+* fixed bwcrop single-axis crops and pixembed edge fill - Sean Morrison
 * fixed isize in non-mm databases - Josh Baker
 * fixed metaball adjust clearing points on errors - Sean Morrison
 * fixed benchmark+pixcmp correctness validation bug - Sean Morrison
-* fixed mged gqa option and value parsing - Sean Morrison
-* improved gist hierarchy overflow rendering - Chris McGregor
-* fixed rtxray -o grayscale/BW output - Sean Morrison
 * fbhelp writes help to stdout and rest to stderr - Sean Morrison
 * FB now supports /dev/stdout and /dev/stderr devices - Sean Morrison
 * fixed mged rt -C + -o background color to file bug - Sean Morrison
@@ -173,7 +131,7 @@ dependencies also receive substantial updates.
 * added rtsil -o support for bw output - Sean Morrison
 * added cutting plane animation generation to rtwizard - Cliff Yapp
 * fixed pixdiff bug on Windows - Sean Morrison
-* fixed mged startup window placement on Windows - Chris McGregor
+* fixed MGED startup window placement on Windows - Chris McGregor
 * fixed pixborder utility right-side edge bug - Sean Morrison
 * fixed 'wavelet' image processing tool on Windows - Sean Morrison
 * quieted mged killtree spurious lookup failure output - Sean Morrison
@@ -184,16 +142,14 @@ dependencies also receive substantial updates.
 * new doc guide, a visual overview of BRL-CAD objects - Sean Morrison
 * improved triangle quality produced by bot remesh - Cliff Yapp
 * updated numerous manpage command references - Josh Baker
+* fixed submodel paths relative to parent databases - Sean Morrison
 * improved performance making huge combinations - Sean Morrison
-* added PNTS get, adjust, form, & g2asc/asc2g support - Chris McGregor
-* fixed ANNOT get, adjust, & g2asc/asc2g support - Chris McGregor
 * added make -t primitive type listing - Chris McGregor
 * added gcv import conversion support for vrml/wrl 1.0 - Sean Morrison
 * fixed gcv appending to existing geometry databases - Sean Morrison
 * added CSS names, hex, HSL/HSV to libbu color parsing - Sean Morrison
 * improved OpenNURBS B-Rep import and adjustment safety - Cliff Yapp
-* improved Debian application menus/MIME associations - Jordi Sayol
-* updated Tcl/Tk dependency to v8.6.18 using bundled Itcl - Cliff Yapp
+* improved Debian application menus and MIME associations - Jordi Sayol
 * updated Clipper2 dependency to v1.5.4 - Cliff Yapp
 * updated Eigen dependency to v5.0.1 - Cliff Yapp
 * updated GDAL dependency to v3.13.3 - Cliff Yapp
@@ -205,6 +161,7 @@ dependencies also receive substantial updates.
 * updated OpenNURBS dependency to v8.34.26223.11001 - Cliff Yapp
 * updated PROJ dependency to v9.8.1 - Cliff Yapp
 * updated PoissonRecon dependency to v18.75 - Cliff Yapp
+* updated Tcl/Tk dependency to v8.6.18, using bundled Itcl - Cliff Yapp
 * updated libjpeg-turbo dependency to v3.2.0 - Cliff Yapp
 * updated libosmesa dependency to SHA1 bfa90de - Cliff Yapp
 * updated libpng dependency to v1.6.58 - Cliff Yapp
@@ -313,7 +270,7 @@ and extensively improving existing docs.
 * added metadata-aware image diff & nirt analysis - Cliff Yapp
 * added basic PNTS point-cloud ray tracing - Sean Morrison
 * improved nirt behavior in Windows bash consoles - Cliff Yapp
-* added g2asc/asc2g support for DSP and EBM - Sean Morrison
+* added g2asc/asc2g support for PNTS, DSP, and EBM - Sean Morrison
 * updated gct mesh decimation dependency to mmesh - Cliff Yapp
 * added procedural wing, propeller & tree generators - Cliff Yapp
 * added ISO tire dimensions & named tread presets - Cliff Yapp
@@ -585,7 +542,6 @@ for more information.
 ----------------------------------------------------------------------
 --- 2023-11-06  Release 7.38.0                                     ---
 ----------------------------------------------------------------------
-
 This major release of BRL-CAD introduces significant improvements to
 the robustness of the facetize command by leveraging the Manifold
 library for boolean evaluation of meshes. Facetization is an operation
@@ -624,7 +580,6 @@ processing high level assemblies in BRL-CAD's example geometry files:
 ----------------------------------------------------------------------
 --- 2023-07-26  Release 7.36.0                                     ---
 ----------------------------------------------------------------------
-
 This major release of BRL-CAD introduces significant updates to
 geometry conversion support. The update includes modernization for
 NURBS importing and interaction, BoT manipulation, and incorporates

README

@@ -1,6 +1,6 @@
-                          BRL-CAD
-                      Release 7.46.0
-                    http://brlcad.org/
+			    BRL-CAD
+			 Release 7.44.2
+		       http://brlcad.org/
 
 BRL-CAD is a powerful cross-platform open source combinatorial
 solid modeling system that includes an interactive 3D solid geometry

annot_viewobj.txt

@@ -0,0 +1,246 @@
+Annotation and view-object rendering notes
+==========================================
+
+Scope
+-----
+
+BRL-CAD has two related, but deliberately distinct, ways to represent visual
+markup:
+
+* An annotation is an ID_ANNOT database primitive.  It is persistent model
+  content with enough semantic and presentation information to be reproduced
+  outside the application that created it.
+* A view object is application-owned scene state.  It is useful for interactive
+  construction and manipulation, but is intentionally not backed by a database
+  object.
+
+They share portions of the interactive display infrastructure, but they do not
+currently have interchangeable ownership, lifetime, or rendering semantics.
+
+
+Annotation capabilities
+-----------------------
+
+An rt_annot_internal stores a model-space anchor, a set of 2D vertices, and
+either a model plane basis or screen-space placement.  Its segment types cover
+lines, circular arcs, NURBS, Beziers, text, and filled polygons.  Segment style
+records include a semantic role, line pattern and width, RGBA color, font or
+symbol, text scaling, and decorative flags.
+
+The libged annotate command can create styled text, leaders, linear/angular/
+radius/diameter/ordinate dimensions, and bounding-box autodim groups.  It can
+regenerate derived annotations from saved attributes and associate them with
+source model objects without making them part of those objects' CSG trees.
+
+Annotations therefore have the properties needed for persistent markup:
+
+* database naming, paths, transforms, and instances;
+* survival across save/reopen and transfer to another application;
+* semantic segment roles and per-segment presentation data;
+* model-space and view-dependent screen-space placement; and
+* high-level generation metadata and source-object associations.
+
+An annotation is informational coverage, not a CSG solid.  The raytracing
+implementation preserves that distinction.  It prepares annotation coverage
+from active database paths and composites it on primary samples after ordinary
+geometry has established color and depth.  Model-space annotation coverage is
+occluded by nearer solid geometry.  Screen-space annotation coverage uses the
+render view and image dimensions.  Annotations do not acquire artificial
+volume, alter model bounds, cast shadows, or participate in reflection and
+refraction rays.
+
+Visible coverage is returned in back-to-front layers.  Model-space layers are
+ordered by ray distance, followed by view-space layers ordered by projected
+anchor depth; segment draw order breaks ties within an annotation.  This lets
+the image renderers apply source-over alpha consistently while retaining the
+single-frontmost-hit query for compatibility.  NIRT continues to report one
+frontmost segment per model-space annotation rather than exposing these
+presentation layers as solid partitions.
+
+Color follows normal database tree state.  An annotation segment with an
+explicit RGBA style wins; otherwise it inherits the resolved combination-path
+color; otherwise it is white.  The legacy region-ID color table is not applied
+to annotation display or coverage, since annotations are informational leaves
+rather than regions.
+
+Prepared coverage uses the existing libbg RTree implementation.  A 3D RTree
+culls model-space annotation instances and a 2D local RTree culls segments
+inside every annotation.  Screen-space instance extents cannot be placed in a
+static world-space index because they change with the view, so those instances
+remain in a compact list while still using their local segment RTree.  This
+avoids a new spatial data structure and keeps prepared scenes immutable and
+safe for concurrent sample queries.
+
+The rt, ART, and ADRT/ISST integrations follow this overlay model.  ERT invokes
+rt with the database, displayed database paths, and view, so database-backed
+annotations can follow that route as well.  ERT does not currently serialize
+or pass application-owned view objects to its subprocess.
+
+NIRT keeps annotation diagnostics opt-in.  `useannot 1` enables model-space
+annotation intersections, and the `fmt a` action reports them.  Default NIRT
+output is unchanged.  Screen-space annotations are deliberately omitted
+because a NIRT ray does not define an image size, pixel footprint, or complete
+model-to-view mapping.
+
+
+View-object capabilities
+------------------------
+
+The core representation is bv_scene_obj.  A view set owns shared view objects,
+and a particular bview can also own local objects.  Independent views force
+local ownership.  bv_obj_put recycles objects when their application lifetime
+ends; the objects are not written to the .g database.
+
+The `view obj` command exposes the principal markup types:
+
+* Lines are model-coordinate polylines represented by vlists.  Generic scene
+  rendering can add arrowheads to segment ends.
+* Labels have a model anchor and optional model-space target/leader.  Their text
+  is projected and drawn as screen-facing HUD text.  If the command omits Z,
+  screen X/Y coordinates are projected onto the current view plane at creation
+  time.
+* Polygons include general polygons, circles, ellipses, rectangles, and
+  squares.  Creation and editing operate on a saved view plane with a view-Z
+  offset.  Selection, movement, open/close, area, overlap, and 2D boolean
+  operations are supported.  A "fill" is currently a set of hatch lines, not
+  an opaque surface.  Polygons can be explicitly imported from or exported to
+  persistent sketch primitives.
+
+Scene objects also provide object color, transparency, pixel line width, dash
+style, and arrow dimensions.  Not all of these properties have a uniform
+command interface for every object type.
+
+There are other visual layers that should not be conflated with these objects.
+Legacy gv_tcl data/sdata arrows, axes, lines, labels, and polygons have a
+separate representation.  Faceplate/HUD elements are screen UI.  Database
+objects temporarily cloned for editing are also transient, but retain database
+primitive semantics and live in the database-object scene containers.
+
+
+Comparison
+----------
+
+Both annotations and view objects eventually produce display-manager pixels,
+often through bv_scene_obj and vlists.  That common final path does not make
+their source representations equivalent.
+
+Annotations are the canonical choice when markup must be saved, named in a
+database, instanced, shared, raytraced from a .g file, or understood by another
+program.  View objects are the right choice for application/session state that
+must be created and manipulated without changing the database.
+
+Some view objects have a direct annotation analogue, but conversion is not
+lossless without type-aware logic:
+
+* A view polyline maps naturally to line segments, including arrowhead style.
+* A polygon outline and its hatch lines can map to annotation coverage, but the
+  existing view "fill" must not silently become an opaque filled polygon.
+* A label is hybrid.  Its target and leader are model-space geometry while its
+  glyphs remain screen-facing at a projected model anchor.
+* Arbitrary vlist extraction is insufficient because display-matrix commands,
+  point and line sizes, child objects, glyph behavior, and semantic roles are
+  otherwise lost or ambiguous.
+
+View-object line width and annotation line width also need an explicit unit
+policy.  View-object widths are display pixels.  Screen-space annotation widths
+are view-dependent coverage widths.  A model-space annotation's line_width is
+used as a scale on the automatically selected model-space stroke width.  In
+other words, removing a separate raytracer "stroke-width" setting does not
+remove width control: the annotation's existing line_width/`--line-width`
+style remains the single user-facing control.  The renderer interprets it in
+the appropriate coordinate system rather than exposing a second competing
+knob.
+
+
+Recommended transient rendering design
+--------------------------------------
+
+Do not make librt or the raytracers depend on bv_scene_obj, libbv application
+containers, or libged.  Do not manufacture tubes, TGCs, meshes, or geometry
+just in front of the view plane.  Artificial geometry would affect bounds,
+depth, shadows, and secondary rays, and would make antialiasing and perspective
+behavior depend on a fabricated distance from the eye.
+
+If transient view objects need to appear in raytraced output, use an explicit
+application-to-renderer snapshot boundary:
+
+1. The GED/application layer snapshots selected live view objects when a render
+   is launched.  This is the only layer that owns their lifetime and semantics.
+2. A type-aware adapter in libged or libbv translates supported objects into a
+   renderer-neutral overlay/coverage description.  It should reuse annotation
+   coverage preparation through an in-memory annotation-instance interface,
+   rather than duplicate line, glyph, and arrow intersection algorithms.
+3. The snapshot travels with the exact view state: model-to-view transform,
+   perspective, image width and height, and sampling convention.  In-process
+   renderers can receive it directly; subprocess rt/ERT needs an explicit
+   serialized side channel.
+4. Raytracers composite this coverage only on primary image samples, after
+   geometry color and depth are known.  It remains outside CSG, model bounds,
+   shadows, reflections, and refractions.
+5. Initial support should be deliberately finite: polyline strokes and
+   arrowheads, label glyphs and leaders, polygon outlines, and polygon hatches.
+   Faceplate elements and legacy gv_tcl overlays should be separately specified
+   opt-in layers rather than being included accidentally.
+6. Snapshot records should preserve object name and type, RGBA, visibility,
+   local/shared ownership, presentation units, and source identity.  Unsupported
+   object types should produce a deliberate diagnostic or documented omission.
+
+Application-side compositing of the interactive display-manager layer over an
+ERT framebuffer is a viable smaller feature if only ERT screen captures need
+view objects.  It would closely reproduce the interactive view, but it would
+not make the objects available to standalone rt, ART, or ISST output and should
+not be presented as the general renderer solution.
+
+Conversion between the two systems should remain explicit.  Automatically
+turning view objects into database annotations would violate their deliberately
+transient contract.  A user command to promote selected view markup to a
+persistent annotation would be reasonable; polygon-to-sketch export is an
+existing precedent.  Conversely, a view-object editor may be used as a front
+end for constructing an annotation, provided the persistence action is clear.
+
+
+Maintenance findings before implementing the adapter
+----------------------------------------------------
+
+Several representation issues should be corrected before type-driven snapshot
+conversion is relied upon:
+
+* Container-selection flags (BV_DB_OBJS, BV_VIEW_OBJS, BV_LOCAL_OBJS, and
+  BV_CHILD_OBJS) reuse numeric bits also used for semantic s_type_flags
+  (BV_DBOBJ_BASED, BV_VIEWONLY, BV_LINES, and others).  bv_obj_create currently
+  assigns its type argument to s_type_flags.  A local view object consequently
+  acquires BV_LINES incidentally, while the line command does not explicitly
+  tag lines.  The namespaces should be separated and every object explicitly
+  typed.
+* Lines, labels, and polygons should have unambiguous type tags and per-type
+  lifecycle/extraction callbacks.  Snapshot conversion must inspect typed data,
+  not infer semantics from a rendered vlist.
+* Label construction should initialize every bv_label field deliberately,
+  including size, line_flag, anchor, arrow, and target.  Current construction
+  does not establish all of those values on every path.
+* The status of legacy gv_tcl overlay data should be decided explicitly:
+  migrate it, support it through a separate adapter, or document its exclusion.
+* Screen-space annotation display and raytrace units need a parity test.  The
+  legacy display-matrix vlist path and the ray coverage path currently express
+  screen offsets through different mechanisms; their apparent agreement must
+  be verified rather than assumed.
+* Per-segment annotation colors exceed the legacy interactive vlist object's
+  single-color model.  Uniform annotate-command colors work today, but true
+  multicolor interactive parity may require split scene objects or a richer
+  display representation.
+
+Tests for a future adapter should cover shared and local ownership, cleanup,
+label anchors and leaders, polygon holes/hatches/view-Z, perspective, pixel
+width, model-space occlusion, screen-space overlay, and serialization across
+the ERT subprocess boundary.
+
+
+Conclusion
+----------
+
+Persistent annotations and transient view objects should continue to coexist.
+The annotation coverage scene is the correct renderer-side abstraction for
+both saved annotations and a future snapshot of transient markup.  The
+application should perform the semantic conversion and transfer a view-complete
+snapshot; the renderer should treat it as primary-ray coverage, never as
+invented scene geometry.

doc/asciidoc/system/man1/nirt.adoc

@@ -132,6 +132,10 @@ Command to set the backout flag. With no option, prints the current value. When
 Sets the _useair_ member of the __rt_i__ structure to the integer _n_. If _n_ is 0, then *nirt* ignores any air in the geometry. If this command is invoked with no arguments, *nirt* prints the current value of _useair_. Default is 0.
 
 
+*useannot [* _n_*]*::
+Enables (_n_=1) or disables (_n_=0) diagnostic annotation intersections. When enabled, *nirt* evaluates all model-space annotations below the active database paths and invokes the *a* output action for each crossing. Screen-space annotations are not evaluated because an individual NIRT ray does not supply a pixel size or model-to-view transform. If this command is invoked with no arguments, *nirt* prints the current value. Default is 0.
+
+
 *overlap_claims [* _n_*]*::
 Specifies how to handle multiple regions' claims to segments of a ray. If _n_ is 0 or "resolve", then all overlaps are resolved in favor of a single region and any other claimants are ignored. If _n_ is 2 or "rebuild_all", then all overlaps are rebuilt, so any overlapping regions along the ray create individual (geometrically intersecting) partitions. If _n_ is 3 or "retain", then all overlaps are retained. In this case, the resulting partitions are always geometrically disjoint, each one is owned by a single region according to the current overlap resolution strategy, but every claimant is recorded. If _n_ is 1 or "rebuild_fastgen", then *nirt* takes on _FASTGEN_ behavior, so overlaps of plate-mode primitives are rebuilt, but other overlaps are retained. This command is useful with the __claimant_count__, __claimant_list__, and __claimant_listn__ output items. Default is "resolve".
 
@@ -197,7 +201,11 @@ Quits *nirt*.
 [[output_statements]]
 == OUTPUT STATEMENTS
 
-*nirt* allows the user a high degree of control, via the _fmt_ command, over what information gets printed out for each ray and in what format. There are six types of output statements, each of which is executed under appropriate circumstances. The types and their use are:
+*nirt* allows the user a high degree of control, via the _fmt_ command, over what information gets printed out for each ray and in what format. There are eight types of output statements, each of which is executed under appropriate circumstances. The types and their use are:
+
+
+*a*::
+Annotation. Executed once for each model-space annotation crossing when *useannot* is enabled. Crossings are reported in ray-distance order after the ordinary solid shotline output. The default format is empty.
 
 
 *r*::
@@ -278,6 +286,41 @@ azimuth of view (i.e., of ray direction).
 elevation of view (i.e., of ray direction).
 
 
+Annotation Information
+
+
+*annot_path*::
+full database path of the annotation primitive.
+
+
+*annot_segment*::
+zero-based annotation segment index.
+
+
+*annot_role*::
+semantic role of the annotation segment, such as `text`, `leader`, `dimension`, or `arrowhead`.
+
+
+*annot_x*::
+_x_ coordinate of the annotation hit point.
+
+
+*annot_y*::
+_y_ coordinate of the annotation hit point.
+
+
+*annot_z*::
+_z_ coordinate of the annotation hit point.
+
+
+*annot_d*::
+_d_ grid-plane coordinate of the annotation hit point.
+
+
+*annot_visible*::
+1 if no non-air solid surface is closer to the ray origin, otherwise 0. Occluded annotation crossings are still reported so they can be inspected during debugging.
+
+
 
 Partition Information
 
@@ -521,6 +564,15 @@ nirt -A material_id -A los_density model.g all
 Each region the ray passes through is reported together with its _material_id_ and _los_density_ attribute values via the *attributes* output item, which is useful for vulnerability and mass-properties analyses.
 ====
 
+.Reporting annotation crossings for debugging
+====
+[source,bash]
+----
+nirt -s -e 'useannot 1;fmt a "%s %d (%g %g %g) visible=%d\n" annot_path annot_segment annot_x annot_y annot_z annot_visible;s' model.g component dimensions
+----
+This opt-in format reports every model-space annotation crossing and whether it is in front of the first solid surface. Ordinary NIRT output is unchanged unless *useannot* is enabled and an *a* format is supplied.
+====
+
 
 [[hints]]
 == HINTS

doc/asciidoc/system/man1/pixdiff.adoc

@@ -4,54 +4,38 @@
 :manmanual: BRL-CAD
 
 == NAME
-pixdiff - highlight differences between two pix or bw images
+pixdiff - highlight samples which differ between two pix or bw files
 
 == SYNOPSIS
 
 [source]
 ----
 pixdiff _file1.pix_ _file2.pix_ > _out.pix_
-pixdiff -b _file1.bw_ _file2.bw_ > _out.bw_
+pixdiff _file1.bw_ _file2.bw_ > _out.bw_
 ----
 
 
 [[description]]
 == DESCRIPTION
 
-By default, *pixdiff* compares raw *pix*(5) images as RGB pixels and writes
-a PIX difference image to standard output. When all three channels of a
-pixel match, it writes a dim grayscale version of that pixel to all three
-output channels. When any channel differs, each matching channel is written
-as zero, each channel differing by one is written as 192, and each channel
-differing by more than one is written as 255.
+*pixdiff* compares two raw *pix*(5) or *bw*(5) images byte by byte and
+writes a difference image of the same size to standard output. Bytes that
+differ by +/-1 are written at intensity 192, while larger differences are
+written at full intensity 255. Matching bytes are written at half their
+original intensity, providing visual context for locating differences.
 
-With *-b*, *pixdiff* compares raw *bw*(5) images one byte at a time and
-writes a BW difference image. Differing samples use the same 192 and 255
-intensities; matching samples are written at half their original intensity.
-
-No image dimensions are needed. The two inputs must be the same size in
-bytes, and PIX input must contain complete three-byte pixels. Otherwise,
-*pixdiff* reports an error. Either input file may be given as
+Because the comparison does not interpret pixel structure, no format option
+is needed and *pixdiff* is independent of image resolution. The two inputs
+must be the same size in bytes; otherwise the program reports the mismatch
+and exits. Either input file may be given as
 *-* to read from standard input. Output must be redirected; *pixdiff*
 refuses to write a raw image stream to a terminal.
 
-Upon completion, *pixdiff* reports matching, off-by-one, and larger-difference
-channel or sample counts on standard error.
+Upon completion, the count of bytes that matched, that were off by 1, and
+that were off by many is reported on standard error.
 
-*bwdiff -n -b* produces a similar BW highlight image, but writes 255 rather
-than 192 for off-by-one differences. The default *bwdiff*(1) output instead
-shows signed differences around a neutral gray.
-
-
-[[options]]
-== OPTIONS
-
-*-b*::
-Treat both inputs and the output as one-byte-per-pixel BW images. Without
-this option, use three-byte-per-pixel PIX images.
-
-*-h*::
-Show usage information.
+Also see *bwdiff*(1), which can also operate on *pix*(5) files and can
+produce other kinds of difference output.
 
 
 [[examples]]
@@ -75,10 +59,10 @@ on standard error give a quick quantitative sense of how much changed.
 ====
 [source,bash]
 ----
-pixdiff -b depth_a.bw depth_b.bw > depth_diff.bw
+pixdiff depth_a.bw depth_b.bw > depth_diff.bw
 ----
-Writes a grayscale highlight image using byte-wise comparison. Matching
-samples retain half their input intensity as context.
+Writes a grayscale highlight image using the same byte comparison and
+statistics as the color example; no format option is needed.
 ====
 
 

doc/asciidoc/system/man3/libanalyze.adoc

@@ -18,8 +18,7 @@ Link with *-lanalyze*.
 The public interface is organized across the headers included by
 _analyze.h_, including _analyze/diff.h_, _analyze/density.h_,
 _analyze/grid.h_, _analyze/heal.h_, _analyze/info.h_, _analyze/pnts.h_,
-_analyze/contour.h_, _analyze/polygonize.h_, _analyze/nirt.h_,
-_analyze/voxelize.h_, and
+_analyze/polygonize.h_, _analyze/nirt.h_, _analyze/voxelize.h_, and
 _analyze/worker.h_.
 
 == DESCRIPTION
@@ -33,7 +32,7 @@ to gather the information they need.
 Facilities provided by the library include model differencing, material
 density and mass property analysis, ray-based gridding and sampling,
 geometry healing, model information queries, point-cloud generation and
-processing, manifold dual contouring and voxelization of geometry, and the
+processing, polygonization and voxelization of geometry, and the
 interactive ray-tracing engine used by the *nirt*(1) command. A parallel
 worker interface supports distributing analysis work across multiple
 threads.

doc/asciidoc/system/man3/librt.adoc

@@ -200,56 +200,6 @@ ____
 
 The LIBRT_DEBUG environment variable may be set to the hexadecimal bit vector as described by the DEBUG_* flags listed in the raytrace.h header. Similarly, many LIBRT applications provide a "-x #" command-line option which will take precedence over the environment variable.
 
-The RT_DEBUG_RTCMP flag (0x10000000) (set with LIBRT_DEBUG=10000000)
-sets up rtcmp capture, but recording starts disabled. Set
-LIBRT_RTCMP_FILE before enabling the flag, then call
-rt_rtcmp_capture_set_enabled(1) from the application or a debugger to
-record selected rays. Call rt_rtcmp_capture_set_enabled(0) to stop
-recording. Enabling returns -1 if the debug flag is unset or the output
-writer failed to initialize; disabling returns 0 even before setup. A
-ray must begin while recording is enabled, and leaving recording
-disabled when it finishes suppresses its record. The debug flag must be
-set before concurrent ray shooting; use the setter to change recording
-state while rays are running.
-
-Captured hits append one high-precision rtcmp JSON record per ray with
-final partitions. Misses produce no record. librt never clears the
-output file; the caller must clear or rotate it before a new comparison
-as needed. Separate processes must not write to the same file
-concurrently. Pending records are flushed within two seconds under
-normal operation and during normal process shutdown; call
-rt_rtcmp_capture_flush() to wait for pending records when using a
-debugger. Disabling recording does not discard records already queued. A
-process terminated abnormally may lose pending records. rtcmp comparison
-files should use the same hit-only capture mode on both sides; its
-generation mode also includes misses. The output can be expanded to also
-include primitive segments if LIBRT_RTCMP_PRIMITIVES is set to 1 before
-capture is set up.
-
-For example, run a single-threaded ray tracing application in a debugger
-and stop at `rt_shootray` for the suspect ray. Capture just that call,
-then compare it with a capture of the same ray from a reference build:
-
-[source,shell]
-----
-$ LIBRT_DEBUG=10000000 LIBRT_RTCMP_FILE=trial.json \
-    gdb --args ./ray_application model.g
-(gdb) break rt_shootray
-(gdb) run
-(gdb) call rt_rtcmp_capture_set_enabled(1)
-(gdb) finish
-(gdb) call rt_rtcmp_capture_set_enabled(0)
-(gdb) call rt_rtcmp_capture_flush()
-(gdb) quit
-$ rtcmp -c reference.json trial.json
-----
-
-Continue to the suspect ray before enabling capture if the breakpoint
-first stops on another ray. Produce `reference.json` with the same ray
-and the same capture mode; set LIBRT_RTCMP_PRIMITIVES=1 for both runs
-when primitive segments are needed to locate a difference before Boolean
-evaluation.
-
 The LIBRT_V4FLIP environment variable may be set to a true or false value to respectively force or disable endianness interpretation of v4 geometry database files. The default behavior is to automatically detect whether flipping the endian interpretation will help LIBRT read a binary-incompatible v4 geometry database file. Setting LIBRT_V4FLIP will override automatic detection.
 
 

doc/asciidoc/system/man5/nirt_formatting.adoc

@@ -25,7 +25,7 @@ where the "char" argument is replaced by a single character identifying the cate
 [[nirt_format_categories]]
 == Output Categories
 
-NIRT offers seven categories of output reporting:
+NIRT offers eight categories of output reporting:
 
 [cols="3*"]
 [%noheader]
@@ -33,6 +33,9 @@ NIRT offers seven categories of output reporting:
 |*Category Name*
 |*Char*
 |*Description*
+|*Annotation*
+|*a*
+|Output for each model-space annotation crossing when annotation reporting is enabled with *useannot 1*. The default annotation format is empty.
 |*Ray*
 |*r*
 |Report of ray information. The formatting associated with r will print once per ray, regardless of whether a region is encountered.
@@ -62,7 +65,7 @@ NIRT offers seven categories of output reporting:
 [[nirt_format_str]]
 == Format Strings
 
-NIRT's format specifications are a subset of those found in printf implementations. Specifically, the d, i, and f types are supported, as well as width and precision specifiers.
+NIRT's format specifications are a subset of those found in printf implementations. The d, i, f, e, E, g, G, and s types are supported, as well as width and precision specifiers.
 
 
 [[nirt_format_keys]]
@@ -100,6 +103,33 @@ Unlike the full C/C++ printf, NIRT accepts only pre-defined string keys that ide
 
 Footnote and Miss categories have no particular keys unique to them, although they (like all categories) can use the Ray keys if desired. Partition, Overlap, and Gap, on the other hand, have keys that are unique to their contexts and valid only within those contexts:
 
+.Annotation Variables
+[cols="2*"]
+[%noheader]
+|===
+|annot_path
+|Full database path of the annotation primitive.
+|annot_segment
+|Zero-based annotation segment index.
+|annot_role
+|Semantic role of the annotation segment, such as `text`, `leader`, `dimension`, or `arrowhead`.
+|annot_x
+|x coordinate of the annotation hit point.
+|annot_y
+|y coordinate of the annotation hit point.
+|annot_z
+|z coordinate of the annotation hit point.
+|annot_d
+|d grid-plane coordinate of the annotation hit point.
+|annot_visible
+|1 if no non-air solid surface is closer to the ray origin, otherwise 0.
+|===
+
+Annotation reporting is diagnostic and opt-in. It returns all model-space
+annotation crossings in ray-distance order, including occluded crossings.
+Screen-space annotations are not reported because a NIRT shot does not define
+the pixel dimensions and model-to-view transform needed to evaluate them.
+
 .Partition Variables
 [cols="2*"]
 [%noheader]

doc/asciidoc/system/mann/annotate.adoc

@@ -40,12 +40,11 @@ Coordinates and explicit sizes use the database's current modeling units.
 Creation commands draw the new annotation immediately.  Use *--no-draw* to
 create it without changing the display.  Names must not already exist.
 
-Model-space annotations also participate in raytracing.  Librt expands their
-strokes and filled areas into transient mesh geometry during preparation; the
-generated geometry is not written to the database.  It uses the normal region
-material and color resolved from the containing combination tree.  Screen-space
-annotations remain display-only because primitive preparation has no render
-view or image dimensions with which to place them.
+Ray-traced annotations are zero-thickness informational coverage rather than
+solid regions.  Model-space annotations are depth-tested against ordinary
+geometry.  Screen-space annotations are generated for the current camera and
+composited in the foreground; they do not participate in shadows, reflections,
+line-of-sight thickness, or Boolean evaluation.
 
 The supported dimensions follow the general presentation conventions of ISO
 129-1.  Lettering defaults to OSIFont, an ISO 3098-style technical font shipped
@@ -209,9 +208,9 @@ Set the persistent annotation color using integer red, green, and blue values
 from 0 through 255.
 
 *--line-width width*::
-Set the display line width.  For a model-space raytrace, this value scales the
-physical stroke width inferred from the annotation's text height or overall
-extent.
+Set the presentation line width.  Display and screen-space ray rendering use
+the value as a pixel width.  Model-space ray rendering uses it to scale the
+automatic stroke width derived from the annotation's text height or extent.
 
 *--line-style style*::
 Set semantic line-pattern metadata to *continuous*, *dashed*, *dotted*,

doc/asciidoc/system/mann/facetize.adoc

@@ -10,7 +10,7 @@ facetize - Creates __new_object__ using __old_object__ as a guide, creating an e
 
 [source]
 ----
-facetize [-h] [-v] [-q] [-n] [-r] [-j jobs] [-s suffix] [-p prefix] [--in-place] [--max-time #] [--max-pnts #] [--resume] [--methods m1,m2,...] [--method-opts METHOD opt1=val opt2=val...] [--no-empty] [--log-file filename] [--nmg-booleval] [--disable-fixup] [--perturb] [--no-perturb] [--perturb-sa-tol #] [--perturb-vol-tol #] [--tolerate-failures] [-B -t threshold] {old_object new_object}
+facetize [-h] [-v] [-q] [-n] [-r] [-s suffix] [-p prefix] [--in-place] [--max-time #] [--max-pnts #] [--resume] [--methods m1,m2,...] [--method-opts METHOD opt1=val opt2=val...] [--no-empty] [--log-file filename] [--nmg-booleval] [--disable-fixup] [--perturb] [--no-perturb] [--perturb-sa-tol #] [--perturb-vol-tol #] [--tolerate-failures] [-B -t threshold] {old_object new_object}
 ----
 
 
@@ -45,8 +45,6 @@ By default, the Manifold BoT path can make use of a coplanarity-avoidance pertur
 
 When the perturbation step is active, per-region *facetize* validates the initial Boolean output against the original CSG geometry using a Crofton stochastic raytrace. The default tolerance thresholds that trigger a perturb retry are 10% for both surface area and volume; these can be adjusted with *--perturb-sa-tol* and *--perturb-vol-tol*. At the end of a per-region facetization run *facetize* prints a concise validation summary listing how many regions passed, triggered a perturb retry, or encountered one of the following noteworthy conditions:
 
-_Unbounded halfspace (validation unavailable)_: Crofton surface-area and volume comparisons require a finite sampling domain. A standalone halfspace, or another CSG result without finite bounds, is therefore not sampled and is counted as having validation unavailable. A bounded object intersected with or subtracted by a halfspace can still be validated because its finite operand supplies the sampling bounds.
-
 _Sub-mm / few-hit bypass (accepted with note)_: When the Crofton sampler finds only a very small number of ray-surface intersections (below an internal threshold of 50 crossings), the sampling noise alone would exceed 14% relative error, making any SA comparison meaningless. Geometry in this category is typically sub-millimetre in cross-section (e.g., a TGC with a radius under 0.2 mm) but the stochastic sampler has confirmed that some geometry is present. The BoT conversion is accepted and the region is counted separately in the summary. Users should verify that such extreme geometry reflects actual modeling intent rather than an accidental near-degenerate primitive.
 
 _Zero-hit replacement_: When the Crofton sampler fires rays but records zero intersections with the CSG geometry while the corresponding BoT is non-empty, *facetize* concludes that the Boolean evaluation most likely produces empty geometry (for example, a subtractor that fully engulfs the base object). In this case the non-empty BoT is discarded and replaced with an empty BoT to match the expected raytrace behavior. If this replacement appears unexpected for a given region, the *lint* command can be used to inspect the result.
@@ -63,21 +61,13 @@ Different tolerance settings will create finer or coarser meshes. See the *tol*
 
 Primitive tessellation is performed in a child process against a private working copy of the database in the BRL-CAD cache.  The original database is not modified as individual primitive meshes are generated.  Only after the conversion succeeds are the requested results transferred from the working database.  This separation also permits an interrupted conversion to retain its intermediate state for use with *--resume*.
 
-The parent command sends primitive names to a small pool of persistent tessellation processes.  Each child reports a result for every name, so a timeout, conversion failure, or child-process crash can be associated directly with the affected primitive.  A child is restarted when necessary and processing continues with the next name; the parent does not subdivide and retry the whole group to discover which object failed.  Unless output is suppressed with *--quiet*, long tessellation batches also report their completed-object count and active-worker state at approximately five-second intervals.
+The parent command sends primitive names to a persistent tessellation process one at a time.  The child reports a result for each name, so a timeout, conversion failure, or child-process crash that occurs before a database write can be associated directly with the affected primitive.  The child is restarted when necessary and processing continues with the next name; the parent does not subdivide and retry the whole group to discover which object failed.  Unless output is suppressed with *--quiet*, long tessellation batches also report their current object and completed-object count at approximately five-second intervals.
 
-By default, *facetize* uses at most two worker processes and only enables the second when both CPU and memory information indicate it is safe.  It reserves the larger of one GiB or one quarter of total physical memory for the rest of the system.  Starting another job also requires at least 512 MiB beyond that reserve, or the largest worker resident-memory measurement seen so far if that is greater.  If the platform cannot report the required memory information, *facetize* falls back to one worker.  The available CPU count is divided among worker processes and is also imposed as the thread limit inside each worker, preventing nested algorithms from independently using all CPUs.  *--jobs* changes the maximum process count; it does not bypass the memory admission checks.
-
-Per-region BoT conversion uses this worker policy for both independent CSG validation and initial Manifold Boolean evaluation.  When more than one worker is available, persistent validation processes each open their own read-only handle to the original database and compute Crofton surface-area and volume references for different regions.  Persistent Boolean processes then evaluate different region trees using independent read-only handles to one shared snapshot of the tessellated working database.  That snapshot also serves as the rollback checkpoint, so increasing *--jobs* does not create additional full database copies.  This process isolation avoids sharing mutable raytrace preparation or database state between concurrent evaluations.  Each result remains associated with the requested region name.  The parent commits completed BoTs through one writer and continues to perform hierarchy changes, validation decisions, and any perturb retries serially.  A one-worker run retains the established in-process order and avoids the snapshot phase.
-
-Validation and Boolean workers report their completed-region count and active-worker state approximately every five seconds.  When *--max-time* is specified, it also bounds each parallel CSG reference or Boolean request; a timed-out process is terminated, and the failure is associated with that region.  Failure or timeout of supplementary CSG validation does not cause unrelated regions to be retried.
-
-Database updates use a separate write handshake.  After preparing a mesh in memory, a child announces that it is ready to write and reports the approximate mesh size and its current resident memory, then waits.  The parent grants write permission, the child acknowledges when it is starting the write, and the child explicitly reports when the update is complete.  Consequently, the normal per-object tessellation limit, including *--max-time* and method-specific limits, applies to mesh preparation and does not interrupt an acknowledged write.
-
-With multiple workers, primitive tessellation children write to small staging databases rather than modifying the shared working database.  Region Boolean children read one immutable cache-local snapshot while completed results are committed to the working database.  The snapshot doubles as the recovery checkpoint and is opened independently by each worker; it is not duplicated per worker.  For very large files, snapshot preparation reports the amount copied and elapsed time at five-second intervals.  A region result is first completed in memory and enters the write handshake only when it is ready to replace the prior staged result.  After a child reports a completed staging write, the parent submits that result to a single persistent writer process.  The writer announces that it is ready, waits for the parent's permission, acknowledges that it is starting the write, transfers the result through one persistent working-database handle, and reports completion.  This permits tessellation and initial region Boolean evaluation to proceed in parallel without concurrent database allocation updates or repeated full directory scans, while allowing the parent to supervise the canonical database write itself.  With one worker, the established direct path avoids the extra staging I/O and database snapshot.
+Database updates use a separate write handshake.  After preparing a mesh in memory, the child announces that it is ready to write and reports the approximate mesh size, then waits.  The parent acknowledges that it is waiting for the write before the child may update the working database.  The child explicitly reports when the update is complete.  Consequently, the normal per-object tessellation limit, including *--max-time* and method-specific limits, applies to mesh preparation and does not interrupt an acknowledged database write.
 
 Before the first primitive conversion, the parent measures write throughput on the same cache filesystem as the working database.  It writes a temporary file in 1 MiB chunks until approximately one second has elapsed or 128 MiB has been written, then removes the file immediately.  Successful mesh writes of at least 1 MiB subsequently refine this profile.  For each database update, the parent allows ten times the projected write time, with a minimum of five seconds and a maximum of five minutes.  This measured allowance protects legitimate long writes on slow storage while still detecting a genuinely stalled update.
 
-If a child fails before acknowledging that its write has started, the working database is known to be unchanged and no rollback is needed.  In parallel mode, a failure during a child write affects only that worker's disposable staging database, so the failed object is reported and the worker can be restarted without replaying unrelated objects.  A serial child failure during a direct write, or a writer-process failure after acknowledging write start, is ambiguous because the working database may contain a partial update.  In those cases *facetize* restores its cache-local checkpoint and replays only the conversions not already tied to a specific failure.  The original database remains untouched during this recovery.
+If a child fails before write permission is granted, the working database is known to be unchanged and no rollback is needed.  A failure, crash, or write-deadline expiration after permission is granted but before completion is ambiguous: the working database may contain a partial update.  In that case *facetize* restores its cache-local checkpoint and replays only the primitive conversions already known to have succeeded.  The original database remains untouched during this recovery.
 
 
 [[facetize_options]]
@@ -104,10 +94,6 @@ Instructs *facetize* to produce NMG solids rather than BoT solids.
 Rather than converting all specified inputs to one top level object, walk specified trees and individually convert regions to facetized objects. This mode is useful when wanting to display models in shaded mode while preserving color and region ID properties. If no regions are present a normal conversion is performed instead.
 
 
-*-j, --jobs #*::
-Set the maximum number of facetize worker processes. A value of 1 forces serial processing. The default value of 0 selects a conservative automatic limit based on available CPUs and memory. Memory checks may delay or reduce parallel work even when a larger value is requested. In per-region BoT mode the same limit controls parallel CSG reference sampling and initial Manifold Boolean evaluation.
-
-
 *-s, --suffix*::
 User specified suffix to use avoiding naming conflicts (naming conflict avoidance is needed when creating objects in *-r* mode.)
 
@@ -133,7 +119,7 @@ Resume an incomplete facetization process (will utilize cached data from a previ
 
 
 *--methods*::
-Accepts a comma separated list of one or more facetization methods to use. The default tries NMG tessellation first, Manifold Dual Contouring (MDC) second, and SPSR as a final fallback. Other listed methods may be selected explicitly if a user wants to try a specific methodology.
+Accepts a comma separated list of one or more facetization methods to use. The default is to try all algorithms available as needed, but if a user knows they want to try a specific methodology they can call it out using this option.
 
 
 *--method-opts*::
@@ -245,15 +231,21 @@ mged>**facetize --tolerate-failures large_model large_model.partial.bot**
 
 ====
 
-.Refine a Manifold Dual Contouring conversion
+.Increase Decimation in CM Mode
 [example]
 ====
-This selects MDC directly and requests a finer cell size with a larger ray
-budget than the default:
+mged>**facetize --methods CM --method-opts feature_size=0.06 d_feature_size=2.0 ellipse.r ell.bot**
 
-mged>**facetize --methods MDC --method-opts "MDC feature_size=0.1 max_rays=400000" ellipse.r ell.bot**
+The default settings on this object in CM mode will produce a little over 8000 faces in the final bot, but the addition of the *d_feature_size* option with a >1.5 value has dramatically reduced the count to 104.
 
-Smaller feature sizes increase mesh detail, ray count, memory use, and runtime.
+====
+
+.Decrease Decimation in CM Mode
+[example]
+====
+mged>**facetize --methods CM --method-opts feature_size=0.06 d_feature_size=0.01 ellipse.r ell.bot**
+
+If the intent is preserve a more detailed mesh, reducing *decimation-feature-size* will keep more of the original triangles - in the above example, almost 250 thousand of the original 6.9 million faces were retained. NOTE: use lower values here cautiously, as extremely large meshes can be very difficult for many tools to deal with.
 
 ====
 

doc/inventory.adoc

@@ -679,21 +679,23 @@ Convert CSG geometry to a 3D voxel grid by:
 fill fraction
 * Level-of-detail parameter controls sub-voxel sampling
 
-==== 7.5 Manifold Dual Contouring
+==== 7.5 Implicit Surface Polygonizer (Marching Cubes / Tetrahedral)
 
-`contour.cpp` evaluates CSG geometry directly with `librt` and constructs
-closed, oriented triangle meshes:
+An implicit surface polygonizer (Bloomenthal tetrahedral marching) for
+converting any inside/outside function to a triangle mesh.
 
-* Recursively refines an implicit octree using cached axis-aligned rays along
-node box edges and interior quarter- and center-lines
-* Reuses complete ray partitions at the conforming finest grid and retains
-evaluated boundary points with analytic surface normals
-* Resolves ambiguous faces consistently and permits multiple vertices in a
-cell to preserve local topology
-* Places vertices with bounded quadratic-error minimization
-* Validates the resulting indexed mesh before returning it
-* Backs the legacy `analyze_polygonize` interface without retaining the old
-  seed-driven marching-tetrahedra implementation
+`polygonizer.c` (based on Jules Bloomenthal’s "`An Implicit Surface
+Polygonizer,`" Graphics Gems IV, 1994) converts any implicit function
+(inside/outside test) to a triangle mesh:
+
+* Maintains a hash table of evaluated lattice corners and edges
+* Expanding front: starts at a seed point and spreads outward, testing
+adjacent cube faces
+* Each transitive cube is optionally decomposed into 6 tetrahedra for
+sub-cube accuracy
+* Edge-surface crossing is found by bisection (`converge`) then the
+vertex normal is estimated from finite differences
+* Used in `libanalyze` for converting CSG regions to faceted form
 
 ==== 7.6 Geometry Differencing (raydiff)
 
@@ -1710,6 +1712,7 @@ import
 runtime
 * *MT19937*: Mersenne Twister PRNG
 * *pbrt-v3* (derived): HLBVH construction algorithm
+* *Bloomenthal polygonizer*: implicit surface polygonization
 * *Screened Poisson Reconstruction* (Michael Kazhdan): surface
 reconstruction from point clouds
 * *Ball Pivoting Algorithm*: surface reconstruction

doc/legal/embedded/CMakeLists.txt

@@ -55,6 +55,7 @@ brlcad_manage_files(pbrt.txt                   ${DOC_DIR}/legal/embedded TRIGGER
 brlcad_manage_files(picohash.txt               ${DOC_DIR}/legal/embedded TRIGGER gist)
 brlcad_manage_files(pixcii.txt                 ${DOC_DIR}/legal/embedded TRIGGER libicv)
 brlcad_manage_files(point_in_polygon.txt       ${DOC_DIR}/legal/embedded TRIGGER libbg)
+brlcad_manage_files(polygonizer.txt            ${DOC_DIR}/legal/embedded TRIGGER libanalyze)
 brlcad_manage_files(predicates.txt             ${DOC_DIR}/legal/embedded TRIGGER libbg)
 brlcad_manage_files(profont.txt                ${DOC_DIR}/legal/embedded TRIGGER libbu)
 brlcad_manage_files(pstdint.txt                ${DOC_DIR}/legal/embedded TRIGGER libbu)

doc/legal/embedded/polygonizer.txt

@@ -0,0 +1,9 @@
+https://www.bloomenthal.com/JBloom/code/polygonizer.c
+
+Authored by Jules Bloomenthal, Xerox PARC.
+Copyright (c) Xerox Corporation, 1991.  All rights reserved.
+Permission is granted to reproduce, use and distribute this code for
+any and all purposes, provided that this notice appears in all copies.
+
+file:/src/libanalyze/polygonizer.c
+

doc/legal/embedded/tcl_code.txt

@@ -61,4 +61,5 @@ Author: Paul Mackerras (paulus@cs.anu.edu.au),
 
 file:/src/libtclcad/tkImgFmtPIX.c
 file:/src/tclscripts/cad_clrpick.tcl
+file:/src/tclscripts/html_library.tcl
 

include/analyze.h

@@ -31,7 +31,6 @@
 #include "common.h"
 
 #include "analyze/defines.h"
-#include "analyze/contour.h"
 #include "analyze/debug.h"
 #include "analyze/diff.h"
 #include "analyze/density.h"

include/analyze/CMakeLists.txt

@@ -1,6 +1,5 @@
 set(
   analyze_headers
-  contour.h
   debug.h
   defines.h
   density.h

include/analyze/contour.h

@@ -1,118 +0,0 @@
-/*                      C O N T O U R . H
- * BRL-CAD
- *
- * Copyright (c) 2026 United States Government as represented by
- * the U.S. Army Research Laboratory.
- *
- * This library is free software; you can redistribute it and/or
- * modify it under the terms of the GNU Lesser General Public License
- * version 2.1 as published by the Free Software Foundation.
- *
- * This library is distributed in the hope that it will be useful, but
- * WITHOUT ANY WARRANTY; without even the implied warranty of
- * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
- * Lesser General Public License for more details.
- *
- * You should have received a copy of the GNU Lesser General Public
- * License along with this file; see the file named COPYING for more
- * information.
- */
-/** @addtogroup libanalyze
- * @{ */
-/** @file analyze/contour.h
- *
- * Ray-driven manifold dual contouring of evaluated BRL-CAD geometry.
- */
-
-#ifndef ANALYZE_CONTOUR_H
-#define ANALYZE_CONTOUR_H
-
-#include "common.h"
-
-#include <stddef.h>
-
-#include "vmath.h"
-#include "analyze/defines.h"
-
-__BEGIN_DECLS
-
-#define ANALYZE_MDC_DEFAULT_MAX_RAYS 0
-#define ANALYZE_MDC_DEFAULT_MINIMUM_FREE_MEM 150000000
-#define ANALYZE_MDC_DEFAULT_MIN_DEPTH 2
-#define ANALYZE_MDC_DEFAULT_MAX_DEPTH 10
-#define ANALYZE_MDC_MAX_DEPTH 20
-#define ANALYZE_MDC_DEFAULT_MAX_TIME 600
-
-struct db_i;
-/** Resource and accuracy controls for analyze_mdc(). */
-struct analyze_mdc_params {
-    /** Requested maximum finest-cell edge length.  Zero selects an
-     * implementation default based on the model bounds. */
-    fastf_t feature_size;
-
-    /** Optional user ceiling for the automatically derived ray budget.
-     * Zero applies only geometry, depth, time, and memory limits. */
-    size_t max_rays;
-
-    /** Abort before available memory drops below this many bytes. */
-    size_t minimum_free_mem;
-
-    /** Minimum and maximum power-of-two grid depths. */
-    int min_depth;
-    int max_depth;
-
-    /** Maximum wall-clock processing time in seconds.  Zero is unlimited. */
-    int max_time;
-
-    /** Emit progress diagnostics when non-zero. */
-    int verbosity;
-};
-
-#define ANALYZE_MDC_PARAMS_DEFAULT { \
-    0.0, \
-    ANALYZE_MDC_DEFAULT_MAX_RAYS, \
-    ANALYZE_MDC_DEFAULT_MINIMUM_FREE_MEM, \
-    ANALYZE_MDC_DEFAULT_MIN_DEPTH, \
-    ANALYZE_MDC_DEFAULT_MAX_DEPTH, \
-    ANALYZE_MDC_DEFAULT_MAX_TIME, \
-    0 \
-}
-
-enum analyze_mdc_status {
-    ANALYZE_MDC_OK = 0,
-    ANALYZE_MDC_INVALID_INPUT = 1,
-    ANALYZE_MDC_PREP_FAILED = 2,
-    ANALYZE_MDC_NO_SURFACE = 3,
-    ANALYZE_MDC_RAY_LIMIT = 4,
-    ANALYZE_MDC_TIMEOUT = 5,
-    ANALYZE_MDC_MEMORY_LIMIT = 6,
-    ANALYZE_MDC_AMBIGUOUS = 7,
-    ANALYZE_MDC_NOT_MANIFOLD = 8
-};
-
-/**
- * Evaluate @p object with librt and construct an indexed triangle mesh.
- *
- * On success, the caller owns @p faces and @p vertices and must release them
- * with bu_free().  Failure leaves all outputs empty.
- */
-ANALYZE_EXPORT extern enum analyze_mdc_status
-analyze_mdc(int **faces, size_t *num_faces, point_t **vertices,
-	    size_t *num_vertices, const char *object, struct db_i *dbip,
-	    const struct analyze_mdc_params *params);
-
-__END_DECLS
-
-#endif /* ANALYZE_CONTOUR_H */
-
-/** @} */
-
-/*
- * Local Variables:
- * mode: C
- * tab-width: 8
- * indent-tabs-mode: t
- * c-file-style: "stroustrup"
- * End:
- * ex: shiftwidth=4 tabstop=8
- */

include/analyze/polygonize.h

@@ -43,29 +43,9 @@ struct analyze_polygonize_params {
     size_t minimum_free_mem;
     int verbosity;
 };
-#define ANALYZE_POLYGONIZE_DEFAULT_MAX_CYCLE_TIME 30
-#define ANALYZE_POLYGONIZE_DEFAULT_MINIMUM_FREE_MEM 150000000
-#define ANALYZE_POLYGONIZE_PARAMS_DEFAULT { 0, 0, \
-    ANALYZE_POLYGONIZE_DEFAULT_MAX_CYCLE_TIME, \
-    ANALYZE_POLYGONIZE_DEFAULT_MINIMUM_FREE_MEM, 0 }
-
-/**
- * Legacy polygonization interface, implemented by analyze_mdc().
- *
- * @p size is the requested maximum cell edge length.  @p p_s is retained for
- * source and binary compatibility but is no longer used as a starting seed.
- * The caller owns successful output arrays and must release them with
- * bu_free().  Failure leaves all outputs empty.
- *
- * @return 0 on success
- * @return -1 on invalid input or polygonization failure
- * @return 2 when a time limit is reached
- * @return 3 when the minimum-free-memory limit is reached
- */
+#define ANALYZE_POLYGONIZE_PARAMS_DEFAULT { 0, 0, 30, 150000000, 0 }
 ANALYZE_EXPORT extern int
-analyze_polygonize(int **faces, int *num_faces, point_t **vertices,
-	int *num_vertices, fastf_t size, point_t p_s, const char *obj,
-	struct db_i *dbip, struct analyze_polygonize_params *p);
+analyze_polygonize(int **faces, int *num_faces, point_t **vertices, int *num_vertices, fastf_t size, point_t p_s, const char *obj, struct db_i *dbip, struct analyze_polygonize_params *p);
 
 __END_DECLS
 

include/bg.h

@@ -55,7 +55,6 @@
 #include "bg/clip.h"
 #include "bg/lseg.h"
 #include "bg/obr.h"
-#include "bg/pca.h"
 #include "bg/plane.h"
 #include "bg/polygon.h"
 #include "bg/sat.h"

include/bg/clip.h

@@ -50,12 +50,12 @@ BG_EXPORT extern int bg_lseg_clip(fastf_t *xp1, fastf_t *yp1, fastf_t *xp2, fast
  * The RPP has faces parallel to the coordinate planes and is defined
  * by a minimum point and a maximum point.
  *
- * Despite its historical name, this function clips the finite line segment
- * from a to b.
+ * FIXME: the function name implies this takes a point,dir for a,b but
+ * it actually takes a line segment going from points a to b!
  *
  * Returns -
- * 0 if the segment does not hit the RPP,
- * !0 if the segment hits the RPP.
+ * 0 if ray does not hit RPP,
+ * !0 if ray hits RPP.
  *
  * Implicit Return -
  * if !0 was returned, "a" and "b" have been clipped to the RPP.

include/bg/obr.h

@@ -126,12 +126,9 @@ BG_EXPORT extern int bg_3d_coplanar_obr(point_t *center,
  * \endverbatim
  *
  *
- * @param[out] pnts     eight non-NULL, caller-provided point_t storage
- *                       locations for the oriented bounding box vertices
- * @param[in] points_3d flat array containing 3 * pnt_cnt coordinate values
- * @param[in] pnt_cnt   number of points in points_3d
- *
- * @return 0 on success and -1 for invalid input or when no box can be fit
+ * @param[out] pnts     eight points of oriented bounding box
+ * @param points_3d	array of 3D points
+ * @param pnt_cnt	number of points in pnts array
  */
 BG_EXPORT extern int bg_3d_obb(point_t **pnts,
 			       const fastf_t *points_3d,

include/bg/pca.h

@@ -21,7 +21,7 @@
 /*----------------------------------------------------------------------*/
 /** @addtogroup bg_pca
  *
- * Principal Component Analysis
+ * Principle Component Analysis
  *
  * Calculates an XYZ coordinate system such that it aligns with the largest
  * variations in the supplied data.  Intuitively, it "aligns" with the
@@ -50,7 +50,7 @@
  * point_t p;
  * for (size_t i = 0; i < pntcnt; i++) {
  *     MAT4X3PNT(p, RT, pnts[i]);
- *     VMOVE(pnts[i], p);
+ *     VMOVE(pnts[i], v);
  * }
  * @endcode
  */
@@ -69,7 +69,7 @@ __BEGIN_DECLS
 
 /**
  * @brief
- * Perform a Principal Component Analysis on a set of points.
+ * Perform a Principle Component Analysis on a set of points.
  *
  * Outputs are a center point and XYZ vectors for the coordinate system.
  *
@@ -88,7 +88,7 @@ BG_EXPORT extern int bg_pca(point_t *c, vect_t *xaxis, vect_t *yaxis, vect_t *za
 
 __END_DECLS
 
-#endif  /* BG_PCA_H */
+#endif  /* BG_PLANE_H */
 /** @} */
 /*
  * Local Variables:

include/bg/plane.h

@@ -1171,10 +1171,9 @@ BG_EXPORT extern int bg_isect_planes(point_t pt,
 
 /**
  * @brief
- * Given an origin and a normal, create a plane_t with a unit normal.
- * The input point and normal are not modified.  A zero-length normal fails.
+ * Given an origin and a normal, create a plane_t.
  */
-BG_EXPORT extern int bg_plane_pt_nrml(plane_t *p, const point_t pt, const vect_t nrml);
+BG_EXPORT extern int bg_plane_pt_nrml(plane_t *p, point_t pt, vect_t nrml);
 
 /**
  * @brief
@@ -1189,15 +1188,13 @@ BG_EXPORT extern int bg_fit_plane(point_t *c, vect_t *n, size_t npnts, point_t *
 
 /**
  * @brief
- * Find the closest U,V point on the plane p to 3d point pt.  The plane normal
- * need not be unit length.
+ * Find the closest U,V point on the plane p to 3d point pt.
  */
 BG_EXPORT extern int bg_plane_closest_pt(fastf_t *u, fastf_t *v, plane_t *p, point_t *pt);
 
 /**
  * @brief
- * Return the 3D point on the plane at parametric coordinates u, v.  The plane
- * normal need not be unit length.
+ * Return the 3D point on the plane at parametric coordinates u, v.
  */
 BG_EXPORT extern int bg_plane_pt_at(point_t *pt, plane_t *p, fastf_t u, fastf_t v);
 

include/bg/polygon.h

@@ -293,7 +293,6 @@ BG_EXPORT extern int bg_3d_polygon_area(fastf_t *area, size_t npts, const point_
  * @param[in] npts Number of point_ts, stored in pts
  * @param[in] pts all points of the polygon, sorted counter-clockwise.
  * The array mustn't contain duplicated points or non-coplanar points.
- * A zero-area polygon fails.
  *
  * @return 0 if calculation was successful
  * @return 1 if calculation failed, e.g. because one in-parameter is a NULL-pointer
@@ -308,7 +307,7 @@ BG_EXPORT extern int bg_3d_polygon_centroid(point_t *cent, size_t npts, const po
  *@param[in] npts Number of points, pts contains
  *@param pts Array of point_ts, building a convex polygon. Duplicated points
  *aren't allowed. The points in the array will be sorted counter-clockwise.
- *@param[in] cmp Plane equation of the polygon.  Its normal must be nonzero.
+ *@param[in] cmp Plane equation of the polygon
  *
  *@return 0 if calculation was successful
  *@return 1 if calculation failed, e.g. because pts is a NULL-pointer

include/bg/sat.h

@@ -39,9 +39,7 @@
 __BEGIN_DECLS
 
 /**
- * Test for an intersection between an infinite line and an Axis-Aligned
- * Bounding Box (AABB).  ldir need not be unit length, but must be nonzero.
- * aabb_extent supplies the positive half extents along the coordinate axes.
+ * Test for an intersection between a line and an Axis-Aligned Bounding Box (AABB).
  *
  * Returns 1 if they intersect, 0 otherwise.
  */
@@ -52,10 +50,7 @@ bg_sat_line_aabb(
 	);
 
 /**
- * Test for an intersection between an infinite line and an Oriented Bounding
- * Box (OBB).  ldir need not be unit length, but must be nonzero.  Each OBB
- * extent vector points from the center to a face center; nonzero vectors must
- * be mutually perpendicular, and a zero vector represents a degenerate axis.
+ * Test for an intersection between a line and an Oriented Bounding Box (OBB).
  *
  * Returns 1 if they intersect, 0 otherwise.
  */
@@ -66,9 +61,7 @@ bg_sat_line_obb(
 	);
 
 /**
- * Test for an intersection between a triangle and an Axis-Aligned Bounding
- * Box (AABB).  extent supplies the positive half extents along the coordinate
- * axes.
+ * Test for an intersection between a triangle and an Axis-Aligned Bounding Box (AABB).
  *
  * Returns 1 if they intersect, 0 otherwise.
  */
@@ -79,10 +72,7 @@ bg_sat_tri_aabb(
 	);
 
 /**
- * Test for an intersection between a triangle and an Oriented Bounding Box
- * (OBB).  Each extent vector points from the center to a face center; nonzero
- * vectors must be mutually perpendicular, and a zero vector represents a
- * degenerate axis.
+ * Test for an intersection between a triangle and an Oriented Bounding Box (OBB).
  *
  * Returns 1 if they intersect, 0 otherwise.
  */
@@ -95,9 +85,8 @@ bg_sat_tri_obb(
 /**
  * Test for an intersection between an Axis-Aligned Bounding Box (AABB) and an
  * Oriented Bounding Box (OBB). The latter is defined by a center point and
- * three extent vectors from the center to the centers of the various faces.
- * Nonzero extent vectors must be mutually perpendicular; a zero vector
- * represents a degenerate axis.
+ * three perpendicular vectors from the center to the centers of the various
+ * faces.
  *
  * Returns 1 if they intersect, 0 otherwise.
  */
@@ -109,9 +98,8 @@ bg_sat_aabb_obb(
 
 /**
  * Test for an intersection between two Oriented Bounding Boxes (OBBs). The
- * boxes are defined by a center point and three extent vectors from the center
- * to the centers of the various faces.  Nonzero extent vectors must be
- * mutually perpendicular; a zero vector represents a degenerate axis.
+ * boxes are defined by a center point and three perpendicular vectors from the
+ * center to the centers of the various faces.
  *
  * Returns 1 if they intersect, 0 otherwise.
  */

include/bg/spsr.h

@@ -67,7 +67,6 @@
 #define BG_SPSR_H
 
 #include "common.h"
-#include <stdint.h>
 #include "vmath.h"
 #include "bg/defines.h"
 
@@ -110,7 +109,7 @@ struct bg_3d_spsr_opts {
 #define BG_3D_SPSR_BOUNDARY_DIRICHLET        3
 
 #define BG_3D_SPSR_DEFAULT_DEGREE            1     /* DEFAULT_FEM_DEGREE */
-#define BG_3D_SPSR_DEFAULT_DEPTH             8
+#define BG_3D_SPSR_DEFAULT_DEPTH             11
 #define BG_3D_SPSR_DEFAULT_KERNELDEPTH       0
 #define BG_3D_SPSR_DEFAULT_ITERATIONS        0
 #define BG_3D_SPSR_DEFAULT_FULL_DEPTH        0
@@ -118,13 +117,13 @@ struct bg_3d_spsr_opts {
 #define BG_3D_SPSR_DEFAULT_BASEVCYCLES       0
 #define BG_3D_SPSR_DEFAULT_MAX_MEM           0
 #define BG_3D_SPSR_DEFAULT_THREADS           0
-#define BG_3D_SPSR_DEFAULT_SAMPLES_PER_NODE  1.5
-#define BG_3D_SPSR_DEFAULT_SCALE             1.1
+#define BG_3D_SPSR_DEFAULT_SAMPLES_PER_NODE  1.1
+#define BG_3D_SPSR_DEFAULT_SCALE             1.0
 #define BG_3D_SPSR_DEFAULT_WIDTH             0.0
 #define BG_3D_SPSR_DEFAULT_CONFIDENCE        0.0
 #define BG_3D_SPSR_DEFAULT_CONFIDENCE_BIAS   0.0
 #define BG_3D_SPSR_DEFAULT_CGSOLVER_ACCURACY 1.0e-3
-#define BG_3D_SPSR_DEFAULT_POINT_WEIGHT      2.0
+#define BG_3D_SPSR_DEFAULT_POINT_WEIGHT      8.0    /* DefaultPointWeightMultiplier * Degree */
 #define BG_3D_SPSR_DEFAULT_NONMANIFOLD       0
 #define BG_3D_SPSR_DEFAULT_LINEARFIT         0
 #define BG_3D_SPSR_DEFAULT_EXACT             1
@@ -168,131 +167,10 @@ struct bg_3d_spsr_opts {
  * @return 0 if successful, else error
  *
  */
-DEPRECATED BG_EXPORT int bg_3d_spsr(int **faces, int *num_faces, point_t **vertices, int *num_vertices,
+BG_EXPORT int bg_3d_spsr(int **faces, int *num_faces, point_t **vertices, int *num_vertices,
 			 const point_t *input_points_3d, const vect_t *input_normals_3d,
 			 int num_input_pnts, struct bg_3d_spsr_opts *opts);
 
-/** Reasons the SPSR solver may request additional source samples. */
-#define BG_3D_SPSR_REFINE_RESIDUAL          0x01u
-#define BG_3D_SPSR_REFINE_NORMAL            0x02u
-#define BG_3D_SPSR_REFINE_DENSITY           0x04u
-#define BG_3D_SPSR_REFINE_SURFACE_VARIATION 0x08u
-
-/** Adaptive SPSR termination states. */
-#define BG_3D_SPSR_COMPLETE                 0
-#define BG_3D_SPSR_CALLBACK_STOP            1
-#define BG_3D_SPSR_POINT_LIMIT              2
-#define BG_3D_SPSR_TIME_LIMIT               3
-#define BG_3D_SPSR_NO_NEW_SAMPLES           4
-#define BG_3D_SPSR_SOLVER_ERROR              5
-
-/**
- * One oriented source sample.  A positive @p thickness records the length of
- * the solid ray partition that supplied the sample.  Non-zero equal pair IDs
- * identify the entry and exit samples of the same partition.
- */
-struct bg_3d_spsr_sample {
-    point_t point;
-    vect_t normal;
-    fastf_t thickness;
-    uint64_t pair_id;
-};
-
-/** A spatial location at which the solver would benefit from more samples. */
-struct bg_3d_spsr_refinement_hint {
-    point_t point;
-    vect_t normal;
-    fastf_t error;
-    unsigned int reasons;
-};
-
-/** Source-to-candidate measurements supplied by an adaptive callback. */
-struct bg_3d_spsr_validation {
-    size_t ray_count;
-    size_t source_hit_rays;
-    size_t partition_mismatch_rays;
-    size_t thin_partition_count;
-    size_t thin_partition_mismatches;
-    fastf_t endpoint_error_p95;
-    fastf_t endpoint_error_p99;
-    fastf_t chord_error_fraction;
-    fastf_t surface_area_error_fraction;
-    fastf_t volume_error_fraction;
-    int passed;
-};
-
-/** Information passed to an adaptive source-sampling callback. */
-struct bg_3d_spsr_refinement_request {
-    size_t pass;
-    fastf_t target_feature_size;
-    double remaining_time;
-    const point_t *vertices;
-    size_t vertex_count;
-    const int *faces;
-    size_t face_count;
-    const struct bg_3d_spsr_refinement_hint *hints;
-    size_t hint_count;
-};
-
-/**
- * Callback response.  Sample storage remains owned by the callback and need
- * only remain valid until the callback returns; libbg copies it immediately.
- */
-struct bg_3d_spsr_refinement_response {
-    const struct bg_3d_spsr_sample *samples;
-    size_t sample_count;
-    struct bg_3d_spsr_validation validation;
-    /** Non-zero asks the adaptive driver to stop after this candidate. */
-    int stop_refinement;
-};
-
-typedef int (*bg_3d_spsr_refinement_func_t)(
-    struct bg_3d_spsr_refinement_response *response,
-    const struct bg_3d_spsr_refinement_request *request,
-    void *client_data);
-
-/** Options controlling a bounded adaptive SPSR reconstruction. */
-struct bg_3d_spsr_adaptive_opts {
-    struct bg_3d_spsr_opts solver;
-    size_t max_refinement_passes;
-    size_t max_points;
-    fastf_t target_feature_size;
-    double max_time;
-};
-
-/** Last-resort solve budget; callbacks may stop earlier when progress stalls. */
-#define BG_3D_SPSR_DEFAULT_REFINEMENT_PASSES 8
-#define BG_3D_SPSR_ADAPTIVE_OPTS_DEFAULT { \
-	BG_3D_SPSR_OPTS_DEFAULT, \
-	BG_3D_SPSR_DEFAULT_REFINEMENT_PASSES, \
-	0, \
-	0.0, \
-	0.0 }
-
-/** Summary of an adaptive SPSR attempt. */
-struct bg_3d_spsr_report {
-    size_t solve_count;
-    size_t initial_sample_count;
-    size_t final_sample_count;
-    size_t requested_sample_count;
-    size_t accepted_sample_count;
-    double elapsed_time;
-    int termination;
-    struct bg_3d_spsr_validation validation;
-};
-
-/**
- * Reconstruct a surface and, when requested, obtain targeted additional
- * samples from @p refine.  A callback is required when refinement passes are
- * enabled.  The final candidate must be marked as passed by the callback.
- */
-BG_EXPORT int bg_3d_spsr_adaptive(int **faces, int *num_faces,
-        point_t **vertices, int *num_vertices,
-        const struct bg_3d_spsr_sample *input_samples, size_t sample_count,
-        const struct bg_3d_spsr_adaptive_opts *opts,
-        bg_3d_spsr_refinement_func_t refine, void *client_data,
-        struct bg_3d_spsr_report *report);
-
 __END_DECLS
 
 #endif  /* BG_SPSR_H */

include/bg/tri_pt.h

@@ -52,7 +52,7 @@ BG_EXPORT extern double bg_tri_closest_pt(
 
 __END_DECLS
 
-#endif  /* BG_TRI_PT_H */
+#endif  /* BG_TRI_TRI_H */
 /** @} */
 /*
  * Local Variables:

include/bg/tri_ray.h

@@ -54,7 +54,7 @@ BG_EXPORT extern int bg_isect_tri_line(point_t *isect,
  * If isect is non-NULL, it is written only when an intersection is found.
  */
 BG_EXPORT extern int bg_isect_tri_ray(const point_t orig,
-				      const vect_t dir,
+				      const point_t dir,
 				      const point_t vert0,
 				      const point_t vert1,
 				      const point_t vert2,

include/bg/trimesh.h

@@ -87,11 +87,7 @@ BG_EXPORT extern int bg_trimesh_oriented(int vcnt, int fcnt, fastf_t *v, int *f)
 /**
  * Check if a mesh is topologically solid. Returns 1 if the mesh is NOT SOLID
  * and 0 if the mesh is SOLID. A SOLID (0) outcome indicates the mesh satisfies
- * all three criteria:  Closed, Manifold, Oriented.  If @p bedges is non-NULL,
- * it is set to NULL when no edge diagnostics are present, or to a bu_calloc()
- * allocated array of vertex-index pairs for unmatched, misoriented, and excess
- * edges.  The caller must bu_free() a non-NULL result.  Use
- * bg_trimesh_solid2() when structured diagnostics and their counts are needed.
+ * all three criteria:  Closed, Manifold, Oriented
  */
 BG_EXPORT extern int bg_trimesh_solid(int vcnt, int fcnt, fastf_t *v, int *f, int **bedges);
 

include/bn/mat.h

@@ -214,8 +214,8 @@ BN_EXPORT extern void bn_mat_ae(mat_t m,
 				double elev);
 
 /**
- * Find the azimuth and elevation angles, in degrees, that correspond
- * to the direction (not including twist) given by a direction vector.
+ * Find the azimuth and elevation angles that correspond to the
+ * direction (not including twist) given by a direction vector.
  */
 BN_EXPORT extern void bn_ae_vec(fastf_t *azp,
 				fastf_t *elp,
@@ -236,16 +236,14 @@ BN_EXPORT extern void bn_aet_vec(fastf_t *az,
 				 fastf_t accuracy);
 
 /**
- * Find a unit vector from the origin given azimuth and elevation in
- * radians.
+ * Find a unit vector from the origin given azimuth and elevation.
  */
 BN_EXPORT extern void bn_vec_ae(vect_t vec,
 				fastf_t az,
 				fastf_t el);
 
 /**
- * Find a vector from the origin given azimuth and elevation in radians,
- * and distance.
+ * Find a vector from the origin given azimuth, elevation, and distance.
  */
 BN_EXPORT extern void bn_vec_aed(vect_t vec,
 				 fastf_t az,
@@ -312,10 +310,8 @@ BN_EXPORT extern void bn_vec_perp(vect_t new_vec,
  * space by the angle between the two.  There are many candidate
  * matrices.
  *
- * The input 'from' and 'to' vectors need not be unit length.  They are
- * normalized internally.  Near-parallel inputs are handled using the
- * angular tolerance in 'tol'.  If either input cannot be normalized,
- * the result is the identity matrix.
+ * The input 'from' and 'to' vectors need not be unit length.
+ * MAT4X3VEC(to, m, from) is the identity that is created.
  *
  */
 BN_EXPORT extern void bn_mat_fromto(mat_t m,
@@ -348,11 +344,12 @@ BN_EXPORT extern void bn_mat_zrot(mat_t m,
 				  double cosz);
 
 /**
- * Given a nonzero direction vector D, produce a matrix which
- * rotates D onto the -Z axis.  D is normalized internally.  Z-axis
- * directions use a deterministic roll, and a zero direction produces
- * the identity matrix.  This matrix will be suitable for use as a
- * model2view matrix.
+ * Given a direction vector D of unit length, product a matrix which
+ * rotates that vector D onto the -Z axis.  This matrix will be
+ * suitable for use as a "model2view" matrix.
+ *
+ * XXX This routine will fail if the vector is already more or less
+ * aligned with the Z axis.
  *
  * This is done in several steps.
  *

include/bn/multipoly.h

@@ -37,42 +37,11 @@ __BEGIN_DECLS
 
 typedef struct bn_multipoly {
     uint32_t magic;
-    /**< Number of coefficients in the first polynomial variable. */
     int dgrs;
-    /**< Number of coefficients in the second polynomial variable. */
     int dgrt;
-    /**< Coefficient of S^s * T^t is cf[s][t]. */
     double **cf;
 }  bn_multipoly_t;
 
-/**
- * Allocate a zero-valued bivariate polynomial with the specified
- * coefficient counts.  Both counts must be positive.
- */
-BN_EXPORT extern struct bn_multipoly *bn_multipoly_new(int dgrs, int dgrt);
-
-/** Release a bivariate polynomial allocated by bn_multipoly_new(). */
-BN_EXPORT extern void bn_multipoly_free(struct bn_multipoly *poly);
-
-/** Grow a polynomial to at least the specified coefficient counts. */
-BN_EXPORT extern struct bn_multipoly *bn_multipoly_grow(struct bn_multipoly *poly,
-							 int dgrs,
-							 int dgrt);
-
-/** Set the coefficient of S^s * T^t, growing the polynomial as needed. */
-BN_EXPORT extern struct bn_multipoly *bn_multipoly_set(struct bn_multipoly *poly,
-							int s,
-							int t,
-							double value);
-
-/** Return the sum of two bivariate polynomials. */
-BN_EXPORT extern struct bn_multipoly *bn_multipoly_add(const struct bn_multipoly *p1,
-							const struct bn_multipoly *p2);
-
-/** Return the product of two bivariate polynomials. */
-BN_EXPORT extern struct bn_multipoly *bn_multipoly_mul(const struct bn_multipoly *p1,
-							const struct bn_multipoly *p2);
-
 __END_DECLS
 
 #endif  /* BN_MULTIPOLY_H */

include/bn/poly.h

@@ -45,10 +45,7 @@ __BEGIN_DECLS
 
 /**
  * Polynomial data type
- *
- * Coefficients are stored from the highest to the lowest degree:
- * `cf[0]` is the coefficient of `X^dgr` and `cf[dgr]` is the
- * constant coefficient.
+ * bn_poly->cf[n] corresponds to X^n
  */
 typedef struct bn_poly {
     uint32_t magic;
@@ -66,9 +63,6 @@ typedef struct bn_poly {
  *
  * @brief
  * multiply two polynomials
- *
- * The output may alias either input.  If the result exceeds
- * BN_MAX_POLY_DEGREE, returns BN_POLY_NULL without modifying product.
  */
 BN_EXPORT extern struct bn_poly *bn_poly_mul(struct bn_poly *product,
 					     const struct bn_poly *m1,
@@ -87,8 +81,6 @@ BN_EXPORT extern struct bn_poly *bn_poly_scale(struct bn_poly *eqn,
  * bn_poly_add
  * @brief
  * add two polynomials
- *
- * The output may alias either input.
  */
 BN_EXPORT extern struct bn_poly *bn_poly_add(struct bn_poly *sum,
 					     const struct bn_poly *poly1,
@@ -98,8 +90,6 @@ BN_EXPORT extern struct bn_poly *bn_poly_add(struct bn_poly *sum,
  * bn_poly_sub
  * @brief
  * subtract two polynomials
- *
- * The output may alias either input.
  */
 BN_EXPORT extern struct bn_poly *bn_poly_sub(struct bn_poly *diff,
 					     const struct bn_poly *poly1,
@@ -108,10 +98,7 @@ BN_EXPORT extern struct bn_poly *bn_poly_sub(struct bn_poly *diff,
 /**
  * @brief
  * Divides any polynomial into any other polynomial using synthetic
- * division.  Both polynomials must have real coefficients and the
- * divisor's leading coefficient must be nonzero.  The quotient and
- * remainder outputs must be distinct; the quotient may alias the
- * dividend.  All other output/input pairs must be distinct.
+ * division.  Both polynomials must have real coefficients.
  */
 BN_EXPORT extern void bn_poly_synthetic_division(struct bn_poly *quo,
 						 struct bn_poly *rem,
@@ -121,8 +108,7 @@ BN_EXPORT extern void bn_poly_synthetic_division(struct bn_poly *quo,
 /**
  *@brief
  * Uses the quadratic formula to find the roots (in `complex' form) of
- * any quadratic equation with real coefficients.  A nonzero leading
- * coefficient need not be one.
+ * any quadratic equation with real coefficients.
  *
  *	@return 1 for success
  *	@return 0 for fail.
@@ -158,7 +144,6 @@ BN_EXPORT extern int bn_poly_quadratic_roots(struct bn_complex roots[],
  * If D < 0, there will be three unequal real roots.
  *
  * Returns 1 for success, 0 for fail.
- * A nonzero leading coefficient need not be one.
  */
 BN_EXPORT extern int bn_poly_cubic_roots(struct bn_complex roots[],
 					 const struct bn_poly *eqn);
@@ -166,8 +151,7 @@ BN_EXPORT extern int bn_poly_cubic_roots(struct bn_complex roots[],
 /**
  *@brief
  * Uses the quartic formula to find the roots (in `complex' form)
- * of any quartic equation with real coefficients.  A nonzero leading
- * coefficient need not be one.
+ * of any quartic equation with real coefficients.
  *
  *	@return 1 for success
  *	@return 0 for fail.

include/bn/qmath.h

@@ -100,8 +100,8 @@ BN_EXPORT extern void quat_double(quat_t qout,
 /**
  *@brief Gives the bisector of quaternions q1 and q2.
  *
- * Inputs must be unit quaternions.  Antipodal inputs use a valid
- * great-circle midpoint.
+ * (Could be done with quat_slerp and factor 0.5)
+ * [I believe they must be unit quaternions this to work]
  */
 BN_EXPORT extern void quat_bisect(quat_t qout,
 				  const quat_t q1,
@@ -161,8 +161,7 @@ BN_EXPORT extern void quat_exp(quat_t out,
 
 /**
  *@brief
- *  Take the natural logarithm of a unit quaternion.  The identity
- *  quaternion maps to the zero quaternion.
+ *  Take the natural logarithm of a unit quaternion.
  *  Code by Ken Shoemake.
  */
 BN_EXPORT extern void quat_log(quat_t out,

include/bu/log.h

@@ -205,30 +205,23 @@ BU_EXPORT extern int bu_sscanf(const char *src, const char *fmt, ...) _BU_ATTR_S
 /** Routines for scanning certain kinds of data. */
 
 /**
- * Scans a sequence of fastf_t numbers from a string or stdin.
+ * Scans a sequence of fastf_t numbers from a string or stdin
  *
  * Scanning fastf_t numbers with bu_sscanf() is difficult, because
  * doing so requires scanning to some intermediate type like double
  * and then assigning to the fastf_t variable to convert the value to
  * whatever type fastf_t really is.  This function makes it possible
- * to scan a series of fastf_t numbers easily, by doing the required
- * conversion internally.  Delimiters form a non-empty character set:
- * one or more delimiter characters must separate values, and repeated
- * delimiters are skipped.  Empty scan fields are therefore not
- * supported (e.g., "0.0,,0.0,1.0" scans as 3 fields, not 4 with the
- * second skipped).  A trailing delimiter is not consumed after the
- * final requested value.
+ * to scan a series of fastf_t numbers separated by some character(s)
+ * easily, by doing the required conversion internally to the
+ * functions.  As series of delimiter characters will be skipped,
+ * empty scan fields are not supported (e.g., "0.0,,0.0,1.0" will scan
+ * as 3 fields, not 4 with the 2nd skipped).
  *
- * @param[out] c Optional destination for the number of input characters
- *     consumed.  It is set to zero before scanning and remains zero for
- *     invalid arguments or when no value is scanned.
+ * @param[out] c Returns number of characters scanned by the function
  * @param[in] src A source string to scan from, or NULL to read from stdin
- * @param[in] delim Non-empty set of delimiter characters to skip between
- *     scan values
+ * @param[in] delim Any delimiter character(s) to skip between scan values
  * @param[in] n Number of fastf_t values to scan from the src input string
- * @param[out] ... Exactly n fastf_t pointers for storing scanned values;
- *     individual pointers may be NULL
- * @return Number of values scanned
+ * @param[out] ... Pointers to fastf_t for storing scanned values (optional)
  *
  */
 BU_EXPORT extern int bu_scan_fastf_t(int *c, const char *src, const char *delim, size_t n, ...);

include/bu/parallel.h

@@ -103,23 +103,10 @@ BU_EXPORT extern void bu_nice_set(int newnice);
 
 /**
  * Return the maximum number of physical CPUs that are considered to
- * be available to this process now, subject to any process-local limit set by
- * bu_avail_cpus_set().
+ * be available to this process now.
  */
 BU_EXPORT extern size_t bu_avail_cpus(void);
 
-/**
- * Limit the CPU count reported by bu_avail_cpus() and honored by subsequent
- * bu_parallel() calls in this process.  A value of zero removes the limit.
- * Call this before starting worker threads; changing the limit does not resize
- * work that is already running.
- *
- * This is useful when several independent processes share one machine and
- * each process contains code that selects its own parallelism with
- * bu_avail_cpus().
- */
-BU_EXPORT extern void bu_avail_cpus_set(size_t ncpu);
-
 
 /**
  * Create parallel threads of execution.

include/conf/MINOR

@@ -1 +1 @@
-47
+45

include/rt/debug.h

@@ -84,7 +84,7 @@ __BEGIN_DECLS
 #define RT_DEBUG_HF		0x04000000	/**< @brief 27 Height Field solids */
 
 #define RT_DEBUG_MESHING	0x08000000	/**< @brief 28 Print meshing/triangulation details */
-#define RT_DEBUG_RTCMP		0x10000000      /**< @brief 29 Set up rtcmp capture; recording uses a runtime toggle.  If LIBRT_RTCMP_PRIMITIVES=1 includes primitive segments. */
+#define RT_DEBUG_UNUSED_3	0x10000000	/**< @brief 29 Unassigned */
 #define RT_DEBUG_UNUSED_4	0x20000000	/**< @brief 30 Unassigned */
 
 /* Options which will cause the library to write binary debugging output */
@@ -96,7 +96,7 @@ __BEGIN_DECLS
     "\040PL_BOX" \
     "\037PL_SOLIDS" \
     "\036UNUSED_4" \
-    "\035RTCMP" \
+    "\035UNUSED_3" \
     "\034UNUSED_2" \
     "\033HF" \
     "\032EBM" \

include/rt/edit.h

@@ -465,17 +465,6 @@ rt_get_solid_keypoint(struct rt_edit *s, point_t *pt, const char **strp, fastf_t
 RT_EXPORT extern void
 rt_edit_set_edflag(struct rt_edit *s, int edflag);
 
-/**
- * Set a parameter-edit translation target expressed in model base units.
- *
- * rt_edit's parameter-edit handlers consume e_para in database local units,
- * while keypoints and edit axes are maintained in model base units.  This
- * helper performs the boundary conversion for callers that have computed a
- * target from a keypoint, axis, or view transform.
- */
-RT_EXPORT extern void
-rt_edit_set_translation_target(struct rt_edit *s, const point_t target);
-
 RT_EXPORT extern int
 rt_edit_knob_cmd_process(
 	struct rt_edit *s,
@@ -510,10 +499,8 @@ rt_knob_edit_sca(
 	int matrix_edit);
 
 /* Equivalent to sedit - run editing logic after input data is set in
- * rt_edit container.  Transient input flags are consumed whether the edit
- * succeeds or fails; parameter values remain available to the caller.
- * Returns BRLCAD_ERROR if the primitive edit fails. */
-RT_EXPORT extern int
+ * rt_edit container */
+RT_EXPORT extern void
 rt_edit_process(struct rt_edit *s);
 
 /**
@@ -534,7 +521,7 @@ rt_edit_snap_point(point2d_t pt, const struct rt_edit *s);
  * restored later with rt_edit_revert().
  *
  * The snapshot is stored inside the rt_edit struct.  Calling this
- * function again replaces the previous snapshot only on success.
+ * function again overwrites any previous snapshot (single-level undo).
  *
  * @return BRLCAD_OK on success, BRLCAD_ERROR if the export failed.
  */
@@ -545,8 +532,8 @@ rt_edit_checkpoint(struct rt_edit *s);
  * Restore primitive parameters from the snapshot saved by
  * rt_edit_checkpoint().
  *
- * The snapshot remains available for repeated restores.  If none has
- * been saved, this function logs a message and returns BRLCAD_ERROR.
+ * If no snapshot has been saved (or the last snapshot was already
+ * consumed) this function logs a message and returns BRLCAD_ERROR.
  *
  * @return BRLCAD_OK on success, BRLCAD_ERROR otherwise.
  */
@@ -710,11 +697,7 @@ enum rt_constraint_edit_op_kind {
     RT_CONSTRAINT_EDIT_OP_SET_BEND,
     RT_CONSTRAINT_EDIT_OP_SCALE_OD,
     RT_CONSTRAINT_EDIT_OP_SCALE_ID,
-    RT_CONSTRAINT_EDIT_OP_SCALE_BEND,
-    /* Set every point when an all-zero dimension has no scale reference. */
-    RT_CONSTRAINT_EDIT_OP_SET_ALL_OD,
-    RT_CONSTRAINT_EDIT_OP_SET_ALL_ID,
-    RT_CONSTRAINT_EDIT_OP_SET_ALL_BEND
+    RT_CONSTRAINT_EDIT_OP_SCALE_BEND
 };
 
 struct rt_constraint_edit_param_ref {

include/rt/primitives/annot.h

@@ -29,10 +29,19 @@
 #include "bu/list.h"
 #include "bu/vls.h"
 #include "bn/tol.h"
+#include "rt/geom.h"
 #include "rt/defines.h"
 
 __BEGIN_DECLS
 
+struct db_i;
+struct xray;
+
+#define RT_ANNOT_SCREEN_DPI 96.0
+#define RT_ANNOT_MILLIMETERS_PER_INCH 25.4
+#define RT_ANNOT_SCREEN_PIXELS_PER_MM \
+    (RT_ANNOT_SCREEN_DPI / RT_ANNOT_MILLIMETERS_PER_INCH)
+
 RT_EXPORT extern struct rt_annot_internal *rt_copy_annot(const struct rt_annot_internal *annot_ip);
 
 /** Validate annotation topology, model-space placement, and optional segment
@@ -40,6 +49,82 @@ RT_EXPORT extern struct rt_annot_internal *rt_copy_annot(const struct rt_annot_i
 RT_EXPORT extern int rt_annot_validate(const struct rt_annot_internal *annot_ip,
 	struct bu_vls *messages);
 
+/** Opaque, immutable collection of prepared annotation instances. */
+struct rt_annot_scene;
+
+/** Camera state needed to evaluate view-space annotations.  Pixel sample
+ * coordinates use a lower-left origin. */
+struct rt_annot_view {
+    mat_t model2view;
+    size_t width;
+    size_t height;
+    fastf_t perspective;
+};
+
+/** A zero-thickness annotation coverage event.  path is owned by the scene. */
+struct rt_annot_hit {
+    const char *path;
+    size_t segment;
+    uint32_t role;
+    point_t point;
+    fastf_t distance;
+    unsigned char color[4];
+    int screen_space;
+    int visible;
+};
+
+/** Prepare annotations referenced by the supplied top-level database paths.
+ * Boolean operations are intentionally ignored: annotations are
+ * informational leaves, not solids. */
+RT_EXPORT extern struct rt_annot_scene *rt_annot_scene_create(
+	struct db_i *dbip, int path_count, const char * const *paths,
+	const struct bg_tess_tol *ttol, const struct bn_tol *tol);
+
+RT_EXPORT extern void rt_annot_scene_destroy(struct rt_annot_scene *scene);
+RT_EXPORT extern size_t rt_annot_scene_count(const struct rt_annot_scene *scene);
+/** Return model-space bounds for prepared annotations.  Screen-space geometry
+ * contributes its model-space anchor because its extent is view-dependent. */
+RT_EXPORT extern int rt_annot_scene_bounds(const struct rt_annot_scene *scene,
+	point_t minimum, point_t maximum);
+
+/** Return the top visible coverage event for a primary sample. */
+RT_EXPORT extern int rt_annot_scene_query(const struct rt_annot_scene *scene,
+	const struct rt_annot_view *view, const struct xray *ray,
+	fastf_t sample_x, fastf_t sample_y, fastf_t scene_distance,
+	struct rt_annot_hit *hit);
+
+/** Return all visible primary-sample coverage events in back-to-front order.
+ * Model-space events precede screen-space events.  When hits is NULL or
+ * capacity is zero, only the required result count is returned. */
+RT_EXPORT extern size_t rt_annot_scene_query_layers(
+	const struct rt_annot_scene *scene, const struct rt_annot_view *view,
+	const struct xray *ray, fastf_t sample_x, fastf_t sample_y,
+	fastf_t scene_distance, struct rt_annot_hit *hits, size_t capacity);
+
+/** Return all model-space coverage events, sorted by ray distance.  When hits
+ * is NULL or capacity is zero, only the required result count is returned. */
+RT_EXPORT extern size_t rt_annot_scene_query_model(const struct rt_annot_scene *scene,
+	const struct xray *ray, fastf_t scene_distance,
+	struct rt_annot_hit *hits, size_t capacity);
+
+/** Alpha-composite a coverage result over a linear RGB sample. */
+RT_EXPORT extern void rt_annot_hit_blend(fastf_t color[3],
+	const struct rt_annot_hit *hit);
+
+/** Alpha-composite all visible primary-sample coverage over color.  If
+ * front_hit is non-NULL, it receives the frontmost event. */
+RT_EXPORT extern int rt_annot_scene_composite(const struct rt_annot_scene *scene,
+	const struct rt_annot_view *view, const struct xray *ray,
+	fastf_t sample_x, fastf_t sample_y, fastf_t scene_distance,
+	fastf_t color[3], struct rt_annot_hit *front_hit);
+
+/** Generate one annotation segment in annotation-local XY coordinates.  This
+ * is shared by plotting and the coverage preparer to keep curve and font
+ * expansion identical. */
+RT_EXPORT extern int rt_annot_segment_vlist(struct bu_list *vlfree,
+	struct bu_list *vhead, const struct bg_tess_tol *ttol,
+	const struct rt_annot_internal *annot_ip, size_t segment);
+
 __END_DECLS
 
 /** @} */

include/rt/primitives/bot.h

@@ -55,11 +55,14 @@ struct clt_tri_specific {
 };
 #endif
 
+#define RT_BOT_SURFNO_INTERIOR (-1)
+
 /* Shared between bot and ars at the moment */
 struct bot_specific {
     unsigned char bot_mode;
     unsigned char bot_orientation;
     unsigned char bot_flags;
+    unsigned char bot_parity_shot;
     size_t bot_ntri;
     fastf_t *bot_thickness;
     struct bu_bitv *bot_facemode;

include/rt/primitives/nmg.h

@@ -41,7 +41,6 @@ struct rt_nmg_edit {
     struct loopuse *lu_copy;    /* copy of loop to be extruded */
     plane_t lu_pl;              /* plane equation for loop to be extruded */
     struct shell *es_s;         /* Shell where extrusion is to end up */
-    int es_s_from_extrusion;    /* Replace this edit's prior preview */
     point_t lu_keypoint;        /* keypoint of lu_copy for extrusion */
     struct vertex *es_v;        /* Currently selected NMG vertex (ECMD_NMG_VPICK) */
     struct faceuse *es_fu;      /* Currently selected NMG faceuse (ECMD_NMG_FPICK) */

include/rt/shoot.h

@@ -55,18 +55,6 @@
 
 __BEGIN_DECLS
 
-#define RT_RTCMP_CAPTURE_RUNTIME_TOGGLE 1
-
-/** Set RT_DEBUG_RTCMP recording: nonzero enables, zero disables.
- * Enabling requires the debug flag and a usable output writer.
- * Return 0 on success and -1 on failure. Disabling succeeds even
- * before setup.
- */
-RT_EXPORT extern int rt_rtcmp_capture_set_enabled(int enabled);
-
-/** Wait until pending RT_DEBUG_RTCMP records have reached the output file. */
-RT_EXPORT extern int rt_rtcmp_capture_flush(void);
-
 /**
  * @brief
  * Shoot a ray

misc/CMake/BRLCAD_EXT_Setup.cmake

@@ -30,7 +30,6 @@
 # SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 
 include(CMakeParseArguments)
-include("${CMAKE_CURRENT_LIST_DIR}/BRLCAD_Repository.cmake")
 
 # Try to locate a usable git executable. On Windows, Visual Studio may
 # provide git in non-standard locations; try common Program Files and
@@ -89,15 +88,6 @@ function(brlcad_bext_init BEXT_SHA1)
     set(BRLCAD_EXT_DIR "${BRLCAD_EXT_DIR_ENV}")
   endif(BRLCAD_EXT_DIR_ENV AND NOT DEFINED BRLCAD_EXT_DIR)
 
-  # Resolve paths before probing them so later subconfigures use the same
-  # directories even when launched from a different working directory.
-  foreach(_brlcad_ext_path_var BRLCAD_EXT_DIR BRLCAD_EXT_SOURCE_DIR)
-    if(DEFINED ${_brlcad_ext_path_var} AND NOT "${${_brlcad_ext_path_var}}" STREQUAL "")
-      get_filename_component(${_brlcad_ext_path_var} "${${_brlcad_ext_path_var}}" ABSOLUTE BASE_DIR "${CMAKE_SOURCE_DIR}")
-      set(${_brlcad_ext_path_var} "${${_brlcad_ext_path_var}}" PARENT_SCOPE)
-    endif()
-  endforeach()
-
   set(_brlcad_ext_dir_specified FALSE)
   set(_brlcad_ext_dir_managed FALSE)
 
@@ -158,14 +148,14 @@ function(brlcad_bext_init BEXT_SHA1)
     elseif(NOT _brlcad_ext_dir_specified)
       message(
 	WARNING
-	"External dependencies will be acquired from ${_BRLCAD_REPOSITORY_ROOT}, configured, and built automatically as-needed. Set BRLCAD_EXT_DIR to specify instead."
+	"External dependencies will be downloaded, configured, and built automatically as-needed. Set BRLCAD_EXT_DIR to specify instead."
 	)
     endif()
     message(
       STATUS
       "BRLCAD_EXT_DIR specifies prebuilt external dependencies directory\n"
       "containing 'install' and 'noinstall' folders, outputs from installing\n"
-      "${_BRLCAD_BEXT_REPOSITORY}.\n"
+      "https://github.com/BRL-CAD/bext.\n"
     )
   endif(NOT BRLCAD_EXT_NOINSTALL_DIR OR NOT BRLCAD_EXT_INSTALL_DIR)
 
@@ -221,24 +211,28 @@ function(brlcad_rel_version BRLCAD_REL)
   set(${BRLCAD_REL} "rel-${VMAJOR}-${VMINOR}-${VPATCH}" PARENT_SCOPE)
 endfunction()
 
-# Branch/tag queries are advisory; an unavailable source cannot validate them.
+# Look up branch or tag SHA1 hashes from primary github.com repo
 function(remote_sha1 SHA1_VAR BRANCH)
-  set(${SHA1_VAR} "" PARENT_SCOPE)
   if(NOT GIT_EXEC)
+    set(${SHA1_VAR} "" PARENT_SCOPE)
     return()
   endif()
-  brlcad_repository_git_command(git_command)
+  # This needs a working internet connection to succeed (and GitHub
+  # must be up and working as well.)
   execute_process(
-    COMMAND "${CMAKE_COMMAND}" -E env GIT_TERMINAL_PROMPT=0 --
-      ${git_command} -c credential.helper= -c credential.interactive=false
-      ls-remote "${_BRLCAD_BEXT_REPOSITORY}" "${BRANCH}"
-    WORKING_DIRECTORY "${CMAKE_SOURCE_DIR}"
-    RESULT_VARIABLE git_result
-    OUTPUT_VARIABLE remote_refs
-    ERROR_QUIET
-  )
-  if(git_result STREQUAL "0" AND remote_refs MATCHES "^([0-9a-f]+)[ \t]")
-    set(${SHA1_VAR} "${CMAKE_MATCH_1}" PARENT_SCOPE)
+    COMMAND ${GIT_EXEC} ls-remote https://github.com/BRL-CAD/bext.git ${BRANCH}
+    WORKING_DIRECTORY ${CMAKE_SOURCE_DIR}
+    RESULT_VARIABLE LS_REMOTE_STATUS
+    OUTPUT_VARIABLE LS_REMOTE_STR
+    ERROR_VARIABLE LS_REMOTE_STR
+    )
+  if(NOT LS_REMOTE_STATUS AND LS_REMOTE_STR)
+    string(REPLACE " " ";" LS_REMOTE_STR "${LS_REMOTE_STR}")
+    string(REPLACE "\t" ";" LS_REMOTE_STR "${LS_REMOTE_STR}")
+    list(GET LS_REMOTE_STR 0 REMOTE_SHA1)
+    set(${SHA1_VAR} "${REMOTE_SHA1}" PARENT_SCOPE)
+  else()
+    set(${SHA1_VAR} "" PARENT_SCOPE)
   endif()
 endfunction()
 
@@ -373,29 +367,24 @@ function(setup_bext_dir)
     # no pre-existing sources, so we're cloning into the build dir
     set(BRLCAD_EXT_SOURCE_DIR "${CMAKE_CURRENT_BINARY_DIR}/bext")
 
-    brlcad_repository_git_command(git_command)
-    message(STATUS "Cloning bext from ${_BRLCAD_BEXT_REPOSITORY}")
-    execute_process(
-      COMMAND ${git_command} clone "${_BRLCAD_BEXT_REPOSITORY}" "${BRLCAD_EXT_SOURCE_DIR}"
-      WORKING_DIRECTORY "${CMAKE_CURRENT_BINARY_DIR}"
-      RESULT_VARIABLE git_result
-    )
-    if(NOT git_result STREQUAL "0")
-      message(FATAL_ERROR "Unable to clone bext from ${_BRLCAD_BEXT_REPOSITORY} (${git_result})")
-    endif()
-    execute_process(
-      COMMAND ${git_command} -c advice.detachedHead=false checkout "${BEXT_SHA1}"
-      WORKING_DIRECTORY "${BRLCAD_EXT_SOURCE_DIR}"
-      RESULT_VARIABLE git_result
-    )
-    if(NOT git_result STREQUAL "0")
-      message(FATAL_ERROR "Unable to check out pinned bext commit ${BEXT_SHA1} from ${_BRLCAD_BEXT_REPOSITORY} (${git_result}); ensure the repository contains this commit")
-    endif()
-
-    # Record the pin only after checkout succeeds.
+    # Write current SHA1 for subsequent processing
     file(WRITE "${CMAKE_CURRENT_BINARY_DIR}/bext.sha1" "${BEXT_SHA1}")
     distclean("${CMAKE_CURRENT_BINARY_DIR}/bext.sha1")
 
+    # Clone
+    message("BRL-CAD bext clone command: ${GIT_EXEC} clone https://github.com/BRL-CAD/bext.git")
+    execute_process(
+      COMMAND ${GIT_EXEC} clone https://github.com/BRL-CAD/bext.git
+      WORKING_DIRECTORY ${CMAKE_CURRENT_BINARY_DIR}
+      COMMAND_ERROR_IS_FATAL ANY
+    )
+    message("BRL-CAD bext checkout command: ${GIT_EXEC} -c advice.detachedHead=false checkout ${BEXT_SHA1}")
+    execute_process(
+      COMMAND ${GIT_EXEC} -c advice.detachedHead=false checkout ${BEXT_SHA1}
+      WORKING_DIRECTORY ${CMAKE_CURRENT_BINARY_DIR}/bext
+      COMMAND_ERROR_IS_FATAL ANY
+    )
+
     # Configure process will need to clean up bext
     if (BRLCAD_BEXT_CLEANUP)
       set(BEXT_SRC_CLEANUP TRUE PARENT_SCOPE)
@@ -436,12 +425,6 @@ function(brlcad_ext_setup)
     message(FATAL_ERROR "Invalid bext directory: ${BRLCAD_EXT_SOURCE_DIR}")
   endif()
 
-  if(NOT _BRLCAD_REPOSITORY_ROOT STREQUAL "${_BRLCAD_REPOSITORY_DEFAULT_ROOT}/"
-      AND NOT EXISTS "${BRLCAD_EXT_SOURCE_DIR}/CMake/BextGit.cmake")
-    message(FATAL_ERROR
-      "bext sources in ${BRLCAD_EXT_SOURCE_DIR} do not support BEXT_REPOSITORY_ROOT; update bext before building with BRLCAD_REPOSITORY_ROOT=${BRLCAD_REPOSITORY_ROOT}")
-  endif()
-
   # We do allow the user to specify a build dir to be reused, however
   # this is unpredictable and may cause unexpected results depending
   # on what is or isn't rebuilt - use with care.
@@ -480,18 +463,11 @@ function(brlcad_ext_setup)
   endif(BRLCAD_ENABLE_APPLESEED)
 
   # Define CMake arguments we'll need regardless of mode
-  set(bext_shallow_clone ON)
-  if(_BRLCAD_REPOSITORY_LOCAL)
-    set(bext_shallow_clone OFF)
-  endif()
   set(CMAKE_CMD_ARGS
-    -S "${BRLCAD_EXT_SOURCE_DIR}" -B "${BRLCAD_EXT_BUILD_DIR}"
-    "-DBEXT_REPOSITORY_ROOT:STRING=${_BRLCAD_REPOSITORY_ROOT}"
-    "-DGIT_SHALLOW_CLONE:BOOL=${bext_shallow_clone}"
-    "-DENABLE_ALL=${BEXT_ENABLE_ALL}" "-DUSE_GDAL=${BEXT_USE_GDAL}"
-    "-DUSE_QT=${BEXT_USE_QT}" "-DUSE_TCL=${BEXT_USE_TCL}"
-    "-DUSE_APPLESEED=${BEXT_USE_APPLESEED}" "-DCMAKE_BUILD_TYPE=${CMAKE_BUILD_TYPE}"
-    "-DCMAKE_INSTALL_PREFIX=${BRLCAD_EXT_INSTALL_DIR}"
+    -S ${BRLCAD_EXT_SOURCE_DIR} -B ${BRLCAD_EXT_BUILD_DIR} -DGIT_SHALLOW_CLONE=ON
+    -DENABLE_ALL=${BEXT_ENABLE_ALL} -DUSE_GDAL=${BEXT_USE_GDAL} -DUSE_QT=${BEXT_USE_QT} -DUSE_TCL=${BEXT_USE_TCL}
+    -DUSE_APPLESEED=${BEXT_USE_APPLESEED} -DCMAKE_BUILD_TYPE=${CMAKE_BUILD_TYPE}
+    -DCMAKE_INSTALL_PREFIX=${BRLCAD_EXT_INSTALL_DIR}
     )
   if(BRLCAD_X11_PROVIDER_RESOLVED STREQUAL "XMIN" AND BRLCAD_XMIN_ROOT)
     list(APPEND CMAKE_CMD_ARGS

misc/CMake/BRLCAD_ExternalDeps.cmake

@@ -1108,9 +1108,9 @@ function(brlcad_install_strclear_replace stamp_dir strclear source install_file
   endif()
 endfunction()
 
-function(brlcad_install_binary_postprocess stamp_dir strclear source install_file install_type mode rpath_tool install_rpath build_lib_path rel_rpath verbose)
+function(brlcad_install_binary_postprocess stamp_dir strclear source install_file install_type mode rpath_tool install_rpath build_lib_path rel_rpath use_selective_rpath verbose)
   set(_brlcad_physical_file "$ENV{DESTDIR}${install_file}")
-  set(_brlcad_signature "binary-postprocess|${install_type}|${mode}|${rpath_tool}|${install_rpath}|${build_lib_path}|${rel_rpath}|${strclear}")
+  set(_brlcad_signature "binary-postprocess|${install_type}|${mode}|${rpath_tool}|${install_rpath}|${build_lib_path}|${rel_rpath}|${use_selective_rpath}|${strclear}")
   _brlcad_postprocess_needed(_brlcad_needed "${stamp_dir}" "${source}" "${_brlcad_physical_file}" "${_brlcad_signature}")
   if(NOT _brlcad_needed)
     return()
@@ -1121,41 +1121,32 @@ function(brlcad_install_binary_postprocess stamp_dir strclear source install_fil
     message(FATAL_ERROR "Post-install copy did not create ${_brlcad_physical_file}")
   endif()
   set(_brlcad_result 0)
-  set(_brlcad_rpath_output)
-  set(_brlcad_rpath_error)
   if("${mode}" STREQUAL "RPATH_TOOL")
-    # Install the exact project-defined search order.  A selective update can
-    # leave an existing relative-only RPATH without the canonical install path.
+    set(_brlcad_selective_rpath_args)
+    if(use_selective_rpath)
+      set(_brlcad_selective_rpath_args --set-rpath-if-needed --set-rpath-if-needed-prepend --stale-rpath-prefix "${build_lib_path}")
+    endif()
     execute_process(
-      COMMAND "${rpath_tool}" --set-rpath "${install_rpath}" "${_brlcad_physical_file}"
+      COMMAND "${rpath_tool}" --set-rpath "${install_rpath}" ${_brlcad_selective_rpath_args} "${_brlcad_physical_file}"
       RESULT_VARIABLE _brlcad_result
-      OUTPUT_VARIABLE _brlcad_rpath_output
-      ERROR_VARIABLE _brlcad_rpath_error
     )
   elseif("${mode}" STREQUAL "APPLE")
     execute_process(
       COMMAND install_name_tool -delete_rpath "${build_lib_path}" "${_brlcad_physical_file}"
       RESULT_VARIABLE _brlcad_result
-      OUTPUT_VARIABLE _brlcad_rpath_output
-      ERROR_VARIABLE _brlcad_rpath_error
+      OUTPUT_VARIABLE _brlcad_output
+      ERROR_VARIABLE _brlcad_error
     )
     if(_brlcad_result EQUAL 0)
       execute_process(
         COMMAND install_name_tool -add_rpath "${rel_rpath}" "${_brlcad_physical_file}"
         RESULT_VARIABLE _brlcad_result
-        OUTPUT_VARIABLE _brlcad_rpath_output
-        ERROR_VARIABLE _brlcad_rpath_error
       )
     endif()
   endif()
 
   if(NOT _brlcad_result EQUAL 0)
-    message(
-      WARNING
-        "Post-install RPATH update failed for ${_brlcad_physical_file}\n"
-        "stdout: ${_brlcad_rpath_output}\n"
-        "stderr: ${_brlcad_rpath_error}"
-    )
+    message(WARNING "Post-install RPATH update failed for ${_brlcad_physical_file}")
     return()
   endif()
 
@@ -1849,9 +1840,13 @@ endfunction()
       elseif(APPLE)
         set(_brlcad_install_postprocess_mode "APPLE")
       endif(P_RPATH_EXECUTABLE)
+      set(_brlcad_install_use_selective_rpath FALSE)
+      if(P_RPATH_SUPPORTS_SET_IF_NEEDED_PREPEND)
+        set(_brlcad_install_use_selective_rpath TRUE)
+      endif(P_RPATH_SUPPORTS_SET_IF_NEEDED_PREPEND)
       install(
         CODE
-          "include(\"${BRLCAD_EXT_INSTALL_POSTPROCESS_SCRIPT}\")\nbrlcad_install_binary_postprocess(\"${BRLCAD_EXT_INSTALL_POSTPROCESS_STAMP_DIR}\" \"${STRCLEAR_EXECUTABLE}\" \"${CMAKE_BINARY_DIR}/${tf}\" \"\${CMAKE_INSTALL_PREFIX}/${tf}\" \"PROGRAM\" \"${_brlcad_install_postprocess_mode}\" \"${_brlcad_install_postprocess_tool}\" \"${_brlcad_install_postprocess_rpath}\" \"${CMAKE_BINARY_DIR}/${LIB_DIR}\" \"${REL_RPATH}\" \"${BRLCAD_VERBOSE}\")"
+          "include(\"${BRLCAD_EXT_INSTALL_POSTPROCESS_SCRIPT}\")\nbrlcad_install_binary_postprocess(\"${BRLCAD_EXT_INSTALL_POSTPROCESS_STAMP_DIR}\" \"${STRCLEAR_EXECUTABLE}\" \"${CMAKE_BINARY_DIR}/${tf}\" \"\${CMAKE_INSTALL_PREFIX}/${tf}\" \"PROGRAM\" \"${_brlcad_install_postprocess_mode}\" \"${_brlcad_install_postprocess_tool}\" \"${_brlcad_install_postprocess_rpath}\" \"${CMAKE_BINARY_DIR}/${LIB_DIR}\" \"${REL_RPATH}\" \"${_brlcad_install_use_selective_rpath}\" \"${BRLCAD_VERBOSE}\")"
       )
       continue()
     endif("${tf}" IN_LIST ALL_BINARY_FILES)
@@ -2205,11 +2200,7 @@ macro(find_package_opencv)
   set(OpenCV_DIR_TMP "${OpenCV_DIR}")
   set(OpenCV_DIR "${CMAKE_BINARY_DIR}/${LIB_DIR}/cmake/opencv4")
   set(OpenCV_ROOT ${CMAKE_BINARY_DIR})
-  find_package(OpenCV CONFIG COMPONENTS core imgproc imgcodecs highgui)
-  if(OpenCV_FOUND AND OpenCV_VERSION_MAJOR GREATER_EQUAL 5)
-    # OpenCV 5 moved boundingRect from imgproc to geometry.
-    find_package(OpenCV CONFIG COMPONENTS geometry)
-  endif()
+  find_package(OpenCV COMPONENTS core features2d imgproc highgui)
   unset(OpenCV_ROOT)
 
   # If no bundled copy, see what the system has

misc/CMake/BRLCAD_Repository.cmake

@@ -1,108 +0,0 @@
-# Copyright (c) 2026 United States Government as represented by
-# the U.S. Army Research Laboratory.
-#
-# Redistribution and use in source and binary forms, with or without
-# modification, are permitted provided that the following conditions
-# are met:
-#
-# 1. Redistributions of source code must retain the above copyright
-# notice, this list of conditions and the following disclaimer.
-#
-# 2. Redistributions in binary form must reproduce the above
-# copyright notice, this list of conditions and the following
-# disclaimer in the documentation and/or other materials provided
-# with the distribution.
-#
-# 3. The name of the author may not be used to endorse or promote
-# products derived from this software without specific prior written
-# permission.
-#
-# THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS
-# OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
-# WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
-# ARE DISCLAIMED. IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY
-# DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
-# DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE
-# GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
-# INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
-# WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
-# NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
-# SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
-
-include_guard(GLOBAL)
-
-set(_BRLCAD_REPOSITORY_DEFAULT_ROOT "https://github.com/BRL-CAD")
-set(BRLCAD_REPOSITORY_ROOT "${_BRLCAD_REPOSITORY_DEFAULT_ROOT}" CACHE STRING
-  "Repository URL or local mirror directory containing bext and its dependency repositories")
-
-function(brlcad_repository_root out_root out_local)
-  set(root "${BRLCAD_REPOSITORY_ROOT}")
-  set(local FALSE)
-  if(root STREQUAL "")
-    message(FATAL_ERROR "BRLCAD_REPOSITORY_ROOT must be a repository URL or local mirror directory")
-  endif()
-
-  # Match bext's URL handling, but resolve relative paths against BRL-CAD.
-  # Drive paths must be checked before scp-style SSH locations.
-  if(root MATCHES "^file://")
-    set(local TRUE)
-    string(REPLACE "=" "%3D" root "${root}")
-  elseif(root MATCHES "^[A-Za-z]:[/\\\\]" OR
-      NOT (root MATCHES "^[A-Za-z][A-Za-z0-9+.-]*://" OR root MATCHES "^[^/]+:.+"))
-    set(local TRUE)
-    string(REPLACE "\\" "/" root "${root}")
-    if(WIN32 OR NOT root MATCHES "^[A-Za-z]:/")
-      get_filename_component(root "${root}" ABSOLUTE BASE_DIR "${CMAKE_SOURCE_DIR}")
-    endif()
-    # Escape literal percent signs first so native paths cannot introduce
-    # URL escapes.  An equals sign also has meaning in Git's -c assignments.
-    string(REPLACE "%" "%25" root "${root}")
-    string(REPLACE " " "%20" root "${root}")
-    string(REPLACE "#" "%23" root "${root}")
-    string(REPLACE "?" "%3F" root "${root}")
-    string(REPLACE "=" "%3D" root "${root}")
-    if(root MATCHES "^[A-Za-z]:/")
-      set(root "file:///${root}")
-    elseif(root MATCHES "^//")
-      set(root "file:${root}")
-    else()
-      set(root "file://${root}")
-    endif()
-  endif()
-
-  if(NOT root MATCHES "^file:///*$")
-    string(REGEX REPLACE "/+$" "" root "${root}")
-    string(APPEND root "/")
-  else()
-    set(root "file:///")
-  endif()
-  set(${out_root} "${root}" PARENT_SCOPE)
-  set(${out_local} "${local}" PARENT_SCOPE)
-endfunction()
-
-function(brlcad_repository_git_command out_command)
-  set(git_command "${GIT_EXEC}")
-  if(WIN32)
-    list(APPEND git_command -c core.longpaths=true)
-  endif()
-  if(_BRLCAD_REPOSITORY_LOCAL)
-    list(APPEND git_command -c protocol.file.allow=always)
-    # This also constrains recursive Git processes and user URL rewrites.
-    list(PREPEND git_command "${CMAKE_COMMAND}" -E env GIT_ALLOW_PROTOCOL=file --)
-  endif()
-  set(${out_command} "${git_command}" PARENT_SCOPE)
-endfunction()
-
-brlcad_repository_root(_BRLCAD_REPOSITORY_ROOT _BRLCAD_REPOSITORY_LOCAL)
-set(_BRLCAD_BEXT_REPOSITORY "${_BRLCAD_REPOSITORY_ROOT}bext.git")
-message(STATUS "BRL-CAD repository root: ${_BRLCAD_REPOSITORY_ROOT}")
-if(_BRLCAD_REPOSITORY_LOCAL)
-  message(STATUS "BRL-CAD Git acquisition: file transport only")
-endif()
-
-# Local Variables:
-# tab-width: 8
-# mode: cmake
-# indent-tabs-mode: t
-# End:
-# ex: shiftwidth=2 tabstop=8

misc/CMake/CMakeLists.txt

@@ -15,7 +15,6 @@ set(
   BRLCAD_Options.cmake
   BRLCAD_Probe_Validation.cmake
   BRLCAD_Regress_Util.cmake
-  BRLCAD_Repository.cmake
   BRLCAD_Summary.cmake
   BRLCAD_Targets.cmake
   BRLCAD_Test_Wrappers.cmake

misc/CMake/Distcheck.cmake

@@ -95,33 +95,6 @@ add_custom_target(
 )
 set_target_properties(distcheck-source_archives PROPERTIES FOLDER "BRL-CAD Distribution Checking")
 
-# Dependency paths must survive CMake, build-tool, and shell parsing.  Load
-# their values from a cache file rather than embedding them in command text.
-set(distcheck_dependency_cache "${CMAKE_CURRENT_BINARY_DIR}/CMakeTmp/distcheck_dependencies.cmake")
-set(distcheck_dependency_vars BRLCAD_REPOSITORY_ROOT BRLCAD_EXT_DIR)
-foreach(distcheck_optional_var BRLCAD_EXT_SOURCE_DIR BRLCAD_EXT_PARALLEL)
-  if(DEFINED ${distcheck_optional_var})
-    list(APPEND distcheck_dependency_vars ${distcheck_optional_var})
-  endif()
-endforeach()
-file(WRITE "${distcheck_dependency_cache}" "")
-foreach(distcheck_var IN LISTS distcheck_dependency_vars)
-  set(distcheck_value "${${distcheck_var}}")
-  if(distcheck_var STREQUAL "BRLCAD_REPOSITORY_ROOT")
-    set(distcheck_value "${_BRLCAD_REPOSITORY_ROOT}")
-  endif()
-  set(distcheck_type STRING)
-  if(distcheck_var MATCHES "_DIR$")
-    set(distcheck_type PATH)
-  endif()
-  string(REPLACE "\\" "\\\\" distcheck_value "${distcheck_value}")
-  string(REPLACE "\"" "\\\"" distcheck_value "${distcheck_value}")
-  string(REPLACE "$" "\\$" distcheck_value "${distcheck_value}")
-  file(APPEND "${distcheck_dependency_cache}"
-    "set(${distcheck_var} \"${distcheck_value}\" CACHE ${distcheck_type} \"Inherited dependency setting\" FORCE)\n")
-endforeach()
-distclean("${distcheck_dependency_cache}")
-
 # Utility function for defining individual distcheck targets
 macro(
   create_distcheck
@@ -147,10 +120,7 @@ macro(
 
     # Need to set these locally so configure_file will pick them up...
     set(TARGET_SUFFIX ${TARGET_SUFFIX})
-    # Keep machine-specific settings out of the options saved for CI above.
-    # The shared cache replaces each variant's unquoted output argument.
-    string(REPLACE "-DBRLCAD_EXT_DIR=${BRLCAD_EXT_DIR}" "" CMAKE_OPTS "${CMAKE_OPTS_IN}")
-    string(APPEND CMAKE_OPTS " -C \"${distcheck_dependency_cache}\"")
+    set(CMAKE_OPTS ${CMAKE_OPTS_IN})
 
     # For configure_file, need to set these as variables not just input parameters
     set(source_dir "${source_dir}")

misc/debian/changelog

@@ -1,9 +1,3 @@
-brlcad (7.46.0-0) unstable; urgency=low
-
-  * update brlcad version
-
- -- Cliff Yapp <cyapp@brlcad.org>  Thu, 07 Oct 2026 20:24:00 +0000
-
 brlcad (7.44.0-0) unstable; urgency=low
 
   * update brlcad version

misc/macosx/Resources/ReadMe.rtfd/TXT.rtf

@@ -11,7 +11,7 @@
 \f0\i\b\fs72 \cf0 README
 \fs36 \
 
-\f1\i0\b0\fs28 This is BRL-CAD 7.46.0 for Mac OS X.  Please send any bugs, comments, questions, or suggestions to the project website.\
+\f1\i0\b0\fs28 This is BRL-CAD 7.44.0 for Mac OS X.  Please send any bugs, comments, questions, or suggestions to the project website.\
 
 \fs22 \
 

misc/macosx/Resources/Welcome.rtfd/TXT.rtf

@@ -4,7 +4,7 @@
 \margl1440\margr1440\vieww8020\viewh6280\viewkind0
 \pard\tx560\tx1120\tx1680\tx2240\tx2800\tx3360\tx3920\tx4480\tx5040\tx5600\tx6160\tx6720\ql\qnatural
 
-\f0\fs22 \cf0 This installer will guide you through the installation of BRL-CAD 7.46.0 for Mac OS X and Darwin.\
+\f0\fs22 \cf0 This installer will guide you through the installation of BRL-CAD 7.44.0 for Mac OS X and Darwin.\
 \
 This distribution of BRL-CAD requires approximately 500 MB of available disk space and does not include source code.\
 \

misc/nirt/default.nrt

@@ -1,5 +1,6 @@
 # default.nrt
 # Description: default   - Standard NIRT Reporting Format
+fmt a ""
 fmt r "Origin (x y z) = (%.8f %.8f %.8f)  (h v d) = (%.4f %.4f %.4f)\nDirection (x y z) = (%.8f %.8f %.8f)  (az el) = (%.8f %.8f)\n" x_orig y_orig z_orig h v d_orig x_dir y_dir z_dir a e
 fmt h "    Region Name               Entry (x y z)              LOS  Obliq_in Attrib\n"
 fmt p "%-20s (%9.4f %9.4f %9.4f) %8.4f %8.4f%s\n" reg_name x_in y_in z_in los obliq_in attributes

regress/fuzz/CMakeLists.txt

@@ -117,20 +117,6 @@ if(TARGET libged)
 
 endif(TARGET libged)
 
-if(TARGET librt AND FUZZER_SUPPORTED)
-  set(RT_EDIT_FUZZ_DIR "${CMAKE_SOURCE_DIR}/src/librt/tests/edit")
-  brlcad_addexec(
-    fuzz_rt_edit_test
-    "${RT_EDIT_FUZZ_DIR}/fuzz.cpp;${RT_EDIT_FUZZ_DIR}/ars.cpp;${RT_EDIT_FUZZ_DIR}/bot.cpp"
-    "libbu;librt;libnmg;libwdb;libbrep"
-    TEST
-  )
-  set_target_properties(fuzz_rt_edit_test PROPERTIES UNITY_BUILD OFF)
-  target_compile_options(fuzz_rt_edit_test PRIVATE ${FUZZER_EXE_COMPILE_OPTIONS})
-  target_link_options(fuzz_rt_edit_test PRIVATE ${FUZZER_LINK_OPTIONS})
-  BRLCAD_ADD_TEST(NAME fuzz_rt_edit_smoke COMMAND fuzz_rt_edit_test -runs=1 -max_total_time=5)
-endif()
-
 cmakefiles(
   CMakeLists.txt
   fuzz_ged.cpp

regress/gcv/CMakeLists.txt

@@ -79,9 +79,7 @@ endif(NOT EXISTS \"\${GCV_EXEC}\")\n"
   # setup conversion tests
   foreach(in_file ${INPUT_FILES})
     get_filename_component(curr_in_basename ${in_file} NAME_WLE)
-    # Tests in the same directory may run concurrently, so their output paths
-    # must not overlap.
-    set(out_path "${CMAKE_CURRENT_BINARY_DIR}/conv-${TARGET_DIRECTION}-${curr_in_basename}")
+    set(out_path "${CMAKE_CURRENT_BINARY_DIR}/conv-${curr_in_basename}")
 
     # need to split suffix so we can distinguish from different converters
     if(${fmt}_EXPORTER)
@@ -122,7 +120,7 @@ endif(NOT EXISTS \"\${GCV_EXEC}\")\n"
         )
       endif(${fmt}_GCV_OPTIONS)
 
-      set(out_distclean ${out_distclean} "${out_path}-${out_fmt}${out_suffix}")
+      set(out_distclean ${out_distclean}} "${out_path}-${out_fmt}${out_suffix}")
       set(expected_output_files "${expected_output_files} \"${out_path}-${out_fmt}${out_suffix}\"")
     endforeach(out_fmt ${${fmt}_CONVERTERS})
 
@@ -164,8 +162,19 @@ execute_process(
 
 file(APPEND \"\${LOGFILE}\" \"\${log_msg}\\n\")
 unset(log_msg)
+  "
+    )
 
-# check output files exist and are a reasonable size
+    set(OUTPUT_FILES ${OUTPUT_FILES} ${out_distclean})
+    unset(expected_output_files)
+    unset(out_distclean)
+  endforeach(in_file ${INPUT_FILES})
+
+  file(
+    APPEND
+    ${cmake_in_path}
+    "
+# check output files exists and are a reasonable size
 foreach(out_file \${ofiles})
   if(NOT EXISTS \"\${out_file}\")
     file(APPEND \"\${LOGFILE}\" \"FAILURE: \${curr_result}, \${out_file} not found.\")
@@ -182,18 +191,13 @@ foreach(out_file \${ofiles})
       file(APPEND \"\${LOGFILE}\" \"FAILURE: \${curr_result} file size is \${file_size} for \${out_file}\\n    Expecting ${EXPECTED_FILE_SIZE} with delta ${FILE_SIZE_DELTA}\")
       message(FATAL_ERROR \"[gcv] Failure while converting ${in_file}, see \${LOGFILE} for more info.\")
     endif(\"\${size_delta}\" GREATER \"${FILE_SIZE_DELTA}\")
-  endif(NOT EXISTS \"\${out_file}\")
+  endif(NOT EXISTS \"\${out_file}\")\n
 
   # cleanup file from build dir
   file(REMOVE \"\${out_file}\")
 endforeach(out_file \${ofiles})
-  "
-    )
-
-    set(OUTPUT_FILES ${OUTPUT_FILES} ${out_distclean})
-    unset(expected_output_files)
-    unset(out_distclean)
-  endforeach(in_file ${INPUT_FILES})
+    "
+  )
 
   # assume the fmt is the gcv dependency unless specified
   if(${fmt}_DEPS)

regress/gcv/gltf/CMakeLists.txt

@@ -2,7 +2,7 @@ if(TARGET gcv-assetimport)
   gcv_regress_util("gltf" "tetra-gltf2.0.gltf;tetra-glb2.0.glb" DEPS "gcv;gcv-assetimport" GCV_OPTIONS "--input-format=assetimport")
   gcv_regress_util("glb" "tetra.g" DEPS "gcv;gcv-assetimport" GCV_OPTIONS "--output-format=assetimport" EXPECTED_FILE_SIZE "1476" EXPORTER)
   gcv_regress_util("gltf" "tetra.g" DEPS "gcv;gcv-assetimport" GCV_OPTIONS "--output-format=assetimport" EXPECTED_FILE_SIZE "3356" EXPORTER)
-  distclean("${CMAKE_CURRENT_BINARY_DIR}/conv-regress-g-gltf-tetra.bin")
+  distclean("${CMAKE_CURRENT_BINARY_DIR}/conv-tetra.bin")
 endif(TARGET gcv-assetimport)
 
 if(TARGET gcv-gltf)

regress/gcv/obj/CMakeLists.txt

@@ -1,7 +1,7 @@
 if(TARGET gcv-assetimport)
   set(dep_str DEPS "gcv;gcv-obj;gcv-assetimport" CONVERTERS "assetimport")
 
-  distclean("${CMAKE_CURRENT_BINARY_DIR}/conv-regress-g-obj-tetra-assetimport.mtl") # hard to track companion files
+  distclean("${CMAKE_CURRENT_BINARY_DIR}/conv-tetra-assetimport.mtl") # hard to track companion files
 else()
   set(dep_str "")
 endif(TARGET gcv-assetimport)
@@ -9,7 +9,7 @@ endif(TARGET gcv-assetimport)
 gcv_regress_util("obj" "tetra.obj" ${dep_str})
 gcv_regress_util("obj_unnamed_group" "tetra-unnamed-group.obj" DEPS "gcv;gcv-obj")
 gcv_regress_util("obj" "tetra.g" ${dep_str} EXPORTER)
-distclean("${CMAKE_CURRENT_BINARY_DIR}/conv-regress-g-obj-tetra.mtl")
+distclean("${CMAKE_CURRENT_BINARY_DIR}/conv-tetra.mtl")
 
 cmakefiles(
   CMakeLists.txt

regress/ged/regress_facetize.cpp

@@ -68,13 +68,6 @@
  *          ("/tc6_reg.c/sub6.s#sub"), so the lookup missed and no variant was
  *          applied, leaving the phantom face.  After the fix the plan is built
  *          from region roots so keys match and the phantom face is eliminated.
- *
- *   TC7  - NMG Boolean isolation and output modes: verifies that an interrupted
- *          worker cannot commit output, then exercises BoT, NMG, region, and
- *          in-place output paths.
- *
- *   TC8  - implicit solid roots survive the working-database close and remain
- *          valid members of the imported output hierarchy.
  */
 
 #include "common.h"
@@ -85,9 +78,7 @@
 #include <string>
 
 #include "bu/app.h"
-#include "bu/env.h"
 #include "bu/file.h"
-#include "bu/log.h"
 #include "bu/malloc.h"
 #include "bu/str.h"
 #include "bu/vls.h"
@@ -151,36 +142,6 @@ write_comb2(struct rt_wdb *wdbp,
    0, 0, 0, 0, 0, 0, 0);
 }
 
-static int
-write_comb1(struct rt_wdb *wdbp, const char *comb_name, const char *member, int is_region)
-{
-    struct bu_list head;
-    BU_LIST_INIT(&head);
-    mk_addmember(member, &head, NULL, WMOP_UNION);
-    return mk_comb(wdbp, comb_name, &head, is_region, NULL, NULL, NULL,
-	    0, 0, 0, 0, 0, 0, 0);
-}
-
-static int
-run_facetize_command(const char *gfile, int argc, const char **argv, int verbose)
-{
-    struct ged *gedp = ged_open("db", gfile, 1);
-    if (!gedp) {
-	bu_log("[regress_facetize] ged_open(%s) failed\n", gfile);
-	return BRLCAD_ERROR;
-    }
-
-    int ret = ged_exec(gedp, argc, argv);
-    if (verbose || ret != BRLCAD_OK) {
-	const char *log = bu_vls_cstr(gedp->ged_result_str);
-	if (log && log[0])
-	    bu_log("[facetize] %s\n", log);
-    }
-
-    ged_close(gedp);
-    return ret;
-}
-
 /**
  * Open @a gfile with ged_open and run "facetize @a input @a output".
  * Returns BRLCAD_OK on success, BRLCAD_ERROR otherwise.
@@ -188,8 +149,23 @@ run_facetize_command(const char *gfile, int argc, const char **argv, int verbose
 static int
 run_facetize(const char *gfile, const char *input, const char *output, int verbose)
 {
+    struct ged *gedp = ged_open("db", gfile, 1);
+    if (!gedp) {
+bu_log("[regress_facetize] ged_open(%s) failed\n", gfile);
+return BRLCAD_ERROR;
+    }
+
     const char *av[4] = {"facetize", input, output, NULL};
-    return run_facetize_command(gfile, 3, av, verbose);
+    int ret = ged_exec(gedp, 3, av);
+
+    if (verbose || ret != BRLCAD_OK) {
+const char *log = bu_vls_cstr(gedp->ged_result_str);
+if (log && log[0])
+    bu_log("[facetize] %s\n", log);
+    }
+
+    ged_close(gedp);
+    return ret;
 }
 
 /**
@@ -236,46 +212,6 @@ bu_log("[regress_facetize] check_bot: '%s' type %d is not a BoT\n",
     return ok;
 }
 
-static int
-check_object_type(const char *gfile, const char *object_name, int expected_type)
-{
-    struct db_i *dbip = db_open(gfile, DB_OPEN_READONLY);
-    if (!dbip || db_dirbuild(dbip) < 0) {
-	if (dbip)
-	    db_close(dbip);
-	return BRLCAD_ERROR;
-    }
-
-    struct directory *dp = db_lookup(dbip, object_name, LOOKUP_QUIET);
-    if (!dp) {
-	bu_log("[regress_facetize] '%s' does not exist\n", object_name);
-	db_close(dbip);
-	return BRLCAD_ERROR;
-    }
-
-    struct rt_db_internal intern;
-    RT_DB_INTERNAL_INIT(&intern);
-    if (rt_db_get_internal(&intern, dp, dbip, NULL) < 0) {
-	db_close(dbip);
-	return BRLCAD_ERROR;
-    }
-    int ret = (intern.idb_minor_type == expected_type) ? BRLCAD_OK : BRLCAD_ERROR;
-    if (ret != BRLCAD_OK)
-	bu_log("[regress_facetize] '%s' type %d, expected %d\n",
-		object_name, intern.idb_minor_type, expected_type);
-    rt_db_free_internal(&intern);
-    db_close(dbip);
-    return ret;
-}
-
-static int
-count_log_hook(void *data, void *UNUSED(message))
-{
-    int *count = (int *)data;
-    (*count)++;
-    return 0;
-}
-
 /**
  * Run rt_crofton_shoot on @a obj_name in @a gfile.
  * Stores SA (mm2) and VOL (mm3) in *sa / *vol.
@@ -1184,183 +1120,6 @@ tc6_havoc_wind9_pattern(const char *tmpdir, int verbose)
     return BRLCAD_OK;
 }
 
-static int
-tc7_nmg_modes(const char *tmpdir, int verbose)
-{
-    bu_log("[regress_facetize] TC7: NMG Boolean and output modes...\n");
-
-    struct bu_vls gpath = BU_VLS_INIT_ZERO;
-    bu_vls_printf(&gpath, "%s/tc7_nmg.g", tmpdir);
-    const char *gfile = bu_vls_cstr(&gpath);
-    if (bu_file_exists(gfile, NULL))
-	bu_file_delete(gfile);
-
-    struct db_i *dbip = db_create(gfile, 5);
-    if (!dbip) {
-	bu_vls_free(&gpath);
-	return BRLCAD_ERROR;
-    }
-    struct rt_wdb *wdbp = wdb_dbopen(dbip, RT_WDB_TYPE_DB_DEFAULT);
-    point_t box_pts[8] = {
-	{0.0,  0.0,  0.0}, {10.0, 0.0,  0.0},
-	{10.0, 10.0, 0.0}, {0.0,  10.0, 0.0},
-	{0.0,  0.0, 10.0}, {10.0, 0.0, 10.0},
-	{10.0, 10.0, 10.0}, {0.0, 10.0, 10.0}
-    };
-    if (write_arb8(wdbp, "box.s", box_pts) < 0 ||
-	    write_comb1(wdbp, "box.r", "box.s", 1) < 0 ||
-	    write_comb1(wdbp, "box.c", "box.s", 0) < 0) {
-	wdb_close(wdbp);
-	bu_vls_free(&gpath);
-	return BRLCAD_ERROR;
-    }
-    wdb_close(wdbp);
-
-    const char *bot_av[] = {"facetize", "--nmg-booleval", "box.r", "box.bot", NULL};
-    const char *fault_av[] = {"facetize", "--nmg-booleval", "box.r", "fault.bot", NULL};
-    const char *nmg_av[] = {"facetize", "-n", "box.r", "box.nmg", NULL};
-    const char *regions_av[] = {"facetize", "-r", "--nmg-booleval", "box.r", "box.regions", NULL};
-    const char *fallback_av[] = {"facetize", "-r", "--nmg-booleval", "box.c", "box-fallback.bot", NULL};
-    const char *object_in_place_av[] = {"facetize", "--in-place", "--nmg-booleval", "box.c", NULL};
-    const char *in_place_av[] = {"facetize", "-r", "--in-place", "--nmg-booleval", "box.r", NULL};
-
-    int ret = BRLCAD_OK;
-    auto check_tc7 = [&ret](bool passed, const char *description) {
-	if (passed)
-	    return;
-	bu_log("[regress_facetize] TC7: FAIL - %s\n", description);
-	ret = BRLCAD_ERROR;
-    };
-
-    struct bu_vls fault_path = BU_VLS_INIT_ZERO;
-    bu_vls_printf(&fault_path, "%s/tc7_nmg_worker.fault", tmpdir);
-    (void)bu_file_delete(bu_vls_cstr(&fault_path));
-    (void)bu_setenv("LIBGED_FACETIZE_TEST_FAULT",
-	    "nmg-before-ready:fault.bot", 1);
-    (void)bu_setenv("LIBGED_FACETIZE_TEST_FAULT_FILE",
-	    bu_vls_cstr(&fault_path), 1);
-    int fault_ret = run_facetize_command(gfile, 4, fault_av, verbose);
-    (void)bu_setenv("LIBGED_FACETIZE_TEST_FAULT", "", 1);
-    (void)bu_setenv("LIBGED_FACETIZE_TEST_FAULT_FILE", "", 1);
-    check_tc7(fault_ret == BRLCAD_ERROR,
-	    "an interrupted NMG Boolean worker did not report failure");
-    check_tc7(bu_file_exists(bu_vls_cstr(&fault_path), NULL),
-	    "NMG Boolean worker fault did not trigger");
-    check_tc7(check_object_type(gfile, "fault.bot", ID_BOT) != BRLCAD_OK,
-	    "an interrupted NMG Boolean worker committed an output");
-    (void)bu_file_delete(bu_vls_cstr(&fault_path));
-    bu_vls_free(&fault_path);
-
-    int hook_calls = 0;
-    bu_log_add_hook(count_log_hook, &hook_calls);
-    int bot_ret = run_facetize_command(gfile, 4, bot_av, verbose);
-    int hook_calls_before_probe = hook_calls;
-    bu_log("facetize NMG hook restoration probe\n");
-    bu_log_delete_hook(count_log_hook, &hook_calls);
-
-    check_tc7(bot_ret == BRLCAD_OK, "NMG Boolean BoT command failed");
-    check_tc7(hook_calls == hook_calls_before_probe + 1,
-	    "temporary log hook was not restored");
-    check_tc7(check_bot_exists(gfile, "box.bot") == BRLCAD_OK,
-	    "NMG Boolean BoT output is missing");
-
-    int nmg_ret = run_facetize_command(gfile, 4, nmg_av, verbose);
-    check_tc7(nmg_ret == BRLCAD_OK, "NMG output command failed");
-    check_tc7(check_object_type(gfile, "box.nmg", ID_NMG) == BRLCAD_OK,
-	    "NMG output has the wrong type");
-
-    int regions_ret = run_facetize_command(gfile, 5, regions_av, verbose);
-    check_tc7(regions_ret == BRLCAD_OK, "NMG Boolean region command failed");
-    check_tc7(check_object_type(gfile, "box.regions", ID_COMBINATION) == BRLCAD_OK,
-	    "NMG Boolean region output is missing");
-
-    int fallback_ret = run_facetize_command(gfile, 5, fallback_av, verbose);
-    check_tc7(fallback_ret == BRLCAD_OK,
-	    "region mode object fallback command failed");
-    check_tc7(check_bot_exists(gfile, "box-fallback.bot") == BRLCAD_OK,
-	    "region mode object fallback output is missing");
-
-    int object_in_place_ret = run_facetize_command(gfile, 4,
-	    object_in_place_av, verbose);
-    check_tc7(object_in_place_ret == BRLCAD_OK,
-	    "in-place NMG Boolean object command failed");
-    check_tc7(check_object_type(gfile, "box.c", ID_BOT) == BRLCAD_OK,
-	    "in-place NMG Boolean object has the wrong type");
-
-    int in_place_ret = run_facetize_command(gfile, 5, in_place_av, verbose);
-    check_tc7(in_place_ret == BRLCAD_OK,
-	    "in-place NMG Boolean region command failed");
-    check_tc7(check_object_type(gfile, "box.r", ID_COMBINATION) == BRLCAD_OK,
-	    "in-place region root has the wrong type");
-    check_tc7(check_bot_exists(gfile, "box.r.bot") == BRLCAD_OK,
-	    "in-place region BoT output is missing");
-
-    bu_log("[regress_facetize] TC7: %s\n", ret == BRLCAD_OK ? "PASS" : "FAIL");
-    bu_vls_free(&gpath);
-    return ret;
-}
-
-/* ------------------------------------------------------------------ */
-/* TC8: implicit root replacement persists across database reopen      */
-/* ------------------------------------------------------------------ */
-static int
-tc8_implicit_root_persistence(const char *tmpdir, int verbose)
-{
-    bu_log("[regress_facetize] TC8: implicit root persistence...\n");
-
-    struct bu_vls gpath = BU_VLS_INIT_ZERO;
-    bu_vls_printf(&gpath, "%s/tc8_implicit_root.g", tmpdir);
-    const char *gfile = bu_vls_cstr(&gpath);
-    if (bu_file_exists(gfile, NULL))
-	bu_file_delete(gfile);
-
-    struct db_i *dbip = db_create(gfile, 5);
-    if (!dbip) {
-	bu_vls_free(&gpath);
-	return BRLCAD_ERROR;
-    }
-    struct rt_wdb *wdbp = wdb_dbopen(dbip, RT_WDB_TYPE_DB_DEFAULT);
-    point_t region_points[8] = {
-	{0.0, 0.0, 0.0}, {10.0, 0.0, 0.0},
-	{10.0, 10.0, 0.0}, {0.0, 10.0, 0.0},
-	{0.0, 0.0, 10.0}, {10.0, 0.0, 10.0},
-	{10.0, 10.0, 10.0}, {0.0, 10.0, 10.0}
-    };
-    point_t loose_points[8] = {
-	{20.0, 0.0, 0.0}, {30.0, 0.0, 0.0},
-	{30.0, 10.0, 0.0}, {20.0, 10.0, 0.0},
-	{20.0, 0.0, 10.0}, {30.0, 0.0, 10.0},
-	{30.0, 10.0, 10.0}, {20.0, 10.0, 10.0}
-    };
-    if (!wdbp || write_arb8(wdbp, "region.s", region_points) < 0 ||
-	    write_arb8(wdbp, "loose.s", loose_points) < 0 ||
-	    write_comb1(wdbp, "box.r", "region.s", 1) < 0 ||
-	    write_comb2(wdbp, "assembly.c", "box.r", 'u',
-		    "loose.s", 'u') < 0) {
-	if (wdbp)
-	    wdb_close(wdbp);
-	bu_vls_free(&gpath);
-	return BRLCAD_ERROR;
-    }
-    wdb_close(wdbp);
-
-    const char *facetize_args[] = {
-	"facetize", "-r", "assembly.c", "assembly.regions", NULL
-    };
-    int ret = run_facetize_command(gfile, 4, facetize_args, verbose);
-    if (ret == BRLCAD_OK &&
-	    (check_object_type(gfile, "assembly.regions", ID_COMBINATION) != BRLCAD_OK ||
-	     check_object_type(gfile, "facetize_assembly.c", ID_COMBINATION) != BRLCAD_OK ||
-	     check_bot_exists(gfile, "box.r.bot") != BRLCAD_OK ||
-	     check_bot_exists(gfile, "facetize_loose.s") != BRLCAD_OK))
-	ret = BRLCAD_ERROR;
-
-    bu_log("[regress_facetize] TC8: %s\n",
-	    ret == BRLCAD_OK ? "PASS" : "FAIL");
-    bu_vls_free(&gpath);
-    return ret;
-}
-
 /* ------------------------------------------------------------------ */
 /* main                                                                 */
 /* ------------------------------------------------------------------ */
@@ -1392,8 +1151,6 @@ return 1;
     if (tc4_havoc_reuse(tmpdir, verbose)         != BRLCAD_OK) ret = 1;
     if (tc5_near_coplanar_subTOL(tmpdir, verbose) != BRLCAD_OK) ret = 1;
     if (tc6_havoc_wind9_pattern(tmpdir, verbose) != BRLCAD_OK) ret = 1;
-    if (tc7_nmg_modes(tmpdir, verbose)             != BRLCAD_OK) ret = 1;
-    if (tc8_implicit_root_persistence(tmpdir, verbose) != BRLCAD_OK) ret = 1;
 
     if (ret == 0)
 bu_log("[regress_facetize] All tests PASSED\n");

regress/gqa/CMakeLists.txt

@@ -23,7 +23,6 @@ set(
   gqa.mass_mix.am.out
   gqa.mass_mix.aw.out
   gqa.mass_mix.u0.out
-  gqa.missing_density.out
   gqa.mged
   gqa.overlaps.plot3
   gqa.sphere_adaptive.out

regress/gqa/gqa.sh

@@ -67,11 +67,6 @@ run_capture ( ) {
     return $ret
 }
 
-expect_failure ( ) {
-    "$@"
-    test $? -ne 0
-}
-
 extract_last_number ( ) {
     label="$1"
     file="$2"
@@ -106,19 +101,6 @@ assert_close ( ) {
     fi
 }
 
-assert_contains ( ) {
-    expected="$1"
-    file="$2"
-    desc="$3"
-
-    if grep -F -q "$expected" "$file" ; then
-	log "PASS: $desc"
-    else
-	log "FAIL: $desc"
-	STATUS="`expr $STATUS + 1`"
-    fi
-}
-
 MGED="`ensearch mged`"
 if test ! -f "$MGED" ; then
     log "Unable to find mged, aborting"
@@ -195,10 +177,6 @@ adjust mass_air.r GIFTmater 2
 g mass_air_only.g mass_air.r
 g mass_air_mix.g mass_solid.r mass_air.r
 
-in missing_density.s rpp 2 3 0 1 0 1
-r missing_density.r u missing_density.s
-adjust missing_density.r GIFTmater 123
-
 in sphere.s sph 0 0 0 1
 r sphere.r u sphere.s
 adjust sphere.r GIFTmater 5
@@ -265,7 +243,6 @@ run_capture gqa.mass_mix.aw.out $GQABIN -u m,m^3,kg -g 250mm-50mm -Aw gqa.g mass
 run_capture gqa.mass_mix.am.out $GQABIN -u m,m^3,kg -g 250mm-50mm -Am gqa.g mass_air_mix.g
 run_capture gqa.mass_mix.u0.out $GQABIN -u m,m^3,kg -g 250mm-50mm -U 0 -Aw gqa.g mass_air_mix.g
 run_capture gqa.sphere_adaptive.out $GQABIN -u m,m^3,kg -q -Av -g 0.05m gqa.g sphere.r
-run_capture gqa.missing_density.out expect_failure $GQABIN -u m,m^3,kg -g 250mm-50mm -f density_table.txt -Aw gqa.g missing_density.r
 
 MASS_AIR_AW="`extract_last_number 'Average total weight:' gqa.mass_air.aw.out`"
 MASS_AIR_AM="`extract_last_number 'Average total weight:' gqa.mass_air.am.out`"
@@ -280,7 +257,6 @@ assert_close 1001 "$MASS_MIX_AW" 0.001 "mixed solid+air modeled weight with -Aw"
 assert_close "$MASS_MIX_AW" "$MASS_MIX_AM" 0.001 "-Aw and -Am agree on mixed modeled air"
 assert_close 1000 "$MASS_MIX_U0" 0.001 "-U 0 excludes modeled air from weight"
 assert_close 4.18879 "$SPHERE_VOLUME" 0.01 "adaptive sphere volume refines beyond symmetric first pass"
-assert_contains "Density 123 on region /missing_density.r is not defined in the density table." gqa.missing_density.out "undefined density produces an error"
 
 
 if [ $STATUS = 0 ] ; then

regress/mged/CMakeLists.txt

@@ -68,7 +68,6 @@ set(
   sca_view.mged
   search_exec.mged
   sed.mged
-  sed_reset.mged
   set_perspective.mged
   setview.mged
   size.mged
@@ -101,7 +100,6 @@ if(TARGET mged)
       continue()
     endif(TNAME STREQUAL "")
 
-    set_property(DIRECTORY APPEND PROPERTY CMAKE_CONFIGURE_DEPENDS "${CMAKE_CURRENT_SOURCE_DIR}/${testfile}")
     file(READ "${CMAKE_CURRENT_SOURCE_DIR}/${testfile}" TEST_CONTENTS)
     if(TNAME STREQUAL "search_exec")
       string(CONFIGURE "${TEST_CONTENTS}" TEST_CONTENTS @ONLY)
@@ -128,12 +126,6 @@ if(TARGET mged)
     elseif(TNAME STREQUAL "arb")
       set(TEST_EXPECT "'arb' testing completed")
       set(TEST_SOURCE "${CMAKE_CURRENT_BINARY_DIR}/regress_${testfile}")
-    elseif(TNAME MATCHES "^(extrude|facedef|mirface|permute)$")
-      set(TEST_EXPECT "'${TNAME}' testing completed")
-      set(TEST_SOURCE "${CMAKE_CURRENT_BINARY_DIR}/regress_${testfile}")
-    elseif(TNAME STREQUAL "edgedir")
-      set(TEST_EXPECT "'edarb edgedir' testing completed")
-      set(TEST_SOURCE "${CMAKE_CURRENT_BINARY_DIR}/regress_${testfile}")
     elseif(TNAME STREQUAL "killtree")
       set(TEST_EXPECT "KILL Solid:  killtree_dup.s")
       set(TEST_REJECT "db_lookup(killtree_dup.s) failed")
@@ -146,10 +138,6 @@ if(TARGET mged)
     elseif(TNAME STREQUAL "prim_edit")
       set(TEST_EXPECT "Primitive parameter editing testing completed")
       set(TEST_REJECT "FAIL:")
-    elseif(TNAME STREQUAL "sed_reset")
-      set(TEST_EXPECT "PASS: 'sed_reset' testing completed")
-      set(TEST_REJECT "FAIL:")
-      set(TEST_SOURCE "${CMAKE_CURRENT_BINARY_DIR}/regress_${testfile}")
     elseif(TNAME STREQUAL "overlap_plotting")
       set(TEST_DB "${CMAKE_SOURCE_DIR}/src/libged/tests/ged_tests.g")
       set(TEST_EXPECT "PASS overlap plotting")
@@ -199,7 +187,6 @@ distclean(${mged_outfiles})
 
 cmakefiles(
   CMakeLists.txt
-  csg_edit.mged
   mged.tcl
   regression_resources.tcl
   search_exec_procs.tcl

regress/mged/arb.mged

@@ -119,35 +119,4 @@ if {![catch {get $invalid_inner}]} {
 }
 
 
-# ARBN uses one independent thickness per plane.  Exercise unequal offsets so
-# the test also verifies the plane ordering and inward normal convention.
-set arbn_outer inside_arbn_outer.s
-set arbn_inner inside_arbn_inner.s
-set arbn_thicknesses {1 2 3 4 1.5 0.5}
-in $arbn_outer arbn 6 \
-    1 0 0 10  -1 0 0 0 \
-    0 1 0 10   0 -1 0 0 \
-    0 0 1 10   0 0 -1 0
-if {[catch {inside $arbn_outer $arbn_inner {*}$arbn_thicknesses} result]} {
-    error "inside failed for ARBN: $result"
-}
-set arbn_definition [get $arbn_inner]
-set expected_offsets {9 -2 7 -4 8.5 -0.5}
-for {set plane 0} {$plane < 6} {incr plane} {
-    set equation [lindex $arbn_definition [expr {4 + 2 * $plane}]]
-    set actual_offset [lindex $equation 3]
-    if {abs($actual_offset - [lindex $expected_offsets $plane]) > $plane_tolerance} {
-	error "ARBN inside plane $plane has offset $actual_offset: $arbn_definition"
-    }
-}
-facetize -n inside_arbn_inner_nmg.s $arbn_inner
-
-set invalid_arbn_inner inside_arbn_invalid.s
-if {![catch {inside $arbn_outer $invalid_arbn_inner {*}[lrepeat 6 6]}]} {
-    error "ARBN inside accepted thicknesses that collapse the solid"
-}
-if {![catch {get $invalid_arbn_inner}]} {
-    error "ARBN inside wrote an invalid collapsed solid"
-}
-
 puts "*** 'arb' testing completed ***"

regress/mged/csg_edit.mged

@@ -962,140 +962,8 @@ foreach test [list \
     csg_matrix_mouse_check {*}$test
 }
 
-puts "--- Repeated graphical object transforms ---"
-proc csg_matrix_repeated_mouse_check {} {
-    global orig_gui transform csg_mouse_dm_pixels csg_mouse_view_size
-    global csg_drag_scale
-    set prim csg_matrix_repeat.s
-    set inner csg_matrix_repeat_inner.c
-    set outer csg_matrix_repeat_outer.c
-    set saved_orig_gui $orig_gui
-    set expected [csg_edit_identity_matrix]
-    set failed [catch {
-        in $prim sph 0 0 0 10
-        comb $inner u $prim
-        comb $outer u $inner
-        set orig_gui 0
-        for {set step 0} {$step < 16} {incr step} {
-            e $outer
-            center 0 0 0
-            size $csg_mouse_view_size
-            setview 0 0 0
-            dm size $csg_mouse_dm_pixels $csg_mouse_dm_pixels
-            _mged_press oill
-            _mged_ill -e -i 1 /$outer/$inner/$prim
-            _mged_matpick 2
-            if {[status state] ne "OBJ EDIT" || $transform ne "e"} {
-                error "step $step did not enter object editing"
-            }
-            set view_before [csg_edit_view_state]
-            if {$step % 2 == 0} {
-                press Rotate
-                dm con r z 100 100
-                set operation [csg_matrix_axis_rotation z 10]
-                set motion {dm motion 140 100}
-            } else {
-                press Scale
-                dm am s 100 100
-                set operation [csg_matrix_scale all $csg_drag_scale]
-                set motion {dm motion 164 100}
-            }
-            set motion_code [catch {uplevel #0 $motion} motion_result]
-            if {$motion_code != 0 && $motion_code != 2} {
-                error "step $step mouse motion failed: $motion_result"
-            }
-            dm idle
-            press accept
-            set expected [csg_matrix_multiply $operation $expected]
-            csg_edit_comb_leaf_matrix_at $inner $prim $expected
-            csg_edit_view_state_at $view_before
-            d $outer
-        }
-    } failure]
-    set orig_gui $saved_orig_gui
-    catch {press reject}
-    catch {Z}
-    foreach object [list $outer $inner $prim] {
-        catch {kill $object}
-    }
-    if {$failed} {
-        puts "  FAIL: repeated object mouse rotations/scales ($failure)"
-    } else {
-        puts "  PASS: repeated object mouse rotations/scales"
-    }
-}
-csg_matrix_repeated_mouse_check
-
 puts "--- All ARB point, edge, and face states ---"
 
-proc csg_edit_sed_point {attribute} {
-    set serialized [get_sed]
-    set index [lsearch -exact $serialized $attribute]
-    if {$index < 0 || $index + 1 >= [llength $serialized]} {
-        error "get_sed has no $attribute: $serialized"
-    }
-    return [lindex $serialized [expr {$index + 1}]]
-}
-
-proc csg_arb_repeated_mouse_check {} {
-    global orig_gui csg_bv_coord_half_range csg_mouse_dm_pixels
-    global csg_mouse_view_size
-    set source csg_lifecycle_arb5.s
-    set prim csg_arb_repeat.s
-    set saved_orig_gui $orig_gui
-    set expected {}
-    set failed [catch {
-        cp $source $prim
-        e $prim
-        center 0 0 0
-        size $csg_mouse_view_size
-        setview 0 0 0
-        dm size $csg_mouse_dm_pixels $csg_mouse_dm_pixels
-        set orig_gui 0
-        sed $prim
-        press {Move Edges}
-        press {Move Point 5}
-        set original_z [lindex [db get $source V5] 2]
-        for {set step 0} {$step < 64} {incr step} {
-            if {$step > 0 && $step % 8 == 0} {
-                press sedit
-                press {Move Edges}
-                press {Move Point 5}
-            }
-            set mouse_x [expr {150 + 2 * ($step % 5)}]
-            set mouse_y [expr {150 + 2 * (($step * 3) % 5)}]
-            M 1 $mouse_x $mouse_y
-            set expected [list \
-                [expr {$mouse_x / $csg_bv_coord_half_range * $csg_mouse_view_size / 2.0}] \
-                [expr {$mouse_y / $csg_bv_coord_half_range * $csg_mouse_view_size / 2.0}] \
-                $original_z]
-            set actual [csg_edit_sed_point V5]
-            if {[csg_edit_point_distance $actual $expected] > 0.000001} {
-                error "step $step expected V5 $expected, got $actual"
-            }
-        }
-        press accept
-        d $prim
-        if {[csg_edit_point_distance [db get $prim V5] $expected] > 0.000001} {
-            error "accepted V5 differs from the final mouse target"
-        }
-        foreach vertex {V1 V2 V3 V4} {
-            csg_edit_point_matches_source $source $vertex $prim
-        }
-        arb_standard_storage_check arb5 $prim
-    } failure]
-    set orig_gui $saved_orig_gui
-    catch {press reject}
-    catch {d $prim}
-    catch {kill $prim}
-    if {$failed} {
-        puts "  FAIL: repeated ARB5 mouse movement ($failure)"
-    } else {
-        puts "  PASS: repeated ARB5 mouse movement"
-    }
-}
-csg_arb_repeated_mouse_check
-
 # Every ARB menu entry is checked independently.  The variants do not share
 # menu indices, so sampling one edge or face would miss precisely the class of
 # integration failure this test is intended to catch.  Move each selection a
@@ -1454,280 +1322,6 @@ csg_edit_sequence_check csg_lifecycle_dsp.s "Toggle Smooth" {
     {press {Toggle Smooth}}
 } [list csg_edit_scalar_at sm 1]
 
-puts "--- Specialized mouse editing and menu inventory ---"
-proc csg_edit_prepare_mouse_view {} {
-    global csg_dispatch_view_size
-    center 0 0 0
-    size $csg_dispatch_view_size
-    setview 0 0 0
-}
-array set csg_mouse_covered {}
-proc csg_mouse_case {source menu_items label commands validator} {
-    global csg_mouse_covered
-    csg_edit_sequence_check $source $label $commands $validator
-    set type [lindex [db get $source] 0]
-    if {![info exists csg_mouse_covered($type)]} {
-        set csg_mouse_covered($type) {}
-    }
-    foreach item $menu_items {
-        if {[lsearch -exact $csg_mouse_covered($type) $item] < 0} {
-            lappend csg_mouse_covered($type) $item
-        }
-    }
-}
-proc csg_edit_menu_inventory_check {source ignored} {
-    global csg_mouse_covered
-    if {[catch {
-        Z
-        e $source
-        sed $source
-        set actual [mmenu_get 0]
-        set type [lindex [db get $source] 0]
-        set covered {}
-        if {[info exists csg_mouse_covered($type)]} {
-            set covered $csg_mouse_covered($type)
-        }
-        foreach item $actual {
-            set tested [expr {[lsearch -exact $covered $item] >= 0}]
-            set exempt [expr {[lsearch -exact $ignored $item] >= 0}]
-            if {$tested == $exempt} {
-                error "menu entry '$item' needs exactly one mouse-test or exclusion"
-            }
-        }
-        foreach item [concat $covered $ignored] {
-            if {[lsearch -exact $actual $item] < 0} {
-                error "recorded menu entry '$item' is absent from $actual"
-            }
-        }
-        press reject
-        d $source
-    } failure]} {
-        puts "  FAIL: menu inventory $source ($failure)"
-        catch {press reject}
-        catch {Z}
-    } else {
-        puts "  PASS: menu inventory $source"
-    }
-}
-
-csg_edit_prepare_mouse_view
-set csg_special_mouse_target [list \
-    [expr {512.0 / $csg_bv_coord_half_range * $csg_dispatch_view_size / 2.0}] \
-    [expr {512.0 / $csg_bv_coord_half_range * $csg_dispatch_view_size / 2.0}] 0]
-set csg_special_mouse_scale [expr {1.0 + 2047.0 / ($csg_bv_coord_half_range * 4.0)}]
-set csg_special_mouse_shrink [expr {
-    1.0 / (1.0 + 1024.0 / ($csg_bv_coord_half_range * 4.0))}]
-csg_mouse_case csg_lifecycle_pipe.s {{Move Point}} "mouse select and move point" {
-    {csg_edit_prepare_mouse_view}
-    {press {Select Point}}
-    {p 0 0 0}
-    {press {Move Point}}
-    {M 1 512 512}
-} [list csg_edit_point_at V0 $csg_special_mouse_target]
-in csg_mouse_pipe_topology.s pipe 4 \
-    0 0 0 1 2 3 \
-    0 0 30 1 2 3 \
-    30 40 80 1 2 3 \
-    80 80 100 1 2 3
-csg_mouse_case csg_mouse_pipe_topology.s {{Select Point} {Move Point}} "mouse select point" {
-    {csg_edit_prepare_mouse_view}
-    {press {Select Point}}
-    {M 1 0 0}
-    {press {Move Point}}
-    {M 1 512 512}
-} [list csg_edit_point_at V0 $csg_special_mouse_target]
-csg_mouse_case csg_mouse_pipe_topology.s {{Append Point}} "mouse append point" {
-    {csg_edit_prepare_mouse_view}
-    {press {Append Point}}
-    {M 1 512 512}
-} [list csg_edit_point_at V4 {5 5 0}]
-csg_mouse_case csg_mouse_pipe_topology.s {{Prepend Point}} "mouse prepend point" {
-    {csg_edit_prepare_mouse_view}
-    {press {Prepend Point}}
-    {M 1 512 512}
-} [list csg_edit_point_at V0 {5 5 0}]
-csg_mouse_case csg_mouse_pipe_topology.s {{Split Segment}} "mouse split segment" {
-    {csg_edit_prepare_mouse_view}
-    {press {Select Point}}
-    {p 0 0 0}
-    {press {Split Segment}}
-    {M 1 512 512}
-} [list csg_edit_point_at V1 {5 5 0}]
-foreach {item attribute initial} {
-    {Set Point OD} O0 5
-    {Set Point ID} I0 3
-    {Set Point Bend} R0 6
-} {
-    csg_mouse_case csg_lifecycle_pipe.s [list $item] "mouse $item" [list \
-        {csg_edit_prepare_mouse_view} \
-        {press {Select Point}} {p 0 0 0} [list press $item] \
-        {M 1 0 -1024}] \
-        [list csg_edit_scalar_at $attribute [expr {$initial * $csg_special_mouse_shrink}]]
-}
-foreach {item attributes initial} {
-    {Set Pipe OD} {O0 O1 O2 O3} {5 5 6 6}
-    {Set Pipe ID} {I0 I1 I2 I3} {3 3 2 0}
-    {Set Pipe Bend} {R0 R1 R2 R3} {6 7 10 8}
-} {
-    set expected {}
-    foreach value $initial {
-        lappend expected [expr {$value * $csg_special_mouse_shrink}]
-    }
-    csg_mouse_case csg_lifecycle_pipe.s [list $item] "mouse $item" [list \
-        {csg_edit_prepare_mouse_view} \
-        [list press $item] {M 1 0 -1024}] \
-        [list csg_edit_scalar_values_at $attributes $expected]
-}
-proc csg_edit_bot_vertices_moved {source moved delta prim} {
-    set vertices [db get $source V]
-    for {set index 0} {$index < [llength $vertices]} {incr index} {
-        set expected [lindex $vertices $index]
-        if {[lsearch -exact $moved $index] >= 0} {
-            set translated {}
-            foreach value $expected shift $delta {
-                lappend translated [expr {$value + $shift}]
-            }
-            set expected $translated
-        }
-        csg_edit_list_value_at V $index $expected $prim
-    }
-}
-csg_mouse_case csg_lifecycle_bot.s {{Move Vertex}} "mouse move vertex" {
-    {csg_edit_prepare_mouse_view}
-    {press {Pick Vertex}}
-    {p 0}
-    {press {Move Vertex}}
-    {M 1 512 512}
-} [list csg_edit_bot_vertices_moved csg_lifecycle_bot.s {0} \
-    $csg_special_mouse_target]
-csg_mouse_case csg_lifecycle_bot.s {{Pick Vertex} {Move Vertex}} "mouse pick and move vertex" {
-    {csg_edit_prepare_mouse_view}
-    {press {Pick Vertex}}
-    {M 1 0 0}
-    {press {Move Vertex}}
-    {M 1 512 512}
-} [list csg_edit_bot_vertices_moved csg_lifecycle_bot.s {0} \
-    $csg_special_mouse_target]
-csg_mouse_case csg_lifecycle_bot.s {{Move Edge}} "mouse move edge" {
-    {csg_edit_prepare_mouse_view}
-    {press {Pick Edge}}
-    {p 0 1}
-    {press {Move Edge}}
-    {M 1 1024 512}
-} [list csg_edit_bot_vertices_moved csg_lifecycle_bot.s {0 1} {10 5 0}]
-csg_mouse_case csg_lifecycle_bot.s {{Pick Edge} {Move Edge}} "mouse pick and move edge" {
-    {csg_edit_prepare_mouse_view}
-    {press {Pick Edge}}
-    {M 1 512 512}
-    {press {Move Edge}}
-    {M 1 1024 512}
-} [list csg_edit_bot_vertices_moved csg_lifecycle_bot.s {0 1} {10 5 0}]
-csg_mouse_case csg_lifecycle_bot.s {{Move Triangle}} "mouse move triangle" {
-    {csg_edit_prepare_mouse_view}
-    {press {Pick Triangle}}
-    {p 0 1 2}
-    {press {Move Triangle}}
-    {M 1 512 512}
-} [list csg_edit_bot_vertices_moved csg_lifecycle_bot.s {0 1 2} \
-    $csg_special_mouse_target]
-csg_mouse_case csg_lifecycle_ebm.s {{Extrude Depth}} "mouse extrusion depth" {
-    {csg_edit_prepare_mouse_view}
-    {press {Extrude Depth}}
-    {M 1 0 2047}
-} [list csg_edit_scalar_at H [expr {4.0 * $csg_special_mouse_scale}]]
-csg_mouse_case csg_lifecycle_vol.s {{Threshold (low)}} "mouse low threshold" {
-    {csg_edit_prepare_mouse_view}
-    {press {Threshold (low)}}
-    {p 10}
-    {M 1 0 2047}
-} [list csg_edit_scalar_at lo [expr {int(10.0 * $csg_special_mouse_scale)}]]
-csg_mouse_case csg_lifecycle_vol.s {{Voxel Size (X Y Z)}} "mouse voxel size" {
-    {csg_edit_prepare_mouse_view}
-    {press {Voxel Size (X Y Z)}}
-    {M 1 0 2047}
-} [list csg_edit_point_at size [list \
-    $csg_special_mouse_scale $csg_special_mouse_scale $csg_special_mouse_scale]]
-csg_mouse_case csg_lifecycle_vol.s {{Threshold (hi)}} "mouse high threshold" {
-    {csg_edit_prepare_mouse_view}
-    {press {Threshold (hi)}}
-    {M 1 0 -1024}
-} [list csg_edit_scalar_at hi [expr {int(255.0 * $csg_special_mouse_shrink)}]]
-csg_mouse_case csg_lifecycle_dsp.s {{Set X}} "mouse X scale" {
-    {csg_edit_prepare_mouse_view}
-    {press {Set X}}
-    {M 1 0 2047}
-} [list csg_edit_matrix_component_at stom 0 $csg_special_mouse_scale]
-foreach {item index} {{Set Y} 5 {Set ALT} 10} {
-    csg_mouse_case csg_lifecycle_dsp.s [list $item] "mouse $item" [list \
-        {csg_edit_prepare_mouse_view} \
-        [list press $item] {M 1 0 2047}] \
-        [list csg_edit_matrix_component_at stom $index $csg_special_mouse_scale]
-}
-foreach type {arbn brep} {
-    set source csg_lifecycle_${type}.s
-    csg_edit_sequence_check $source "generic mouse translation" {
-        {csg_edit_prepare_mouse_view}
-        {keypoint 0 0 0}
-        {press sxy}
-        {M 1 512 512}
-    } [list csg_edit_bbox_translated $source $csg_special_mouse_target]
-}
-
-# Every specialized menu entry must be accounted for.  Exclusions are menu
-# headings, non-positional actions/settings, or the BOT triangle picker (tested
-# under Xmin because it needs Tk).  A new entry without a mouse case or
-# explicit exclusion fails here.
-foreach {source ignored} {
-    csg_lifecycle_pipe.s {
-        {PIPE MENU} {Next Point} {Previous Point} {Delete Point}
-    }
-    csg_lifecycle_bot.s {
-        {BOT MENU} {Pick Triangle} {Move Vertex List}
-        {Split Edge} {Split Face} {Fuse Vertices} {Fuse Faces}
-        {Delete Triangle} {Select Mode} {Select Orientation}
-        {Set flags} {Set Face Thickness} {Set Face Mode}
-    }
-    csg_lifecycle_dsp.s {
-        {DSP MENU} Name {Toggle Smooth} {Set Data Source}
-    }
-    csg_lifecycle_ebm.s {
-        {EBM MENU} {File Name} {File Size (W N)}
-    }
-    csg_lifecycle_vol.s {
-        {VOL MENU} {File Name} {File Size (X Y Z)}
-    }
-    csg_lifecycle_arbn.s {}
-    csg_lifecycle_brep.s {}
-} {
-    csg_edit_menu_inventory_check $source $ignored
-}
-
-# A view-space request must remain view-space even when the editor's current
-# coordinate system is object-relative.  This is the fast path underlying the
-# corresponding real mouse drag in the Xmin GUI test.
-set csg_saved_coords [rset var coords]
-set csg_saved_transform [rset var transform]
-if {[catch {
-    Z
-    e csg_lifecycle_bot.s
-    sed csg_lifecycle_bot.s
-    center 0 0 0
-    rset var coords o
-    rset var transform v
-    knob -i -v aX 1
-    if {[csg_edit_point_distance [center] {-1 0 0}] > 0.000001} {
-        error "explicit view knob left center at [center]"
-    }
-} csg_view_override_error]} {
-    puts "  FAIL: view-space translation with object coordinates ($csg_view_override_error)"
-} else {
-    puts "  PASS: view-space translation with object coordinates"
-}
-catch {press reject}
-catch {Z}
-rset var coords $csg_saved_coords
-rset var transform $csg_saved_transform
-
 file delete -force $csg_data_file
 file delete -force $csg_dsp_file
 

regress/mged/extrude.mged

@@ -43,37 +43,4 @@ foreach {arb_type face expected_type} $arb_extrude_cases {
     arb_standard_storage_check $expected_type $prim
 }
 
-# A one-inch extrusion must extend a face by 25.4 base millimeters through
-# both the direct libged command and MGED's interactive solid editor.
-set inch_in_mm 25.4
-set tolerance 0.001
-foreach mode {direct interactive} {
-    set suffix "_inch_$mode"
-    in_arb8 extrude $suffix
-    set prim "extrude_arb8${suffix}.s"
-    set base_x [lindex [db get $prim V1] 0]
-
-    units in
-    if {$mode eq "direct"} {
-	if {[catch {edarb extrude $prim 1234 1} result]} {
-	    error "inch edarb extrusion failed: $result"
-	}
-    } else {
-	Z
-	e $prim
-	sed $prim
-	if {[catch {extrude 1234 1} result]} {
-	    error "inch MGED extrusion failed: $result"
-	}
-	accept
-	Z
-    }
-    units mm
-
-    set new_x [lindex [db get $prim V5] 0]
-    if {abs($new_x - $base_x - $inch_in_mm) > $tolerance} {
-	error "inch $mode extrusion moved V5 by [expr {$new_x - $base_x}] mm"
-    }
-}
-
 puts "*** 'extrude' testing completed ***"

regress/mged/prim_edit.mged

@@ -31,7 +31,6 @@ in tor_pe.s tor 0 0 0  0 0 1  20 2
 prim_edit_check      tor_pe.s "Set Radius 1" 30 r_a 30 -1
 prim_edit_check      tor_pe.s "Set Radius 2"  5 r_h  5 -1
 prim_edit_mouse_check tor_pe.s "Set Radius 1" r_a 30 -1
-prim_edit_mouse_check tor_pe.s "Set Radius 2" r_h 5 -1
 kill tor_pe.s
 
 puts "--- ELL (ellipsoid: A, B, C) ---"
@@ -40,8 +39,6 @@ prim_edit_check      ell_pe.s "Set A" 8 A 8 0
 prim_edit_check      ell_pe.s "Set B" 6 B 6 1
 prim_edit_check      ell_pe.s "Set C" 4 C 4 2
 prim_edit_mouse_check ell_pe.s "Set A" A 8 0
-prim_edit_mouse_check ell_pe.s "Set B" B 6 1
-prim_edit_mouse_check ell_pe.s "Set C" C 4 2
 kill ell_pe.s
 
 puts "--- RPC (right parabolic cylinder: H, r) ---"
@@ -49,7 +46,6 @@ in rpc_pe.s rpc 0 0 0  0 0 10  0 5 0  3
 prim_edit_check      rpc_pe.s "Set H" 20 H 20 2
 prim_edit_check      rpc_pe.s "Set r"  6 r  6 -1
 prim_edit_mouse_check rpc_pe.s "Set H" H 20 2
-prim_edit_mouse_check rpc_pe.s "Set r" r 6 -1
 kill rpc_pe.s
 
 puts "--- RHC (right hyperbolic cylinder: H, r, c) ---"
@@ -58,8 +54,6 @@ prim_edit_check      rhc_pe.s "Set H" 20 H 20 2
 prim_edit_check      rhc_pe.s "Set r"  6 r  6 -1
 prim_edit_check      rhc_pe.s "Set c"  3 c  3 -1
 prim_edit_mouse_check rhc_pe.s "Set r" r 6 -1
-prim_edit_mouse_check rhc_pe.s "Set H" H 20 2
-prim_edit_mouse_check rhc_pe.s "Set c" c 3 -1
 kill rhc_pe.s
 
 puts "--- EPA (elliptical paraboloid: H, r_1, r_2) ---"
@@ -68,8 +62,6 @@ prim_edit_check      epa_pe.s "Set H" 20 H   20 2
 prim_edit_check      epa_pe.s "Set A"  8 r_1  8 -1
 prim_edit_check      epa_pe.s "Set B"  4 r_2  4 -1
 prim_edit_mouse_check epa_pe.s "Set H" H 20 2
-prim_edit_mouse_check epa_pe.s "Set A" r_1 8 -1
-prim_edit_mouse_check epa_pe.s "Set B" r_2 4 -1
 kill epa_pe.s
 
 puts "--- EHY (elliptical hyperboloid: H, r_1, c) ---"
@@ -78,8 +70,6 @@ prim_edit_check      ehy_pe.s "Set H" 20 H   20 2
 prim_edit_check      ehy_pe.s "Set A"  8 r_1  8 -1
 prim_edit_check      ehy_pe.s "Set c"  3 c    3 -1
 prim_edit_mouse_check ehy_pe.s "Set A" r_1 8 -1
-prim_edit_mouse_check ehy_pe.s "Set H" H 20 2
-prim_edit_mouse_check ehy_pe.s "Set c" c 3 -1
 kill ehy_pe.s
 
 puts "--- PART (particle: H, r_v, r_h) ---"
@@ -88,8 +78,6 @@ prim_edit_check      part_pe.s "Set H" 20 H   20 2
 prim_edit_check      part_pe.s "Set v"  4 r_v  4 -1
 prim_edit_check      part_pe.s "Set h"  3 r_h  3 -1
 prim_edit_mouse_check part_pe.s "Set v" r_v 4 -1
-prim_edit_mouse_check part_pe.s "Set H" H 20 2
-prim_edit_mouse_check part_pe.s "Set h" r_h 3 -1
 kill part_pe.s
 
 puts "--- TGC (truncated general cone: H, A, B, C, D) ---"
@@ -103,10 +91,6 @@ prim_edit_check      tgc_pe.s "Set B"  3 B   3 1
 prim_edit_check      tgc_pe.s "Set C"  3 C   3 0
 prim_edit_check      tgc_pe.s "Set D"  4 D   4 1
 prim_edit_mouse_check tgc_pe.s "Set H"  H 8 2
-prim_edit_mouse_check tgc_pe.s "Set A"  A 4 0
-prim_edit_mouse_check tgc_pe.s "Set B"  B 3 1
-prim_edit_mouse_check tgc_pe.s "Set C"  C 3 0
-prim_edit_mouse_check tgc_pe.s "Set D"  D 4 1
 kill tgc_pe.s
 
 puts "--- ETO (elliptical torus: r, C, D) ---"
@@ -116,7 +100,6 @@ prim_edit_check      eto_pe.s "Set r"  8   r   8   -1
 prim_edit_check      eto_pe.s "Set D"  1.5 r_d 1.5 -1
 prim_edit_mouse_check eto_pe.s "Set r"  r   8   -1
 prim_edit_mouse_check eto_pe.s "Set C"  C   5   0
-prim_edit_mouse_check eto_pe.s "Set D"  r_d 1.5 -1
 kill eto_pe.s
 
 puts "--- HYP (one-sheet hyperboloid: H, A, B, c) ---"
@@ -128,9 +111,6 @@ prim_edit_check      hyp_pe.s "Set A"  2   A   10   1
 prim_edit_check      hyp_pe.s "Set B"  2   b    4   -1
 prim_edit_check      hyp_pe.s "Set c"  0.5 bnr  0.25 -1
 prim_edit_mouse_check hyp_pe.s "Set H"  H  20   2
-prim_edit_mouse_check hyp_pe.s "Set A"  A  10   1
-prim_edit_mouse_check hyp_pe.s "Set B"  b  4    -1
-prim_edit_mouse_check hyp_pe.s "Set c"  bnr 0.25 -1
 kill hyp_pe.s
 
 puts "--- SUPERELL (superellipsoid: A, B, C) ---"
@@ -139,8 +119,6 @@ prim_edit_check      superell_pe.s "Set A" 5 A 5 0
 prim_edit_check      superell_pe.s "Set B" 4 B 4 1
 prim_edit_check      superell_pe.s "Set C" 3 C 3 2
 prim_edit_mouse_check superell_pe.s "Set A" A 5 0
-prim_edit_mouse_check superell_pe.s "Set B" B 4 1
-prim_edit_mouse_check superell_pe.s "Set C" C 3 2
 kill superell_pe.s
 
 puts "--- CLINE (cylinder line: H, R, plate thickness) ---"
@@ -151,7 +129,6 @@ prim_edit_check      cline_pe.s "Set R"              4   R 4   -1
 prim_edit_check      cline_pe.s "Set plate thickness" 0.4 T 0.4 -1
 prim_edit_mouse_check cline_pe.s "Set H" H 8   2
 prim_edit_mouse_check cline_pe.s "Set R" R 4   -1
-prim_edit_mouse_check cline_pe.s "Set plate thickness" T 0.4 -1
 kill cline_pe.s
 
 puts "--- NMG (non-manifold geometry: Pick Edge, Move Edge, Next EU) ---"

regress/mged/sed_reset.mged

@@ -1,64 +0,0 @@
-puts "*** Testing 'sed_reset' command ***"
-
-proc sed_reset_check {primitive edit_commands} {
-    e $primitive
-    sed $primitive
-    if {[status state] ne "SOL EDIT"} {
-	error "$primitive entered [status state], not SOL EDIT"
-    }
-
-    set original [get_sed]
-    foreach command $edit_commands {
-	uplevel #0 $command
-    }
-    if {[get_sed] eq $original} {
-	error "$primitive edit did not change the working geometry"
-    }
-
-    sed_reset
-    if {[status state] ne "SOL EDIT"} {
-	error "$primitive reset entered [status state], not SOL EDIT"
-    }
-    if {[get_sed] ne $original} {
-	error "$primitive reset did not restore the original geometry"
-    }
-
-    press reject
-    if {[status state] ne "VIEWING"} {
-	error "$primitive reject entered [status state], not VIEWING"
-    }
-    d $primitive
-}
-
-if {[catch {
-    in sed_reset_arb8.s arb8 \
-	0 0 0  10 0 0  10 10 0  0 10 0 \
-	0 0 10  10 0 10  10 10 10  0 10 10
-    sed_reset_check sed_reset_arb8.s {
-	{press {Move Faces}}
-	{press {Move Face 1234}}
-	{p 0 0 -1}
-    }
-} failure]} {
-    puts "FAIL: ARB8 face move reset ($failure)"
-    catch {press reject}
-    catch {Z}
-} else {
-    puts "PASS: ARB8 face move reset"
-}
-
-if {[catch {
-    make sed_reset_cline.s cline
-    sed_reset_check sed_reset_cline.s {
-	{press {Move End H}}
-	{p 3 4 12}
-    }
-} failure]} {
-    puts "FAIL: CLINE edit reset ($failure)"
-    catch {press reject}
-    catch {Z}
-} else {
-    puts "PASS: CLINE edit reset"
-}
-
-puts "PASS: 'sed_reset' testing completed"

regress/xmin/CMakeLists.txt

@@ -93,7 +93,6 @@ if(BRLCAD_X11_PROVIDER_RESOLVED STREQUAL "XMIN" AND
       "PLOT3_ASC_BIN=$<TARGET_FILE:plot3-asc>"
     )
     mged_add_xmin_gui_fixture(help "mged;dm-tkswrast;m35.g")
-    mged_add_xmin_gui_fixture(edit-interaction "mged;dm-tkswrast;m35.g")
     mged_add_xmin_gui_fixture(bot "mged;dm-tkswrast;m35.g")
     mged_add_xmin_gui_fixture(editors "mged;dm-tkswrast;m35.g")
     mged_add_xmin_gui_fixture(query-pattern "mged;dm-tkswrast;m35.g")
@@ -144,18 +143,6 @@ if(BRLCAD_X11_PROVIDER_RESOLVED STREQUAL "XMIN" AND
       )
       brlcad_regression_test(regress-archer-embedded-xmin-gui
         "archer;bwish" STAND_ALONE TEST_DEFINED)
-
-      brlcad_add_gui_test(
-        NAME regress-archer-geometry-xmin-gui
-        TIMEOUT 120
-        ENVIRONMENT ${_gui_smoke_environment}
-          "GUI_TEST_LIBRARY=${_gui_test_library}"
-        COMMAND "${SH_EXEC}" "${_gui_tcl_runner}" archer-geometry
-          $<TARGET_FILE:bwish>
-          "${CMAKE_CURRENT_SOURCE_DIR}/archer_geometry_xmin_gui.tcl"
-      )
-      brlcad_regression_test(regress-archer-geometry-xmin-gui
-        "archer;bwish" STAND_ALONE TEST_DEFINED)
     endif()
   endif()
 
@@ -213,7 +200,6 @@ cmakefiles(
   archer_action.sh
   archer_smoke.sh
   archer_embedded_xmin_gui.tcl
-  archer_geometry_xmin_gui.tcl
   archer_menu_manifest.txt
   archer_xmin_gui.tcl
   qged_action.sh
@@ -226,19 +212,18 @@ cmakefiles(
   mged/mged_xmin_gui.sh
   mged/xmin_gui_bot.tcl
   mged/xmin_gui_controls.tcl
-  mged/xmin_gui_edit.tcl
-  mged/xmin_gui_edit_interaction.tcl
   mged/xmin_gui_editors.tcl
   mged/xmin_gui_help.tcl
   mged/xmin_gui_lod.tcl
   mged/xmin_gui_overlaps.tcl
-  mged/xmin_gui_pipe_oracle.tcl
   mged/xmin_gui_query_pattern.tcl
   mged/xmin_gui_render_outputs.tcl
   mged/xmin_gui_search_exec.tcl
   mged/xmin_gui_shotvis.tcl
   mged/xmin_gui_test.tcl
   mged/xmin_gui_workflows.tcl
+  mged/xmin_gui_edit.tcl
+  mged/xmin_gui_pipe_oracle.tcl
 )
 
 # Local Variables:

regress/xmin/archer_geometry_xmin_gui.tcl

@@ -1,134 +0,0 @@
-#    A R C H E R _ G E O M E T R Y _ X M I N _ G U I . T C L
-# BRL-CAD
-#
-# Copyright (c) 2026 United States Government as represented by
-# the U.S. Army Research Laboratory.
-#
-# This library is free software; you can redistribute it and/or
-# modify it under the terms of the GNU Lesser General Public License
-# version 2.1 as published by the Free Software Foundation.
-#
-# This library is distributed in the hope that it will be useful, but
-# WITHOUT ANY WARRANTY; without even the implied warranty of
-# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
-# Lesser General Public License for more details.
-#
-# You should have received a copy of the GNU Lesser General Public
-# License along with this file; see the file named COPYING for more
-# information.
-#
-###
-if {![info exists ::env(GUI_TEST_LIBRARY)]} {
-    puts stderr "GUI_TEST_LIBRARY is required"
-    exit 2
-}
-source $::env(GUI_TEST_LIBRARY)
-set ::env(ARCHER_PREFS_FILE) [file join $::env(GUI_TEST_DIR) .archerrc]
-
-proc archer_geometry_preferences {{geometry ""}} {
-    set channel [open $::env(ARCHER_PREFS_FILE) w]
-    if {$geometry ne ""} {
-	puts $channel [list set mWindowGeometry $geometry]
-    }
-    close $channel
-}
-
-proc archer_geometry_saved {} {
-    set channel [open $::env(ARCHER_PREFS_FILE) r]
-    set contents [read $channel]
-    close $channel
-    if {[regexp -line {^set mWindowGeometry (.+)$} $contents -> geometry]} {
-	return $geometry
-    }
-    return ""
-}
-
-proc archer_geometry_settle {} {
-    # Allow MapNotify and the deferred geometry restore to reach the event loop.
-    set settle_delay_ms 100
-    after $settle_delay_ms [list set ::archer_geometry_settled 1]
-    vwait ::archer_geometry_settled
-    update
-    ::gui::test::require_no_background_errors
-}
-
-if {[catch {
-    package require Archer 1.0
-    ::gui::test::capture_background_errors
-    wm withdraw .
-
-    set saved_geometry 900x650+120+90
-    set startup_geometry 1000x700+20+30
-    set updated_geometry 950x680+60+70
-    set default_geometry 920x660+40+50
-
-    archer_geometry_preferences $saved_geometry
-    Archer .unmapped
-    ::gui::test::require {[wm state .unmapped] eq "withdrawn"} \
-	"startup fixture unexpectedly showed its window"
-    ::itcl::delete object .unmapped
-    archer_geometry_settle
-    ::gui::test::require {[archer_geometry_saved] eq $saved_geometry} \
-	"closing an unmapped Archer replaced its saved geometry"
-
-    archer_geometry_preferences
-    Archer .unmapped_new
-    ::itcl::delete object .unmapped_new
-    archer_geometry_settle
-    ::gui::test::require {[archer_geometry_saved] eq ""} \
-	"closing an unmapped Archer saved an unfinished layout size"
-
-    archer_geometry_preferences $saved_geometry
-    Archer .restore
-    .restore configure -geometry $startup_geometry
-    wm geometry .restore $startup_geometry
-    event generate [.restore component hpane] <Map>
-    archer_geometry_settle
-    ::gui::test::require {[wm geometry .restore] eq $startup_geometry} \
-	"a child Map event restored the main window's geometry prematurely"
-
-    wm deiconify .restore
-    archer_geometry_settle
-    ::gui::test::require {[wm geometry .restore] eq $saved_geometry} \
-	"mapping the main window did not restore its saved geometry"
-
-    wm geometry .restore $updated_geometry
-    event generate [.restore component hpane] <Map>
-    archer_geometry_settle
-    ::gui::test::require {[wm geometry .restore] eq $updated_geometry} \
-	"a later child Map event reset the resized main window"
-    wm withdraw .restore
-    ::itcl::delete object .restore
-    archer_geometry_settle
-    ::gui::test::require {[archer_geometry_saved] eq $updated_geometry} \
-	"closing a previously mapped window did not save its latest geometry"
-
-    archer_geometry_preferences
-    Archer .default
-    wm geometry .default $default_geometry
-    wm deiconify .default
-    archer_geometry_settle
-    ::gui::test::require {[wm geometry .default] eq $default_geometry} \
-	"an empty geometry option reset the window to its layout size"
-    wm withdraw .default
-    ::itcl::delete object .default
-    archer_geometry_settle
-} message options]} {
-    puts stderr $message
-    if {[dict exists $options -errorinfo]} {
-	puts stderr [dict get $options -errorinfo]
-    }
-    puts "FAIL: Archer window geometry regression"
-    exit 1
-}
-
-puts "PASS: Archer window geometry restoration and startup cancellation"
-exit 0
-
-# Local Variables:
-# tab-width: 8
-# mode: Tcl
-# tcl-indent-level: 4
-# indent-tabs-mode: t
-# End:
-# ex: shiftwidth=4 tabstop=8

regress/xmin/mged/xmin_gui_help.tcl

@@ -50,18 +50,6 @@ proc ::mged::xmin::help::exercise_command_manual {top} {
     ::gui::test::require {
 	[llength $selection] == 1 && [$toc get $selection] eq "rt"
     } "man rt did not select and load the rt manual page"
-
-    event generate $toc <<ListboxSelect>>
-    ::mged::gui::test::settle
-
-    man search
-    set html_view [$browser component browser].htmlview
-    set html_widget [$html_view html]
-    $html_view goto {#hierarchy_performance}
-    ::mged::gui::test::settle
-    ::gui::test::require {[lindex [$html_widget yview] 0] > 0.0} \
-	"manual browser did not follow an in-page hyperlink"
-
     ::gui::test::require {
 	[catch {man xmin_page_that_does_not_exist} message] &&
 	[string first "couldn't find manual page" $message] >= 0
@@ -163,38 +151,10 @@ proc ::mged::xmin::help::exercise_html_manual {top} {
     set mged_browser [file join $::env(GUI_TEST_DIR) missing-browser]
     ::mged::gui::test::invoke $top {Help Manual}
     set fallback $top.man
-    set viewer $fallback.html
     ::gui::test::require {
 	[winfo exists $fallback] && [winfo ismapped $fallback] &&
-	[winfo exists $viewer] && [$viewer title] eq "Table of Contents"
+	[string length [string trim [$fallback.text get 1.0 end]]] > 100
     } "Manual did not load MGED's internal HTML viewer"
-
-    set special_path [file join [pwd] {manual path} {contents#1.html}]
-    ::gui::test::require {
-	[::mged::manual::local_path \
-	    [::mged::manual::file_uri $special_path]] eq \
-	[file normalize $special_path]
-    } "internal Manual viewer did not preserve special path characters"
-
-    set html [$viewer html]
-    set contents [$html text text]
-    ::gui::test::require {
-	[string first "Preface" $contents] >= 0 &&
-	[string first "background-color" $contents] < 0
-    } "internal Manual viewer did not render its table of contents"
-
-    $viewer goto mged.html#csg
-    ::mged::gui::test::settle
-    set contents [$html text text]
-    ::gui::test::require {
-	[$viewer title] eq "Mged User's Manual" &&
-	[string first "Modeling With CSG" $contents] >= 0
-    } "internal Manual viewer did not follow its document links"
-
-    $fallback.f.back invoke
-    ::mged::gui::test::settle
-    ::gui::test::require {[$viewer title] eq "Table of Contents"} \
-	"internal Manual viewer Back button did not restore the prior page"
     destroy $fallback
 }
 

src/adrt/CMakeLists.txt

@@ -53,7 +53,7 @@ set(
 )
 
 brlcad_addlib(librender "${LIBRENDER_SRCS}" "" "${ADRT_LOCAL_INCLUDE_DIRS}"
-  PRIVATE_LIBS libgcv ZLIB::ZLIB ${M_LIBRARY}
+  PRIVATE_LIBS libgcv librt ZLIB::ZLIB ${M_LIBRARY}
 )
 set_target_properties(librender PROPERTIES VERSION 20.0.1 SOVERSION 20)
 add_dependencies(librender gcv_plugins)
@@ -109,8 +109,9 @@ cmakefiles(isst.c)
 
 add_subdirectory(master)
 add_subdirectory(slave)
+add_subdirectory(tests)
 
-set(adrt_ignore_files load.h tienet.h)
+set(adrt_ignore_files load.h load_g.h tienet.h)
 cmakefiles(${adrt_ignore_files})
 cmakefiles(CMakeLists.txt)
 

src/adrt/adrt.h

@@ -90,11 +90,19 @@ enum
 
 #define ADRT_NAME_SIZE 256
 
+/* TIE geometry is normalized from millimeters to meters while loading. */
+#define ADRT_MODEL_UNITS_PER_TIE_UNIT 1000.0
+
 #define ADRT_LOAD_FORMAT_G 42
 #define ADRT_LOAD_FORMAT_REG 43
 #define ADRT_LOAD_FORMAT_KDTREE 44
 
+struct rt_annot_scene;
+
 RENDER_EXPORT extern int load_g(struct tie_s *tie, const char *db, int argc, const char **argv, struct adrt_mesh_s **);
+RENDER_EXPORT extern int load_g_annotations(struct tie_s *tie, const char *db,
+	int argc, const char **argv, struct adrt_mesh_s **,
+	struct rt_annot_scene **annotations);
 
 #endif
 

src/adrt/isst.c

@@ -59,6 +59,7 @@ struct isst_s {
     struct render_camera_s camera;
     struct camera_tile_s tile;
     struct adrt_mesh_s *meshes;
+    struct rt_annot_scene *annotations;
     tienet_buffer_t buffer_image;
     int ogl, sflags, w, h, gs, ui;
     double dt, fps, uic;
@@ -67,6 +68,7 @@ struct isst_s {
     void *texdata;
     vect_t camera_pos_init;
     vect_t camera_focus_init;
+    fastf_t scene_radius;
     int64_t t1;
     int64_t t2;
     int dirty;
@@ -119,13 +121,7 @@ resize_isst(struct isst_s *isstp)
     glTexParameteri (GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
     glTexParameteri (GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
     glTexEnvf (GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_MODULATE);
-    {
-	/* realloc through a temporary so a failure does not leak (and null out)
-	 * the existing buffer. */
-	void *tmp = realloc(isstp->texdata, isstp->camera.w * isstp->camera.h * 3);
-	if (tmp)
-	    isstp->texdata = tmp;
-    }
+    isstp->texdata = realloc(isstp->texdata, isstp->camera.w * isstp->camera.h * 3);
     glTexImage2D (GL_TEXTURE_2D, 0, GL_RGB, isstp->camera.w, isstp->camera.h, 0, GL_RGB, GL_UNSIGNED_BYTE, isstp->texdata);
     glDisable(GL_LIGHTING);
 
@@ -158,11 +154,14 @@ isst_load_g(ClientData UNUSED(clientData), Tcl_Interp *interp, int objc,
     argv = (char **)malloc(sizeof(char *) * (strlen(Tcl_GetString(objv[3])) + 1));	/* allocate way too much. */
     argc = (int)bu_argv_from_string(argv, strlen(Tcl_GetString(objv[3])), Tcl_GetString(objv[3]));
 
-    load_g(isst->tie, Tcl_GetString(objv[2]), argc, (const char **)argv, &(isst->meshes));
+    rt_annot_scene_destroy(isst->annotations);
+    isst->annotations = NULL;
+    load_g_annotations(isst->tie, Tcl_GetString(objv[2]), argc,
+	(const char **)argv, &(isst->meshes), &isst->annotations);
     free(argv);
 
-    VSETALL(isst->camera.pos, isst->tie->radius);
-    VMOVE(isst->camera.focus, isst->tie->mid);
+    isst->scene_radius = render_camera_fit_scene(&isst->camera, isst->tie,
+	isst->annotations, ADRT_MODEL_UNITS_PER_TIE_UNIT);
     VMOVE(isst->camera_pos_init, isst->camera.pos);
     VMOVE(isst->camera_focus_init, isst->camera.focus);
 
@@ -239,7 +238,9 @@ paint_window(ClientData UNUSED(clientData), Tcl_Interp *interp, int objc, Tcl_Ob
 	isst->buffer_image.ind = 0;
 
 	render_camera_prep(&isst->camera);
-	render_camera_render(&isst->camera, isst->tie, &isst->tile, &isst->buffer_image);
+	render_camera_render_annotations(&isst->camera, isst->tie, &isst->tile,
+	    &isst->buffer_image, isst->annotations,
+	    ADRT_MODEL_UNITS_PER_TIE_UNIT);
 
 	isst->t1 = bu_gettime();
 
@@ -282,7 +283,6 @@ set_resolution(ClientData UNUSED(clientData), Tcl_Interp *interp, int objc, Tcl_
     }
 
     CLAMP(resolution, 1, 20);
-    /* gs is the render grid size; 0 means render at full window resolution */
     if (resolution == 20)
 	isst->gs = 0;
     else
@@ -326,13 +326,10 @@ isst_zap(ClientData UNUSED(clientData), Tcl_Interp *interp, int objc, Tcl_Obj *c
 	return TCL_ERROR;
     }
 
-    /* Release the buffers owned directly by the workspace before the struct
-     * itself.  (The tie and mesh list have a more involved lifecycle and are
-     * left for a dedicated teardown.) */
     TIENET_BUFFER_FREE(isst->buffer_image);
     if (isst->texdata)
 	free(isst->texdata);
-
+    rt_annot_scene_destroy(isst->annotations);
     bu_free(isst, "isst free");
     isst = NULL;
 
@@ -354,7 +351,7 @@ render_mode(ClientData UNUSED(clientData), Tcl_Interp *interp, int objc, Tcl_Obj
     if (objc == 4)
 	buf = Tcl_GetString(objv[3]);
 
-    /* pack the 'rest' into buf - probably should use a vls for this */
+    /* pack the 'rest' into buf - probably should use a vls for this*/
     if ( strlen(mode) == 3 && bu_strncmp("cut", mode, 3) == 0 ) {
 	struct adrt_mesh_s *mesh;
 
@@ -391,26 +388,22 @@ zero_view(ClientData UNUSED(clientData), Tcl_Interp *UNUSED(interp), int UNUSED(
 
 
 static int
-move_walk(ClientData UNUSED(clientData), Tcl_Interp *interp, int objc, Tcl_Obj *const *objv)
+move_walk(ClientData UNUSED(clientData), Tcl_Interp *interp, int UNUSED(objc), Tcl_Obj *const *objv)
 {
     vect_t vec;
     int flag;
 
-    if (objc < 3) {
-	Tcl_WrongNumArgs(interp, 1, objv, "pathName flag");
-	return TCL_ERROR;
-    }
     if (Tcl_GetIntFromObj(interp, objv[2], &flag) != TCL_OK)
 	return TCL_ERROR;
 
     if (flag >= 0) {
 	VSUB2(vec, isst->camera.focus, isst->camera.pos);
-	VSCALE(vec, vec, 0.1 * isst->tie->radius);
+	VSCALE(vec, vec, 0.1 * isst->scene_radius);
 	VADD2(isst->camera.pos, isst->camera.pos, vec);
 	VADD2(isst->camera.focus, isst->camera.focus, vec);
     } else {
 	VSUB2(vec, isst->camera.pos, isst->camera.focus);
-	VSCALE(vec, vec, 0.1 * isst->tie->radius);
+	VSCALE(vec, vec, 0.1 * isst->scene_radius);
 	VADD2(isst->camera.pos, isst->camera.pos, vec);
 	VADD2(isst->camera.focus, isst->camera.focus, vec);
     }
@@ -419,16 +412,12 @@ move_walk(ClientData UNUSED(clientData), Tcl_Interp *interp, int objc, Tcl_Obj *
 }
 
 static int
-move_strafe(ClientData UNUSED(clientData), Tcl_Interp *interp, int objc, Tcl_Obj *const *objv)
+move_strafe(ClientData UNUSED(clientData), Tcl_Interp *interp, int UNUSED(objc), Tcl_Obj *const *objv)
 {
     vect_t vec, dir, up;
 
     int flag;
 
-    if (objc < 3) {
-	Tcl_WrongNumArgs(interp, 1, objv, "pathName flag");
-	return TCL_ERROR;
-    }
     if (Tcl_GetIntFromObj(interp, objv[2], &flag) != TCL_OK)
 	return TCL_ERROR;
 
@@ -437,13 +426,13 @@ move_strafe(ClientData UNUSED(clientData), Tcl_Interp *interp, int objc, Tcl_Obj
     if (flag >= 0) {
 	VSUB2(dir, isst->camera.focus, isst->camera.pos);
 	VCROSS(vec, dir, up);
-	VSCALE(vec, vec, 0.1 * isst->tie->radius);
+	VSCALE(vec, vec, 0.1 * isst->scene_radius);
 	VADD2(isst->camera.pos, isst->camera.pos, vec);
 	VADD2(isst->camera.focus, isst->camera.pos, dir);
     } else {
 	VSUB2(dir, isst->camera.focus, isst->camera.pos);
 	VCROSS(vec, dir, up);
-	VSCALE(vec, vec, -0.1 * isst->tie->radius);
+	VSCALE(vec, vec, -0.1 * isst->scene_radius);
 	VADD2(isst->camera.pos, isst->camera.pos, vec);
 	VADD2(isst->camera.focus, isst->camera.pos, dir);
     }
@@ -551,14 +540,10 @@ aerotate(ClientData UNUSED(clientData), Tcl_Interp *interp, int objc, Tcl_Obj *c
 	VSCALE(vecdfoc, vecdfoc, mag_focus);
 	VADD2(isst->camera.focus, isst->camera_focus_init, vecdfoc);
     }
-    /* Update the tcl copies of the az/el vars.  Store them in the same units
-     * as isst_load_g (radians via -DEG2RAD); previously aerotate wrote the raw
-     * degrees from bn_ae_vec, leaving the Tcl vars inconsistent after a drag. */
+    /* Update the tcl copies of the az/el vars */
     VSUB2(vec, isst->camera.focus, isst->camera.pos);
     VUNITIZE(vec);
     bn_ae_vec(&az, &el, vec);
-    az = az * -DEG2RAD;
-    el = el * -DEG2RAD;
     bu_vls_sprintf(&tclstr, "%f", az);
     Tcl_SetVar(interp, "az", bu_vls_addr(&tclstr), 0);
     bu_vls_sprintf(&tclstr, "%f", el);
@@ -576,7 +561,7 @@ open_dm(ClientData UNUSED(cdata), Tcl_Interp *interp, int UNUSED(objc), Tcl_Obj
     dmp = dm_open(NULL, (void *)interp, dm_default_type(), sizeof(av)/sizeof(void*)-1, (const char **)av);
 
     if (dmp == DM_NULL) {
-	printf("Failed to open display manager\n");
+	printf("dm failed?\n");
 	return TCL_ERROR;
     }
 
@@ -654,7 +639,7 @@ const char *fullname;
     argv = __argv;
 #endif
 
-    /* initialize progname for run-time resource finding */
+    /* initialize progname for run-tim resource finding */
     bu_setprogname(argv[0]);
 
 #ifdef HAVE_WINDOWS_H

src/adrt/librender/camera.c

@@ -29,6 +29,7 @@
 #include "bu/dylib.h"
 #include "bu/log.h"
 #include "bu/str.h"
+#include "rt/primitives/annot.h"
 
 #include "./camera.h"
 
@@ -43,13 +44,27 @@ struct render_shader_s {
 
 static struct render_shader_s *shaders = NULL;
 
-/* number of leading characters compared when matching a shader by name */
-#define RENDER_SHADER_NAME_CMP_LEN	8
+struct render_camera_annotation_data {
+    const struct rt_annot_scene *scene;
+    fastf_t model_units_per_tie_unit;
+    struct rt_annot_view view;
+};
+
+
+struct render_camera_parallel_data {
+    render_camera_thread_data_t thread;
+    struct render_camera_annotation_data annotation;
+};
+
+static const fastf_t render_camera_default_scene_radius = 1.0;
+
 
 void render_camera_render_thread(int cpu, void *ptr);	/* for bu_parallel */
 static void render_camera_prep_ortho(render_camera_t *camera);
 static void render_camera_prep_persp(render_camera_t *camera);
 static void render_camera_prep_persp_dof(render_camera_t *camera);
+static void render_camera_annotation_view(const render_camera_t *camera,
+	fastf_t model_units_per_tie_unit, struct rt_annot_view *view);
 
 static struct render_shader_s *render_shader_register (const char *name, int (*init)(render_t *, const char *));
 
@@ -393,20 +408,147 @@ render_camera_prep(render_camera_t *camera)
 }
 
 
-void
-render_camera_render_thread(int UNUSED(cpu), void *ptr)
+fastf_t
+render_camera_fit_scene(render_camera_t *camera, const struct tie_s *tie,
+	const struct rt_annot_scene *annotations,
+	fastf_t model_units_per_tie_unit)
+{
+    point_t annotation_min, annotation_max;
+    point_t scene_min, scene_max;
+    vect_t offset;
+    fastf_t scene_radius;
+
+    if (!camera || !tie)
+	return 0.0;
+
+    if (!annotations || !rt_annot_scene_bounds(annotations,
+	    annotation_min, annotation_max)) {
+	VSETALL(camera->pos, tie->radius);
+	VMOVE(camera->focus, tie->mid);
+	return tie->radius;
+    }
+
+    if (!isfinite(model_units_per_tie_unit) ||
+	    model_units_per_tie_unit <= 0.0)
+	model_units_per_tie_unit = 1.0;
+    VSCALE(annotation_min, annotation_min,
+	1.0 / model_units_per_tie_unit);
+    VSCALE(annotation_max, annotation_max,
+	1.0 / model_units_per_tie_unit);
+    VMOVE(scene_min, annotation_min);
+    VMOVE(scene_max, annotation_max);
+    if (tie->tri_num) {
+	VMINMAX(scene_min, scene_max, tie->amin);
+	VMINMAX(scene_min, scene_max, tie->amax);
+    }
+
+    VADD2SCALE(camera->focus, scene_min, scene_max, 0.5);
+    VSUB2(offset, scene_max, scene_min);
+    scene_radius = 0.5 * MAGNITUDE(offset);
+    if (!isfinite(scene_radius) || scene_radius <= SMALL_FASTF)
+	scene_radius = render_camera_default_scene_radius;
+
+    /* Keep the established isometric initial direction while positioning it
+     * relative to the combined scene center. */
+    VSETALL(offset, scene_radius);
+    VADD2(camera->pos, camera->focus, offset);
+    return scene_radius;
+}
+
+
+static void
+render_camera_annotation_view(const render_camera_t *camera,
+	fastf_t model_units_per_tie_unit, struct rt_annot_view *view)
+{
+    vect_t look, up, side, temp, right;
+    point_t model_position;
+    fastf_t angle, sine, cosine, xy_scale = 1.0;
+    mat_t *matrix = &view->model2view;
+
+    VSUB2(look, camera->focus, camera->pos);
+    VUNITIZE(look);
+    VSET(up, 0.0, 0.0, 1.0);
+    VMOVE(temp, look);
+    angle = VDOT(up, temp);
+    VSCALE(temp, temp, angle);
+    VSUB2(up, up, temp);
+    VUNITIZE(up);
+    VCROSS(side, up, look);
+    sine = sin(-camera->tilt * DEG2RAD);
+    cosine = cos(-camera->tilt * DEG2RAD);
+    VSCALE(up, up, cosine);
+    VSCALE(side, side, sine);
+    VADD2(up, up, side);
+    VCROSS(side, up, look);
+    VREVERSE(right, side);
+
+    MAT_IDN(*matrix);
+    if (camera->type == RENDER_CAMERA_ORTHOGRAPHIC &&
+	    camera->aspect > 0.0 && camera->gridsize > 0.0) {
+	xy_scale = 2.0 / (camera->aspect * camera->gridsize *
+		model_units_per_tie_unit);
+    } else if (camera->aspect > 0.0) {
+	xy_scale = 1.0 / camera->aspect;
+    }
+    for (size_t axis = 0; axis < 3; ++axis) {
+	(*matrix)[axis] = right[axis] * xy_scale;
+	(*matrix)[4 + axis] = up[axis] * xy_scale;
+	(*matrix)[8 + axis] = -look[axis];
+    }
+    VSCALE(model_position, camera->pos, model_units_per_tie_unit);
+    (*matrix)[3] = -VDOT(&(*matrix)[0], model_position);
+    (*matrix)[7] = -VDOT(&(*matrix)[4], model_position);
+    (*matrix)[11] = -VDOT(&(*matrix)[8], model_position);
+    view->width = camera->w;
+    view->height = camera->h;
+    view->perspective =
+	camera->type == RENDER_CAMERA_PERSPECTIVE ? 2.0 * camera->fov : 0.0;
+}
+
+
+static void
+render_camera_composite_annotation(render_camera_thread_data_t *td,
+	const struct render_camera_annotation_data *annotation,
+	const struct tie_ray_s *tie_ray, fastf_t sample_x, fastf_t sample_y,
+	vect_t *pixel)
+{
+    struct rt_annot_hit back_layer = {0};
+    struct xray ray;
+    fastf_t scene_distance = INFINITY;
+
+    if (!annotation || !annotation->scene)
+	return;
+    VSCALE(ray.r_pt, tie_ray->pos, annotation->model_units_per_tie_unit);
+    VMOVE(ray.r_dir, tie_ray->dir);
+    if (!rt_annot_scene_query_layers(annotation->scene, &annotation->view,
+	    &ray, sample_x, sample_y, INFINITY, &back_layer, 1))
+	return;
+
+    if (!back_layer.screen_space) {
+	struct tie_id_s id;
+	struct tie_ray_s depth_ray = *tie_ray;
+	if (TIE_WORK(td->tie, &depth_ray, &id, render_hit, NULL))
+	    scene_distance = id.dist * annotation->model_units_per_tie_unit;
+    }
+    (void)rt_annot_scene_composite(annotation->scene, &annotation->view,
+	&ray, sample_x, sample_y, scene_distance, *pixel, NULL);
+}
+
+
+static void
+render_camera_render_thread_impl(int UNUSED(cpu), render_camera_thread_data_t *td,
+	const struct render_camera_annotation_data *annotation)
 {
-    render_camera_thread_data_t *td;
     int d, n, res_ind, scanline, v_scanline;
     vect_t pixel, accum;
     vect_t v1 = VINIT_ZERO;
     vect_t v2 = VINIT_ZERO;
     struct tie_ray_s ray;
+    struct tie_ray_s annotation_ray;
     fastf_t view_inv;
 
     VSETALL(v1, 0);
 
-    td = (render_camera_thread_data_t *)ptr;
     view_inv = 1.0 / td->camera->view_num;
 
     td->camera->render.tie = td->tie;
@@ -457,9 +599,15 @@ render_camera_render_thread(int UNUSED(cpu), void *ptr)
 		    VMOVE(ray.pos, td->camera->view_list[d].pos);
 		    ray.depth = 0;
 		    VUNITIZE(ray.dir);
+		    annotation_ray = ray;
 
 		    /* Compute pixel value using this ray */
 		    td->camera->render.work(&td->camera->render, td->tie, &ray, &pixel);
+		    render_camera_composite_annotation(td, annotation,
+			&annotation_ray,
+			(fastf_t)n + 0.5,
+			(fastf_t)td->camera->h - (fastf_t)v_scanline - 0.5,
+			&pixel);
 
 		    VADD2(accum, accum, pixel);
 		}
@@ -476,6 +624,7 @@ render_camera_render_thread(int UNUSED(cpu), void *ptr)
 		    VMOVE(ray.pos, td->camera->view_list[0].pos);
 		    ray.depth = 0;
 		    VUNITIZE(ray.dir);
+		    annotation_ray = ray;
 
 		    /* Compute pixel value using this ray */
 		    td->camera->render.work(&td->camera->render, td->tie, &ray, &pixel);
@@ -490,12 +639,19 @@ render_camera_render_thread(int UNUSED(cpu), void *ptr)
 
 		    VSET(pixel, (TFLOAT)RENDER_CAMERA_BGR, (TFLOAT)RENDER_CAMERA_BGG, (TFLOAT)RENDER_CAMERA_BGB);
 		    ray.depth = 0;
+		    annotation_ray = ray;
 
 		    /* Compute pixel value using this ray */
 		    td->camera->render.work(&td->camera->render, td->tie, &ray, &pixel);
 		}
 	    }
 
+	    if (td->camera->view_num == 1)
+		render_camera_composite_annotation(td, annotation,
+		    &annotation_ray, (fastf_t)n + 0.5,
+		    (fastf_t)td->camera->h - (fastf_t)v_scanline - 0.5,
+		    &pixel);
+
 
 	    if (td->tile->format == RENDER_CAMERA_BIT_DEPTH_24) {
 		V_MIN(pixel[0], 1);
@@ -523,9 +679,30 @@ render_camera_render_thread(int UNUSED(cpu), void *ptr)
 
 
 void
-render_camera_render(render_camera_t *camera, struct tie_s *tie, camera_tile_t *tile, tienet_buffer_t *result)
+render_camera_render_thread(int cpu, void *ptr)
 {
-    render_camera_thread_data_t td;
+    render_camera_render_thread_impl(cpu, (render_camera_thread_data_t *)ptr,
+	NULL);
+}
+
+
+static void
+render_camera_render_annotations_thread(int cpu, void *ptr)
+{
+    struct render_camera_parallel_data *data =
+	(struct render_camera_parallel_data *)ptr;
+
+    render_camera_render_thread_impl(cpu, &data->thread, &data->annotation);
+}
+
+
+static void
+render_camera_render_internal(render_camera_t *camera, struct tie_s *tie,
+	camera_tile_t *tile, tienet_buffer_t *result,
+	const struct rt_annot_scene *annotations,
+	fastf_t model_units_per_tie_unit)
+{
+    struct render_camera_parallel_data data = {0};
     unsigned int scanline;
     uint32_t ind;
 
@@ -543,19 +720,47 @@ render_camera_render(render_camera_t *camera, struct tie_s *tie, camera_tile_t *
     TCOPY(camera_tile_t, tile, 0, result->data, result->ind);
     result->ind += sizeof(camera_tile_t);
 
-    td.tie = tie;
-    td.camera = camera;
-    td.tile = tile;
-    td.res_buf = &((char *)result->data)[result->ind];
+    data.thread.tie = tie;
+    data.thread.camera = camera;
+    data.thread.tile = tile;
+    data.thread.res_buf = &((char *)result->data)[result->ind];
     scanline = 0;
-    td.scanline = &scanline;
-    td.sem_tie_worker = bu_semaphore_register("sem_tie_worker");
+    data.thread.scanline = &scanline;
+    data.thread.sem_tie_worker = bu_semaphore_register("sem_tie_worker");
+    data.annotation.scene = annotations;
+    data.annotation.model_units_per_tie_unit =
+	model_units_per_tie_unit > 0.0 ? model_units_per_tie_unit : 1.0;
+    if (annotations)
+	render_camera_annotation_view(camera,
+	    data.annotation.model_units_per_tie_unit, &data.annotation.view);
 
-    bu_parallel(render_camera_render_thread, camera->thread_num, &td);
+    if (annotations)
+	bu_parallel(render_camera_render_annotations_thread,
+	    camera->thread_num, &data);
+    else
+	bu_parallel(render_camera_render_thread, camera->thread_num,
+	    &data.thread);
 
     result->ind = ind;
+}
 
-    return;
+
+void
+render_camera_render(render_camera_t *camera, struct tie_s *tie,
+	camera_tile_t *tile, tienet_buffer_t *result)
+{
+    render_camera_render_internal(camera, tie, tile, result, NULL, 1.0);
+}
+
+
+void
+render_camera_render_annotations(render_camera_t *camera, struct tie_s *tie,
+	camera_tile_t *tile, tienet_buffer_t *result,
+	const struct rt_annot_scene *annotations,
+	fastf_t model_units_per_tie_unit)
+{
+    render_camera_render_internal(camera, tie, tile, result, annotations,
+	model_units_per_tie_unit);
 }
 
 
@@ -620,9 +825,9 @@ render_shader_unload_plugin(render_t *r, const char *name)
 {
 #ifdef HAVE_DLFCN_H
     struct render_shader_s *t, *s = shaders, *meh;
-    if (!bu_strncmp(s->name, name, RENDER_SHADER_NAME_CMP_LEN)) {
+    if (!bu_strncmp(s->name, name, 8)) {
 	t = s->next;
-	if (r && r->shader && !bu_strncmp(r->shader, name, RENDER_SHADER_NAME_CMP_LEN)) {
+	if (r && r->shader && !bu_strncmp(r->shader, name, 8)) {
 	    meh = s->next;
 	    while (meh) {
 		if (render_shader_init(r, meh->name, NULL) != -1)
@@ -641,7 +846,7 @@ LOADED:
     }
 
     while (s->next) {
-	if (!bu_strncmp(s->next->name, name, RENDER_SHADER_NAME_CMP_LEN)) {
+	if (!bu_strncmp(s->next->name, name, 8)) {
 	    if (r)
 		render_shader_init(r, s->name, NULL);
 	    if (s->next->dlh)
@@ -651,7 +856,6 @@ LOADED:
 	    bu_free(t, "unload shader");
 	    return 0;
 	}
-	s = s->next;
     }
 
     bu_log("Could not find shader \"%s\"\n", name);
@@ -667,7 +871,7 @@ render_shader_init(render_t *r, const char *name, const char *buf)
 {
     struct render_shader_s *s = shaders;
     while (s) {
-	if (!bu_strncmp(s->name, name, RENDER_SHADER_NAME_CMP_LEN)) {
+	if (!bu_strncmp(s->name, name, 8)) {
 	    s->init(r, buf);
 	    r->shader = s->name;
 	    return 0;

src/adrt/librender/camera.h

@@ -99,7 +99,18 @@ __BEGIN_DECLS
 RENDER_EXPORT extern void render_camera_init(render_camera_t *camera, size_t threads);
 RENDER_EXPORT extern void render_camera_free(render_camera_t *camera);
 RENDER_EXPORT extern void render_camera_prep(render_camera_t *camera);
+/** Center a camera on TIE geometry and any prepared annotations.  Returns the
+ * effective scene radius in TIE coordinates. */
+RENDER_EXPORT extern fastf_t render_camera_fit_scene(render_camera_t *camera,
+	const struct tie_s *tie, const struct rt_annot_scene *annotations,
+	fastf_t model_units_per_tie_unit);
 RENDER_EXPORT extern void render_camera_render(render_camera_t *camera, struct tie_s *tie, camera_tile_t *tile, tienet_buffer_t *result);
+/** Render with an externally owned annotation scene.  The coordinate scale
+ * converts TIE coordinates to database model coordinates. */
+RENDER_EXPORT extern void render_camera_render_annotations(render_camera_t *camera,
+	struct tie_s *tie, camera_tile_t *tile, tienet_buffer_t *result,
+	const struct rt_annot_scene *annotations,
+	fastf_t model_units_per_tie_unit);
 
 RENDER_EXPORT extern int render_shader_init(render_t *, const char *name, const char *buf);
 RENDER_EXPORT extern const char *render_shader_load_plugin(const char *filename);

src/adrt/load_g.c

@@ -40,11 +40,13 @@
 #include "vmath.h"
 #include "nmg.h"
 #include "rt/geom.h"
+#include "rt/search.h"
 #include "raytrace.h"
 
 #include "rt/tie.h"
 #include "adrt.h"
 #include "adrt_struct.h"
+#include "load_g.h"
 
 
 static struct bn_tol tol;		/* calculation tolerance */
@@ -104,6 +106,42 @@ struct gcv_data {
 static struct gcv_data gcvwriter = {{load_nmg_to_adrt_gcvwrite, NULL, NULL}, NULL};
 
 
+int
+adrt_path_has_geometry(struct db_i *database, const char *path)
+{
+    static const char *geometry_search = "! -type comb ! -type annot";
+    struct directory *root;
+
+    if (!database || !path)
+	return 0;
+    root = db_lookup(database, path, LOOKUP_QUIET);
+
+    /* Let db_walk_tree report a missing path rather than hiding the error. */
+    if (!root)
+	return 1;
+
+    return db_search(NULL, DB_SEARCH_QUIET, geometry_search, 1, &root,
+	database, NULL, NULL, NULL) != 0;
+}
+
+
+union tree *
+adrt_leaf_tess(struct db_tree_state *tsp, const struct db_full_path *pathp,
+	struct rt_db_internal *ip, void *client_data)
+{
+    union tree *result;
+
+    if (!ip || ip->idb_type != ID_ANNOT)
+	return rt_booltree_leaf_tess(tsp, pathp, ip, client_data);
+
+    /* Annotations are rendered directly and have no NMG representation. */
+    BU_GET(result, union tree);
+    RT_TREE_INIT(result);
+    result->tr_op = OP_NOP;
+    return result;
+}
+
+
 /* load the region into the tie image */
 static void
 load_nmg_to_adrt_internal(struct adrt_mesh_s *mesh, struct nmgregion *r)
@@ -153,7 +191,8 @@ load_nmg_to_adrt_internal(struct adrt_mesh_s *mesh, struct nmgregion *r)
 		    NMG_CK_VERTEX(v);
 
 		    /* convert mm to m */
-		    VSCALE((*tribuf[vert_count]).v, v->vg_p->coord, 1.0/1000.0);
+		    VSCALE((*tribuf[vert_count]).v, v->vg_p->coord,
+			1.0 / ADRT_MODEL_UNITS_PER_TIE_UNIT);
 		    vert_count++;
 		}
 		if (vert_count > 3)
@@ -248,9 +287,12 @@ load_nmg_to_adrt_regstart(struct db_tree_state *ts, const struct db_full_path *p
 
 	for (i=0;i<bot->num_faces;i++)
 	{
-	    VSCALE((*tribuf[0]).v, (bot->vertices+3*bot->faces[3*i+0]), 1.0/1000.0);
-	    VSCALE((*tribuf[1]).v, (bot->vertices+3*bot->faces[3*i+1]), 1.0/1000.0);
-	    VSCALE((*tribuf[2]).v, (bot->vertices+3*bot->faces[3*i+2]), 1.0/1000.0);
+	    VSCALE((*tribuf[0]).v, (bot->vertices+3*bot->faces[3*i+0]),
+		1.0 / ADRT_MODEL_UNITS_PER_TIE_UNIT);
+	    VSCALE((*tribuf[1]).v, (bot->vertices+3*bot->faces[3*i+1]),
+		1.0 / ADRT_MODEL_UNITS_PER_TIE_UNIT);
+	    VSCALE((*tribuf[2]).v, (bot->vertices+3*bot->faces[3*i+2]),
+		1.0 / ADRT_MODEL_UNITS_PER_TIE_UNIT);
 
 	    TIE_VAL(tie_push)(cur_tie, tribuf, 1, mesh, 0);
 	}
@@ -298,11 +340,16 @@ load_nmg_to_adrt_gcvwrite(struct nmgregion *r, const struct db_full_path *pathp,
 
 
 int
-load_g(struct tie_s *tie, const char *db, int argc, const char **argv, struct adrt_mesh_s **meshes)
+load_g_annotations(struct tie_s *tie, const char *db, int argc,
+	const char **argv, struct adrt_mesh_s **meshes,
+	struct rt_annot_scene **annotations)
 {
     struct model *the_model;
     struct bg_tess_tol ttol;		/* tessellation tolerance in mm */
     struct db_tree_state tree_state;	/* includes tol & model */
+    const char **geometry_paths = NULL;
+    int geometry_count = 0;
+    int walk_result = 0;
 
     cur_tie = tie;	/* NOTE: cur_tie is a file-scope global; this is not
 			 * reentrant and would need locking to be thread-safe. */
@@ -338,13 +385,21 @@ load_g(struct tie_s *tie, const char *db, int argc, const char **argv, struct ad
     if ((dbip = db_open(db, DB_OPEN_READONLY)) == DBI_NULL) {
 	perror(db);
 	bu_log("Unable to open geometry database file (%s)\n", db);
+	nmg_km(the_model);
 	return -1;
     }
     if (db_dirbuild(dbip)) {
 	bu_log("ERROR: db_dirbuild failed\n");
+	db_close(dbip);
+	dbip = DBI_NULL;
+	nmg_km(the_model);
 	return -1;
     }
 
+    if (annotations) {
+	*annotations = rt_annot_scene_create(dbip, argc, argv, &ttol, &tol);
+    }
+
     BN_CK_TOL(tree_state.ts_tol);
     BG_CK_TESS_TOL(tree_state.ts_ttol);
 
@@ -362,15 +417,26 @@ load_g(struct tie_s *tie, const char *db, int argc, const char **argv, struct ad
     tribuf[1] = (TIE_3 *)bu_malloc(sizeof(TIE_3) * 3, "triangle tribuffer");
     tribuf[2] = (TIE_3 *)bu_malloc(sizeof(TIE_3) * 3, "triangle tribuffer");
 
-    (void) db_walk_tree(dbip,
-			argc,			/* number of toplevel regions */
-			argv,			/* region names */
+    if (argc > 0) {
+	geometry_paths = (const char **)bu_calloc((size_t)argc,
+	    sizeof(const char *), "ADRT geometry paths");
+	for (int i = 0; i < argc; ++i) {
+	    if (adrt_path_has_geometry(dbip, argv[i]))
+		geometry_paths[geometry_count++] = argv[i];
+	}
+    }
+
+    if (geometry_count) {
+	walk_result = db_walk_tree(dbip,
+			geometry_count,		/* number of toplevel regions */
+			geometry_paths,		/* region names */
 			1,			/* ncpu */
 			&tree_state,		/* initial tree state */
 			load_nmg_to_adrt_regstart,	/* region start function */
 			gcv_region_end,		/* region end function */
-			rt_booltree_leaf_tess,	/* leaf func */
+			adrt_leaf_tess,		/* leaf func */
 			(void *)&gcvwriter);	/* client data */
+    }
 
     /* Release dynamic storage */
     nmg_km(the_model);
@@ -380,10 +446,20 @@ load_g(struct tie_s *tie, const char *db, int argc, const char **argv, struct ad
     bu_free(tribuf[1], "vert");
     bu_free(tribuf[2], "vert");
     bu_free(tribuf, "tri");
+    if (geometry_paths)
+	bu_free(geometry_paths, "ADRT geometry paths");
 
     TIE_VAL(tie_prep)(cur_tie);
 
-    return 0;
+    return walk_result;
+}
+
+
+int
+load_g(struct tie_s *tie, const char *db, int argc, const char **argv,
+	struct adrt_mesh_s **meshes)
+{
+    return load_g_annotations(tie, db, argc, argv, meshes, NULL);
 }
 
 /*

src/adrt/load_g.h

@@ -0,0 +1,43 @@
+/*                      L O A D _ G . H
+ * BRL-CAD / ADRT
+ *
+ * Copyright (c) 2026 United States Government as represented by
+ * the U.S. Army Research Laboratory.
+ *
+ * SPDX-License-Identifier: LGPL-2.1-or-later
+ */
+
+#ifndef ADRT_LOAD_G_H
+#define ADRT_LOAD_G_H
+
+#include "common.h"
+
+#include "adrt.h"
+
+struct db_i;
+struct db_tree_state;
+struct db_full_path;
+struct rt_db_internal;
+union tree;
+
+__BEGIN_DECLS
+
+RENDER_EXPORT extern int adrt_path_has_geometry(struct db_i *dbip,
+	const char *path);
+RENDER_EXPORT extern union tree *adrt_leaf_tess(struct db_tree_state *tsp,
+	const struct db_full_path *pathp, struct rt_db_internal *ip,
+	void *client_data);
+
+__END_DECLS
+
+#endif /* ADRT_LOAD_G_H */
+
+/*
+ * Local Variables:
+ * tab-width: 8
+ * mode: C
+ * indent-tabs-mode: t
+ * c-file-style: "stroustrup"
+ * End:
+ * ex: shiftwidth=4 tabstop=8
+ */

src/adrt/tests/CMakeLists.txt

@@ -0,0 +1,16 @@
+brlcad_addexec(adrt_annotation_test annotation.c
+  "librender;libwdb;librt;libbu;${M_LIBRARY}" TEST)
+target_include_directories(adrt_annotation_test PRIVATE
+  ${BRLCAD_SOURCE_DIR}/src/adrt
+  ${BRLCAD_SOURCE_DIR}/src/adrt/librender
+)
+brlcad_add_test(NAME adrt_annotation_test COMMAND adrt_annotation_test)
+
+cmakefiles(annotation.c CMakeLists.txt)
+
+# Local Variables:
+# tab-width: 8
+# mode: cmake
+# indent-tabs-mode: t
+# End:
+# ex: shiftwidth=2 tabstop=8

src/adrt/tests/annotation.c

@@ -0,0 +1,179 @@
+/*                 A N N O T A T I O N . C
+ * BRL-CAD
+ *
+ * Copyright (c) 2026 United States Government as represented by
+ * the U.S. Army Research Laboratory.
+ *
+ * SPDX-License-Identifier: LGPL-2.1-or-later
+ */
+
+#include "common.h"
+
+#include <math.h>
+#include <stdio.h>
+
+#include "bu/app.h"
+#include "bu/file.h"
+#include "raytrace.h"
+#include "rt/primitives/annot.h"
+#include "wdb.h"
+
+#include "load_g.h"
+#include "librender/camera.h"
+
+
+static void
+write_test_database(struct db_i *dbip)
+{
+    struct rt_wdb *wdbp = wdb_dbopen(dbip, RT_WDB_TYPE_DB_DISK);
+    struct rt_annot_internal annot = {0};
+    point2d_t vertices[2];
+    struct line_seg line = {0};
+    void *segments[1];
+    int reverse[1] = {0};
+    point_t center = VINIT_ZERO;
+    struct wmember members;
+
+    if (wdbp == RT_WDB_NULL)
+	bu_exit(1, "unable to open ADRT annotation test database\n");
+    annot.magic = RT_ANNOT_INTERNAL_MAGIC;
+    VSET(annot.V, 10.0, 0.0, 0.0);
+    annot.flags = RT_ANNOT_MODEL_SPACE;
+    VSET(annot.u_vec, 0.0, 1.0, 0.0);
+    VSET(annot.v_vec, 0.0, 0.0, 1.0);
+    annot.vert_count = 2;
+    annot.verts = vertices;
+    V2SET(vertices[0], -1.0, 0.0);
+    V2SET(vertices[1], 1.0, 0.0);
+    annot.ant.count = 1;
+    annot.ant.reverse = reverse;
+    annot.ant.segments = segments;
+    line.magic = CURVE_LSEG_MAGIC;
+    line.start = 0;
+    line.end = 1;
+    segments[0] = &line;
+    if (rt_annot_validate(&annot, NULL) ||
+	    mk_annot(wdbp, "test.annot", &annot) < 0 ||
+	    mk_sph(wdbp, "test.sph", center, 1.0) < 0)
+	bu_exit(1, "unable to write ADRT annotation test objects\n");
+
+    BU_LIST_INIT(&members.l);
+    if (!mk_addmember("test.annot", &members.l, NULL, WMOP_UNION) ||
+	    mk_lcomb(wdbp, "annotation-only.c", &members, 0, NULL, NULL,
+		NULL, 0))
+	bu_exit(1, "unable to write annotation-only test combination\n");
+
+    BU_LIST_INIT(&members.l);
+    if (!mk_addmember("test.annot", &members.l, NULL, WMOP_UNION) ||
+	    !mk_addmember("test.sph", &members.l, NULL, WMOP_UNION) ||
+	    mk_lcomb(wdbp, "mixed.c", &members, 0, NULL, NULL, NULL, 0))
+	bu_exit(1, "unable to write mixed test combination\n");
+}
+
+
+static void
+check_geometry_filter(struct db_i *dbip)
+{
+    struct rt_db_internal intern;
+    union tree *leaf;
+
+    if (adrt_path_has_geometry(dbip, "test.annot") ||
+	    adrt_path_has_geometry(dbip, "annotation-only.c") ||
+	    !adrt_path_has_geometry(dbip, "test.sph") ||
+	    !adrt_path_has_geometry(dbip, "mixed.c"))
+	bu_exit(1, "ADRT annotation geometry filtering failed\n");
+
+    RT_DB_INTERNAL_INIT(&intern);
+    intern.idb_type = ID_ANNOT;
+    leaf = adrt_leaf_tess(NULL, NULL, &intern, NULL);
+    if (!leaf || leaf->tr_op != OP_NOP)
+	bu_exit(1, "ADRT annotation leaf did not become a no-op\n");
+    db_free_tree(leaf);
+}
+
+
+static void
+check_camera_fit(struct db_i *dbip)
+{
+    const char *paths[] = {"test.annot"};
+    struct bg_tess_tol ttol = BG_TESS_TOL_INIT_TOL;
+    struct bn_tol tol = BN_TOL_INIT_TOL;
+    struct rt_annot_scene *scene = rt_annot_scene_create(dbip, 1, paths,
+	&ttol, &tol);
+    render_camera_t camera = {0};
+    struct tie_s tie = {0};
+    fastf_t radius;
+    const fastf_t combined_radius = 0.5 * sqrt(41.0);
+
+    if (!scene || rt_annot_scene_count(scene) != 1)
+	bu_exit(1, "unable to prepare ADRT annotation test scene\n");
+    tie.radius = 2.0;
+    VSET(tie.mid, 1.0, 2.0, 3.0);
+    radius = render_camera_fit_scene(&camera, &tie, NULL, 1.0);
+    if (!NEAR_EQUAL(radius, 2.0, BN_TOL_DIST) ||
+	    !VNEAR_EQUAL(camera.focus, tie.mid, BN_TOL_DIST) ||
+	    !NEAR_EQUAL(camera.pos[X], 2.0, BN_TOL_DIST) ||
+	    !NEAR_EQUAL(camera.pos[Y], 2.0, BN_TOL_DIST) ||
+	    !NEAR_EQUAL(camera.pos[Z], 2.0, BN_TOL_DIST))
+	bu_exit(1, "ADRT geometry-only camera framing changed\n");
+
+    radius = render_camera_fit_scene(&camera, &tie, scene, 10.0);
+    if (radius <= SMALL_FASTF ||
+	    !NEAR_EQUAL(camera.focus[X], 1.0, BN_TOL_DIST) ||
+	    !NEAR_ZERO(camera.focus[Y], BN_TOL_DIST) ||
+	    !NEAR_ZERO(camera.focus[Z], BN_TOL_DIST) ||
+	    !NEAR_EQUAL(camera.pos[X] - camera.focus[X], radius, BN_TOL_DIST) ||
+	    !NEAR_EQUAL(camera.pos[Y] - camera.focus[Y], radius, BN_TOL_DIST) ||
+	    !NEAR_EQUAL(camera.pos[Z] - camera.focus[Z], radius, BN_TOL_DIST))
+	bu_exit(1, "ADRT annotation-only camera framing failed: radius=%g "
+	    "focus=(%g %g %g) pos=(%g %g %g)\n", radius,
+	    V3ARGS(camera.focus), V3ARGS(camera.pos));
+
+    tie.tri_num = 1;
+    VSET(tie.amin, -2.0, -2.0, -2.0);
+    VSET(tie.amax, 0.0, 2.0, 2.0);
+    radius = render_camera_fit_scene(&camera, &tie, scene, 10.0);
+    if (!NEAR_EQUAL(radius, combined_radius, BN_TOL_DIST) ||
+	    !NEAR_EQUAL(camera.focus[X], -0.5, BN_TOL_DIST) ||
+	    !NEAR_ZERO(camera.focus[Y], BN_TOL_DIST) ||
+	    !NEAR_ZERO(camera.focus[Z], BN_TOL_DIST))
+	bu_exit(1, "ADRT combined geometry and annotation framing failed\n");
+    rt_annot_scene_destroy(scene);
+}
+
+
+int
+main(int argc, char **argv)
+{
+    char path[MAXPATHLEN];
+    struct db_i *dbip;
+
+    bu_setprogname(argv[0]);
+    if (argc != 1)
+	return 1;
+    {
+	FILE *fp = bu_temp_file(path, MAXPATHLEN);
+	if (!fp)
+	    bu_exit(1, "unable to create ADRT annotation test path\n");
+	fclose(fp);
+    }
+    dbip = db_create(path, 5);
+    if (dbip == DBI_NULL)
+	bu_exit(1, "unable to create ADRT annotation test database\n");
+    write_test_database(dbip);
+    check_geometry_filter(dbip);
+    check_camera_fit(dbip);
+    db_close(dbip);
+    bu_file_delete(path);
+    return 0;
+}
+
+/*
+ * Local Variables:
+ * tab-width: 8
+ * mode: C
+ * indent-tabs-mode: t
+ * c-file-style: "stroustrup"
+ * End:
+ * ex: shiftwidth=4 tabstop=8
+ */

src/archer/CMakeLists.txt

@@ -53,41 +53,30 @@ endif(BRLCAD_ENABLE_TCL)
 # Archer plugins
 set(
   archer_utility_FILES
-  plugins/Utility/README
   plugins/Utility/attrGroupsDisplayUtilityP.tcl
-  plugins/Utility/attrGroupsDisplayUtilityP/AttrGroupsDisplayUtilityP.tcl
   plugins/Utility/botUtilityP.tcl
-  plugins/Utility/botUtilityP/BotUtilityP.tcl
   plugins/Utility/lodUtilityP.tcl
-  plugins/Utility/lodUtilityP/LODUtilityP.tcl
-)
-set(archer_wizard_FILES
-  plugins/Wizards/humanwizard.tcl
-  plugins/Wizards/humanwizard/HumanWizard.tcl
-  plugins/Wizards/tankwizard.tcl
-  plugins/Wizards/tankwizard/TankWizard.tcl
-  plugins/Wizards/tankwizard/images/tank.png
-  plugins/Wizards/tirewizard.tcl
-  plugins/Wizards/tirewizard/TireWizard.tcl
+  plugins/Utility/README
 )
+set(archer_wizard_FILES plugins/Wizards/humanwizard.tcl plugins/Wizards/tankwizard.tcl plugins/Wizards/tirewizard.tcl)
 
 if(BRLCAD_ENABLE_TCL)
-  # Archer scans only the plugin category directories for registration scripts.
-  # Keep class definitions in their subdirectories so they are not sourced twice.
-  foreach(plugin_file IN LISTS archer_utility_FILES archer_wizard_FILES)
-    get_filename_component(plugin_dir "${plugin_file}" DIRECTORY)
-    file(RELATIVE_PATH plugin_dir "${CMAKE_CURRENT_SOURCE_DIR}/plugins" "${CMAKE_CURRENT_SOURCE_DIR}/${plugin_dir}")
-    brlcad_adddata("${plugin_file}" "plugins/archer/${plugin_dir}")
-  endforeach()
+  brlcad_adddata(archer_utility_FILES plugins/archer/Utility)
+  brlcad_adddata(archer_wizard_FILES plugins/archer/Wizards)
   brlcad_adddata(plugins/Commands/README plugins/archer/Command)
   brlcad_adddata(plugins/Core/README plugins/archer/Core)
+  brlcad_adddata(plugins/Utility/attrGroupsDisplayUtilityP/AttrGroupsDisplayUtilityP.tcl plugins/archer/Utility/attrGroupsDisplayUtilityP)
+  brlcad_adddata(plugins/Utility/botUtilityP/BotUtilityP.tcl plugins/archer/Utility/botUtilityP)
+  brlcad_adddata(plugins/Utility/lodUtilityP/LODUtilityP.tcl plugins/archer/Utility/lodUtilityP)
+  brlcad_adddata(plugins/Wizards/humanwizard/HumanWizard.tcl plugins/archer/Wizards/humanwizard)
+  brlcad_adddata(plugins/Wizards/tankwizard/TankWizard.tcl plugins/archer/Wizards/tankwizard)
+  brlcad_adddata(plugins/Wizards/tankwizard/images/tank.png plugins/archer/Wizards/tankwizardIA/images)
+  brlcad_adddata(plugins/Wizards/tirewizard/TireWizard.tcl plugins/archer/Wizards/tirewizard)
 endif(BRLCAD_ENABLE_TCL)
 
 set(
   archer_ignore_files
   CMakeLists.txt
-  plugins/Commands/README
-  plugins/Core/README
   ${archer_utility_FILES}
   ${archer_wizard_FILES}
   ${archer_txt}

src/art/tile.cpp

@@ -64,16 +64,17 @@ ArtTileCallback::on_tile_end(const renderer::Frame* frame, const size_t tile_x,
 {
     foundation::Tile& t = frame->image().tile(tile_x, tile_y);
     const foundation::Tile rgb(t, PixelFormatUInt8);
+    const foundation::CanvasProperties& properties = frame->image().properties();
     // fb_write(fbp, tile_x, tile_y, rgb.get_storage(), rgb.get_size());
     // printf("yay!\n");
     // printf("%lu %lu \n", tile_x, tile_y);
     // printf("%lu %lu \n", rgb.get_width(), rgb.get_height());
     for (size_t y = 0; y < rgb.get_height(); y++) {
 	for (size_t x = 0; x < rgb.get_width(); x++) {
-	    size_t x_coord = tile_x * rgb.get_width() + x;
-	    size_t y_coord = tile_y * rgb.get_height() + y;
-	    size_t img_h = frame->image().properties().m_canvas_height;
-	    fb_write(fbp, (int)x_coord, (int)(img_h - y_coord), rgb.get_storage()+((y * rgb.get_width() * 4) + (x * 4)), 1);
+	    size_t x_coord = tile_x * properties.m_tile_width + x;
+	    size_t y_coord = tile_y * properties.m_tile_height + y;
+	    size_t img_h = properties.m_canvas_height;
+	    fb_write(fbp, (int)x_coord, (int)(img_h - y_coord - 1), rgb.get_storage()+((y * rgb.get_width() * 4) + (x * 4)), 1);
 	}
     }
 }

src/external/Creo/part.cpp

@@ -752,7 +752,7 @@ tessellate_part(struct creo_conv_info *cinfo, ProMdl model, struct bu_vls **snam
             }
 
         /* Check solidity */
-        success = !bg_trimesh_solid((int)vert_tree->curr_vert, (size_t)faces.size() / 3,
+        success = bg_trimesh_solid((int)vert_tree->curr_vert, (size_t)faces.size() / 3,
                                         vert_tree->the_array, &faces[0], NULL);
 
         /* If it's not solid and we're testing solidity, keep trying... */

src/gtools/gist/InformationGatherer.cpp

@@ -23,20 +23,10 @@
 #include "InformationGatherer.h"
 #include "bu/log.h"
 
-// flush results, if we have any, then exit
-static void reportAndExit(const std::string &message, const std::string &result = "")
-{
-    if (!result.empty())
-	bu_log("%s%s", result.c_str(), result.back() == '\n' ? "" : "\n");
-
-    bu_exit(BRLCAD_ERROR, "%s", message.c_str());
-}
-
-
 std::string getCmdPath(std::string exeDir, const char* cmd) {
     char buf[MAXPATHLEN] = {0};
     if (!bu_dir(buf, MAXPATHLEN, exeDir.c_str(), cmd, BU_DIR_EXT, NULL))
-	reportAndExit("Couldn't find " + std::string(cmd) + ", aborting.\n");
+        bu_exit(BRLCAD_ERROR, "Couldn't find %s, aborting.\n", cmd);
 
     return std::string(buf);
 }
@@ -68,13 +58,13 @@ getSurfaceArea(Options* opt, std::map<std::string, std::string> UNUSED(map), std
     bu_process_create(&p, rtarea_av, BU_PROCESS_HIDE_WINDOW | BU_PROCESS_OUT_EQ_ERR);
 
     if (bu_process_pid(p) <= 0) {
-	reportAndExit("Problem with getSurfaceArea, aborting\n");
+        bu_exit(BRLCAD_ERROR, "Problem with getSurfaceArea, aborting\n");
     }
 
     char buffer[128];
-    std::string result;
+    std::string result = "";
     int read_cnt = 0;
-    while ((read_cnt = bu_process_read_n(p, BU_PROCESS_STDOUT, sizeof(buffer) - 1, buffer)) > 0) {
+    while ((read_cnt = bu_process_read_n(p, BU_PROCESS_STDOUT, 128-1, buffer)) > 0) {
         /* NOTE: read does not ensure null-termination, thus buffersize-1 */
         buffer[read_cnt] = '\0';
         result += buffer;
@@ -90,7 +80,7 @@ getSurfaceArea(Options* opt, std::map<std::string, std::string> UNUSED(map), std
         try {
             surfArea += stod(result);
         } catch (const std::invalid_argument& ia) {
-	    reportAndExit("Invalid argument for getSurfaceArea. Got (" + std::string(ia.what()) + "), aborting.\n", result);
+            bu_exit(BRLCAD_ERROR, "Invalid argument for getSurfaceArea. Got (%s), aborting.\n", ia.what());
         }
     }
 }
@@ -210,7 +200,7 @@ getVerificationData(struct ged* UNUSED(g), Options* opt, std::map<std::string, s
     struct bu_process* p;
     int read_cnt = 0;
     char buffer[1024] = {0};
-    std::string result;
+    std::string result = "";
     std::string ncpu(std::to_string(opt->getNCPU()));
     std::string dFile(opt->getDensityFile());
 
@@ -243,23 +233,21 @@ getVerificationData(struct ged* UNUSED(g), Options* opt, std::map<std::string, s
     bu_process_create(&p, const_cast<const char**>(gqa_av_vec.data()), BU_PROCESS_HIDE_WINDOW | BU_PROCESS_OUT_EQ_ERR);
 
     if (bu_process_pid(p) <= 0) {
-	reportAndExit("Problem in getVerificationData gqa process creation, aborting\n");
+        bu_exit(BRLCAD_ERROR, "Problem in getVerificationData gqa process creation, aborting\n");
     }
 
-    while ((read_cnt = bu_process_read_n(p, BU_PROCESS_STDOUT, sizeof(buffer) - 1, buffer)) > 0) {
+    result = "";
+    while ((read_cnt = bu_process_read_n(p, BU_PROCESS_STDOUT, 1024-1, buffer)) > 0) {
         buffer[read_cnt] = '\0';
         result += buffer;
     }
 
-    const int gqa_status = bu_process_wait_n(&p, 0);
-
-    // sanity - we should have clean status code and not the "no density found" message
-    const bool no_density_data = result.find(no_density_msg) != std::string::npos;
-    if (gqa_status != 0 && (!no_density_data || !dFile.empty())) {
-	reportAndExit("Problem collecting gqa volume and mass, aborting\n", result);
+    if (bu_process_wait_n(&p, 0) != 0) {
+        bu_exit(BRLCAD_ERROR, "Problem collecting gqa volume and mass, aborting\n");
     }
 
-    if (!no_density_data) {
+    // make sure we did not get 'no density file' error message
+    if (result.find(no_density_msg) == std::string::npos) {
         // NOTE: tmp_val is zeroed at each parseLabeledNumber call
         if (parseLabeledNumber(tmp_val, result, volume_label))
             volume += tmp_val;
@@ -275,17 +263,17 @@ getVerificationData(struct ged* UNUSED(g), Options* opt, std::map<std::string, s
     bu_process_create(&p, const_cast<const char**>(gqa_av_vec.data()), BU_PROCESS_HIDE_WINDOW | BU_PROCESS_OUT_EQ_ERR);
 
     if (bu_process_pid(p) <= 0) {
-	reportAndExit("Problem in getVerificationData gqa process creation, aborting\n");
+        bu_exit(BRLCAD_ERROR, "Problem in getVerificationData gqa process creation, aborting\n");
     }
 
-    result.clear();
-    while ((read_cnt = bu_process_read_n(p, BU_PROCESS_STDOUT, sizeof(buffer) - 1, buffer)) > 0) {
+    result = "";
+    while ((read_cnt = bu_process_read_n(p, BU_PROCESS_STDOUT, 1024-1, buffer)) > 0) {
         buffer[read_cnt] = '\0';
         result += buffer;
     }
 
     if (bu_process_wait_n(&p, 0) != 0) {
-	reportAndExit("Problem collecting gqa volume, aborting\n", result);
+        bu_exit(BRLCAD_ERROR, "Problem collecting gqa volume, aborting\n");
     }
 
     if (parseLabeledNumber(tmp_val, result, volume_label))
@@ -393,7 +381,7 @@ InformationGatherer::getMainComp()
         ged_exec_exists(g, 2, cmd);
         std::string res = bu_vls_addr(g->ged_result_str);
         if (res != "1") {
-	    reportAndExit("Could not find component (" + topcomp + "), aborting.\n");
+            bu_exit(BRLCAD_ERROR, "Could not find component (%s), aborting.\n", topname);
         }
 
         int entities = getNumEntities(opt->getTopComp());

src/gtools/gist/pch.h

@@ -34,9 +34,6 @@
 
 // OpenCV header files
 #include <opencv2/opencv.hpp>
-#if CV_VERSION_MAJOR >= 5
-#include <opencv2/geometry.hpp>
-#endif
 #include <opencv2/core/core.hpp>
 #include "opencv2/highgui/highgui.hpp"
 

src/gtools/gqa.c

@@ -43,7 +43,6 @@ main(int argc, char *argv[])
     int i, j;
     int db_index;
     int c;
-    int ret;
     const char **av;
     struct ged *gedp;
 
@@ -104,14 +103,14 @@ main(int argc, char *argv[])
 	bu_exit(1, usage, argv[0]);
     }
 
-    ret = ged_exec_gqa(gedp, j, av);
+    (void)ged_exec_gqa(gedp, j, av);
     if (bu_vls_strlen(gedp->ged_result_str) > 0)
 	bu_log("%s", bu_vls_addr(gedp->ged_result_str));
     ged_close(gedp);
 
     bu_free((void *)av, "av");
 
-    return ret;
+    return 0;
 }
 
 

src/isst/CMakeLists.txt

@@ -42,7 +42,6 @@ if(BRLCAD_ENABLE_QT AND OPENGL_TARGETS)
       Qt6::Widgets
       Qt6::OpenGL
       Qt6::OpenGLWidgets
-      ${OPENGL_TARGETS}
     )
   else(Qt6Widgets_FOUND)
     qt5_wrap_cpp(moc_srcs ${isst_headers})

src/isst/gfile.cpp

@@ -26,6 +26,7 @@
 #include <QFile>
 #include <QPlainTextEdit>
 #include <QTextStream>
+#include <vector>
 #include "gfile.h"
 
 #include "bu/malloc.h"
@@ -33,6 +34,7 @@
 #include "gcv.h"
 #include "nmg.h"
 #include "rt/geom.h"
+#include "load_g.h"
 
 /* Replace the load_g globals from adrt with a struct that is
  * passed through the callers. */
@@ -102,7 +104,8 @@ nmg_to_adrt_internal(TIE_3 **tribuf, struct tie_s *cur_tie, struct adrt_mesh_s *
 		    NMG_CK_VERTEX(v);
 
 		    /* convert mm to m */
-		    VSCALE((*tribuf[vert_count]).v, v->vg_p->coord, 1.0/1000.0);
+		    VSCALE((*tribuf[vert_count]).v, v->vg_p->coord,
+			1.0 / ADRT_MODEL_UNITS_PER_TIE_UNIT);
 		    vert_count++;
 		}
 		if (vert_count > 3)
@@ -185,9 +188,12 @@ nmg_to_adrt_regstart(struct db_tree_state *ts, const struct db_full_path *path,
 
 	for (i=0;i<bot->num_faces;i++)
 	{
-	    VSCALE((*d->tribuf[0]).v, (bot->vertices+3*bot->faces[3*i+0]), 1.0/1000.0);
-	    VSCALE((*d->tribuf[1]).v, (bot->vertices+3*bot->faces[3*i+1]), 1.0/1000.0);
-	    VSCALE((*d->tribuf[2]).v, (bot->vertices+3*bot->faces[3*i+2]), 1.0/1000.0);
+	    VSCALE((*d->tribuf[0]).v, (bot->vertices+3*bot->faces[3*i+0]),
+		1.0 / ADRT_MODEL_UNITS_PER_TIE_UNIT);
+	    VSCALE((*d->tribuf[1]).v, (bot->vertices+3*bot->faces[3*i+1]),
+		1.0 / ADRT_MODEL_UNITS_PER_TIE_UNIT);
+	    VSCALE((*d->tribuf[2]).v, (bot->vertices+3*bot->faces[3*i+2]),
+		1.0 / ADRT_MODEL_UNITS_PER_TIE_UNIT);
 
 	    TIE_VAL(tie_push)(d->cur_tie, d->tribuf, 1, mesh, 0);
 	}
@@ -233,6 +239,13 @@ nmg_to_adrt_gcvwrite(struct nmgregion *r, const struct db_full_path *pathp, stru
     nmg_to_adrt_internal(d->tribuf, d->cur_tie, mesh, r);
 }
 
+
+GFile::~GFile()
+{
+    closedb();
+}
+
+
 int
 GFile::load_g(const char *filename, int argc, const char *argv[])
 {
@@ -241,6 +254,9 @@ GFile::load_g(const char *filename, int argc, const char *argv[])
     struct bg_tess_tol ttol;		/* tessellation tolerance in mm */
     struct db_tree_state tree_state;	/* includes tol & model */
     struct isst_nmg_data d;
+    std::vector<const char *> geometry_paths;
+
+    closedb();
 
     RT_DBTS_INIT(&tree_state);
     tree_state.ts_tol = &tol;
@@ -278,6 +294,8 @@ GFile::load_g(const char *filename, int argc, const char *argv[])
     }
     if (db_dirbuild(dbip)) {
 	bu_log("ERROR: db_dirbuild failed\n");
+	db_close(dbip);
+	dbip = NULL;
 	return -1;
     }
     d.dbip = dbip;
@@ -304,20 +322,30 @@ GFile::load_g(const char *filename, int argc, const char *argv[])
     tribuf[2] = (TIE_3 *)bu_malloc(sizeof(TIE_3) * 3, "triangle tribuffer");
     d.tribuf = this->tribuf;
 
-    (void) db_walk_tree(dbip,
-			argc,			/* number of toplevel regions */
-			argv,			/* region names */
-			1,			/* ncpu */
-			&tree_state,		/* initial tree state */
-			nmg_to_adrt_regstart,	/* region start function */
-			gcv_region_end,		/* region end function */
-			rt_booltree_leaf_tess,	/* leaf func */
-			(void *)&gcvwriter);	/* client data */
+    annotations = rt_annot_scene_create(dbip, argc, argv, &ttol, &tol);
+
+    for (int i = 0; i < argc; ++i) {
+	if (adrt_path_has_geometry(dbip, argv[i]))
+	    geometry_paths.push_back(argv[i]);
+    }
+
+    if (!geometry_paths.empty()) {
+	(void) db_walk_tree(dbip,
+		(int)geometry_paths.size(), /* number of toplevel regions */
+		geometry_paths.data(),	/* region names */
+		1,			/* ncpu */
+		&tree_state,		/* initial tree state */
+		nmg_to_adrt_regstart,	/* region start function */
+		gcv_region_end,		/* region end function */
+		adrt_leaf_tess,		/* leaf func */
+		(void *)&gcvwriter);	/* client data */
+    }
 
     /* Release dynamic storage */
     nmg_km(the_model);
     rt_vlist_cleanup();
     db_close(dbip);
+    dbip = NULL;
     bu_free(tribuf[0], "vert");
     bu_free(tribuf[1], "vert");
     bu_free(tribuf[2], "vert");
@@ -331,7 +359,17 @@ GFile::load_g(const char *filename, int argc, const char *argv[])
 void
 GFile::closedb()
 {
-    bu_free(tie, "free tie");
+    rt_annot_scene_destroy(annotations);
+    annotations = NULL;
+    if (tie) {
+	TIE_FREE(tie);
+	bu_free(tie, "free tie");
+	tie = NULL;
+    }
+    if (dbip) {
+	db_close(dbip);
+	dbip = NULL;
+    }
     current_file.clear();
 }
 
@@ -345,4 +383,3 @@ GFile::closedb()
  * End:
  * ex: shiftwidth=4 tabstop=8
  */
-

src/isst/gfile.h

@@ -32,6 +32,7 @@
 extern "C" {
 #include "bn/tol.h"
 #include "raytrace.h"
+#include "rt/primitives/annot.h"
 
 #include "rt/tie.h"
 #include "adrt.h"
@@ -43,13 +44,16 @@ class GFile : public QObject
     Q_OBJECT
 
     public:
+	~GFile() override;
+
 	int load_g(const char *filename, int argc, const char **argv);
 	void closedb();
 
-	struct tie_s *tie;
-	struct adrt_mesh_s *meshes;
-	struct db_i *dbip;
-	TIE_3 **tribuf;
+	struct tie_s *tie = NULL;
+	struct adrt_mesh_s *meshes = NULL;
+	struct db_i *dbip = NULL;
+	struct rt_annot_scene *annotations = NULL;
+	TIE_3 **tribuf = NULL;
 
 	QString current_file;
 };
@@ -66,4 +70,3 @@ class GFile : public QObject
  * End:
  * ex: shiftwidth=4 tabstop=8
  */
-

src/isst/isstapp.cpp

@@ -28,6 +28,9 @@
 int
 ISSTApp::load_g(const char *filename, int argc, const char *argv[])
 {
+    /* Synchronize with the render thread before GFile releases the old
+     * annotation scene and TIE data. */
+    w.canvas->clear_tie();
     if (g.load_g(filename, argc, argv)) {
 	w.statusBar()->showMessage("open failed");
 	return -1;
@@ -35,7 +38,7 @@ ISSTApp::load_g(const char *filename, int argc, const char *argv[])
 
     // The OpenGL widget manages the rendering, so let it know about the
     // TIE data structure associated with the current model
-    w.canvas->set_tie(g.tie);
+    w.canvas->set_tie(g.tie, g.annotations);
 
     return 0;
 }
@@ -51,4 +54,3 @@ ISSTApp::load_g(const char *filename, int argc, const char *argv[])
  * End:
  * ex: shiftwidth=4 tabstop=8
  */
-

src/isst/isstapp.h

@@ -36,16 +36,14 @@
 
 class ISSTApp : public QApplication
 {
-    Q_OBJECT
-
     public:
 	ISSTApp(int &argc, char *argv[]) :QApplication(argc, argv) {};
 	~ISSTApp() {};
 
 	int load_g(const char *filename, int argc, const char **argv);
 
-	ISST_MainWindow w;
 	GFile g;
+	ISST_MainWindow w;
 };
 
 #endif // ISSTAPP_H
@@ -60,4 +58,3 @@ class ISSTApp : public QApplication
  * End:
  * ex: shiftwidth=4 tabstop=8
  */
-

src/isst/isstgl.cpp

@@ -159,8 +159,13 @@ void TIERenderer::render()
     m_grabCond.wait(&m_grabMutex);
     QMutexLocker lock(&m_renderMutex);
     m_grabMutex.unlock();
-    if (m_exiting)
+
+    /* The GUI may clear a scene while a queued render is waiting for the
+     * context.  Return the context before leaving that stale request. */
+    if (m_exiting || !tie) {
+	ctx->moveToThread(qGuiApp->thread());
 	return;
+    }
     Q_ASSERT(ctx->thread() == QThread::currentThread());
 
 
@@ -183,7 +188,8 @@ void TIERenderer::render()
 
     // Core TIE render
     render_camera_prep(&camera);
-    render_camera_render(&camera, tie, &tile, &buffer_image);
+    render_camera_render_annotations(&camera, tie, &tile, &buffer_image,
+	annotations, ADRT_MODEL_UNITS_PER_TIE_UNIT);
 
     glDisable(GL_LIGHTING);
 
@@ -297,17 +303,31 @@ void isstGL::grabContext()
 }
 
 void
-isstGL::set_tie(struct tie_s *in_tie)
+isstGL::clear_tie()
 {
-    m_renderer->tie = in_tie;
+    m_renderer->lockRenderer();
+    m_renderer->tie = NULL;
+    m_renderer->annotations = NULL;
+    m_renderer->changed = false;
+    m_renderer->unlockRenderer();
+}
 
-    // Initialize the camera position
-    VSETALL(m_renderer->camera.pos, m_renderer->tie->radius);
-    VMOVE(m_renderer->camera.focus, m_renderer->tie->mid);
+
+void
+isstGL::set_tie(struct tie_s *in_tie, struct rt_annot_scene *annotations)
+{
+    m_renderer->lockRenderer();
+    m_renderer->tie = in_tie;
+    m_renderer->annotations = annotations;
+
+    render_camera_fit_scene(&m_renderer->camera, m_renderer->tie, annotations,
+	ADRT_MODEL_UNITS_PER_TIE_UNIT);
 
     // Record the initial settings for use in subsequent calculations
-    VSETALL(m_renderer->camera_pos_init, m_renderer->tie->radius);
-    VMOVE(m_renderer->camera_focus_init, m_renderer->tie->mid);
+    VMOVE(m_renderer->camera_pos_init, m_renderer->camera.pos);
+    VMOVE(m_renderer->camera_focus_init, m_renderer->camera.focus);
+    m_renderer->changed = true;
+    m_renderer->unlockRenderer();
 
     // Having just loaded a new TIE scene,
     // we need a new image
@@ -369,4 +389,3 @@ void isstGL::save_image() {
 // c-file-style: "stroustrup"
 // End:
 // ex: shiftwidth=4 tabstop=8
-

src/isst/isstgl.h

@@ -77,6 +77,7 @@ class TIERenderer : public QObject, protected QOpenGLFunctions
 
     public:
 	struct tie_s *tie = NULL; // From parent app
+	struct rt_annot_scene *annotations = NULL; // Owned by the parent app
 	struct render_camera_s camera;
 	vect_t camera_pos_init;
 	vect_t camera_focus_init;
@@ -113,7 +114,8 @@ class isstGL : public QOpenGLWidget
 	explicit isstGL(QWidget *parent = nullptr);
 	~isstGL();
 
-	void set_tie(struct tie_s *in_tie);
+	void clear_tie();
+	void set_tie(struct tie_s *in_tie, struct rt_annot_scene *annotations);
 
 	void save_image();
 
@@ -153,4 +155,3 @@ class isstGL : public QOpenGLWidget
 // c-file-style: "stroustrup"
 // End:
 // ex: shiftwidth=4 tabstop=8
-

src/libanalyze/CMakeLists.txt

@@ -29,7 +29,6 @@ set(
   api.c
   centroid.c
   check_options.c
-  contour.cpp
   density.cpp
   globals.c
   inside.c
@@ -69,7 +68,7 @@ set(
 cmakefiles(${libanalyze_ignored_files})
 
 # Note - libanalyze_deps is defined by ${BRLCAD_SOURCE_DIR}/src/source_dirs.cmake
-brlcad_addlib(libanalyze "${LIBANALYZE_SOURCES}" "" "${CMAKE_CURRENT_BINARY_DIR};${EIGEN3_INCLUDE_DIR}"
+brlcad_addlib(libanalyze "${LIBANALYZE_SOURCES}" "" "${CMAKE_CURRENT_BINARY_DIR}"
   PUBLIC_LIBS ${libanalyze_deps}
 )
 set_target_properties(libanalyze PROPERTIES VERSION 20.0.1 SOVERSION 20)

src/libanalyze/nirt/nirt.cpp

@@ -480,6 +480,22 @@ _nirt_get_resource(struct nirt_state *nss)
 }
 
 
+static void
+_nirt_annotation_result_clear(struct nirt_state *nss)
+{
+    delete nss->i->vals->annotation;
+    nss->i->vals->annotation = NULL;
+}
+
+
+static void
+_nirt_annotation_scene_clear(struct nirt_state *nss)
+{
+    rt_annot_scene_destroy(nss->i->annotation_scene);
+    nss->i->annotation_scene = NULL;
+}
+
+
 int
 _nirt_raytrace_prep(struct nirt_state *nss)
 {
@@ -517,11 +533,17 @@ _nirt_raytrace_prep(struct nirt_state *nss)
 	bu_free(attrs, "objs");
 	return -1;
     }
-    bu_free(objs, "objs");
     bu_free(attrs, "objs");
 
     nmsg(nss, "Prepping the geometry...\n");
     rt_prep(nss->i->ap->a_rt_i);
+    _nirt_annotation_scene_clear(nss);
+    if (nss->i->use_annotations) {
+	nss->i->annotation_scene = rt_annot_scene_create(nss->i->dbip,
+	    ocnt, objs, &nss->i->ap->a_rt_i->rti_ttol,
+	    &nss->i->ap->a_rt_i->rti_tol);
+    }
+    bu_free(objs, "objs");
     nmsg(nss, "%s ", (nss->i->active_paths.size() == 1) ? "Object" : "Objects");
     for (size_t i = 0; i < nss->i->active_paths.size(); i++) {
 	if (i == 0) {
@@ -960,6 +982,17 @@ _nirt_print_fmt_substr(struct nirt_state *nss, struct bu_vls *ostr, const char *
     nirt_print_key("a", r->a);
     nirt_print_key("e", r->e);
 
+    if (r->annotation) {
+	nirt_print_key("annot_path", r->annotation->path.c_str());
+	nirt_print_key("annot_segment", r->annotation->segment);
+	nirt_print_key("annot_role", r->annotation->role.c_str());
+	nirt_print_key("annot_x", r->annotation->point[X] * base2local);
+	nirt_print_key("annot_y", r->annotation->point[Y] * base2local);
+	nirt_print_key("annot_z", r->annotation->point[Z] * base2local);
+	nirt_print_key("annot_d", r->annotation->d * base2local);
+	nirt_print_key("annot_visible", r->annotation->visible);
+    }
+
     if (!r->seg || r->seg->type == 0) return;
 
     if (r->seg->type == NIRT_PARTITION_SEG || r->seg->type == NIRT_ALL_SEG) {
@@ -1056,6 +1089,9 @@ _nirt_report(struct nirt_state *nss, char type, struct nirt_output_record *r)
     std::vector<std::pair<std::string,std::string> >::iterator f_it;
 
     switch (type) {
+	case 'a':
+	    fmt_vect = &nss->i->fmt.annotation;
+	    break;
 	case 'r':
 	    fmt_vect = &nss->i->fmt.ray;
 	    break;
@@ -1091,6 +1127,75 @@ _nirt_report(struct nirt_state *nss, char type, struct nirt_output_record *r)
 }
 
 
+static const char *
+_nirt_annotation_role(uint32_t role)
+{
+    switch (role) {
+	case RT_ANNOT_ROLE_UNSPECIFIED:
+	    return "unspecified";
+	case RT_ANNOT_ROLE_GEOMETRY:
+	    return "geometry";
+	case RT_ANNOT_ROLE_TEXT:
+	    return "text";
+	case RT_ANNOT_ROLE_LEADER:
+	    return "leader";
+	case RT_ANNOT_ROLE_EXTENSION:
+	    return "extension";
+	case RT_ANNOT_ROLE_DIMENSION:
+	    return "dimension";
+	case RT_ANNOT_ROLE_ARROWHEAD:
+	    return "arrowhead";
+	case RT_ANNOT_ROLE_SYMBOL:
+	    return "symbol";
+	case RT_ANNOT_ROLE_MASK:
+	    return "mask";
+	case RT_ANNOT_ROLE_CENTERMARK:
+	    return "centermark";
+	case RT_ANNOT_ROLE_PLACEHOLDER:
+	    return "placeholder";
+	case RT_ANNOT_ROLE_TEXT_DECORATION:
+	    return "text-decoration";
+	default:
+	    return "unknown";
+    }
+}
+
+
+static void
+_nirt_report_annotations(struct nirt_state *nss)
+{
+    const size_t count = rt_annot_scene_query_model(nss->i->annotation_scene,
+	&nss->i->ap->a_ray, nss->i->first_solid_distance, NULL, 0);
+    if (!count)
+	return;
+
+    std::vector<struct rt_annot_hit> hits(count);
+    const size_t found = rt_annot_scene_query_model(nss->i->annotation_scene,
+	&nss->i->ap->a_ray, nss->i->first_solid_distance, hits.data(),
+	hits.size());
+    nirt_seg *saved_segment = nss->i->vals->seg;
+    nss->i->vals->seg = NULL;
+    for (size_t i = 0; i < std::min(found, hits.size()); ++i) {
+	const struct rt_annot_hit &hit = hits[i];
+	if (hit.segment > (size_t)std::numeric_limits<int>::max()) {
+	    nerr(nss, "Annotation segment index is too large for NIRT reporting\n");
+	    continue;
+	}
+	if (!nss->i->vals->annotation)
+	    nss->i->vals->annotation = new nirt_annotation_record;
+	nirt_annotation_record *record = nss->i->vals->annotation;
+	record->path = hit.path ? hit.path : "";
+	record->segment = (int)hit.segment;
+	record->role = _nirt_annotation_role(hit.role);
+	VMOVE(record->point, hit.point);
+	record->d = d_calc(nss, record->point);
+	record->visible = hit.visible;
+	_nirt_report(nss, 'a', nss->i->vals);
+    }
+    nss->i->vals->seg = saved_segment;
+}
+
+
 /************************
  * Raytracing Callbacks *
  ************************/
@@ -1213,6 +1318,21 @@ _nirt_if_hit(struct application *ap, struct partition *part_head, struct seg *UN
     int ev_odd = 1; /* first partition is colored as "odd" */
     point_t out_old = VINIT_ZERO;
     double d_out_old = 0.0;
+
+    nss->i->first_solid_distance = INFINITY;
+    for (part = part_head->pt_forw; part != part_head; part = part->pt_forw) {
+	if (part->pt_regionp->reg_aircode)
+	    continue;
+	if (part->pt_inhit->hit_dist >= 0.0) {
+	    nss->i->first_solid_distance = part->pt_inhit->hit_dist;
+	    break;
+	}
+	if (part->pt_outhit->hit_dist >= 0.0) {
+	    nss->i->first_solid_distance = part->pt_outhit->hit_dist;
+	    break;
+	}
+    }
+
     nirt_seg *s = new nirt_seg;
     if (vals->seg) {
 	delete vals->seg;
@@ -1549,9 +1669,10 @@ static const struct nirt_cmd_desc nirt_descs[] = {
     { "s",              "shoot a ray at the target",                     NULL },
     { "backout",        "back out of model",                             NULL },
     { "useair",         "set/query use of air",                          "<0|1|2|...>" },
+    { "useannot",       "set/query annotation reporting",                "<0|1>" },
     { "units",          "set/query local units",                         "<mm|cm|m|in|ft>" },
     { "overlap_claims", "set/query overlap rebuilding/retention",        "<0|1|2|3>" },
-    { "fmt",            "set/query output formats",                      "{rhpfmog} format item item ..." },
+    { "fmt",            "set/query output formats",                      "{arhpfmog} format item item ..." },
     { "plotfile",       "designate an output file to hold plot data",    "file" },
     { "print",          "query an output item",                          "item" },
     { "debug",          "set/query nirt debug flags",                    "[-h] [-l [lib]] [-C [lib]] [-V [lib] [val]] [lib [flag]]" },
@@ -1828,6 +1949,9 @@ _nirt_cmd_shoot(void *ns, int argc, const char **UNUSED(argv))
 	return -1;
     }
 
+    _nirt_annotation_result_clear(nss);
+    nss->i->first_solid_distance = INFINITY;
+
     /* If we have no active rtip, we don't need to prep or shoot */
     if (!_nirt_get_rtip(nss)) {
 	return 0;
@@ -1860,9 +1984,10 @@ _nirt_cmd_shoot(void *ns, int argc, const char **UNUSED(argv))
 	    nss->i->ap->a_ray.r_pt[2] * nss->i->base2local,
 	    V3ARGS(nss->i->ap->a_ray.r_dir));
 
-    // TODO - any necessary initialization for data collection by callbacks
     _nirt_init_ovlp(nss);
     (void)rt_shootray(nss->i->ap);
+    if (nss->i->use_annotations)
+	_nirt_report_annotations(nss);
 
     // Undo backout
     for (i = 0; i < 3; ++i) {
@@ -1965,6 +2090,45 @@ _nirt_cmd_use_air(void *ns, int argc, const char *argv[])
 }
 
 
+extern "C" int
+_nirt_cmd_use_annotations(void *ns, int argc, const char *argv[])
+{
+    struct nirt_state *nss = (struct nirt_state *)ns;
+    struct bu_vls optparse_msg = BU_VLS_INIT_ZERO;
+    int enabled = 0;
+
+    if (!ns) return -1;
+    if (argc == 1) {
+	nout(nss, "use_annotations = %d\n", nss->i->use_annotations);
+	return 0;
+    }
+    if (argc > 2) {
+	nerr(nss, "Usage:  useannot %s\n", _nirt_get_desc_args("useannot"));
+	return -1;
+    }
+    if (bu_opt_int(&optparse_msg, 1, (const char **)&argv[1],
+	    (void *)&enabled) == -1) {
+	nerr(nss, "Error: bu_opt value read failure: %s\n\nUsage:  useannot %s\n",
+	    bu_vls_cstr(&optparse_msg), _nirt_get_desc_args("useannot"));
+	bu_vls_free(&optparse_msg);
+	return -1;
+    }
+    bu_vls_free(&optparse_msg);
+    if (enabled != 0 && enabled != 1) {
+	nerr(nss, "useannot must be set to 0 (off) or 1 (on)\n");
+	return -1;
+    }
+    if (enabled != nss->i->use_annotations) {
+	nss->i->use_annotations = enabled;
+	_nirt_annotation_scene_clear(nss);
+	_nirt_annotation_result_clear(nss);
+	if (enabled)
+	    nss->i->need_reprep = 1;
+    }
+    return 0;
+}
+
+
 extern "C" int
 _nirt_cmd_units(void *ns, int argc, const char *argv[])
 {
@@ -2085,6 +2249,9 @@ _nirt_cmd_format_output(void *ns, int argc, const char **argv)
 	// print out current fmt str for the specified type
 	struct bu_vls ostr = BU_VLS_INIT_ZERO;
 	switch (type) {
+	    case 'a':
+		_nirt_print_fmt_str(&ostr, nss->i->fmt.annotation);
+		break;
 	    case 'r':
 		_nirt_print_fmt_str(&ostr, nss->i->fmt.ray);
 		break;
@@ -2191,6 +2358,10 @@ _nirt_cmd_format_output(void *ns, int argc, const char **argv)
 	 * formatting instructions with the new. */
 set_fmt:
 	switch (type) {
+	    case 'a':
+		nss->i->fmt.annotation.clear();
+		nss->i->fmt.annotation = fmt_tmp;
+		break;
 	    case 'r':
 		nss->i->fmt.ray.clear();
 		nss->i->fmt.ray = fmt_tmp;
@@ -2543,6 +2714,7 @@ _nirt_cmd_state(void *ns, int argc, const char *argv[])
 	bu_vls_printf(&dumpstr, "xyz %g %g %g\n", V3ARGS(nss->i->vals->orig));
 	bu_vls_printf(&dumpstr, "dir %g %g %g\n", V3ARGS(nss->i->vals->dir));
 	bu_vls_printf(&dumpstr, "useair %d\n", nss->i->use_air);
+	bu_vls_printf(&dumpstr, "useannot %d\n", nss->i->use_annotations);
 	bu_vls_printf(&dumpstr, "units %s\n", bu_units_string(nss->i->local2base));
 	bu_vls_printf(&dumpstr, "overlap_claims ");
 	switch (nss->i->overlap_claims) {
@@ -2563,6 +2735,8 @@ _nirt_cmd_state(void *ns, int argc, const char *argv[])
 		break;
 	}
 	struct bu_vls ostr = BU_VLS_INIT_ZERO;
+	_nirt_print_fmt_cmd(&ostr, 'a', nss->i->fmt.annotation);
+	bu_vls_printf(&dumpstr, "%s", bu_vls_cstr(&ostr));
 	_nirt_print_fmt_cmd(&ostr, 'r', nss->i->fmt.ray);
 	bu_vls_printf(&dumpstr, "%s", bu_vls_cstr(&ostr));
 	_nirt_print_fmt_cmd(&ostr, 'h', nss->i->fmt.head);
@@ -2617,6 +2791,7 @@ const struct bu_cmdtab _libanalyze_nirt_cmds[] = {
     { "state",          _nirt_cmd_state},
     { "units",          _nirt_cmd_units},
     { "useair",         _nirt_cmd_use_air},
+    { "useannot",       _nirt_cmd_use_annotations},
     { "xyz",            _nirt_cmd_target_coor},
     { (char *)NULL,     NULL}
 };
@@ -2783,6 +2958,14 @@ static const char *nirt_cmd_tbl_defs =
 "z_dir,          FNOUNIT,       Ray direction unit vector z component,"
 "a,              FNOUNIT,       Azimuth,"
 "e,              FNOUNIT,       Elevation,"
+"annot_path,     STRING,        Annotation database path,"
+"annot_segment,  INT,           Annotation segment index,"
+"annot_role,     STRING,        Annotation segment role,"
+"annot_x,        FLOAT,         Annotation hit X coordinate,"
+"annot_y,        FLOAT,         Annotation hit Y coordinate,"
+"annot_z,        FLOAT,         Annotation hit Z coordinate,"
+"annot_d,        FLOAT,         Annotation hit gridplane distance,"
+"annot_visible,  INT,           Annotation hit visibility,"
 "x_in,           FLOAT,         ,"
 "y_in,           FLOAT,         ,"
 "z_in,           FLOAT,         ,"
@@ -2895,6 +3078,7 @@ nirt_init(struct nirt_state *ns)
     n->base2local = 0.0;
     n->local2base = 0.0;
     n->use_air = 0;
+    n->use_annotations = 0;
     n->backout = 1;
 
     n->b_state = false;
@@ -2909,6 +3093,8 @@ nirt_init(struct nirt_state *ns)
     BU_GET(n->res_air, struct resource);
     n->rtip = RTI_NULL;
     n->rtip_air = RTI_NULL;
+    n->annotation_scene = NULL;
+    n->first_solid_distance = INFINITY;
     n->need_reprep = 1;
 
     BU_GET(n->val_types, struct bu_attribute_value_set);
@@ -2922,6 +3108,7 @@ nirt_init(struct nirt_state *ns)
     n->vals->h = 0.0;
     n->vals->v = 0.0;
     n->vals->d_orig = 0.0;
+    n->vals->annotation = NULL;
     n->vals->seg = NULL;
 
     _nirt_diff_create(ns);
@@ -3025,6 +3212,8 @@ nirt_clear_dbip(struct nirt_state *ns)
 {
     if (!ns) return -1;
 
+    _nirt_annotation_scene_clear(ns);
+    _nirt_annotation_result_clear(ns);
     db_close(ns->i->dbip);
     ns->i->dbip = NULL;
 
@@ -3045,6 +3234,9 @@ void
 nirt_destroy(struct nirt_state *ns)
 {
     if (!ns) return;
+
+    _nirt_annotation_scene_clear(ns);
+    _nirt_annotation_result_clear(ns);
     bu_vls_free(&ns->nirt_cmd);
     bu_vls_free(ns->i->err);
     bu_vls_free(ns->i->msg);
@@ -3057,6 +3249,7 @@ nirt_destroy(struct nirt_state *ns)
 
     db_close(ns->i->dbip);
 
+    delete ns->i->vals->seg;
     BU_PUT(ns->i->vals, struct nirt_output_record);
 
     BU_PUT(ns->i->res, struct resource);

src/libanalyze/nirt/nirt.h

@@ -83,7 +83,7 @@ extern "C" b_off_t ftello(FILE *);
 #define NIRT_OVLP_RETAIN             3
 
 #define NIRT_PRINTF_SPECIFIERS "difeEgGs"
-#define NIRT_OUTPUT_TYPE_SPECIFIERS "rhpfmog"
+#define NIRT_OUTPUT_TYPE_SPECIFIERS "arhpfmog"
 
 #define HELPFLAG "--print-help"
 #define PURPOSEFLAG "--print-purpose"
@@ -274,6 +274,16 @@ class nirt_seg {
 	}
 };
 
+class nirt_annotation_record {
+    public:
+	std::string path;
+	int segment;
+	std::string role;
+	point_t point;
+	fastf_t d;
+	int visible;
+};
+
 /* This record structure doesn't have to correspond exactly to the above list
  * of available values, but it needs to retain sufficient information to
  * support the ability to generate all of them upon request. */
@@ -286,6 +296,7 @@ struct nirt_output_record {
     fastf_t a;
     fastf_t e;
     struct nirt_overlap ovlp_list;
+    nirt_annotation_record *annotation;
     nirt_seg *seg;
 };
 
@@ -294,6 +305,7 @@ struct nirt_diff_state;
 
 class nirt_fmt_state {
     public:
+	std::vector<std::pair<std::string,std::string> > annotation;
 	std::vector<std::pair<std::string,std::string> > ray;
 	std::vector<std::pair<std::string,std::string> > head;
 	std::vector<std::pair<std::string,std::string> > part;
@@ -303,6 +315,7 @@ class nirt_fmt_state {
 	std::vector<std::pair<std::string,std::string> > gap;
 
 	void clear() {
+	    annotation.clear();
 	    ray.clear();
 	    head.clear();
 	    part.clear();
@@ -347,6 +360,7 @@ struct nirt_state_impl {
     int backout;
     int overlap_claims;
     int use_air;
+    int use_annotations;
     std::set<std::string> attrs;        // active attributes
     std::vector<std::string> active_paths; // active paths for raytracer
     struct nirt_output_record *vals;
@@ -368,6 +382,8 @@ struct nirt_state_impl {
     struct resource *res;
     struct rt_i *rtip_air;
     struct resource *res_air;
+    struct rt_annot_scene *annotation_scene;
+    fastf_t first_solid_distance;
     int need_reprep;
 
     /* internal format specifier arrays */