Diff of the two buildlogs:

--
--- b1/build.log	2024-04-16 21:11:17.974478981 +0000
+++ b2/build.log	2024-04-16 21:13:53.367845489 +0000
@@ -1,6 +1,6 @@
 I: pbuilder: network access will be disabled during build
-I: Current time: Mon May 19 15:32:44 -12 2025
-I: pbuilder-time-stamp: 1747711964
+I: Current time: Wed Apr 17 11:11:20 +14 2024
+I: pbuilder-time-stamp: 1713301880
 I: Building the build Environment
 I: extracting base tarball [/var/cache/pbuilder/trixie-reproducible-base.tgz]
 I: copying local configuration
@@ -27,54 +27,86 @@
 dpkg-source: info: unpacking octave-nurbs_1.4.3-3.debian.tar.xz
 I: Not using root during the build.
 I: Installing the build-deps
-I: user script /srv/workspace/pbuilder/116699/tmp/hooks/D02_print_environment starting
+I: user script /srv/workspace/pbuilder/20177/tmp/hooks/D01_modify_environment starting
+debug: Running on ionos2-i386.
+I: Changing host+domainname to test build reproducibility
+I: Adding a custom variable just for the fun of it...
+I: Changing /bin/sh to bash
+'/bin/sh' -> '/bin/bash'
+lrwxrwxrwx 1 root root 9 Apr 16 21:11 /bin/sh -> /bin/bash
+I: Setting pbuilder2's login shell to /bin/bash
+I: Setting pbuilder2's GECOS to second user,second room,second work-phone,second home-phone,second other
+I: user script /srv/workspace/pbuilder/20177/tmp/hooks/D01_modify_environment finished
+I: user script /srv/workspace/pbuilder/20177/tmp/hooks/D02_print_environment starting
 I: set
-  BUILDDIR='/build/reproducible-path'
-  BUILDUSERGECOS='first user,first room,first work-phone,first home-phone,first other'
-  BUILDUSERNAME='pbuilder1'
-  BUILD_ARCH='i386'
-  DEBIAN_FRONTEND='noninteractive'
-  DEB_BUILD_OPTIONS='buildinfo=+all reproducible=+all parallel=22 '
-  DISTRIBUTION='trixie'
-  HOME='/root'
-  HOST_ARCH='i386'
+  BASH=/bin/sh
+  BASHOPTS=checkwinsize:cmdhist:complete_fullquote:extquote:force_fignore:globasciiranges:globskipdots:hostcomplete:interactive_comments:patsub_replacement:progcomp:promptvars:sourcepath
+  BASH_ALIASES=()
+  BASH_ARGC=()
+  BASH_ARGV=()
+  BASH_CMDS=()
+  BASH_LINENO=([0]="12" [1]="0")
+  BASH_LOADABLES_PATH=/usr/local/lib/bash:/usr/lib/bash:/opt/local/lib/bash:/usr/pkg/lib/bash:/opt/pkg/lib/bash:.
+  BASH_SOURCE=([0]="/tmp/hooks/D02_print_environment" [1]="/tmp/hooks/D02_print_environment")
+  BASH_VERSINFO=([0]="5" [1]="2" [2]="21" [3]="1" [4]="release" [5]="i686-pc-linux-gnu")
+  BASH_VERSION='5.2.21(1)-release'
+  BUILDDIR=/build/reproducible-path
+  BUILDUSERGECOS='second user,second room,second work-phone,second home-phone,second other'
+  BUILDUSERNAME=pbuilder2
+  BUILD_ARCH=i386
+  DEBIAN_FRONTEND=noninteractive
+  DEB_BUILD_OPTIONS='buildinfo=+all reproducible=+all parallel=10 '
+  DIRSTACK=()
+  DISTRIBUTION=trixie
+  EUID=0
+  FUNCNAME=([0]="Echo" [1]="main")
+  GROUPS=()
+  HOME=/root
+  HOSTNAME=i-capture-the-hostname
+  HOSTTYPE=i686
+  HOST_ARCH=i386
   IFS=' 	
   '
-  INVOCATION_ID='7643a6c130554ed5a50f8a034c3546f1'
-  LANG='C'
-  LANGUAGE='en_US:en'
-  LC_ALL='C'
-  LD_LIBRARY_PATH='/usr/lib/libeatmydata'
-  LD_PRELOAD='libeatmydata.so'
-  MAIL='/var/mail/root'
-  OPTIND='1'
-  PATH='/usr/sbin:/usr/bin:/sbin:/bin:/usr/games'
-  PBCURRENTCOMMANDLINEOPERATION='build'
-  PBUILDER_OPERATION='build'
-  PBUILDER_PKGDATADIR='/usr/share/pbuilder'
-  PBUILDER_PKGLIBDIR='/usr/lib/pbuilder'
-  PBUILDER_SYSCONFDIR='/etc'
-  PPID='116699'
-  PS1='# '
-  PS2='> '
+  INVOCATION_ID=d45f655a224143fa8fa7d25f1931eb44
+  LANG=C
+  LANGUAGE=de_CH:de
+  LC_ALL=C
+  LD_LIBRARY_PATH=/usr/lib/libeatmydata
+  LD_PRELOAD=libeatmydata.so
+  MACHTYPE=i686-pc-linux-gnu
+  MAIL=/var/mail/root
+  OPTERR=1
+  OPTIND=1
+  OSTYPE=linux-gnu
+  PATH=/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path
+  PBCURRENTCOMMANDLINEOPERATION=build
+  PBUILDER_OPERATION=build
+  PBUILDER_PKGDATADIR=/usr/share/pbuilder
+  PBUILDER_PKGLIBDIR=/usr/lib/pbuilder
+  PBUILDER_SYSCONFDIR=/etc
+  PIPESTATUS=([0]="0")
+  POSIXLY_CORRECT=y
+  PPID=20177
   PS4='+ '
-  PWD='/'
-  SHELL='/bin/bash'
-  SHLVL='2'
-  SUDO_COMMAND='/usr/bin/timeout -k 18.1h 18h /usr/bin/ionice -c 3 /usr/bin/nice /usr/sbin/pbuilder --build --configfile /srv/reproducible-results/rbuild-debian/r-b-build.73SuaaUa/pbuilderrc_VhZb --distribution trixie --hookdir /etc/pbuilder/first-build-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/trixie-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.73SuaaUa/b1 --logfile b1/build.log octave-nurbs_1.4.3-3.dsc'
-  SUDO_GID='112'
-  SUDO_UID='107'
-  SUDO_USER='jenkins'
-  TERM='unknown'
-  TZ='/usr/share/zoneinfo/Etc/GMT+12'
-  USER='root'
-  _='/usr/bin/systemd-run'
-  http_proxy='http://85.184.249.68:3128'
+  PWD=/
+  SHELL=/bin/bash
+  SHELLOPTS=braceexpand:errexit:hashall:interactive-comments:posix
+  SHLVL=3
+  SUDO_COMMAND='/usr/bin/timeout -k 24.1h 24h /usr/bin/ionice -c 3 /usr/bin/nice -n 11 /usr/bin/unshare --uts -- /usr/sbin/pbuilder --build --configfile /srv/reproducible-results/rbuild-debian/r-b-build.73SuaaUa/pbuilderrc_w7Ho --distribution trixie --hookdir /etc/pbuilder/rebuild-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/trixie-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.73SuaaUa/b2 --logfile b2/build.log octave-nurbs_1.4.3-3.dsc'
+  SUDO_GID=112
+  SUDO_UID=107
+  SUDO_USER=jenkins
+  TERM=unknown
+  TZ=/usr/share/zoneinfo/Etc/GMT-14
+  UID=0
+  USER=root
+  _='I: set'
+  http_proxy=http://78.137.99.97:3128
 I: uname -a
-  Linux ionos16-i386 6.1.0-20-amd64 #1 SMP PREEMPT_DYNAMIC Debian 6.1.85-1 (2024-04-11) x86_64 GNU/Linux
+  Linux i-capture-the-hostname 6.1.0-20-amd64 #1 SMP PREEMPT_DYNAMIC Debian 6.1.85-1 (2024-04-11) x86_64 GNU/Linux
 I: ls -l /bin
-  lrwxrwxrwx 1 root root 7 May 18 17:46 /bin -> usr/bin
-I: user script /srv/workspace/pbuilder/116699/tmp/hooks/D02_print_environment finished
+  lrwxrwxrwx 1 root root 7 Apr 16 11:24 /bin -> usr/bin
+I: user script /srv/workspace/pbuilder/20177/tmp/hooks/D02_print_environment finished
  -> Attempting to satisfy build-dependencies
  -> Creating pbuilder-satisfydepends-dummy package
 Package: pbuilder-satisfydepends-dummy
@@ -620,7 +652,7 @@
 Get: 506 http://deb.debian.org/debian trixie/main i386 gfortran i386 4:13.2.0-7 [1432 B]
 Get: 507 http://deb.debian.org/debian trixie/main i386 octave-dev i386 8.4.0-1+b1 [960 kB]
 Get: 508 http://deb.debian.org/debian trixie/main i386 dh-octave all 1.8.0 [22.7 kB]
-Fetched 212 MB in 6s (32.8 MB/s)
+Fetched 212 MB in 3s (64.0 MB/s)
 debconf: delaying package configuration, since apt-utils is not installed
 Selecting previously unselected package netbase.
 (Reading database ... 
(Reading database ... 5%
(Reading database ... 10%
(Reading database ... 15%
(Reading database ... 20%
(Reading database ... 25%
(Reading database ... 30%
(Reading database ... 35%
(Reading database ... 40%
(Reading database ... 45%
(Reading database ... 50%
(Reading database ... 55%
(Reading database ... 60%
(Reading database ... 65%
(Reading database ... 70%
(Reading database ... 75%
(Reading database ... 80%
(Reading database ... 85%
(Reading database ... 90%
(Reading database ... 95%
(Reading database ... 100%
(Reading database ... 19874 files and directories currently installed.)
@@ -2705,7 +2737,11 @@
 Building tag database...
  -> Finished parsing the build-deps
 I: Building the package
-I: Running cd /build/reproducible-path/octave-nurbs-1.4.3/ && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games" HOME="/nonexistent/first-build" dpkg-buildpackage -us -uc -b && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games" HOME="/nonexistent/first-build" dpkg-genchanges -S  > ../octave-nurbs_1.4.3-3_source.changes
+I: user script /srv/workspace/pbuilder/20177/tmp/hooks/A99_set_merged_usr starting
+Not re-configuring usrmerge for trixie
+I: user script /srv/workspace/pbuilder/20177/tmp/hooks/A99_set_merged_usr finished
+hostname: Name or service not known
+I: Running cd /build/reproducible-path/octave-nurbs-1.4.3/ && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path" HOME="/nonexistent/second-build" dpkg-buildpackage -us -uc -b && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path" HOME="/nonexistent/second-build" dpkg-genchanges -S  > ../octave-nurbs_1.4.3-3_source.changes
 dpkg-buildpackage: info: source package octave-nurbs
 dpkg-buildpackage: info: source version 1.4.3-3
 dpkg-buildpackage: info: source distribution unstable
@@ -2744,38 +2780,38 @@
 /usr/bin/mkoctfile --verbose -c low_level_functions.cc
 g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    low_level_functions.cc -o low_level_functions.o
 /usr/bin/mkoctfile --verbose basisfun.cc low_level_functions.o
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    basisfun.cc -o /tmp/oct-8vkCYU.o
 /usr/bin/mkoctfile --verbose basisfunder.cc low_level_functions.o
 /usr/bin/mkoctfile --verbose bspderiv.cc low_level_functions.o
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    basisfunder.cc -o /tmp/oct-9kukhJ.o
 /usr/bin/mkoctfile --verbose bspeval.cc low_level_functions.o
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    bspderiv.cc -o /tmp/oct-w4VyYs.o
 /usr/bin/mkoctfile --verbose curvederivcpts.cc low_level_functions.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    basisfun.cc -o /tmp/oct-DdldKy.o
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    bspeval.cc -o /tmp/oct-52ISI0.o
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    curvederivcpts.cc -o /tmp/oct-Gkcuyg.o
 /usr/bin/mkoctfile --verbose nrb_srf_basisfun__.cc low_level_functions.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    basisfunder.cc -o /tmp/oct-ZIbIYv.o
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    nrb_srf_basisfun__.cc -o /tmp/oct-0g4eSK.o
 /usr/bin/mkoctfile --verbose nrb_srf_basisfun_der__.cc low_level_functions.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    bspderiv.cc -o /tmp/oct-WZzx9X.o
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    nrb_srf_basisfun_der__.cc -o /tmp/oct-Fv8kMX.o
 /usr/bin/mkoctfile --verbose nrbsurfderiveval.cc low_level_functions.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    bspeval.cc -o /tmp/oct-XyGoqp.o
 /usr/bin/mkoctfile --verbose surfderivcpts.cc low_level_functions.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    curvederivcpts.cc -o /tmp/oct-jXM5Iq.o
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    nrbsurfderiveval.cc -o /tmp/oct-BfsaPW.o
 /usr/bin/mkoctfile --verbose surfderiveval.cc low_level_functions.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    nrb_srf_basisfun__.cc -o /tmp/oct-BWUC5i.o
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    surfderivcpts.cc -o /tmp/oct-CoKrJO.o
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    surfderiveval.cc -o /tmp/oct-nuHEv8.o
 /usr/bin/mkoctfile --verbose tbasisfun.cc low_level_functions.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    nrb_srf_basisfun_der__.cc -o /tmp/oct-mARFAi.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    nrbsurfderiveval.cc -o /tmp/oct-CreADb.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    surfderivcpts.cc -o /tmp/oct-23wW1C.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    surfderiveval.cc -o /tmp/oct-PyDCfo.o
-g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    tbasisfun.cc -o /tmp/oct-gq0prV.o
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o surfderivcpts.oct  low_level_functions.o /tmp/oct-23wW1C.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o curvederivcpts.oct  low_level_functions.o /tmp/oct-jXM5Iq.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o basisfun.oct  low_level_functions.o /tmp/oct-DdldKy.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o basisfunder.oct  low_level_functions.o /tmp/oct-ZIbIYv.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o bspderiv.oct  low_level_functions.o /tmp/oct-WZzx9X.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o surfderiveval.oct  low_level_functions.o /tmp/oct-PyDCfo.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o bspeval.oct  low_level_functions.o /tmp/oct-XyGoqp.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o nrbsurfderiveval.oct  low_level_functions.o /tmp/oct-CreADb.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o tbasisfun.oct  low_level_functions.o /tmp/oct-gq0prV.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o nrb_srf_basisfun__.oct  low_level_functions.o /tmp/oct-BWUC5i.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
-g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o nrb_srf_basisfun_der__.oct  low_level_functions.o /tmp/oct-mARFAi.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security    tbasisfun.cc -o /tmp/oct-NDNR1w.o
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o bspeval.oct  low_level_functions.o /tmp/oct-52ISI0.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o curvederivcpts.oct  low_level_functions.o /tmp/oct-Gkcuyg.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o surfderivcpts.oct  low_level_functions.o /tmp/oct-CoKrJO.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o bspderiv.oct  low_level_functions.o /tmp/oct-w4VyYs.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o basisfunder.oct  low_level_functions.o /tmp/oct-9kukhJ.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o surfderiveval.oct  low_level_functions.o /tmp/oct-nuHEv8.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o basisfun.oct  low_level_functions.o /tmp/oct-8vkCYU.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o nrbsurfderiveval.oct  low_level_functions.o /tmp/oct-BfsaPW.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o tbasisfun.oct  low_level_functions.o /tmp/oct-NDNR1w.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o nrb_srf_basisfun__.oct  low_level_functions.o /tmp/oct-0g4eSK.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
+g++ -I/usr/include/octave-8.4.0/octave/.. -I/usr/include/octave-8.4.0/octave  -pthread -fopenmp -mieee-fp -g -O2 -ffile-prefix-map=/build/reproducible-path/octave-nurbs-1.4.3=. -fstack-protector-strong -Wformat -Werror=format-security  -o nrb_srf_basisfun_der__.oct  low_level_functions.o /tmp/oct-Fv8kMX.o   -shared -Wl,-Bsymbolic -Wl,-z,relro  -L/usr/lib/i386-linux-gnu -shared -Wl,-Bsymbolic  -flto=auto -ffat-lto-objects -Wl,-z,relro 
 make[1]: Leaving directory '/build/reproducible-path/octave-nurbs-1.4.3/src'
 copyfile /build/reproducible-path/octave-nurbs-1.4.3/./src/basisfun.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/basisfunder.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/bspderiv.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/bspeval.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/curvederivcpts.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/nrb_srf_basisfun__.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/nrb_srf_basisfun_der__.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/nrbsurfderiveval.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/surfderivcpts.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/surfderiveval.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/tbasisfun.oct /build/reproducible-path/octave-nurbs-1.4.3/./src/basisfun.cc-tst /build/reproducible-path/octave-nurbs-1.4.3/./src/curvederivcpts.cc-tst /build/reproducible-path/octave-nurbs-1.4.3/./src/nrbsurfderiveval.cc-tst /build/reproducible-path/octave-nurbs-1.4.3/./src/surfderivcpts.cc-tst /build/reproducible-path/octave-nurbs-1.4.3/./src/surfderiveval.cc-tst /build/reproducible-path/octave-nurbs-1.4.3/./src/tbasisfun.cc-tst /build/reproducible-path/octave-nurbs-1.4.3/./inst/i686-pc-linux-gnu-api-v58
 For information about changes from previous versions of the nurbs package, run 'news nurbs'.
@@ -2783,68 +2819,372 @@
 Checking package...
 Run the unit tests...
 Checking m files ...
-[inst/nrbeval.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbeval.m
+[inst/nrbbasisfunder.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbbasisfunder.m
 ***** demo
- srf = nrbtestsrf;
- p = nrbeval(srf,{linspace(0.0,1.0,20) linspace(0.0,1.0,20)});
- h = surf(squeeze(p(1,:,:)),squeeze(p(2,:,:)),squeeze(p(3,:,:)));
- title('Test surface.');
+ U = [0 0 0 0 1 1 1 1];
+ x = [0 1/3 2/3 1] ;
+ y = [0 0 0 0];
+ w = [1 1 1 1];
+ nrb = nrbmak ([x;y;y;w], U);
+ u = linspace(0, 1, 30);
+ [Bu, id] = nrbbasisfunder (u, nrb);
+ plot(u, Bu)
+ title('Derivatives of the cubic Bernstein polynomials')
  hold off
 ***** test
- knots{1} = [0 0 0 1 1 1];
- knots{2} = [0 0 0 .5 1 1 1];
- knots{3} = [0 0 0 0 1 1 1 1];
- cx = [0 0.5 1]; nx = length(cx);
- cy = [0 0.25 0.75 1]; ny = length(cy);
- cz = [0 1/3 2/3 1]; nz = length(cz);
- coefs(1,:,:,:) = repmat(reshape(cx,nx,1,1),[1 ny nz]);
- coefs(2,:,:,:) = repmat(reshape(cy,1,ny,1),[nx 1 nz]);
- coefs(3,:,:,:) = repmat(reshape(cz,1,1,nz),[nx ny 1]);
- coefs(4,:,:,:) = 1;
- nurbs = nrbmak(coefs, knots);
- x = rand(5,1); y = rand(5,1); z = rand(5,1);
- tt = [x y z]';
- points = nrbeval(nurbs,tt);
-
- assert(points,tt,1e-10)
+ U = [0 0 0 0 1 1 1 1];
+ x = [0 1/3 2/3 1] ;
+ y = [0 0 0 0];
+ w = rand(1,4);
+ nrb = nrbmak ([x;y;y;w], U);
+ u = linspace(0, 1, 30);
+ [Bu, id] = nrbbasisfunder (u, nrb);
+ #plot(u, Bu)
+ assert (sum(Bu, 2), zeros(numel(u), 1), 1e-10), 
+***** test
+ U = [0 0 0 0 1/2 1 1 1 1];
+ x = [0 1/4 1/2 3/4 1] ;
+ y = [0 0 0 0 0];
+ w = rand(1,5);
+ nrb = nrbmak ([x;y;y;w], U);
+ u = linspace(0, 1, 300); 
+ [Bu, id] = nrbbasisfunder (u, nrb); 
+ assert (sum(Bu, 2), zeros(numel(u), 1), 1e-10)
+***** test
+ p = 2;   q = 3;   m = 4; n = 5;
+ Lx  = 1; Ly  = 1; 
+ nrb = nrb4surf   ([0 0], [1 0], [0 1], [1 1]);
+ nrb = nrbdegelev (nrb, [p-1, q-1]);
+ aux1 = linspace(0,1,m); aux2 = linspace(0,1,n);
+ nrb = nrbkntins  (nrb, {aux1(2:end-1), aux2(2:end-1)});
+ nrb.coefs (4,:,:) = nrb.coefs(4,:,:) + rand (size (nrb.coefs (4,:,:)));
+ [Bu, Bv, N] = nrbbasisfunder ({rand(1, 20), rand(1, 20)}, nrb);
+ #plot3(squeeze(u(1,:,:)), squeeze(u(2,:,:)), reshape(Bu(:,10), 20, 20),'o')
+ assert (sum (Bu, 2), zeros(20^2, 1), 1e-10)
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/findspan.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/findspan.m
+***** test
+  n = 3; 
+  U = [0 0 0 1/2 1 1 1]; 
+  p = 2; 
+  u = linspace(0, 1, 10);  
+  s = findspan (n, p, u, U);
+  assert (s, [2*ones(1, 5) 3*ones(1, 5)]);
+***** test
+ p = 2; m = 7; n = m - p - 1;
+ U = [zeros(1,p)  linspace(0,1,m+1-2*p) ones(1,p)];
+ u = [ 0   0.11880   0.55118   0.93141   0.40068   0.35492 0.44392   0.88360   0.35414   0.92186   0.83085   1];
+ s = [2   2   3   4   3   3   3   4   3   4   4   4];
+ assert (findspan (n, p, u, U), s, 1e-10);
+2 tests, 2 passed, 0 known failure, 0 skipped
+[inst/nrbcrvderiveval.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbcrvderiveval.m
+***** test
+ knots = [0 0 0 1 1 1];
+ coefs(:,1) = [0; 0; 0; 1];
+ coefs(:,2) = [1; 0; 1; 1];
+ coefs(:,3) = [1; 1; 1; 2];
+ crv = nrbmak (coefs, knots);
+ u = linspace (0, 1, 100);
+ ck = nrbcrvderiveval (crv, u, 2); 
+ w  = @(x) 1 + x.^2;
+ dw = @(x) 2*x;
+ F1 = @(x) (2*x - x.^2)./w(x);
+ F2 = @(x) x.^2./w(x);
+ F3 = @(x) (2*x - x.^2)./w(x);
+ dF1 = @(x) (2 - 2*x)./w(x) - 2*(2*x - x.^2).*x./w(x).^2;
+ dF2 = @(x) 2*x./w(x) - 2*x.^3./w(x).^2;
+ dF3 = @(x) (2 - 2*x)./w(x) - 2*(2*x - x.^2).*x./w(x).^2;
+ d2F1 = @(x) -2./w(x) - 2*x.*(2-2*x)./w(x).^2 - (8*x-6*x.^2)./w(x).^2 + 8*x.^2.*(2*x-x.^2)./w(x).^3;
+ d2F2 = @(x) 2./w(x) - 4*x.^2./w(x).^2 - 6*x.^2./w(x).^2 + 8*x.^4./w(x).^3;
+ d2F3 = @(x) -2./w(x) - 2*x.*(2-2*x)./w(x).^2 - (8*x-6*x.^2)./w(x).^2 + 8*x.^2.*(2*x-x.^2)./w(x).^3;
+ assert ([F1(u); F2(u); F3(u)], squeeze(ck(:, 1, :)), 1e2*eps);
+ assert ([dF1(u); dF2(u); dF3(u)], squeeze(ck(:, 2, :)), 1e2*eps);
+ assert ([d2F1(u); d2F2(u); d2F3(u)], squeeze(ck(:, 3, :)), 1e2*eps);                                   
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/crvkntremove.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/crvkntremove.m
+***** test 
+ crv  = nrbdegelev (nrbline (), 3);
+ acrv = nrbkntins (crv, [.11 .11 .11]);
+ [rcrv, t] = crvkntremove (acrv, .11, 8, 3, 3, 1e-10);
+ assert (crv.knots, rcrv.knots, 1e-10);
+ assert (t, 3);
+***** test 
+ crv  = nrbcirc ();
+ acrv = nrbkntins (crv, [.3 .3]);
+ [rcrv, t] = crvkntremove (acrv, .3, 7, 2, 2, 1e-10);
+ assert (crv.knots, rcrv.knots, 1e-10);
+ assert (t, 2);
+2 tests, 2 passed, 0 known failure, 0 skipped
+[inst/kntbrkdegmult.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/kntbrkdegmult.m
+***** test
+ breaks = [0 1 2 3 4];
+ degree = 3;
+ knots = kntbrkdegmult (breaks, degree);
+ assert (knots, [0 0 0 0 1 2 3 4 4 4 4])
+***** test
+ breaks = [0 1 2 3 4];
+ degree = 3;
+ mult   = 2;
+ knots = kntbrkdegmult (breaks, degree, mult);
+ assert (knots, [0 0 0 0 1 1 2 2 3 3 4 4 4 4])
+***** test
+ breaks = [0 1 2 3 4];
+ degree = 3;
+ mult   = [1 2 3];
+ knots = kntbrkdegmult (breaks, degree, mult);
+ assert (knots, [0 0 0 0 1 2 2 3 3 3 4 4 4 4])
+***** test
+ breaks = {[0 1 2 3 4] [0 1 2 3]};
+ degree = [3 2];
+ mult   = {[1 2 3] 2};
+ knots = kntbrkdegmult (breaks, degree, mult);
+ assert (knots, {[0 0 0 0 1 2 2 3 3 3 4 4 4 4] [0 0 0 1 1 2 2 3 3 3]})
+4 tests, 4 passed, 0 known failure, 0 skipped
+[inst/nrbeval_der_w.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbeval_der_w.m
+***** test % 1D
+ nrb = nrbkntins (nrbcirc (1, [0 0], 0, pi/2), .5);
+ u = linspace (0, 1, 11);
+ delta_w = .01;
+ n = nrb.number;
+ der_ex = zeros (3, numel (u), n);
+ der_fd = zeros (3, numel (u), n);
+ for iw = 1:n
+   new_w1   = nrb.coefs (4, iw) + delta_w;
+   new_w2   = nrb.coefs (4, iw) - delta_w;
+   nrb1 = nrbmodw (nrb, new_w1, iw);
+   nrb2 = nrbmodw (nrb, new_w2, iw);
+   der_ex(:,:,iw) = nrbeval_der_w (nrb, iw, u);
+   p2 = nrbeval (nrb2, u);
+   p1 = nrbeval (nrb1, u);
+   der_fd(:,:,iw) = -(p2 - p1) ./ (2*delta_w);  
+ end
+ error = max (abs (der_ex(:) - der_fd(:)));
+ assert (error < 1.e-4)
+
+***** test %2D
+ crv = nrbline([1 0], [2 0]);
+ nrb = nrbtransp (nrbrevolve (crv, [], [0 0 1], pi/2));
+ new_knots = linspace (1/9, 8/9, 8);
+ nrb = nrbkntins (nrb, {new_knots, new_knots});
+ u = linspace (0, 1, 5);
+ v = u;
+ delta_w = .01;
+ n = nrb.number(1) * nrb.number(2);
+ der_ex = zeros (3, numel(u)* numel(v), n);
+ der_fd = zeros (3, numel(u)* numel(v), n);
+ for iw = 1:prod(nrb.number)
+   new_w1   = nrb.coefs (4, iw) + delta_w;
+   new_w2   = nrb.coefs (4, iw) - delta_w;
+   nrb1 = nrbmodw (nrb, new_w1, iw);
+   nrb2 = nrbmodw (nrb, new_w2, iw);
+   der_ex(:,:,iw) = nrbeval_der_w (nrb, iw, {u v});
+   p2 = nrbeval (nrb2, {u v});
+   p1 = nrbeval (nrb1, {u v});
+   der_fd(:,:,iw) = reshape (-(p2 - p1) ./ (2*delta_w), 3, []);
+ end
+ error = max (abs (der_ex(:) - der_fd(:)));
+ assert (error < 1.e-5)
+ 
+***** test % 3D
+ crv = nrbline([1 0], [2 0]);
+ nrb = nrbtransp (nrbrevolve (crv, [], [0 0 1], pi/2));
+ nrb = nrbextrude (nrb, [0 0 1]);
+ u = 0:.33:.99;
+ v = 0:.1:.9;
+ w = [.25 .5 .75];
+ delta_w = .01;
+ n = nrb.number(1) * nrb.number(2) * nrb.number(3);
+ der_ex = zeros (3, numel(u)*numel(v)*numel(w), n);
+ der_fd = zeros (3, numel(u)*numel(v)*numel(w), n);
+ for iw = 1:prod(nrb.number)
+   new_w1   = nrb.coefs (4, iw) + delta_w;
+   new_w2   = nrb.coefs (4, iw) - delta_w;
+   nrb1 = nrbmodw (nrb, new_w1, iw);
+   nrb2 = nrbmodw (nrb, new_w2, iw);
+   der_ex(:,:,iw) = nrbeval_der_w (nrb, iw, {u v w});
+   p2 = nrbeval (nrb2, {u v w});
+   p1 = nrbeval (nrb1, {u v w});
+   der_fd(:,:,iw) = reshape (-(p2 - p1) ./ (2*delta_w), 3, []);
+ end
+ error = max (max (squeeze (max (abs (der_ex - der_fd)))));
+ assert (error < 1.e-4)
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/nrbtransp.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbtransp.m
+***** demo
+ srf = nrb4surf([0 0 0], [1 0 1], [0 1 1], [1 1 2]);
+ nrbplot(srf,[20 5]);
+ title('Plane surface and its transposed (translated)')
+ hold on
+ srf.coefs(3,:,:) = srf.coefs(3,:,:) + 10;
+ srf = nrbtransp(srf);
+ nrbplot(srf,[20 5]);
+ hold off
+***** test
+ srf  = nrbrevolve(nrbline([1 0],[2 0]), [0 0 0], [0 0 1], pi/2);
+ srft = nrbtransp(srf);
+ assert (srf.number, fliplr(srft.number));
+ assert (srf.order, fliplr(srft.order));
+ assert (srf.knots, fliplr(srft.knots));
+ assert (srf.coefs, permute(srft.coefs, [1 3 2]));
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/aveknt.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/aveknt.m
+***** test
+ knt = [0 0 0 0.5 1 1 1];
+ pts = aveknt (knt, 3);
+ assert (pts - [0 1/4 3/4 1] < 1e-14)
 
 ***** test
- knots{1} = [0 0 0 1 1 1];
- knots{2} = [0 0 0 0 1 1 1 1];
- knots{3} = [0 0 1 1];
- cx = [0 0 1]; nx = length(cx);
- cy = [0 0 0 1]; ny = length(cy);
- cz = [0 1]; nz = length(cz);
- coefs(1,:,:,:) = repmat(reshape(cx,nx,1,1),[1 ny nz]);
- coefs(2,:,:,:) = repmat(reshape(cy,1,ny,1),[nx 1 nz]);
- coefs(3,:,:,:) = repmat(reshape(cz,1,1,nz),[nx ny 1]);
- coefs(4,:,:,:) = 1;
- nurbs = nrbmak(coefs, knots);
- x = rand(5,1); y = rand(5,1); z = rand(5,1);
- tt = [x y z]';
- points = nrbeval(nurbs,tt);
- assert(points,[x.^2 y.^3 z]',1e-10);
+ knt = {[0 0 0 0.5 1 1 1] [0 0 0 0 1/3 2/3 1 1 1 1]};
+ pts = aveknt (knt, [3 4]);
+ assert (pts{1} - [0 1/4 3/4 1] < 1e-14);
+ assert (pts{2} - [0 1/9 1/3 2/3 8/9 1] < 1e-14);
 
 ***** test
- knots{1} = [0 0 0 1 1 1];
- knots{2} = [0 0 0 0 1 1 1 1];
- knots{3} = [0 0 1 1];
- cx = [0 0 1]; nx = length(cx);
- cy = [0 0 0 1]; ny = length(cy);
- cz = [0 1]; nz = length(cz);
- coefs(1,:,:,:) = repmat(reshape(cx,nx,1,1),[1 ny nz]);
- coefs(2,:,:,:) = repmat(reshape(cy,1,ny,1),[nx 1 nz]);
- coefs(3,:,:,:) = repmat(reshape(cz,1,1,nz),[nx ny 1]);
- coefs(4,:,:,:) = 1;
- coefs = coefs([2 1 3 4],:,:,:);
- nurbs = nrbmak(coefs, knots);
- x = rand(5,1); y = rand(5,1); z = rand(5,1);
- tt = [x y z]';
- points = nrbeval(nurbs,tt);
- [y.^3 x.^2 z]';
- assert(points,[y.^3 x.^2 z]',1e-10);
+ nrb = nrb4surf([0 0], [1 0], [0 1], [1 1]);
+ nrb = nrbkntins (nrbdegelev (nrb, [1 2]), {[1/2] [1/3 2/3]});
+ pts = aveknt (nrb);
+ assert (pts{1} - [0 1/4 3/4 1] < 1e-14);
+ assert (pts{2} - [0 1/9 1/3 2/3 8/9 1] < 1e-14);
 3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/curvederiveval.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/curvederiveval.m
+***** test
+ k = [0 0 0 1 1 1];
+ coefs(:,1) = [0;0;0;1];
+ coefs(:,2) = [1;0;1;1];
+ coefs(:,3) = [1;1;1;1];
+ crv = nrbmak (coefs, k);
+ ck = curvederiveval (crv.number-1, crv.order-1, crv.knots, squeeze (crv.coefs(1,:,:)), 0.5, 2);
+ assert(ck, [0.75; 1; -2]);
+ ck = curvederiveval (crv.number-1, crv.order-1, crv.knots, squeeze (crv.coefs(2,:,:)), 0.5, 2);
+ assert(ck, [0.25; 1; 2]);
+ ck = curvederiveval (crv.number-1, crv.order-1, crv.knots, squeeze (crv.coefs(3,:,:)), 0.5, 2);
+ assert(ck, [0.75; 1; -2]);
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/nrbinverse.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbinverse.m
+***** test
+ nrb = nrb4surf ([0 0], [1 0], [2 3], [5 4]);
+ p = nrbeval (nrb, {.25 .75});
+ u = nrbinverse (nrb, p, 'Display', false);
+ assert (norm (u - [.25; .75]) < 1e-8);
+
+***** test
+ nrb = nrb4surf ([0 0], [1 0], [2 3], [5 4]);
+ nrb = nrbdegelev (nrbextrude (nrb, [0 2 1]), [3 3 3]);
+ p = nrbeval (nrb, {.25 .75 .05});
+ u = nrbinverse (nrb, p, 'Display', false, 'TolX', 1e-12, 'TolFun', 1e-10);
+ assert (norm (u - [.25; .75; .05]) < 1e-8);
+
+2 tests, 2 passed, 0 known failure, 0 skipped
+[inst/kntbrkdegreg.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/kntbrkdegreg.m
+***** test
+ breaks = [0 1 2 3 4];
+ degree = 3;
+ knots = kntbrkdegreg (breaks, degree);
+ assert (knots, [0 0 0 0 1 2 3 4 4 4 4])
+***** test
+ breaks = [0 1 2 3 4];
+ degree = 3;
+ reg    = 1;
+ knots = kntbrkdegreg (breaks, degree, reg);
+ assert (knots, [0 0 0 0 1 1 2 2 3 3 4 4 4 4])
+***** test
+ breaks = [0 1 2 3 4];
+ degree = 3;
+ reg    = [0 1 2];
+ knots = kntbrkdegreg (breaks, degree, reg);
+ assert (knots, [0 0 0 0 1 1 1 2 2 3 4 4 4 4])
+***** test
+ breaks = {[0 1 2 3 4] [0 1 2 3]};
+ degree = [3 2];
+ reg    = {[0 1 2] 0};
+ knots = kntbrkdegreg (breaks, degree, reg);
+ assert (knots, {[0 0 0 0 1 1 1 2 2 3 4 4 4 4] [0 0 0 1 1 2 2 3 3 3]})
+4 tests, 4 passed, 0 known failure, 0 skipped
+[inst/surfderivcpts.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/surfderivcpts.m
+***** test
+ coefs = cat(3,[0 0; 0 1],[1 1; 0 1]);
+ knots = {[0 0 1 1]  [0 0 1 1]};
+ plane = nrbmak(coefs,knots);
+ pkl = surfderivcpts (plane.number(1)-1, plane.order(1)-1,...
+                       plane.knots{1}, plane.number(2)-1,...
+                       plane.order(2)-1, plane.knots{2}, ...
+                       squeeze (plane.coefs(1,:,:)), 1);
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/numbasisfun.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/numbasisfun.m
+***** test
+  n = 3; 
+  U = [0 0 0 1/2 1 1 1]; 
+  p = 2; 
+  u = linspace (0, 1, 10);  
+  s = findspan (n, p, u, U); 
+  Bref = [0   0   0   0   0   1   1   1   1   1; ...
+          1   1   1   1   1   2   2   2   2   2; ...
+          2   2   2   2   2   3   3   3   3   3].';
+  B = numbasisfun (s, u, p, U);
+  assert (B, Bref)
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/nrbclamp.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbclamp.m
+***** test
+ crv = nrbdegelev (nrbcirc (1, [], 0, pi/2), 2);
+ crv = nrbunclamp (crv, 3);
+ xx = linspace (0, 1, 20);
+ crv1 = nrbclamp (crv);
+ assert (crv1.knots, [0 0 0 0 0 1 1 1 1 1])
+ assert (nrbeval(crv, xx), nrbeval(crv1, xx), 1e-14)
+ crv1 = nrbclamp (crv, 2);
+ assert (crv1.knots, [-3 -2 -1 0 0 1 1 2 3 4])
+ assert (nrbeval(crv, xx), nrbeval(crv1, xx), 1e-14)
+***** test
+ crv1 = nrbcirc(1,[],0,pi/4);
+ crv2 = nrbcirc(2,[],0,pi/4);
+ srf = nrbkntins (nrbdegelev (nrbruled(crv1, crv2), [3 2]), {0.25 []});
+ srf = nrbunclamp (srf, [4 2]);
+ srf1 = nrbclamp (srf);
+ xx = linspace(0,1,20);
+ assert(srf1.knots, {[0 0 0 0 0 0 0.2500 1 1 1 1 1 1]  [0 0 0 0 1 1 1 1]})
+ assert (nrbeval(srf, {xx xx}), nrbeval(srf1, {xx xx}), 1e-14);
+ srf1 = nrbclamp (srf, [3 1]);
+ assert (srf1.knots, {[-2 -1.75 -1 -0.75 0 0 0.25 1 1 1.25 2 2.25 3], [-2 -1 0 0 1 1 2 3]})
+ assert (nrbeval(srf, {xx xx}), nrbeval(srf1, {xx xx}), 1e-14);
+ srf1 = nrbclamp (srf, [], 2);
+ assert(srf1.knots, {[-2.75 -2 -1.75 -1 -0.75 0 0.25 1 1.25 2 2.25 3 3.25] [0 0 0 0 1 1 1 1]})
+ assert (nrbeval(srf, {xx xx}), nrbeval(srf1, {xx xx}), 1e-14);
+2 tests, 2 passed, 0 known failure, 0 skipped
+[inst/basisfun.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/basisfun.m
+***** test
+  n = 3; 
+  U = [0 0 0 1/2 1 1 1]; 
+  p = 2; 
+  u = linspace (0, 1, 10);  
+  s = findspan (n, p, u, U);  
+  Bref = [1.00000   0.00000   0.00000
+          0.60494   0.37037   0.02469
+          0.30864   0.59259   0.09877
+          0.11111   0.66667   0.22222
+          0.01235   0.59259   0.39506
+          0.39506   0.59259   0.01235
+          0.22222   0.66667   0.11111
+          0.09877   0.59259   0.30864
+          0.02469   0.37037   0.60494
+          0.00000   0.00000   1.00000];
+  B = basisfun (s, u, p, U);
+  assert (B, Bref, 1e-5);
+1 test, 1 passed, 0 known failure, 0 skipped
 [inst/nrbpermute.m]
 >>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbpermute.m
 ***** demo
@@ -2871,63 +3211,196 @@
  assert ({vol.knots{perm2}}, vol3.knots)
  assert (permute(vol.coefs, [1, perm2+1]), vol3.coefs)
 1 test, 1 passed, 0 known failure, 0 skipped
-[inst/nrbeval_der_p.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbeval_der_p.m
-***** test %% 1D
- nrb =   nrbkntins (nrbcirc (1, [0 0], 0, pi/2), .5);
- u = 0:.1:.9;
- index = 1:nrb.number;
- e = zeros (numel (u), numel (index), 1);
- for jj = 1:numel (index)
-   deltap = .1 * rand (3, 1);
-   nrb2 = nrbmodp (nrb, deltap, index(jj));
-   der_ex = nrbeval_der_p (nrb, index(jj), u);
-   p2 = nrbeval (nrb2, u);
-   p1 = nrbeval (nrb, u);
-   der_fd = (p2 - p1) ./ deltap;
-   e(:,jj) = sqrt (sum ((repmat (der_ex, 3, 1) - der_fd).^2, 1));
- end
- assert (max(e(:)) < 1.e-8);
- 
-***** test %% 2D
- crv = nrbline([1 0], [2 0]);
- nrb = nrbtransp (nrbrevolve (crv, [], [0 0 1], pi/2));
- new_knots = linspace (1/9, 8/9, 8);
- nrb = nrbkntins (nrb, {new_knots, new_knots});
- u = 0:.1:.9;
- v = u;
- e = zeros (nrb.number(1) * nrb.number(2), numel (u), numel (v));
- for index = 1:prod(nrb.number)
-   deltap = .1 * rand (3, 1);
-   nrb2 = nrbmodp (nrb, deltap, index);
-   der_ex = nrbeval_der_p (nrb, index, {u v});
-   p2 = nrbeval (nrb2, {u v});
-   p1 = nrbeval (nrb, {u v});
-   der_fd = (p2 - p1) ./ deltap;
-   der_ex = reshape (repmat (der_ex, 3, 1), size(der_fd));
-   e(index,:,:) = sqrt (sum ((der_ex - der_fd).^2, 1));
- end
- assert (max(e(:)) < 1.e-8)
- 
-***** test %% 3D
- crv = nrbline([1 0], [2 0]);
- nrb = nrbtransp (nrbrevolve (crv, [], [0 0 1], pi/2));
- nrb = nrbextrude (nrb, [0 0 1]);
- u = 0:.1:.9;
- v = u;
- w = u;
- e = zeros (nrb.number(1) * nrb.number(2) * nrb.number(3), numel(u), numel(v), numel(w));
- for index = 1:prod(nrb.number)
-   deltap = .1 * rand (3, 1);
-   nrb2 = nrbmodp (nrb, deltap, index);
-   der_ex = nrbeval_der_p (nrb, index, {u v w});
-   p2 = nrbeval (nrb2, {u v w});
-   p1 = nrbeval (nrb, {u v w});
-   der_fd = (p2 - p1) ./ deltap;
-   der_ex = reshape (repmat (der_ex, 3, 1), size (der_fd));
-   e(index,:,:,:) = sqrt (sum ((der_ex - der_fd).^2, 1));
- end
- assert (max (e(:)) < 1.e-8);
+[inst/nrbreverse.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbreverse.m
+***** demo
+ pnts = [0.5 1.5 3.0 7.5 8.5;
+         3.0 5.5 1.5 4.0 4.5;
+         0.0 0.0 0.0 0.0 0.0];
+ crv1 = nrbmak(pnts,[0 0 0 1/2 3/4 1 1 1]);
+ crv2 = nrbreverse(crv1);
+ fprintf('Knots of the original curve\n')
+ disp(crv1.knots)
+ fprintf('Knots of the reversed curve\n')
+ disp(crv2.knots)
+ fprintf('Control points of the original curve\n')
+ disp(crv1.coefs(1:2,:))
+ fprintf('Control points of the reversed curve\n')
+ disp(crv2.coefs(1:2,:))
+ nrbplot(crv1,100)
+ hold on
+ nrbplot(crv2,100)
+ title('The curve and its reverse are the same')
+ hold off
+***** test
+ srf  = nrbrevolve(nrbline([1 0],[2 0]), [0 0 0], [0 0 1], pi/2);
+ srf  = nrbkntins (srf, {0.3, 0.6});
+ srf2 = nrbreverse (srf);
+ assert (srf.knots, cellfun(@(x) sort(1-x), srf2.knots, 'UniformOutput', false), 1e-15)
+ assert (srf.coefs, srf2.coefs(:,end:-1:1,end:-1:1))
+***** test
+ srf  = nrbrevolve(nrbline([1 0],[2 0]), [0 0 0], [0 0 1], pi/2);
+ srf  = nrbkntins (srf, {0.3, 0.6});
+ srf2 = nrbreverse (srf, 1);
+ knt{1} = sort(1-srf2.knots{1}); knt{2} = srf2.knots{2};
+ assert (srf.knots, knt, 1e-15)
+ assert (srf.coefs, srf2.coefs(:,end:-1:1,:))
+ srf2 = nrbreverse (srf, 2);
+ knt{1} = srf2.knots{1}; knt{2} = sort(1-srf2.knots{2});
+ assert (srf.knots, knt, 1e-15)
+ assert (srf.coefs, srf2.coefs(:,:,end:-1:1))
+2 tests, 2 passed, 0 known failure, 0 skipped
+[inst/nrbnumbasisfun.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbnumbasisfun.m
+***** test
+ p = 2;   q = 3;   m = 4; n = 5;
+ Lx  = 1; Ly  = 1; 
+ nrb = nrb4surf   ([0 0], [1 0], [0 1], [1 1]);
+ nrb = nrbdegelev (nrb, [p-1, q-1]);
+ ikx = linspace(0,1,m); iky = linspace(0,1,n);
+ nrb = nrbkntins  (nrb, {ikx(2:end-1), iky(2:end-1)});
+ nrb.coefs (4,:,:) = nrb.coefs (4,:,:) + rand (size (nrb.coefs (4,:,:)));
+ u = rand (1, 30); v = rand (1, 10);
+ u = (u-min (u))/max (u-min (u));
+ v = (v-min (v))/max (v-min (v));
+ N = nrbnumbasisfun ({u, v}, nrb);
+ assert (all (all (N>0)), true)
+ assert (all (all (N <= prod (nrb.number))), true)
+ assert (max (max (N)), prod (nrb.number))
+ assert (min (min (N)), 1)
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/tbasisfun.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/tbasisfun.m
+***** demo
+ U = {[0 0 1/2 1 1], [0 0 0 1 1]};
+ p = [3, 3];
+ [X, Y] = meshgrid (linspace(0, 1, 30));
+ u = [X(:), Y(:)]';
+ N = tbasisfun (u, p, U);
+ surf (X, Y, reshape (N, size(X)))
+ title('Basis function associated to a local knot vector')
+ hold off
+***** test
+ U = [0 1/2 1];
+ p = 1;
+ u = [0.3 0.4 0.6 0.7];
+ [N, Nder] = tbasisfun (u, p, U);
+ assert (N, [0.6 0.8 0.8 0.6], 1e-12);
+ assert (Nder, [2 2 -2 -2], 1e-12);
+***** test
+ U = {[0 1/2 1] [0 1/2 1]};
+ p = [1 1];
+ u = [0.3 0.4 0.6 0.7; 0.3 0.4 0.6 0.7];
+ [N, Nder] = tbasisfun (u, p, U);
+ assert (N, [0.36 0.64 0.64 0.36], 1e-12);
+ assert (Nder, [1.2 1.6 -1.6 -1.2; 1.2 1.6 -1.6 -1.2], 1e-12);
+***** test
+ U = {[0 1/2 1] [0 1/2 1] [0 1/2 1]};
+ p = [1 1 1];
+ u = [0.4 0.4 0.6 0.6; 0.4 0.4 0.6 0.6; 0.4 0.6 0.4 0.6];
+ [N, Nder] = tbasisfun (u, p, U);
+ assert (N, [0.512 0.512 0.512 0.512], 1e-12);
+ assert (Nder, [1.28 1.28 -1.28 -1.28; 1.28 1.28 -1.28 -1.28; 1.28 -1.28 1.28 -1.28], 1e-12);
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/surfderiveval.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/surfderiveval.m
+***** shared srf
+***** test
+ k = [0 0 0 1 1 1];
+ c = [0 1/2 1];
+ [coef(2,:,:), coef(1,:,:)] = meshgrid (c, c);
+ srf = nrbmak (coef, {k, k});
+ skl = surfderiveval (srf.number(1)-1, ...
+                      srf.order(1)-1, ...
+                      srf.knots{1}, ...
+                      srf.number(2)-1, ...
+                      srf.order(2)-1, ...
+                      srf.knots{2},...
+                      squeeze(srf.coefs(1,:,:)), .5, .5, 1) ;
+ assert (skl, [.5 0; 1 0])
+***** test
+ srf = nrbkntins (srf, {[], rand(1,2)});
+ skl = surfderiveval (srf.number(1)-1,... 
+                      srf.order(1)-1, ...
+                      srf.knots{1},...
+                      srf.number(2)-1,... 
+                      srf.order(2)-1, ...
+                      srf.knots{2},...
+                      squeeze(srf.coefs(1,:,:)), .5, .5, 1) ;
+ assert (skl, [.5 0; 1 0], 100*eps)
+2 tests, 2 passed, 0 known failure, 0 skipped
+[inst/kntrefine.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/kntrefine.m
+***** shared nrbs
+***** test
+ knots = {[0 0 1 1] [0 0 0 1 1 1]};
+ coefs(1,:,:) = [1 sqrt(2)/2 0; 2 sqrt(2) 0];
+ coefs(2,:,:) = [0 sqrt(2)/2 1; 0 sqrt(2) 2];
+ coefs(4,:,:) = [1 sqrt(2)/2 1; 1 sqrt(2)/2 1];
+ nrbs = nrbmak (coefs, knots);
+ nrbs = nrbkntins (nrbs, {[] [0.5 0.6 0.6]});
+ nrbs = nrbdegelev (nrbs, [0 1]);
+ nrbs = nrbkntins (nrbs, {[] [0.4]});
+ rknots = kntrefine (nrbs.knots, [1 1], [1 1], [0 0]);
+ assert (rknots{1} == [0 0 0.5 1 1]);
+ assert (rknots{2} == [0 0 0.2 0.4 0.45 0.5 0.55 0.6 0.8 1 1]);
+
+***** test
+ rknots = kntrefine (nrbs.knots, [1 1], [3 3], [0 0]);
+ assert (rknots{1}, [0 0 0 0 0.5 0.5 0.5 1 1 1 1]);
+ assert (rknots{2}, [0 0 0 0 0.2 0.2 0.2 0.4 0.4 0.4 0.45 0.45 0.45 0.5 0.5 0.5 0.55 0.55 0.55 0.6 0.6 0.6 0.8 0.8 0.8 1 1 1 1]);
+
+***** test
+ rknots = kntrefine (nrbs.knots, [1 1], [3 3], [2 2]);
+ assert (rknots{1}, [0 0 0 0 0.5 1 1 1 1]);
+ assert (rknots{2}, [0 0 0 0 0.2 0.4 0.45 0.5 0.5 0.55 0.6 0.6 0.6 0.8 1 1 1 1]);
+
+***** test
+ rknots = kntrefine (nrbs.knots, [1 1], [4 4], [0 0]);
+ assert (rknots{1}, [0 0 0 0 0 0.5 0.5 0.5 0.5 1 1 1 1 1]);
+ assert (rknots{2}, [0 0 0 0 0 0.2 0.2 0.2 0.2 0.4 0.4 0.4 0.4 0.45 0.45 0.45 0.45 0.5 0.5 0.5 0.5 0.55 0.55 0.55 0.55 0.6 0.6 0.6 0.6 0.8 0.8 0.8 0.8 1 1 1 1 1]);
+
+***** test
+ rknots = kntrefine (nrbs.knots, [1 1], [4 4], [3 3]);
+ assert (rknots{1}, [0 0 0 0 0 0.5 1 1 1 1 1]);
+ assert (rknots{2}, [0 0 0 0 0 0.2 0.4 0.4 0.45 0.5 0.5 0.5 0.55 0.6 0.6 0.6 0.6 0.8 1 1 1 1 1]);
+
+***** test
+ knots = [0 0 0 0 0.4 0.5 0.5 0.6 0.6 0.6 1 1 1 1];
+ rknots = kntrefine (knots, 1, 4, 3);
+ assert (rknots, [0 0 0 0 0 0.2 0.4 0.4 0.45 0.5 0.5 0.5 0.55 0.6 0.6 0.6 0.6 0.8 1 1 1 1 1]);
+6 tests, 6 passed, 0 known failure, 0 skipped
+[inst/curvederivcpts.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/curvederivcpts.m
+***** test
+ line = nrbmak([0.0 1.5; 0.0 3.0],[0.0 0.0 1.0 1.0]);
+ pk   = curvederivcpts (line.number-1, line.order-1, line.knots,...
+                        line.coefs(1,:), 2);
+ assert (pk, [0 3/2; 3/2 0], 100*eps);
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/nrbmeasure.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbmeasure.m
+***** test 
+ c = nrbcirc (1, [0 0], 0,  pi/3);
+ l = nrbmeasure(c, 0, 1, 1e-7);
+ assert (l, pi/3, 1e-7)
+***** test 
+ c = nrbcirc (1, [0 0], 0,  pi/2);
+ s = zeros (1, 100); e = linspace (0, 1, 100);
+ for ii = 1:100
+   l(ii) = nrbmeasure (c, s(ii), e(ii), 1e-7);
+ endfor
+ xx = nrbeval (c, e);
+ theta = atan2 (xx(2,:), xx(1,:));
+ assert (l, theta, 1e-7)
+***** test 
+ c = nrbcirc (1, [0 0], 0,  pi/2);
+ s = 0; e = linspace (0, 1, 100);
+ for ii = 1:100
+   l(ii) = nrbmeasure (c, s, e(ii), 1e-7);
+ endfor
+ l2 = nrbmeasure (c, s, e, 1e-7);
+ assert (l, l2, eps)
 3 tests, 3 passed, 0 known failure, 0 skipped
 [inst/nrbbasisfun.m]
 >>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbbasisfun.m
@@ -3018,330 +3491,6 @@
  assert (max (max (N)) == prod (nrb.number))
  assert (min (min (N))== 1)
 3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/nrbsquare.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbsquare.m
-***** test
- srf = nrbsquare ([], 1, 2, 2, 4);
- assert (srf.order, [3 3]);
- knt = [0 0 0 1/4 1/2 3/4 1 1 1];
- assert (srf.knots, {knt knt})
- x = linspace (0, 1, 100);
- [X,Y] = ndgrid (x, x);
- vals = nrbeval (srf, {x x});
- assert (squeeze(vals(1,:,:)), X, 1e-15);
- assert (squeeze(vals(2,:,:)), 2*Y, 1e-15);
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/crvkntremove.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/crvkntremove.m
-***** test 
- crv  = nrbdegelev (nrbline (), 3);
- acrv = nrbkntins (crv, [.11 .11 .11]);
- [rcrv, t] = crvkntremove (acrv, .11, 8, 3, 3, 1e-10);
- assert (crv.knots, rcrv.knots, 1e-10);
- assert (t, 3);
-***** test 
- crv  = nrbcirc ();
- acrv = nrbkntins (crv, [.3 .3]);
- [rcrv, t] = crvkntremove (acrv, .3, 7, 2, 2, 1e-10);
- assert (crv.knots, rcrv.knots, 1e-10);
- assert (t, 2);
-2 tests, 2 passed, 0 known failure, 0 skipped
-[inst/nrbcrvderiveval.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbcrvderiveval.m
-***** test
- knots = [0 0 0 1 1 1];
- coefs(:,1) = [0; 0; 0; 1];
- coefs(:,2) = [1; 0; 1; 1];
- coefs(:,3) = [1; 1; 1; 2];
- crv = nrbmak (coefs, knots);
- u = linspace (0, 1, 100);
- ck = nrbcrvderiveval (crv, u, 2); 
- w  = @(x) 1 + x.^2;
- dw = @(x) 2*x;
- F1 = @(x) (2*x - x.^2)./w(x);
- F2 = @(x) x.^2./w(x);
- F3 = @(x) (2*x - x.^2)./w(x);
- dF1 = @(x) (2 - 2*x)./w(x) - 2*(2*x - x.^2).*x./w(x).^2;
- dF2 = @(x) 2*x./w(x) - 2*x.^3./w(x).^2;
- dF3 = @(x) (2 - 2*x)./w(x) - 2*(2*x - x.^2).*x./w(x).^2;
- d2F1 = @(x) -2./w(x) - 2*x.*(2-2*x)./w(x).^2 - (8*x-6*x.^2)./w(x).^2 + 8*x.^2.*(2*x-x.^2)./w(x).^3;
- d2F2 = @(x) 2./w(x) - 4*x.^2./w(x).^2 - 6*x.^2./w(x).^2 + 8*x.^4./w(x).^3;
- d2F3 = @(x) -2./w(x) - 2*x.*(2-2*x)./w(x).^2 - (8*x-6*x.^2)./w(x).^2 + 8*x.^2.*(2*x-x.^2)./w(x).^3;
- assert ([F1(u); F2(u); F3(u)], squeeze(ck(:, 1, :)), 1e2*eps);
- assert ([dF1(u); dF2(u); dF3(u)], squeeze(ck(:, 2, :)), 1e2*eps);
- assert ([d2F1(u); d2F2(u); d2F3(u)], squeeze(ck(:, 3, :)), 1e2*eps);                                   
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/kntbrkdegreg.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/kntbrkdegreg.m
-***** test
- breaks = [0 1 2 3 4];
- degree = 3;
- knots = kntbrkdegreg (breaks, degree);
- assert (knots, [0 0 0 0 1 2 3 4 4 4 4])
-***** test
- breaks = [0 1 2 3 4];
- degree = 3;
- reg    = 1;
- knots = kntbrkdegreg (breaks, degree, reg);
- assert (knots, [0 0 0 0 1 1 2 2 3 3 4 4 4 4])
-***** test
- breaks = [0 1 2 3 4];
- degree = 3;
- reg    = [0 1 2];
- knots = kntbrkdegreg (breaks, degree, reg);
- assert (knots, [0 0 0 0 1 1 1 2 2 3 4 4 4 4])
-***** test
- breaks = {[0 1 2 3 4] [0 1 2 3]};
- degree = [3 2];
- reg    = {[0 1 2] 0};
- knots = kntbrkdegreg (breaks, degree, reg);
- assert (knots, {[0 0 0 0 1 1 1 2 2 3 4 4 4 4] [0 0 0 1 1 2 2 3 3 3]})
-4 tests, 4 passed, 0 known failure, 0 skipped
-[inst/nrbsurfderiveval.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbsurfderiveval.m
-***** test
- k = [0 0 1 1];
- c = [0 1];
- [coef(2,:,:), coef(1,:,:)] = meshgrid (c, c);
- coef(3,:,:) = coef(1,:,:);
- srf = nrbmak (coef, {k, k});
- [u, v] = meshgrid (linspace(0,1,11));
- uv = [u(:)';v(:)'];
- skl = nrbsurfderiveval (srf, uv, 0);
- aux = nrbeval(srf,uv);
- assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
-***** test
- k = [0 0 1 1];
- c = [0 1];
- [coef(2,:,:), coef(1,:,:)] = meshgrid (c, c);
- coef(3,:,:) = coef(1,:,:);
- srf = nrbmak (coef, {k, k});
- srf = nrbkntins (srf, {[], rand(2,1)});
- [u, v] = meshgrid (linspace(0,1,11));
- uv = [u(:)';v(:)'];
- skl = nrbsurfderiveval (srf, uv, 0);
- aux = nrbeval(srf,uv);
- assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
-***** shared srf, uv
-***** test 
- k = [0 0 0 1 1 1];
- c = [0 1/2 1];
- [coef(1,:,:), coef(2,:,:)] = meshgrid (c, c);
- coef(3,:,:) = coef(1,:,:);
- srf = nrbmak (coef, {k, k});
- ders= nrbderiv (srf);
- [u, v] = meshgrid (linspace(0,1,11));
- uv = [u(:)';v(:)'];
- skl = nrbsurfderiveval (srf, uv, 1);
- [fun, der] = nrbdeval (srf, ders, uv);
- assert (squeeze (skl (1:2,1,1,:)), fun(1:2,:), 1e3*eps)
- assert (squeeze (skl (1:2,2,1,:)), der{1}(1:2,:), 1e3*eps)
- assert (squeeze (skl (1:2,1,2,:)), der{2}(1:2,:), 1e3*eps)
-
-***** test 
- srf = nrbdegelev (srf, [3, 1]);
- ders= nrbderiv (srf);
- [fun, der] = nrbdeval (srf, ders, uv);
- skl = nrbsurfderiveval (srf, uv, 1);
- assert (squeeze (skl (1:2,1,1,:)), fun(1:2,:), 1e3*eps)
- assert (squeeze (skl (1:2,2,1,:)), der{1}(1:2,:), 1e3*eps)
- assert (squeeze (skl (1:2,1,2,:)), der{2}(1:2,:), 1e3*eps)
-***** shared uv
-***** test 
- k = [0 0 0 1 1 1];
- c = [0 1/2 1];
- [coef(2,:,:), coef(1,:,:)] = meshgrid (c, c);
- coef(3,:,:) = coef(1,:,:);
- srf = nrbmak (coef, {k, k});
- ders= nrbderiv (srf);
- [u, v] = meshgrid (linspace(0,1,11));
- uv = [u(:)';v(:)'];
- skl = nrbsurfderiveval (srf, uv, 1);
- [fun, der] = nrbdeval (srf, ders, uv);
- assert (squeeze (skl (1:2,1,1,:)), fun(1:2,:), 1e3*eps)
- assert (squeeze (skl (1:2,2,1,:)), der{1}(1:2,:), 1e3*eps)
- assert (squeeze (skl (1:2,1,2,:)), der{2}(1:2,:), 1e3*eps)
-
-***** test 
- p = 3;  q = 3;
- mcp = 5; ncp = 5;
- Lx  = 10*rand(1); Ly  = Lx;
- srf = nrbdegelev (nrb4surf ([0 0], [Lx, 0], [0 Ly], [Lx Ly]), [p-1, q-1]);
- %%srf = nrbkntins (srf, {linspace(0,1,mcp-p+2)(2:end-1), linspace(0,1,ncp-q+2)(2:end-1)});
- %%srf.coefs = permute (srf.coefs, [1 3 2]);
- ders= nrbderiv (srf);
- [fun, der] = nrbdeval (srf, ders, uv);
- skl = nrbsurfderiveval (srf, uv, 1);
- assert (squeeze (skl (1:2,1,1,:)), fun(1:2,:), 1e3*eps)
- assert (squeeze (skl (1:2,2,1,:)), der{1}(1:2,:), 1e3*eps)
- assert (squeeze (skl (1:2,1,2,:)), der{2}(1:2,:), 1e3*eps)
-***** shared srf, uv, P, dPdx, d2Pdx2, c1, c2
-***** test
- [u, v] = meshgrid (linspace(0,1,10));
- uv = [u(:)';v(:)'];
- c1 = nrbmak([0 1/2 1; 0 1 0],[0 0 0 1 1 1]);
- c1 = nrbtform (c1, vecrotx (pi/2));
- c2  = nrbtform(c1, vectrans([0 1 0]));
- srf = nrbdegelev (nrbruled (c1, c2), [3, 1]);
- skl = nrbsurfderiveval (srf, uv, 2);
- P = squeeze(skl(:,1,1,:));
- dPdx = squeeze(skl(:,2,1,:));
- d2Pdx2 = squeeze(skl(:,3,1,:));
-***** assert(P(3,:), 2*(P(1,:)-P(1,:).^2),100*eps)
-***** assert(dPdx(3,:), 2-4*P(1,:), 100*eps)
-***** assert(d2Pdx2(3,:), -4+0*P(1,:), 100*eps)
- srf = nrbdegelev (nrbruled (c1, c2), [5, 6]);
- skl = nrbsurfderiveval (srf, uv, 2);
- P = squeeze(skl(:,1,1,:));
- dPdx = squeeze(skl(:,2,1,:));
- d2Pdx2 = squeeze(skl(:,3,1,:));
- aux = nrbeval(srf,uv);
- assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
-***** assert(P(3,:), 2*(P(1,:)-P(1,:).^2),100*eps)
-***** assert(dPdx(3,:), 2-4*P(1,:), 100*eps)
-***** assert(d2Pdx2(3,:), -4+0*P(1,:), 100*eps)
-
-***** test
- skl = nrbsurfderiveval (srf, uv, 0);
- aux = nrbeval (srf, uv);
- assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
-***** shared dPdu, d2Pdu2, P, srf, uv
-***** test
- [u, v] = meshgrid (linspace(0,1,10));
- uv = [u(:)';v(:)'];
- c1 = nrbmak([0 1/2 1; 0.1 1.6 1.1; 0 0 0],[0 0 0 1 1 1]);
- c2 = nrbmak([0 1/2 1; 0.1 1.6 1.1; 1 1 1],[0 0 0 1 1 1]);
- srf = nrbdegelev (nrbruled (c1, c2), [0, 1]);
- skl = nrbsurfderiveval (srf, uv, 2);
- P = squeeze(skl(:,1,1,:));
- dPdu = squeeze(skl(:,2,1,:));
- dPdv = squeeze(skl(:,1,2,:));
- d2Pdu2 = squeeze(skl(:,3,1,:));
- aux = nrbeval(srf,uv);
- assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
-***** assert(dPdu(2,:), 3-4*P(1,:),100*eps)
-***** assert(d2Pdu2(2,:), -4+0*P(1,:),100*eps)
-
-***** test
- skl = nrbsurfderiveval (srf, uv, 0);
- aux = nrbeval(srf,uv);
- assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
-***** test
- srf = nrb4surf([0 0], [1 0], [0 1], [1 1]);
- geo = nrbdegelev (srf, [3 3]);
- geo = nrbkntins (geo, {[.1:.1:.9], [.2:.2:.8]});
- [u, v] = meshgrid (linspace(0,1,10));
- uv = [u(:)';v(:)'];
- skl = nrbsurfderiveval (geo, uv, 2);
- dgeo = nrbderiv (geo);
- [pnts, ders] = nrbdeval (geo, dgeo, uv);
- assert (ders{1}, squeeze(skl(:,2,1,:)), 1e-9)
- assert (ders{2}, squeeze(skl(:,1,2,:)), 1e-9)
-***** test
- crv = nrbline ([1 0], [2 0]);
- srf = nrbrevolve (crv, [0 0 0], [0 0 1], pi/2);
- srf = nrbtransp (srf);
- [v, u] = meshgrid (linspace (0, 1, 11));
- uv = [u(:)'; v(:)'];
- skl = nrbsurfderiveval (srf, uv, 2);
- c = sqrt(2);
- w      = @(x, y) (2 - c)*y.^2 + (c-2)*y + 1;
- dwdy   = @(x, y) 2*(2-c)*y + c - 2;
- d2wdy2 = @(x, y) 2*(2-c);
- F1 = @(x, y) (x+1) .* ((1-y).^2 + c*y.*(1-y)) ./ w(x,y);
- F2 = @(x, y) (x+1) .* (y.^2 + c*y.*(1-y)) ./ w(x,y);
- dF1dx = @(x, y) ((1-y).^2 + c*y.*(1-y)) ./ w(x,y);
- dF2dx = @(x, y) (y.^2 + c*y.*(1-y)) ./ w(x,y);
- dF1dy = @(x, y) (x+1) .* ((2 - 2*c)*y + c - 2) ./ w(x,y) - (x+1) .* ((1-y).^2 + c*y.*(1-y)) .* dwdy(x,y) ./ w(x,y).^2;
- dF2dy = @(x, y) (x+1) .* ((2 - 2*c)*y + c) ./ w(x,y) - (x+1) .* (y.^2 + c*y.*(1-y)) .* dwdy(x,y) ./ w(x,y).^2;
- d2F1dx2 = @(x, y) zeros (size (x));
- d2F2dx2 = @(x, y) zeros (size (x));
- d2F1dxdy = @(x, y) ((2 - 2*c)*y + c - 2) ./ w(x,y) - ((1-y).^2 + c*y.*(1-y)) .* dwdy(x,y) ./ w(x,y).^2;
- d2F2dxdy = @(x, y) ((2 - 2*c)*y + c) ./ w(x,y) - (y.^2 + c*y.*(1-y)) .* dwdy(x,y) ./ w(x,y).^2;
- d2F1dy2  = @(x, y) (x+1)*(2 - 2*c) ./ w(x,y) - 2*(x+1) .* ((2 - 2*c)*y + c - 2) .* dwdy(x,y) ./ w(x,y).^2 - ...
-                    (x+1) .* ((1-y).^2 + c*y.*(1-y)) * d2wdy2(x,y) ./ w(x,y).^2 + ...
-                    2 * (x+1) .* ((1-y).^2 + c*y.*(1-y)) .* w(x,y) .*dwdy(x,y).^2 ./ w(x,y).^4;
- d2F2dy2  = @(x, y) (x+1)*(2 - 2*c) ./ w(x,y) - 2*(x+1) .* ((2 - 2*c)*y + c) .* dwdy(x,y) ./ w(x,y).^2 - ...
-                    (x+1) .* (y.^2 + c*y.*(1-y)) * d2wdy2(x,y) ./ w(x,y).^2 + ...
-                    2 * (x+1) .* (y.^2 + c*y.*(1-y)) .* w(x,y) .*dwdy(x,y).^2 ./ w(x,y).^4;
- assert ([F1(u(:),v(:)), F2(u(:),v(:))], squeeze(skl(1:2,1,1,:))', 1e2*eps);
- assert ([dF1dx(u(:),v(:)), dF2dx(u(:),v(:))], squeeze(skl(1:2,2,1,:))', 1e2*eps);
- assert ([dF1dy(u(:),v(:)), dF2dy(u(:),v(:))], squeeze(skl(1:2,1,2,:))', 1e2*eps);
- assert ([d2F1dx2(u(:),v(:)), d2F2dx2(u(:),v(:))], squeeze(skl(1:2,3,1,:))', 1e2*eps);
- assert ([d2F1dxdy(u(:),v(:)), d2F2dxdy(u(:),v(:))], squeeze(skl(1:2,2,2,:))', 1e2*eps);
- assert ([d2F1dy2(u(:),v(:)), d2F2dy2(u(:),v(:))], squeeze(skl(1:2,1,3,:))', 1e2*eps);
-***** test
- knots = {[0 0 1 1] [0 0 1 1]};
- coefs(:,1,1) = [0;0;0;1];
- coefs(:,2,1) = [1;0;0;1];
- coefs(:,1,2) = [0;1;0;1];
- coefs(:,2,2) = [1;1;1;2];
- srf = nrbmak (coefs, knots);
- [v, u] = meshgrid (linspace (0, 1, 3));
- uv = [u(:)'; v(:)'];
- skl = nrbsurfderiveval (srf, uv, 2);
- w = @(x, y) x.*y + 1;
- F1 = @(x, y) x ./ w(x,y);
- F2 = @(x, y) y ./ w(x,y);
- F3 = @(x, y) x .* y ./ w(x,y);
- dF1dx = @(x, y) 1./w(x,y) - x.*y./w(x,y).^2;
- dF1dy = @(x, y)  - x.^2./w(x,y).^2;
- dF2dx = @(x, y)  - y.^2./w(x,y).^2;
- dF2dy = @(x, y) 1./w(x,y) - x.*y./w(x,y).^2;
- dF3dx = @(x, y) y./w(x,y) - x.*(y./w(x,y)).^2;
- dF3dy = @(x, y) x./w(x,y) - y.*(x./w(x,y)).^2;
- d2F1dx2  = @(x, y) -2*y./w(x,y).^2 + 2*x.*y.^2./w(x,y).^3;
- d2F1dy2  = @(x, y) 2*x.^3./w(x,y).^3;
- d2F1dxdy = @(x, y) -x./w(x,y).^2 - x./w(x,y).^2 + 2*x.^2.*y./w(x,y).^3;
- d2F2dx2  = @(x, y) 2*y.^3./w(x,y).^3;
- d2F2dy2  = @(x, y) -2*x./w(x,y).^2 + 2*y.*x.^2./w(x,y).^3;
- d2F2dxdy = @(x, y) -y./w(x,y).^2 - y./w(x,y).^2 + 2*y.^2.*x./w(x,y).^3;
- d2F3dx2  = @(x, y) -2*y.^2./w(x,y).^2 + 2*x.*y.^3./w(x,y).^3;
- d2F3dy2  = @(x, y) -2*x.^2./w(x,y).^2 + 2*y.*x.^3./w(x,y).^3;
- d2F3dxdy = @(x, y) 1./w(x,y) - 3*x.*y./w(x,y).^2 + 2*(x.*y).^2./w(x,y).^3;
- assert ([F1(u(:),v(:)), F2(u(:),v(:)), F3(u(:),v(:))], squeeze(skl(1:3,1,1,:))', 1e2*eps);
- assert ([dF1dx(u(:),v(:)), dF2dx(u(:),v(:)), dF3dx(u(:),v(:))], squeeze(skl(1:3,2,1,:))', 1e2*eps);
- assert ([dF1dy(u(:),v(:)), dF2dy(u(:),v(:)), dF3dy(u(:),v(:))], squeeze(skl(1:3,1,2,:))', 1e2*eps);
- assert ([d2F1dx2(u(:),v(:)), d2F2dx2(u(:),v(:)), d2F3dx2(u(:),v(:))], squeeze(skl(1:3,3,1,:))', 1e2*eps);
- assert ([d2F1dy2(u(:),v(:)), d2F2dy2(u(:),v(:)), d2F3dy2(u(:),v(:))], squeeze(skl(1:3,1,3,:))', 1e2*eps);
- assert ([d2F1dxdy(u(:),v(:)), d2F2dxdy(u(:),v(:)), d2F3dxdy(u(:),v(:))], squeeze(skl(1:3,2,2,:))', 1e2*eps);
-21 tests, 21 passed, 0 known failure, 0 skipped
-[inst/nrbinverse.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbinverse.m
-***** test
- nrb = nrb4surf ([0 0], [1 0], [2 3], [5 4]);
- p = nrbeval (nrb, {.25 .75});
- u = nrbinverse (nrb, p, 'Display', false);
- assert (norm (u - [.25; .75]) < 1e-8);
-
-***** test
- nrb = nrb4surf ([0 0], [1 0], [2 3], [5 4]);
- nrb = nrbdegelev (nrbextrude (nrb, [0 2 1]), [3 3 3]);
- p = nrbeval (nrb, {.25 .75 .05});
- u = nrbinverse (nrb, p, 'Display', false, 'TolX', 1e-12, 'TolFun', 1e-10);
- assert (norm (u - [.25; .75; .05]) < 1e-8);
-
-2 tests, 2 passed, 0 known failure, 0 skipped
-[inst/aveknt.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/aveknt.m
-***** test
- knt = [0 0 0 0.5 1 1 1];
- pts = aveknt (knt, 3);
- assert (pts - [0 1/4 3/4 1] < 1e-14)
-
-***** test
- knt = {[0 0 0 0.5 1 1 1] [0 0 0 0 1/3 2/3 1 1 1 1]};
- pts = aveknt (knt, [3 4]);
- assert (pts{1} - [0 1/4 3/4 1] < 1e-14);
- assert (pts{2} - [0 1/9 1/3 2/3 8/9 1] < 1e-14);
-
-***** test
- nrb = nrb4surf([0 0], [1 0], [0 1], [1 1]);
- nrb = nrbkntins (nrbdegelev (nrb, [1 2]), {[1/2] [1/3 2/3]});
- pts = aveknt (nrb);
- assert (pts{1} - [0 1/4 3/4 1] < 1e-14);
- assert (pts{2} - [0 1/9 1/3 2/3 8/9 1] < 1e-14);
-3 tests, 3 passed, 0 known failure, 0 skipped
 [inst/nrbderiv.m]
 >>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbderiv.m
 ***** demo
@@ -3851,6 +4000,277 @@
  d2Fdww(2,:,:,:) = (40.*x.*y.^2.*(x - 1).*(y - 1).*(4.*x.^2.*y.^2 - 4.*x.^2.*y - 4.*x.*y.^2 + 4.*x.*y + 5) + 40.*x.*y.^2.*z.^2.*(x - 1).*(y - 1).*(12.*x.^2.*y.^2 - 12.*x.^2.*y - 12.*x.*y.^2 + 12.*x.*y) - 40.*x.*y.^2.*z.*(x - 1).*(y - 1).*(12.*x.^2.*y.^2 - 12.*x.^2.*y - 12.*x.*y.^2 + 12.*x.*y))./(- 4.*x.^2.*y.^2.*z.^2 + 4.*x.^2.*y.^2.*z + 4.*x.^2.*y.*z.^2 - 4.*x.^2.*y.*z + 4.*x.*y.^2.*z.^2 - 4.*x.*y.^2.*z - 4.*x.*y.*z.^2 + 4.*x.*y.*z + 5).^3;
  d2Fdww(3,:,:,:) = (40.*x.*y.*(x - 1).*(y - 1).*(4.*x.^2.*y.^2.*z.^3 - 4.*x.^2.*y.*z.^3 - 4.*x.*y.^2.*z.^3 + 4.*x.*y.*z.^3 + 15.*z - 5))./(- 4.*x.^2.*y.^2.*z.^2 + 4.*x.^2.*y.^2.*z + 4.*x.^2.*y.*z.^2 - 4.*x.^2.*y.*z + 4.*x.*y.^2.*z.^2 - 4.*x.*y.^2.*z - 4.*x.*y.*z.^2 + 4.*x.*y.*z + 5).^3;
 8 tests, 8 passed, 0 known failure, 0 skipped
+[inst/nrbeval.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbeval.m
+***** demo
+ srf = nrbtestsrf;
+ p = nrbeval(srf,{linspace(0.0,1.0,20) linspace(0.0,1.0,20)});
+ h = surf(squeeze(p(1,:,:)),squeeze(p(2,:,:)),squeeze(p(3,:,:)));
+ title('Test surface.');
+ hold off
+***** test
+ knots{1} = [0 0 0 1 1 1];
+ knots{2} = [0 0 0 .5 1 1 1];
+ knots{3} = [0 0 0 0 1 1 1 1];
+ cx = [0 0.5 1]; nx = length(cx);
+ cy = [0 0.25 0.75 1]; ny = length(cy);
+ cz = [0 1/3 2/3 1]; nz = length(cz);
+ coefs(1,:,:,:) = repmat(reshape(cx,nx,1,1),[1 ny nz]);
+ coefs(2,:,:,:) = repmat(reshape(cy,1,ny,1),[nx 1 nz]);
+ coefs(3,:,:,:) = repmat(reshape(cz,1,1,nz),[nx ny 1]);
+ coefs(4,:,:,:) = 1;
+ nurbs = nrbmak(coefs, knots);
+ x = rand(5,1); y = rand(5,1); z = rand(5,1);
+ tt = [x y z]';
+ points = nrbeval(nurbs,tt);
+
+ assert(points,tt,1e-10)
+
+***** test
+ knots{1} = [0 0 0 1 1 1];
+ knots{2} = [0 0 0 0 1 1 1 1];
+ knots{3} = [0 0 1 1];
+ cx = [0 0 1]; nx = length(cx);
+ cy = [0 0 0 1]; ny = length(cy);
+ cz = [0 1]; nz = length(cz);
+ coefs(1,:,:,:) = repmat(reshape(cx,nx,1,1),[1 ny nz]);
+ coefs(2,:,:,:) = repmat(reshape(cy,1,ny,1),[nx 1 nz]);
+ coefs(3,:,:,:) = repmat(reshape(cz,1,1,nz),[nx ny 1]);
+ coefs(4,:,:,:) = 1;
+ nurbs = nrbmak(coefs, knots);
+ x = rand(5,1); y = rand(5,1); z = rand(5,1);
+ tt = [x y z]';
+ points = nrbeval(nurbs,tt);
+ assert(points,[x.^2 y.^3 z]',1e-10);
+
+***** test
+ knots{1} = [0 0 0 1 1 1];
+ knots{2} = [0 0 0 0 1 1 1 1];
+ knots{3} = [0 0 1 1];
+ cx = [0 0 1]; nx = length(cx);
+ cy = [0 0 0 1]; ny = length(cy);
+ cz = [0 1]; nz = length(cz);
+ coefs(1,:,:,:) = repmat(reshape(cx,nx,1,1),[1 ny nz]);
+ coefs(2,:,:,:) = repmat(reshape(cy,1,ny,1),[nx 1 nz]);
+ coefs(3,:,:,:) = repmat(reshape(cz,1,1,nz),[nx ny 1]);
+ coefs(4,:,:,:) = 1;
+ coefs = coefs([2 1 3 4],:,:,:);
+ nurbs = nrbmak(coefs, knots);
+ x = rand(5,1); y = rand(5,1); z = rand(5,1);
+ tt = [x y z]';
+ points = nrbeval(nurbs,tt);
+ [y.^3 x.^2 z]';
+ assert(points,[y.^3 x.^2 z]',1e-10);
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/nrbsurfderiveval.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbsurfderiveval.m
+***** test
+ k = [0 0 1 1];
+ c = [0 1];
+ [coef(2,:,:), coef(1,:,:)] = meshgrid (c, c);
+ coef(3,:,:) = coef(1,:,:);
+ srf = nrbmak (coef, {k, k});
+ [u, v] = meshgrid (linspace(0,1,11));
+ uv = [u(:)';v(:)'];
+ skl = nrbsurfderiveval (srf, uv, 0);
+ aux = nrbeval(srf,uv);
+ assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
+***** test
+ k = [0 0 1 1];
+ c = [0 1];
+ [coef(2,:,:), coef(1,:,:)] = meshgrid (c, c);
+ coef(3,:,:) = coef(1,:,:);
+ srf = nrbmak (coef, {k, k});
+ srf = nrbkntins (srf, {[], rand(2,1)});
+ [u, v] = meshgrid (linspace(0,1,11));
+ uv = [u(:)';v(:)'];
+ skl = nrbsurfderiveval (srf, uv, 0);
+ aux = nrbeval(srf,uv);
+ assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
+***** shared srf, uv
+***** test 
+ k = [0 0 0 1 1 1];
+ c = [0 1/2 1];
+ [coef(1,:,:), coef(2,:,:)] = meshgrid (c, c);
+ coef(3,:,:) = coef(1,:,:);
+ srf = nrbmak (coef, {k, k});
+ ders= nrbderiv (srf);
+ [u, v] = meshgrid (linspace(0,1,11));
+ uv = [u(:)';v(:)'];
+ skl = nrbsurfderiveval (srf, uv, 1);
+ [fun, der] = nrbdeval (srf, ders, uv);
+ assert (squeeze (skl (1:2,1,1,:)), fun(1:2,:), 1e3*eps)
+ assert (squeeze (skl (1:2,2,1,:)), der{1}(1:2,:), 1e3*eps)
+ assert (squeeze (skl (1:2,1,2,:)), der{2}(1:2,:), 1e3*eps)
+
+***** test 
+ srf = nrbdegelev (srf, [3, 1]);
+ ders= nrbderiv (srf);
+ [fun, der] = nrbdeval (srf, ders, uv);
+ skl = nrbsurfderiveval (srf, uv, 1);
+ assert (squeeze (skl (1:2,1,1,:)), fun(1:2,:), 1e3*eps)
+ assert (squeeze (skl (1:2,2,1,:)), der{1}(1:2,:), 1e3*eps)
+ assert (squeeze (skl (1:2,1,2,:)), der{2}(1:2,:), 1e3*eps)
+***** shared uv
+***** test 
+ k = [0 0 0 1 1 1];
+ c = [0 1/2 1];
+ [coef(2,:,:), coef(1,:,:)] = meshgrid (c, c);
+ coef(3,:,:) = coef(1,:,:);
+ srf = nrbmak (coef, {k, k});
+ ders= nrbderiv (srf);
+ [u, v] = meshgrid (linspace(0,1,11));
+ uv = [u(:)';v(:)'];
+ skl = nrbsurfderiveval (srf, uv, 1);
+ [fun, der] = nrbdeval (srf, ders, uv);
+ assert (squeeze (skl (1:2,1,1,:)), fun(1:2,:), 1e3*eps)
+ assert (squeeze (skl (1:2,2,1,:)), der{1}(1:2,:), 1e3*eps)
+ assert (squeeze (skl (1:2,1,2,:)), der{2}(1:2,:), 1e3*eps)
+
+***** test 
+ p = 3;  q = 3;
+ mcp = 5; ncp = 5;
+ Lx  = 10*rand(1); Ly  = Lx;
+ srf = nrbdegelev (nrb4surf ([0 0], [Lx, 0], [0 Ly], [Lx Ly]), [p-1, q-1]);
+ %%srf = nrbkntins (srf, {linspace(0,1,mcp-p+2)(2:end-1), linspace(0,1,ncp-q+2)(2:end-1)});
+ %%srf.coefs = permute (srf.coefs, [1 3 2]);
+ ders= nrbderiv (srf);
+ [fun, der] = nrbdeval (srf, ders, uv);
+ skl = nrbsurfderiveval (srf, uv, 1);
+ assert (squeeze (skl (1:2,1,1,:)), fun(1:2,:), 1e3*eps)
+ assert (squeeze (skl (1:2,2,1,:)), der{1}(1:2,:), 1e3*eps)
+ assert (squeeze (skl (1:2,1,2,:)), der{2}(1:2,:), 1e3*eps)
+***** shared srf, uv, P, dPdx, d2Pdx2, c1, c2
+***** test
+ [u, v] = meshgrid (linspace(0,1,10));
+ uv = [u(:)';v(:)'];
+ c1 = nrbmak([0 1/2 1; 0 1 0],[0 0 0 1 1 1]);
+ c1 = nrbtform (c1, vecrotx (pi/2));
+ c2  = nrbtform(c1, vectrans([0 1 0]));
+ srf = nrbdegelev (nrbruled (c1, c2), [3, 1]);
+ skl = nrbsurfderiveval (srf, uv, 2);
+ P = squeeze(skl(:,1,1,:));
+ dPdx = squeeze(skl(:,2,1,:));
+ d2Pdx2 = squeeze(skl(:,3,1,:));
+***** assert(P(3,:), 2*(P(1,:)-P(1,:).^2),100*eps)
+***** assert(dPdx(3,:), 2-4*P(1,:), 100*eps)
+***** assert(d2Pdx2(3,:), -4+0*P(1,:), 100*eps)
+ srf = nrbdegelev (nrbruled (c1, c2), [5, 6]);
+ skl = nrbsurfderiveval (srf, uv, 2);
+ P = squeeze(skl(:,1,1,:));
+ dPdx = squeeze(skl(:,2,1,:));
+ d2Pdx2 = squeeze(skl(:,3,1,:));
+ aux = nrbeval(srf,uv);
+ assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
+***** assert(P(3,:), 2*(P(1,:)-P(1,:).^2),100*eps)
+***** assert(dPdx(3,:), 2-4*P(1,:), 100*eps)
+***** assert(d2Pdx2(3,:), -4+0*P(1,:), 100*eps)
+
+***** test
+ skl = nrbsurfderiveval (srf, uv, 0);
+ aux = nrbeval (srf, uv);
+ assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
+***** shared dPdu, d2Pdu2, P, srf, uv
+***** test
+ [u, v] = meshgrid (linspace(0,1,10));
+ uv = [u(:)';v(:)'];
+ c1 = nrbmak([0 1/2 1; 0.1 1.6 1.1; 0 0 0],[0 0 0 1 1 1]);
+ c2 = nrbmak([0 1/2 1; 0.1 1.6 1.1; 1 1 1],[0 0 0 1 1 1]);
+ srf = nrbdegelev (nrbruled (c1, c2), [0, 1]);
+ skl = nrbsurfderiveval (srf, uv, 2);
+ P = squeeze(skl(:,1,1,:));
+ dPdu = squeeze(skl(:,2,1,:));
+ dPdv = squeeze(skl(:,1,2,:));
+ d2Pdu2 = squeeze(skl(:,3,1,:));
+ aux = nrbeval(srf,uv);
+ assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
+***** assert(dPdu(2,:), 3-4*P(1,:),100*eps)
+***** assert(d2Pdu2(2,:), -4+0*P(1,:),100*eps)
+
+***** test
+ skl = nrbsurfderiveval (srf, uv, 0);
+ aux = nrbeval(srf,uv);
+ assert (squeeze (skl (1:2,1,1,:)), aux(1:2,:), 1e3*eps)
+***** test
+ srf = nrb4surf([0 0], [1 0], [0 1], [1 1]);
+ geo = nrbdegelev (srf, [3 3]);
+ geo = nrbkntins (geo, {[.1:.1:.9], [.2:.2:.8]});
+ [u, v] = meshgrid (linspace(0,1,10));
+ uv = [u(:)';v(:)'];
+ skl = nrbsurfderiveval (geo, uv, 2);
+ dgeo = nrbderiv (geo);
+ [pnts, ders] = nrbdeval (geo, dgeo, uv);
+ assert (ders{1}, squeeze(skl(:,2,1,:)), 1e-9)
+ assert (ders{2}, squeeze(skl(:,1,2,:)), 1e-9)
+***** test
+ crv = nrbline ([1 0], [2 0]);
+ srf = nrbrevolve (crv, [0 0 0], [0 0 1], pi/2);
+ srf = nrbtransp (srf);
+ [v, u] = meshgrid (linspace (0, 1, 11));
+ uv = [u(:)'; v(:)'];
+ skl = nrbsurfderiveval (srf, uv, 2);
+ c = sqrt(2);
+ w      = @(x, y) (2 - c)*y.^2 + (c-2)*y + 1;
+ dwdy   = @(x, y) 2*(2-c)*y + c - 2;
+ d2wdy2 = @(x, y) 2*(2-c);
+ F1 = @(x, y) (x+1) .* ((1-y).^2 + c*y.*(1-y)) ./ w(x,y);
+ F2 = @(x, y) (x+1) .* (y.^2 + c*y.*(1-y)) ./ w(x,y);
+ dF1dx = @(x, y) ((1-y).^2 + c*y.*(1-y)) ./ w(x,y);
+ dF2dx = @(x, y) (y.^2 + c*y.*(1-y)) ./ w(x,y);
+ dF1dy = @(x, y) (x+1) .* ((2 - 2*c)*y + c - 2) ./ w(x,y) - (x+1) .* ((1-y).^2 + c*y.*(1-y)) .* dwdy(x,y) ./ w(x,y).^2;
+ dF2dy = @(x, y) (x+1) .* ((2 - 2*c)*y + c) ./ w(x,y) - (x+1) .* (y.^2 + c*y.*(1-y)) .* dwdy(x,y) ./ w(x,y).^2;
+ d2F1dx2 = @(x, y) zeros (size (x));
+ d2F2dx2 = @(x, y) zeros (size (x));
+ d2F1dxdy = @(x, y) ((2 - 2*c)*y + c - 2) ./ w(x,y) - ((1-y).^2 + c*y.*(1-y)) .* dwdy(x,y) ./ w(x,y).^2;
+ d2F2dxdy = @(x, y) ((2 - 2*c)*y + c) ./ w(x,y) - (y.^2 + c*y.*(1-y)) .* dwdy(x,y) ./ w(x,y).^2;
+ d2F1dy2  = @(x, y) (x+1)*(2 - 2*c) ./ w(x,y) - 2*(x+1) .* ((2 - 2*c)*y + c - 2) .* dwdy(x,y) ./ w(x,y).^2 - ...
+                    (x+1) .* ((1-y).^2 + c*y.*(1-y)) * d2wdy2(x,y) ./ w(x,y).^2 + ...
+                    2 * (x+1) .* ((1-y).^2 + c*y.*(1-y)) .* w(x,y) .*dwdy(x,y).^2 ./ w(x,y).^4;
+ d2F2dy2  = @(x, y) (x+1)*(2 - 2*c) ./ w(x,y) - 2*(x+1) .* ((2 - 2*c)*y + c) .* dwdy(x,y) ./ w(x,y).^2 - ...
+                    (x+1) .* (y.^2 + c*y.*(1-y)) * d2wdy2(x,y) ./ w(x,y).^2 + ...
+                    2 * (x+1) .* (y.^2 + c*y.*(1-y)) .* w(x,y) .*dwdy(x,y).^2 ./ w(x,y).^4;
+ assert ([F1(u(:),v(:)), F2(u(:),v(:))], squeeze(skl(1:2,1,1,:))', 1e2*eps);
+ assert ([dF1dx(u(:),v(:)), dF2dx(u(:),v(:))], squeeze(skl(1:2,2,1,:))', 1e2*eps);
+ assert ([dF1dy(u(:),v(:)), dF2dy(u(:),v(:))], squeeze(skl(1:2,1,2,:))', 1e2*eps);
+ assert ([d2F1dx2(u(:),v(:)), d2F2dx2(u(:),v(:))], squeeze(skl(1:2,3,1,:))', 1e2*eps);
+ assert ([d2F1dxdy(u(:),v(:)), d2F2dxdy(u(:),v(:))], squeeze(skl(1:2,2,2,:))', 1e2*eps);
+ assert ([d2F1dy2(u(:),v(:)), d2F2dy2(u(:),v(:))], squeeze(skl(1:2,1,3,:))', 1e2*eps);
+***** test
+ knots = {[0 0 1 1] [0 0 1 1]};
+ coefs(:,1,1) = [0;0;0;1];
+ coefs(:,2,1) = [1;0;0;1];
+ coefs(:,1,2) = [0;1;0;1];
+ coefs(:,2,2) = [1;1;1;2];
+ srf = nrbmak (coefs, knots);
+ [v, u] = meshgrid (linspace (0, 1, 3));
+ uv = [u(:)'; v(:)'];
+ skl = nrbsurfderiveval (srf, uv, 2);
+ w = @(x, y) x.*y + 1;
+ F1 = @(x, y) x ./ w(x,y);
+ F2 = @(x, y) y ./ w(x,y);
+ F3 = @(x, y) x .* y ./ w(x,y);
+ dF1dx = @(x, y) 1./w(x,y) - x.*y./w(x,y).^2;
+ dF1dy = @(x, y)  - x.^2./w(x,y).^2;
+ dF2dx = @(x, y)  - y.^2./w(x,y).^2;
+ dF2dy = @(x, y) 1./w(x,y) - x.*y./w(x,y).^2;
+ dF3dx = @(x, y) y./w(x,y) - x.*(y./w(x,y)).^2;
+ dF3dy = @(x, y) x./w(x,y) - y.*(x./w(x,y)).^2;
+ d2F1dx2  = @(x, y) -2*y./w(x,y).^2 + 2*x.*y.^2./w(x,y).^3;
+ d2F1dy2  = @(x, y) 2*x.^3./w(x,y).^3;
+ d2F1dxdy = @(x, y) -x./w(x,y).^2 - x./w(x,y).^2 + 2*x.^2.*y./w(x,y).^3;
+ d2F2dx2  = @(x, y) 2*y.^3./w(x,y).^3;
+ d2F2dy2  = @(x, y) -2*x./w(x,y).^2 + 2*y.*x.^2./w(x,y).^3;
+ d2F2dxdy = @(x, y) -y./w(x,y).^2 - y./w(x,y).^2 + 2*y.^2.*x./w(x,y).^3;
+ d2F3dx2  = @(x, y) -2*y.^2./w(x,y).^2 + 2*x.*y.^3./w(x,y).^3;
+ d2F3dy2  = @(x, y) -2*x.^2./w(x,y).^2 + 2*y.*x.^3./w(x,y).^3;
+ d2F3dxdy = @(x, y) 1./w(x,y) - 3*x.*y./w(x,y).^2 + 2*(x.*y).^2./w(x,y).^3;
+ assert ([F1(u(:),v(:)), F2(u(:),v(:)), F3(u(:),v(:))], squeeze(skl(1:3,1,1,:))', 1e2*eps);
+ assert ([dF1dx(u(:),v(:)), dF2dx(u(:),v(:)), dF3dx(u(:),v(:))], squeeze(skl(1:3,2,1,:))', 1e2*eps);
+ assert ([dF1dy(u(:),v(:)), dF2dy(u(:),v(:)), dF3dy(u(:),v(:))], squeeze(skl(1:3,1,2,:))', 1e2*eps);
+ assert ([d2F1dx2(u(:),v(:)), d2F2dx2(u(:),v(:)), d2F3dx2(u(:),v(:))], squeeze(skl(1:3,3,1,:))', 1e2*eps);
+ assert ([d2F1dy2(u(:),v(:)), d2F2dy2(u(:),v(:)), d2F3dy2(u(:),v(:))], squeeze(skl(1:3,1,3,:))', 1e2*eps);
+ assert ([d2F1dxdy(u(:),v(:)), d2F2dxdy(u(:),v(:)), d2F3dxdy(u(:),v(:))], squeeze(skl(1:3,2,2,:))', 1e2*eps);
+21 tests, 21 passed, 0 known failure, 0 skipped
 [inst/nrbunclamp.m]
 >>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbunclamp.m
 ***** demo
@@ -3875,435 +4295,64 @@
  ucrv = nrbunclamp (crv, 3);
  assert (F, nrbeval(ucrv, x), 1e-14);
 1 test, 1 passed, 0 known failure, 0 skipped
-[inst/surfderivcpts.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/surfderivcpts.m
-***** test
- coefs = cat(3,[0 0; 0 1],[1 1; 0 1]);
- knots = {[0 0 1 1]  [0 0 1 1]};
- plane = nrbmak(coefs,knots);
- pkl = surfderivcpts (plane.number(1)-1, plane.order(1)-1,...
-                       plane.knots{1}, plane.number(2)-1,...
-                       plane.order(2)-1, plane.knots{2}, ...
-                       squeeze (plane.coefs(1,:,:)), 1);
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/nrbnumbasisfun.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbnumbasisfun.m
-***** test
- p = 2;   q = 3;   m = 4; n = 5;
- Lx  = 1; Ly  = 1; 
- nrb = nrb4surf   ([0 0], [1 0], [0 1], [1 1]);
- nrb = nrbdegelev (nrb, [p-1, q-1]);
- ikx = linspace(0,1,m); iky = linspace(0,1,n);
- nrb = nrbkntins  (nrb, {ikx(2:end-1), iky(2:end-1)});
- nrb.coefs (4,:,:) = nrb.coefs (4,:,:) + rand (size (nrb.coefs (4,:,:)));
- u = rand (1, 30); v = rand (1, 10);
- u = (u-min (u))/max (u-min (u));
- v = (v-min (v))/max (v-min (v));
- N = nrbnumbasisfun ({u, v}, nrb);
- assert (all (all (N>0)), true)
- assert (all (all (N <= prod (nrb.number))), true)
- assert (max (max (N)), prod (nrb.number))
- assert (min (min (N)), 1)
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/surfderiveval.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/surfderiveval.m
-***** shared srf
-***** test
- k = [0 0 0 1 1 1];
- c = [0 1/2 1];
- [coef(2,:,:), coef(1,:,:)] = meshgrid (c, c);
- srf = nrbmak (coef, {k, k});
- skl = surfderiveval (srf.number(1)-1, ...
-                      srf.order(1)-1, ...
-                      srf.knots{1}, ...
-                      srf.number(2)-1, ...
-                      srf.order(2)-1, ...
-                      srf.knots{2},...
-                      squeeze(srf.coefs(1,:,:)), .5, .5, 1) ;
- assert (skl, [.5 0; 1 0])
-***** test
- srf = nrbkntins (srf, {[], rand(1,2)});
- skl = surfderiveval (srf.number(1)-1,... 
-                      srf.order(1)-1, ...
-                      srf.knots{1},...
-                      srf.number(2)-1,... 
-                      srf.order(2)-1, ...
-                      srf.knots{2},...
-                      squeeze(srf.coefs(1,:,:)), .5, .5, 1) ;
- assert (skl, [.5 0; 1 0], 100*eps)
-2 tests, 2 passed, 0 known failure, 0 skipped
-[inst/tbasisfun.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/tbasisfun.m
-***** demo
- U = {[0 0 1/2 1 1], [0 0 0 1 1]};
- p = [3, 3];
- [X, Y] = meshgrid (linspace(0, 1, 30));
- u = [X(:), Y(:)]';
- N = tbasisfun (u, p, U);
- surf (X, Y, reshape (N, size(X)))
- title('Basis function associated to a local knot vector')
- hold off
-***** test
- U = [0 1/2 1];
- p = 1;
- u = [0.3 0.4 0.6 0.7];
- [N, Nder] = tbasisfun (u, p, U);
- assert (N, [0.6 0.8 0.8 0.6], 1e-12);
- assert (Nder, [2 2 -2 -2], 1e-12);
-***** test
- U = {[0 1/2 1] [0 1/2 1]};
- p = [1 1];
- u = [0.3 0.4 0.6 0.7; 0.3 0.4 0.6 0.7];
- [N, Nder] = tbasisfun (u, p, U);
- assert (N, [0.36 0.64 0.64 0.36], 1e-12);
- assert (Nder, [1.2 1.6 -1.6 -1.2; 1.2 1.6 -1.6 -1.2], 1e-12);
-***** test
- U = {[0 1/2 1] [0 1/2 1] [0 1/2 1]};
- p = [1 1 1];
- u = [0.4 0.4 0.6 0.6; 0.4 0.4 0.6 0.6; 0.4 0.6 0.4 0.6];
- [N, Nder] = tbasisfun (u, p, U);
- assert (N, [0.512 0.512 0.512 0.512], 1e-12);
- assert (Nder, [1.28 1.28 -1.28 -1.28; 1.28 1.28 -1.28 -1.28; 1.28 -1.28 1.28 -1.28], 1e-12);
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/basisfun.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/basisfun.m
-***** test
-  n = 3; 
-  U = [0 0 0 1/2 1 1 1]; 
-  p = 2; 
-  u = linspace (0, 1, 10);  
-  s = findspan (n, p, u, U);  
-  Bref = [1.00000   0.00000   0.00000
-          0.60494   0.37037   0.02469
-          0.30864   0.59259   0.09877
-          0.11111   0.66667   0.22222
-          0.01235   0.59259   0.39506
-          0.39506   0.59259   0.01235
-          0.22222   0.66667   0.11111
-          0.09877   0.59259   0.30864
-          0.02469   0.37037   0.60494
-          0.00000   0.00000   1.00000];
-  B = basisfun (s, u, p, U);
-  assert (B, Bref, 1e-5);
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/findspan.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/findspan.m
-***** test
-  n = 3; 
-  U = [0 0 0 1/2 1 1 1]; 
-  p = 2; 
-  u = linspace(0, 1, 10);  
-  s = findspan (n, p, u, U);
-  assert (s, [2*ones(1, 5) 3*ones(1, 5)]);
-***** test
- p = 2; m = 7; n = m - p - 1;
- U = [zeros(1,p)  linspace(0,1,m+1-2*p) ones(1,p)];
- u = [ 0   0.11880   0.55118   0.93141   0.40068   0.35492 0.44392   0.88360   0.35414   0.92186   0.83085   1];
- s = [2   2   3   4   3   3   3   4   3   4   4   4];
- assert (findspan (n, p, u, U), s, 1e-10);
-2 tests, 2 passed, 0 known failure, 0 skipped
-[inst/kntrefine.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/kntrefine.m
-***** shared nrbs
-***** test
- knots = {[0 0 1 1] [0 0 0 1 1 1]};
- coefs(1,:,:) = [1 sqrt(2)/2 0; 2 sqrt(2) 0];
- coefs(2,:,:) = [0 sqrt(2)/2 1; 0 sqrt(2) 2];
- coefs(4,:,:) = [1 sqrt(2)/2 1; 1 sqrt(2)/2 1];
- nrbs = nrbmak (coefs, knots);
- nrbs = nrbkntins (nrbs, {[] [0.5 0.6 0.6]});
- nrbs = nrbdegelev (nrbs, [0 1]);
- nrbs = nrbkntins (nrbs, {[] [0.4]});
- rknots = kntrefine (nrbs.knots, [1 1], [1 1], [0 0]);
- assert (rknots{1} == [0 0 0.5 1 1]);
- assert (rknots{2} == [0 0 0.2 0.4 0.45 0.5 0.55 0.6 0.8 1 1]);
-
-***** test
- rknots = kntrefine (nrbs.knots, [1 1], [3 3], [0 0]);
- assert (rknots{1}, [0 0 0 0 0.5 0.5 0.5 1 1 1 1]);
- assert (rknots{2}, [0 0 0 0 0.2 0.2 0.2 0.4 0.4 0.4 0.45 0.45 0.45 0.5 0.5 0.5 0.55 0.55 0.55 0.6 0.6 0.6 0.8 0.8 0.8 1 1 1 1]);
-
-***** test
- rknots = kntrefine (nrbs.knots, [1 1], [3 3], [2 2]);
- assert (rknots{1}, [0 0 0 0 0.5 1 1 1 1]);
- assert (rknots{2}, [0 0 0 0 0.2 0.4 0.45 0.5 0.5 0.55 0.6 0.6 0.6 0.8 1 1 1 1]);
-
-***** test
- rknots = kntrefine (nrbs.knots, [1 1], [4 4], [0 0]);
- assert (rknots{1}, [0 0 0 0 0 0.5 0.5 0.5 0.5 1 1 1 1 1]);
- assert (rknots{2}, [0 0 0 0 0 0.2 0.2 0.2 0.2 0.4 0.4 0.4 0.4 0.45 0.45 0.45 0.45 0.5 0.5 0.5 0.5 0.55 0.55 0.55 0.55 0.6 0.6 0.6 0.6 0.8 0.8 0.8 0.8 1 1 1 1 1]);
-
-***** test
- rknots = kntrefine (nrbs.knots, [1 1], [4 4], [3 3]);
- assert (rknots{1}, [0 0 0 0 0 0.5 1 1 1 1 1]);
- assert (rknots{2}, [0 0 0 0 0 0.2 0.4 0.4 0.45 0.5 0.5 0.5 0.55 0.6 0.6 0.6 0.6 0.8 1 1 1 1 1]);
-
-***** test
- knots = [0 0 0 0 0.4 0.5 0.5 0.6 0.6 0.6 1 1 1 1];
- rknots = kntrefine (knots, 1, 4, 3);
- assert (rknots, [0 0 0 0 0 0.2 0.4 0.4 0.45 0.5 0.5 0.5 0.55 0.6 0.6 0.6 0.6 0.8 1 1 1 1 1]);
-6 tests, 6 passed, 0 known failure, 0 skipped
-[inst/curvederiveval.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/curvederiveval.m
-***** test
- k = [0 0 0 1 1 1];
- coefs(:,1) = [0;0;0;1];
- coefs(:,2) = [1;0;1;1];
- coefs(:,3) = [1;1;1;1];
- crv = nrbmak (coefs, k);
- ck = curvederiveval (crv.number-1, crv.order-1, crv.knots, squeeze (crv.coefs(1,:,:)), 0.5, 2);
- assert(ck, [0.75; 1; -2]);
- ck = curvederiveval (crv.number-1, crv.order-1, crv.knots, squeeze (crv.coefs(2,:,:)), 0.5, 2);
- assert(ck, [0.25; 1; 2]);
- ck = curvederiveval (crv.number-1, crv.order-1, crv.knots, squeeze (crv.coefs(3,:,:)), 0.5, 2);
- assert(ck, [0.75; 1; -2]);
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/nrbmeasure.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbmeasure.m
-***** test 
- c = nrbcirc (1, [0 0], 0,  pi/3);
- l = nrbmeasure(c, 0, 1, 1e-7);
- assert (l, pi/3, 1e-7)
-***** test 
- c = nrbcirc (1, [0 0], 0,  pi/2);
- s = zeros (1, 100); e = linspace (0, 1, 100);
- for ii = 1:100
-   l(ii) = nrbmeasure (c, s(ii), e(ii), 1e-7);
- endfor
- xx = nrbeval (c, e);
- theta = atan2 (xx(2,:), xx(1,:));
- assert (l, theta, 1e-7)
-***** test 
- c = nrbcirc (1, [0 0], 0,  pi/2);
- s = 0; e = linspace (0, 1, 100);
- for ii = 1:100
-   l(ii) = nrbmeasure (c, s, e(ii), 1e-7);
- endfor
- l2 = nrbmeasure (c, s, e, 1e-7);
- assert (l, l2, eps)
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/numbasisfun.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/numbasisfun.m
-***** test
-  n = 3; 
-  U = [0 0 0 1/2 1 1 1]; 
-  p = 2; 
-  u = linspace (0, 1, 10);  
-  s = findspan (n, p, u, U); 
-  Bref = [0   0   0   0   0   1   1   1   1   1; ...
-          1   1   1   1   1   2   2   2   2   2; ...
-          2   2   2   2   2   3   3   3   3   3].';
-  B = numbasisfun (s, u, p, U);
-  assert (B, Bref)
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/nrbclamp.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbclamp.m
-***** test
- crv = nrbdegelev (nrbcirc (1, [], 0, pi/2), 2);
- crv = nrbunclamp (crv, 3);
- xx = linspace (0, 1, 20);
- crv1 = nrbclamp (crv);
- assert (crv1.knots, [0 0 0 0 0 1 1 1 1 1])
- assert (nrbeval(crv, xx), nrbeval(crv1, xx), 1e-14)
- crv1 = nrbclamp (crv, 2);
- assert (crv1.knots, [-3 -2 -1 0 0 1 1 2 3 4])
- assert (nrbeval(crv, xx), nrbeval(crv1, xx), 1e-14)
-***** test
- crv1 = nrbcirc(1,[],0,pi/4);
- crv2 = nrbcirc(2,[],0,pi/4);
- srf = nrbkntins (nrbdegelev (nrbruled(crv1, crv2), [3 2]), {0.25 []});
- srf = nrbunclamp (srf, [4 2]);
- srf1 = nrbclamp (srf);
- xx = linspace(0,1,20);
- assert(srf1.knots, {[0 0 0 0 0 0 0.2500 1 1 1 1 1 1]  [0 0 0 0 1 1 1 1]})
- assert (nrbeval(srf, {xx xx}), nrbeval(srf1, {xx xx}), 1e-14);
- srf1 = nrbclamp (srf, [3 1]);
- assert (srf1.knots, {[-2 -1.75 -1 -0.75 0 0 0.25 1 1 1.25 2 2.25 3], [-2 -1 0 0 1 1 2 3]})
- assert (nrbeval(srf, {xx xx}), nrbeval(srf1, {xx xx}), 1e-14);
- srf1 = nrbclamp (srf, [], 2);
- assert(srf1.knots, {[-2.75 -2 -1.75 -1 -0.75 0 0.25 1 1.25 2 2.25 3 3.25] [0 0 0 0 1 1 1 1]})
- assert (nrbeval(srf, {xx xx}), nrbeval(srf1, {xx xx}), 1e-14);
-2 tests, 2 passed, 0 known failure, 0 skipped
-[inst/nrbeval_der_w.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbeval_der_w.m
-***** test % 1D
- nrb = nrbkntins (nrbcirc (1, [0 0], 0, pi/2), .5);
- u = linspace (0, 1, 11);
- delta_w = .01;
- n = nrb.number;
- der_ex = zeros (3, numel (u), n);
- der_fd = zeros (3, numel (u), n);
- for iw = 1:n
-   new_w1   = nrb.coefs (4, iw) + delta_w;
-   new_w2   = nrb.coefs (4, iw) - delta_w;
-   nrb1 = nrbmodw (nrb, new_w1, iw);
-   nrb2 = nrbmodw (nrb, new_w2, iw);
-   der_ex(:,:,iw) = nrbeval_der_w (nrb, iw, u);
+[inst/nrbeval_der_p.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbeval_der_p.m
+***** test %% 1D
+ nrb =   nrbkntins (nrbcirc (1, [0 0], 0, pi/2), .5);
+ u = 0:.1:.9;
+ index = 1:nrb.number;
+ e = zeros (numel (u), numel (index), 1);
+ for jj = 1:numel (index)
+   deltap = .1 * rand (3, 1);
+   nrb2 = nrbmodp (nrb, deltap, index(jj));
+   der_ex = nrbeval_der_p (nrb, index(jj), u);
    p2 = nrbeval (nrb2, u);
-   p1 = nrbeval (nrb1, u);
-   der_fd(:,:,iw) = -(p2 - p1) ./ (2*delta_w);  
+   p1 = nrbeval (nrb, u);
+   der_fd = (p2 - p1) ./ deltap;
+   e(:,jj) = sqrt (sum ((repmat (der_ex, 3, 1) - der_fd).^2, 1));
  end
- error = max (abs (der_ex(:) - der_fd(:)));
- assert (error < 1.e-4)
-
-***** test %2D
+ assert (max(e(:)) < 1.e-8);
+ 
+***** test %% 2D
  crv = nrbline([1 0], [2 0]);
  nrb = nrbtransp (nrbrevolve (crv, [], [0 0 1], pi/2));
  new_knots = linspace (1/9, 8/9, 8);
  nrb = nrbkntins (nrb, {new_knots, new_knots});
- u = linspace (0, 1, 5);
+ u = 0:.1:.9;
  v = u;
- delta_w = .01;
- n = nrb.number(1) * nrb.number(2);
- der_ex = zeros (3, numel(u)* numel(v), n);
- der_fd = zeros (3, numel(u)* numel(v), n);
- for iw = 1:prod(nrb.number)
-   new_w1   = nrb.coefs (4, iw) + delta_w;
-   new_w2   = nrb.coefs (4, iw) - delta_w;
-   nrb1 = nrbmodw (nrb, new_w1, iw);
-   nrb2 = nrbmodw (nrb, new_w2, iw);
-   der_ex(:,:,iw) = nrbeval_der_w (nrb, iw, {u v});
+ e = zeros (nrb.number(1) * nrb.number(2), numel (u), numel (v));
+ for index = 1:prod(nrb.number)
+   deltap = .1 * rand (3, 1);
+   nrb2 = nrbmodp (nrb, deltap, index);
+   der_ex = nrbeval_der_p (nrb, index, {u v});
    p2 = nrbeval (nrb2, {u v});
-   p1 = nrbeval (nrb1, {u v});
-   der_fd(:,:,iw) = reshape (-(p2 - p1) ./ (2*delta_w), 3, []);
+   p1 = nrbeval (nrb, {u v});
+   der_fd = (p2 - p1) ./ deltap;
+   der_ex = reshape (repmat (der_ex, 3, 1), size(der_fd));
+   e(index,:,:) = sqrt (sum ((der_ex - der_fd).^2, 1));
  end
- error = max (abs (der_ex(:) - der_fd(:)));
- assert (error < 1.e-5)
+ assert (max(e(:)) < 1.e-8)
  
-***** test % 3D
+***** test %% 3D
  crv = nrbline([1 0], [2 0]);
  nrb = nrbtransp (nrbrevolve (crv, [], [0 0 1], pi/2));
  nrb = nrbextrude (nrb, [0 0 1]);
- u = 0:.33:.99;
- v = 0:.1:.9;
- w = [.25 .5 .75];
- delta_w = .01;
- n = nrb.number(1) * nrb.number(2) * nrb.number(3);
- der_ex = zeros (3, numel(u)*numel(v)*numel(w), n);
- der_fd = zeros (3, numel(u)*numel(v)*numel(w), n);
- for iw = 1:prod(nrb.number)
-   new_w1   = nrb.coefs (4, iw) + delta_w;
-   new_w2   = nrb.coefs (4, iw) - delta_w;
-   nrb1 = nrbmodw (nrb, new_w1, iw);
-   nrb2 = nrbmodw (nrb, new_w2, iw);
-   der_ex(:,:,iw) = nrbeval_der_w (nrb, iw, {u v w});
+ u = 0:.1:.9;
+ v = u;
+ w = u;
+ e = zeros (nrb.number(1) * nrb.number(2) * nrb.number(3), numel(u), numel(v), numel(w));
+ for index = 1:prod(nrb.number)
+   deltap = .1 * rand (3, 1);
+   nrb2 = nrbmodp (nrb, deltap, index);
+   der_ex = nrbeval_der_p (nrb, index, {u v w});
    p2 = nrbeval (nrb2, {u v w});
-   p1 = nrbeval (nrb1, {u v w});
-   der_fd(:,:,iw) = reshape (-(p2 - p1) ./ (2*delta_w), 3, []);
+   p1 = nrbeval (nrb, {u v w});
+   der_fd = (p2 - p1) ./ deltap;
+   der_ex = reshape (repmat (der_ex, 3, 1), size (der_fd));
+   e(index,:,:,:) = sqrt (sum ((der_ex - der_fd).^2, 1));
  end
- error = max (max (squeeze (max (abs (der_ex - der_fd)))));
- assert (error < 1.e-4)
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/nrbreverse.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbreverse.m
-***** demo
- pnts = [0.5 1.5 3.0 7.5 8.5;
-         3.0 5.5 1.5 4.0 4.5;
-         0.0 0.0 0.0 0.0 0.0];
- crv1 = nrbmak(pnts,[0 0 0 1/2 3/4 1 1 1]);
- crv2 = nrbreverse(crv1);
- fprintf('Knots of the original curve\n')
- disp(crv1.knots)
- fprintf('Knots of the reversed curve\n')
- disp(crv2.knots)
- fprintf('Control points of the original curve\n')
- disp(crv1.coefs(1:2,:))
- fprintf('Control points of the reversed curve\n')
- disp(crv2.coefs(1:2,:))
- nrbplot(crv1,100)
- hold on
- nrbplot(crv2,100)
- title('The curve and its reverse are the same')
- hold off
-***** test
- srf  = nrbrevolve(nrbline([1 0],[2 0]), [0 0 0], [0 0 1], pi/2);
- srf  = nrbkntins (srf, {0.3, 0.6});
- srf2 = nrbreverse (srf);
- assert (srf.knots, cellfun(@(x) sort(1-x), srf2.knots, 'UniformOutput', false), 1e-15)
- assert (srf.coefs, srf2.coefs(:,end:-1:1,end:-1:1))
-***** test
- srf  = nrbrevolve(nrbline([1 0],[2 0]), [0 0 0], [0 0 1], pi/2);
- srf  = nrbkntins (srf, {0.3, 0.6});
- srf2 = nrbreverse (srf, 1);
- knt{1} = sort(1-srf2.knots{1}); knt{2} = srf2.knots{2};
- assert (srf.knots, knt, 1e-15)
- assert (srf.coefs, srf2.coefs(:,end:-1:1,:))
- srf2 = nrbreverse (srf, 2);
- knt{1} = srf2.knots{1}; knt{2} = sort(1-srf2.knots{2});
- assert (srf.knots, knt, 1e-15)
- assert (srf.coefs, srf2.coefs(:,:,end:-1:1))
-2 tests, 2 passed, 0 known failure, 0 skipped
-[inst/nrbbasisfunder.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbbasisfunder.m
-***** demo
- U = [0 0 0 0 1 1 1 1];
- x = [0 1/3 2/3 1] ;
- y = [0 0 0 0];
- w = [1 1 1 1];
- nrb = nrbmak ([x;y;y;w], U);
- u = linspace(0, 1, 30);
- [Bu, id] = nrbbasisfunder (u, nrb);
- plot(u, Bu)
- title('Derivatives of the cubic Bernstein polynomials')
- hold off
-***** test
- U = [0 0 0 0 1 1 1 1];
- x = [0 1/3 2/3 1] ;
- y = [0 0 0 0];
- w = rand(1,4);
- nrb = nrbmak ([x;y;y;w], U);
- u = linspace(0, 1, 30);
- [Bu, id] = nrbbasisfunder (u, nrb);
- #plot(u, Bu)
- assert (sum(Bu, 2), zeros(numel(u), 1), 1e-10), 
-***** test
- U = [0 0 0 0 1/2 1 1 1 1];
- x = [0 1/4 1/2 3/4 1] ;
- y = [0 0 0 0 0];
- w = rand(1,5);
- nrb = nrbmak ([x;y;y;w], U);
- u = linspace(0, 1, 300); 
- [Bu, id] = nrbbasisfunder (u, nrb); 
- assert (sum(Bu, 2), zeros(numel(u), 1), 1e-10)
-***** test
- p = 2;   q = 3;   m = 4; n = 5;
- Lx  = 1; Ly  = 1; 
- nrb = nrb4surf   ([0 0], [1 0], [0 1], [1 1]);
- nrb = nrbdegelev (nrb, [p-1, q-1]);
- aux1 = linspace(0,1,m); aux2 = linspace(0,1,n);
- nrb = nrbkntins  (nrb, {aux1(2:end-1), aux2(2:end-1)});
- nrb.coefs (4,:,:) = nrb.coefs(4,:,:) + rand (size (nrb.coefs (4,:,:)));
- [Bu, Bv, N] = nrbbasisfunder ({rand(1, 20), rand(1, 20)}, nrb);
- #plot3(squeeze(u(1,:,:)), squeeze(u(2,:,:)), reshape(Bu(:,10), 20, 20),'o')
- assert (sum (Bu, 2), zeros(20^2, 1), 1e-10)
+ assert (max (e(:)) < 1.e-8);
 3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/nrbtransp.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbtransp.m
-***** demo
- srf = nrb4surf([0 0 0], [1 0 1], [0 1 1], [1 1 2]);
- nrbplot(srf,[20 5]);
- title('Plane surface and its transposed (translated)')
- hold on
- srf.coefs(3,:,:) = srf.coefs(3,:,:) + 10;
- srf = nrbtransp(srf);
- nrbplot(srf,[20 5]);
- hold off
-***** test
- srf  = nrbrevolve(nrbline([1 0],[2 0]), [0 0 0], [0 0 1], pi/2);
- srft = nrbtransp(srf);
- assert (srf.number, fliplr(srft.number));
- assert (srf.order, fliplr(srft.order));
- assert (srf.knots, fliplr(srft.knots));
- assert (srf.coefs, permute(srft.coefs, [1 3 2]));
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/curvederivcpts.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/curvederivcpts.m
-***** test
- line = nrbmak([0.0 1.5; 0.0 3.0],[0.0 0.0 1.0 1.0]);
- pk   = curvederivcpts (line.number-1, line.order-1, line.knots,...
-                        line.coefs(1,:), 2);
- assert (pk, [0 3/2; 3/2 0], 100*eps);
-1 test, 1 passed, 0 known failure, 0 skipped
 [inst/nrbdeval.m]
 >>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbdeval.m
 ***** demo
@@ -4391,6 +4440,19 @@
  assert(jac{1}(2,:,:),2*x',1e-12)
  assert(jac{3}(3,:,:),ones(size(z')),1e-12)
 2 tests, 2 passed, 0 known failure, 0 skipped
+[inst/nrbsquare.m]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/nrbsquare.m
+***** test
+ srf = nrbsquare ([], 1, 2, 2, 4);
+ assert (srf.order, [3 3]);
+ knt = [0 0 0 1/4 1/2 3/4 1 1 1];
+ assert (srf.knots, {knt knt})
+ x = linspace (0, 1, 100);
+ [X,Y] = ndgrid (x, x);
+ vals = nrbeval (srf, {x x});
+ assert (squeeze(vals(1,:,:)), X, 1e-15);
+ assert (squeeze(vals(2,:,:)), 2*Y, 1e-15);
+1 test, 1 passed, 0 known failure, 0 skipped
 [inst/basisfunder.m]
 >>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/basisfunder.m
 ***** test
@@ -4425,32 +4487,6 @@
  assert (ders(:, (p+2):end, :), zeros(numel(u), 8-p-1, p+1), 1e-13)
  assert (all(all(ders(:, 1, :) <= 1)), true)
 3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/kntbrkdegmult.m]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/inst/kntbrkdegmult.m
-***** test
- breaks = [0 1 2 3 4];
- degree = 3;
- knots = kntbrkdegmult (breaks, degree);
- assert (knots, [0 0 0 0 1 2 3 4 4 4 4])
-***** test
- breaks = [0 1 2 3 4];
- degree = 3;
- mult   = 2;
- knots = kntbrkdegmult (breaks, degree, mult);
- assert (knots, [0 0 0 0 1 1 2 2 3 3 4 4 4 4])
-***** test
- breaks = [0 1 2 3 4];
- degree = 3;
- mult   = [1 2 3];
- knots = kntbrkdegmult (breaks, degree, mult);
- assert (knots, [0 0 0 0 1 2 2 3 3 3 4 4 4 4])
-***** test
- breaks = {[0 1 2 3 4] [0 1 2 3]};
- degree = [3 2];
- mult   = {[1 2 3] 2};
- knots = kntbrkdegmult (breaks, degree, mult);
- assert (knots, {[0 0 0 0 1 2 2 3 3 3 4 4 4 4] [0 0 0 1 1 2 2 3 3 3]})
-4 tests, 4 passed, 0 known failure, 0 skipped
 Checking C++ files ...
 [src/nrbsurfderiveval.cc]
 >>>>> /build/reproducible-path/octave-nurbs-1.4.3/src/nrbsurfderiveval.cc
@@ -4707,71 +4743,6 @@
  assert ([d2F1dy2(u(:),v(:)), d2F2dy2(u(:),v(:)), d2F3dy2(u(:),v(:))], squeeze(skl(1:3,1,3,:))', 1e2*eps);
  assert ([d2F1dxdy(u(:),v(:)), d2F2dxdy(u(:),v(:)), d2F3dxdy(u(:),v(:))], squeeze(skl(1:3,2,2,:))', 1e2*eps);
 23 tests, 23 passed, 0 known failure, 0 skipped
-[src/tbasisfun.cc]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/src/tbasisfun.cc
-***** demo
- U = {[0 0 1/2 1 1], [0 0 0 1 1]};
- p = [3, 3];
- [X, Y] = meshgrid (linspace(0, 1, 30));
- u = [X(:), Y(:)]';
- N = tbasisfun (u, p, U);
- surf (X, Y, reshape (N, size(X)))
- title('Basis function associated to a local knot vector')
- hold off
-***** test
- U = [0 1/2 1];
- p = 1;
- u = [0.3 0.4 0.6 0.7];
- [N, Nder] = tbasisfun (u, p, U);
- assert (N, [0.6 0.8 0.8 0.6], 1e-12);
- assert (Nder, [2 2 -2 -2], 1e-12);
-***** test
- U = {[0 1/2 1] [0 1/2 1]};
- p = [1 1];
- u = [0.3 0.4 0.6 0.7; 0.3 0.4 0.6 0.7];
- [N, Nder] = tbasisfun (u, p, U);
- assert (N, [0.36 0.64 0.64 0.36], 1e-12);
- assert (Nder, [1.2 1.6 -1.6 -1.2; 1.2 1.6 -1.6 -1.2], 1e-12);
-***** test
- U = {[0 1/2 1] [0 1/2 1] [0 1/2 1]};
- p = [1 1 1];
- u = [0.4 0.4 0.6 0.6; 0.4 0.4 0.6 0.6; 0.4 0.6 0.4 0.6];
- [N, Nder] = tbasisfun (u, p, U);
- assert (N, [0.512 0.512 0.512 0.512], 1e-12);
- assert (Nder, [1.28 1.28 -1.28 -1.28; 1.28 1.28 -1.28 -1.28; 1.28 -1.28 1.28 -1.28], 1e-12);
-3 tests, 3 passed, 0 known failure, 0 skipped
-[src/curvederivcpts.cc]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/src/curvederivcpts.cc
-***** test
- line = nrbmak([0.0 1.5; 0.0 3.0],[0.0 0.0 1.0 1.0]);
- pk   = curvederivcpts (line.number-1, line.order-1, line.knots,
-                        line.coefs(1,:), 2);
- assert (pk, [0 3/2; 3/2 0], 100*eps);
-1 test, 1 passed, 0 known failure, 0 skipped
-[src/basisfun.cc]
->>>>> /build/reproducible-path/octave-nurbs-1.4.3/src/basisfun.cc
-***** shared n, U, p, u, s
-***** test
-  n = 3; 
-  U = [0 0 0 1/2 1 1 1]; 
-  p = 2; 
-  u = linspace(0, 1, 10);  
-  s = findspan(n, p, u, U); 
-  assert (s, [2*ones(1, 5) 3*ones(1, 5)]);
-***** test
-  Bref = [1.00000   0.00000   0.00000
-          0.60494   0.37037   0.02469
-          0.30864   0.59259   0.09877
-          0.11111   0.66667   0.22222
-          0.01235   0.59259   0.39506
-          0.39506   0.59259   0.01235
-          0.22222   0.66667   0.11111
-          0.09877   0.59259   0.30864
-          0.02469   0.37037   0.60494
-          0.00000   0.00000   1.00000];
-  B = basisfun(s, u, p, U);
-  assert (B, Bref, 1e-5);
-2 tests, 2 passed, 0 known failure, 0 skipped
 [src/surfderiveval.cc]
 >>>>> /build/reproducible-path/octave-nurbs-1.4.3/src/surfderiveval.cc
 ***** shared srf
@@ -4819,6 +4790,71 @@
  
  assert (pkl(:),pkl2);
 1 test, 1 passed, 0 known failure, 0 skipped
+[src/basisfun.cc]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/src/basisfun.cc
+***** shared n, U, p, u, s
+***** test
+  n = 3; 
+  U = [0 0 0 1/2 1 1 1]; 
+  p = 2; 
+  u = linspace(0, 1, 10);  
+  s = findspan(n, p, u, U); 
+  assert (s, [2*ones(1, 5) 3*ones(1, 5)]);
+***** test
+  Bref = [1.00000   0.00000   0.00000
+          0.60494   0.37037   0.02469
+          0.30864   0.59259   0.09877
+          0.11111   0.66667   0.22222
+          0.01235   0.59259   0.39506
+          0.39506   0.59259   0.01235
+          0.22222   0.66667   0.11111
+          0.09877   0.59259   0.30864
+          0.02469   0.37037   0.60494
+          0.00000   0.00000   1.00000];
+  B = basisfun(s, u, p, U);
+  assert (B, Bref, 1e-5);
+2 tests, 2 passed, 0 known failure, 0 skipped
+[src/tbasisfun.cc]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/src/tbasisfun.cc
+***** demo
+ U = {[0 0 1/2 1 1], [0 0 0 1 1]};
+ p = [3, 3];
+ [X, Y] = meshgrid (linspace(0, 1, 30));
+ u = [X(:), Y(:)]';
+ N = tbasisfun (u, p, U);
+ surf (X, Y, reshape (N, size(X)))
+ title('Basis function associated to a local knot vector')
+ hold off
+***** test
+ U = [0 1/2 1];
+ p = 1;
+ u = [0.3 0.4 0.6 0.7];
+ [N, Nder] = tbasisfun (u, p, U);
+ assert (N, [0.6 0.8 0.8 0.6], 1e-12);
+ assert (Nder, [2 2 -2 -2], 1e-12);
+***** test
+ U = {[0 1/2 1] [0 1/2 1]};
+ p = [1 1];
+ u = [0.3 0.4 0.6 0.7; 0.3 0.4 0.6 0.7];
+ [N, Nder] = tbasisfun (u, p, U);
+ assert (N, [0.36 0.64 0.64 0.36], 1e-12);
+ assert (Nder, [1.2 1.6 -1.6 -1.2; 1.2 1.6 -1.6 -1.2], 1e-12);
+***** test
+ U = {[0 1/2 1] [0 1/2 1] [0 1/2 1]};
+ p = [1 1 1];
+ u = [0.4 0.4 0.6 0.6; 0.4 0.4 0.6 0.6; 0.4 0.6 0.4 0.6];
+ [N, Nder] = tbasisfun (u, p, U);
+ assert (N, [0.512 0.512 0.512 0.512], 1e-12);
+ assert (Nder, [1.28 1.28 -1.28 -1.28; 1.28 1.28 -1.28 -1.28; 1.28 -1.28 1.28 -1.28], 1e-12);
+3 tests, 3 passed, 0 known failure, 0 skipped
+[src/curvederivcpts.cc]
+>>>>> /build/reproducible-path/octave-nurbs-1.4.3/src/curvederivcpts.cc
+***** test
+ line = nrbmak([0.0 1.5; 0.0 3.0],[0.0 0.0 1.0 1.0]);
+ pk   = curvederivcpts (line.number-1, line.order-1, line.knots,
+                        line.coefs(1,:), 2);
+ assert (pk, [0 3/2; 3/2 0], 100*eps);
+1 test, 1 passed, 0 known failure, 0 skipped
 Done running the unit tests.
 Summary: 127 tests, 127 passed, 0 known failures, 0 skipped
    dh_installdocs -O--buildsystem=octave
@@ -4852,12 +4888,14 @@
 dpkg-buildpackage: info: binary-only upload (no source included)
 dpkg-genchanges: info: not including original source code in upload
 I: copying local configuration
+I: user script /srv/workspace/pbuilder/20177/tmp/hooks/B01_cleanup starting
+I: user script /srv/workspace/pbuilder/20177/tmp/hooks/B01_cleanup finished
 I: unmounting dev/ptmx filesystem
 I: unmounting dev/pts filesystem
 I: unmounting dev/shm filesystem
 I: unmounting proc filesystem
 I: unmounting sys filesystem
 I: cleaning the build env 
-I: removing directory /srv/workspace/pbuilder/116699 and its subdirectories
-I: Current time: Mon May 19 15:34:17 -12 2025
-I: pbuilder-time-stamp: 1747712057
+I: removing directory /srv/workspace/pbuilder/20177 and its subdirectories
+I: Current time: Wed Apr 17 11:13:52 +14 2024
+I: pbuilder-time-stamp: 1713302032