Diff of the two buildlogs:

--
--- b1/build.log	2020-08-12 12:29:38.993756432 +0000
+++ b2/build.log	2020-08-12 12:36:12.789540194 +0000
@@ -1,6 +1,6 @@
 I: pbuilder: network access will be disabled during build
-I: Current time: Wed Aug 12 00:18:19 -12 2020
-I: pbuilder-time-stamp: 1597234699
+I: Current time: Thu Aug 13 02:29:59 +14 2020
+I: pbuilder-time-stamp: 1597235399
 I: Building the build Environment
 I: extracting base tarball [/var/cache/pbuilder/buster-reproducible-base.tgz]
 I: copying local configuration
@@ -17,7 +17,7 @@
 I: Extracting source
 gpgv: unknown type of key resource 'trustedkeys.kbx'
 gpgv: keyblock resource '/root/.gnupg/trustedkeys.kbx': General error
-gpgv: Signature made Fri Jan  4 20:40:30 2019 -12
+gpgv: Signature made Sat Jan  5 22:40:30 2019 +14
 gpgv:                using RSA key 3F464391498FE874BDB5D98F2124AA1983785C90
 gpgv: Can't check signature: No public key
 dpkg-source: warning: failed to verify signature on ./octave-nurbs_1.3.13-5.dsc
@@ -26,136 +26,170 @@
 dpkg-source: info: unpacking octave-nurbs_1.3.13-5.debian.tar.xz
 I: Not using root during the build.
 I: Installing the build-deps
-I: user script /srv/workspace/pbuilder/27930/tmp/hooks/D02_print_environment starting
+I: user script /srv/workspace/pbuilder/7978/tmp/hooks/D01_modify_environment starting
+debug: Running on ff2a.
+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
+Removing 'diversion of /bin/sh to /bin/sh.distrib by dash'
+Adding 'diversion of /bin/sh to /bin/sh.distrib by bash'
+Removing 'diversion of /usr/share/man/man1/sh.1.gz to /usr/share/man/man1/sh.distrib.1.gz by dash'
+Adding 'diversion of /usr/share/man/man1/sh.1.gz to /usr/share/man/man1/sh.distrib.1.gz by 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/7978/tmp/hooks/D01_modify_environment finished
+I: user script /srv/workspace/pbuilder/7978/tmp/hooks/D02_print_environment starting
 I: set
-  BUILDDIR='/build'
-  BUILDUSERGECOS='first user,first room,first work-phone,first home-phone,first other'
-  BUILDUSERNAME='pbuilder1'
-  BUILD_ARCH='armhf'
-  DEBIAN_FRONTEND='noninteractive'
-  DEB_BUILD_OPTIONS='buildinfo=+all reproducible=+all parallel=3'
-  DISTRIBUTION=''
-  HOME='/root'
-  HOST_ARCH='armhf'
+  BASH=/bin/sh
+  BASHOPTS=checkwinsize:cmdhist:complete_fullquote:extquote:force_fignore:globasciiranges:hostcomplete:interactive_comments:progcomp:promptvars:sourcepath
+  BASH_ALIASES=()
+  BASH_ARGC=()
+  BASH_ARGV=()
+  BASH_CMDS=()
+  BASH_LINENO=([0]="12" [1]="0")
+  BASH_SOURCE=([0]="/tmp/hooks/D02_print_environment" [1]="/tmp/hooks/D02_print_environment")
+  BASH_VERSINFO=([0]="5" [1]="0" [2]="3" [3]="1" [4]="release" [5]="arm-unknown-linux-gnueabihf")
+  BASH_VERSION='5.0.3(1)-release'
+  BUILDDIR=/build
+  BUILDUSERGECOS='second user,second room,second work-phone,second home-phone,second other'
+  BUILDUSERNAME=pbuilder2
+  BUILD_ARCH=armhf
+  DEBIAN_FRONTEND=noninteractive
+  DEB_BUILD_OPTIONS='buildinfo=+all reproducible=+all parallel=4'
+  DIRSTACK=()
+  DISTRIBUTION=
+  EUID=0
+  FUNCNAME=([0]="Echo" [1]="main")
+  GROUPS=()
+  HOME=/root
+  HOSTNAME=i-capture-the-hostname
+  HOSTTYPE=arm
+  HOST_ARCH=armhf
   IFS=' 	
   '
-  INVOCATION_ID='0624ddb18ec34dda9c376f17dbe3ad9b'
-  LANG='C'
-  LANGUAGE='en_US:en'
-  LC_ALL='C'
-  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='27930'
-  PS1='# '
-  PS2='> '
+  INVOCATION_ID=6297e4e6b7764aab9b4c93e43fc9079e
+  LANG=C
+  LANGUAGE=it_CH:it
+  LC_ALL=C
+  MACHTYPE=arm-unknown-linux-gnueabihf
+  MAIL=/var/mail/root
+  OPTERR=1
+  OPTIND=1
+  OSTYPE=linux-gnueabihf
+  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=7978
   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/tmp.wUclLRGWgJ/pbuilderrc_WPdh --hookdir /etc/pbuilder/first-build-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/buster-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/tmp.wUclLRGWgJ/b1 --logfile b1/build.log octave-nurbs_1.3.13-5.dsc'
-  SUDO_GID='115'
-  SUDO_UID='110'
-  SUDO_USER='jenkins'
-  TERM='unknown'
-  TZ='/usr/share/zoneinfo/Etc/GMT+12'
-  USER='root'
-  _='/usr/bin/systemd-run'
-  http_proxy='http://10.0.0.15:8000/'
+  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/tmp.wUclLRGWgJ/pbuilderrc_TzgO --hookdir /etc/pbuilder/rebuild-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/buster-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/tmp.wUclLRGWgJ/b2 --logfile b2/build.log octave-nurbs_1.3.13-5.dsc'
+  SUDO_GID=114
+  SUDO_UID=109
+  SUDO_USER=jenkins
+  TERM=unknown
+  TZ=/usr/share/zoneinfo/Etc/GMT-14
+  UID=0
+  USER=root
+  _='I: set'
+  http_proxy=http://10.0.0.15:8000/
 I: uname -a
-  Linux p64b 4.19.0-10-arm64 #1 SMP Debian 4.19.132-1 (2020-07-24) aarch64 GNU/Linux
+  Linux i-capture-the-hostname 4.19.0-10-armmp-lpae #1 SMP Debian 4.19.132-1 (2020-07-24) armv7l GNU/Linux
 I: ls -l /bin
   total 3328
-  -rwxr-xr-x 1 root root 767656 Apr 17  2019 bash
-  -rwxr-xr-x 3 root root  26052 Jul 10  2019 bunzip2
-  -rwxr-xr-x 3 root root  26052 Jul 10  2019 bzcat
-  lrwxrwxrwx 1 root root      6 Jul 10  2019 bzcmp -> bzdiff
-  -rwxr-xr-x 1 root root   2227 Jul 10  2019 bzdiff
-  lrwxrwxrwx 1 root root      6 Jul 10  2019 bzegrep -> bzgrep
-  -rwxr-xr-x 1 root root   4877 Jun 24  2019 bzexe
-  lrwxrwxrwx 1 root root      6 Jul 10  2019 bzfgrep -> bzgrep
-  -rwxr-xr-x 1 root root   3641 Jul 10  2019 bzgrep
-  -rwxr-xr-x 3 root root  26052 Jul 10  2019 bzip2
-  -rwxr-xr-x 1 root root   9636 Jul 10  2019 bzip2recover
-  lrwxrwxrwx 1 root root      6 Jul 10  2019 bzless -> bzmore
-  -rwxr-xr-x 1 root root   1297 Jul 10  2019 bzmore
-  -rwxr-xr-x 1 root root  22432 Feb 28  2019 cat
-  -rwxr-xr-x 1 root root  38868 Feb 28  2019 chgrp
-  -rwxr-xr-x 1 root root  38836 Feb 28  2019 chmod
-  -rwxr-xr-x 1 root root  42972 Feb 28  2019 chown
-  -rwxr-xr-x 1 root root  88376 Feb 28  2019 cp
-  -rwxr-xr-x 1 root root  75516 Jan 17  2019 dash
-  -rwxr-xr-x 1 root root  71648 Feb 28  2019 date
-  -rwxr-xr-x 1 root root  51212 Feb 28  2019 dd
-  -rwxr-xr-x 1 root root  55672 Feb 28  2019 df
-  -rwxr-xr-x 1 root root  88444 Feb 28  2019 dir
-  -rwxr-xr-x 1 root root  54872 Jan  9  2019 dmesg
-  lrwxrwxrwx 1 root root      8 Sep 26  2018 dnsdomainname -> hostname
-  lrwxrwxrwx 1 root root      8 Sep 26  2018 domainname -> hostname
-  -rwxr-xr-x 1 root root  22364 Feb 28  2019 echo
-  -rwxr-xr-x 1 root root     28 Jan  7  2019 egrep
-  -rwxr-xr-x 1 root root  18260 Feb 28  2019 false
-  -rwxr-xr-x 1 root root     28 Jan  7  2019 fgrep
-  -rwxr-xr-x 1 root root  47356 Jan  9  2019 findmnt
-  -rwsr-xr-x 1 root root  21980 Apr 22 07:38 fusermount
-  -rwxr-xr-x 1 root root 124508 Jan  7  2019 grep
-  -rwxr-xr-x 2 root root   2345 Jan  5  2019 gunzip
-  -rwxr-xr-x 1 root root   6375 Jan  5  2019 gzexe
-  -rwxr-xr-x 1 root root  64232 Jan  5  2019 gzip
-  -rwxr-xr-x 1 root root  13784 Sep 26  2018 hostname
-  -rwxr-xr-x 1 root root  43044 Feb 28  2019 ln
-  -rwxr-xr-x 1 root root  34932 Jul 26  2018 login
-  -rwxr-xr-x 1 root root  88444 Feb 28  2019 ls
-  -rwxr-xr-x 1 root root  67036 Jan  9  2019 lsblk
-  -rwxr-xr-x 1 root root  47168 Feb 28  2019 mkdir
-  -rwxr-xr-x 1 root root  43040 Feb 28  2019 mknod
-  -rwxr-xr-x 1 root root  26552 Feb 28  2019 mktemp
-  -rwxr-xr-x 1 root root  26024 Jan  9  2019 more
-  -rwsr-xr-x 1 root root  34268 Jan  9  2019 mount
-  -rwxr-xr-x 1 root root   9688 Jan  9  2019 mountpoint
-  -rwxr-xr-x 1 root root  84284 Feb 28  2019 mv
-  lrwxrwxrwx 1 root root      8 Sep 26  2018 nisdomainname -> hostname
-  lrwxrwxrwx 1 root root     14 Feb 14  2019 pidof -> /sbin/killall5
-  -rwxr-xr-x 1 root root  22416 Feb 28  2019 pwd
-  lrwxrwxrwx 1 root root      4 Apr 17  2019 rbash -> bash
-  -rwxr-xr-x 1 root root  26504 Feb 28  2019 readlink
-  -rwxr-xr-x 1 root root  42968 Feb 28  2019 rm
-  -rwxr-xr-x 1 root root  26496 Feb 28  2019 rmdir
-  -rwxr-xr-x 1 root root  14136 Jan 21  2019 run-parts
-  -rwxr-xr-x 1 root root  76012 Dec 22  2018 sed
-  lrwxrwxrwx 1 root root      4 Aug 10 20:27 sh -> dash
-  -rwxr-xr-x 1 root root  22384 Feb 28  2019 sleep
-  -rwxr-xr-x 1 root root  51124 Feb 28  2019 stty
-  -rwsr-xr-x 1 root root  42472 Jan  9  2019 su
-  -rwxr-xr-x 1 root root  22392 Feb 28  2019 sync
-  -rwxr-xr-x 1 root root 283324 Apr 23  2019 tar
-  -rwxr-xr-x 1 root root   9808 Jan 21  2019 tempfile
-  -rwxr-xr-x 1 root root  63464 Feb 28  2019 touch
-  -rwxr-xr-x 1 root root  18260 Feb 28  2019 true
-  -rwxr-xr-x 1 root root   9636 Apr 22 07:38 ulockmgr_server
-  -rwsr-xr-x 1 root root  21976 Jan  9  2019 umount
-  -rwxr-xr-x 1 root root  22380 Feb 28  2019 uname
-  -rwxr-xr-x 2 root root   2345 Jan  5  2019 uncompress
-  -rwxr-xr-x 1 root root  88444 Feb 28  2019 vdir
-  -rwxr-xr-x 1 root root  21980 Jan  9  2019 wdctl
-  -rwxr-xr-x 1 root root    946 Jan 21  2019 which
-  lrwxrwxrwx 1 root root      8 Sep 26  2018 ypdomainname -> hostname
-  -rwxr-xr-x 1 root root   1983 Jan  5  2019 zcat
-  -rwxr-xr-x 1 root root   1677 Jan  5  2019 zcmp
-  -rwxr-xr-x 1 root root   5879 Jan  5  2019 zdiff
-  -rwxr-xr-x 1 root root     29 Jan  5  2019 zegrep
-  -rwxr-xr-x 1 root root     29 Jan  5  2019 zfgrep
-  -rwxr-xr-x 1 root root   2080 Jan  5  2019 zforce
-  -rwxr-xr-x 1 root root   7584 Jan  5  2019 zgrep
-  -rwxr-xr-x 1 root root   2205 Jan  5  2019 zless
-  -rwxr-xr-x 1 root root   1841 Jan  5  2019 zmore
-  -rwxr-xr-x 1 root root   4552 Jan  5  2019 znew
-I: user script /srv/workspace/pbuilder/27930/tmp/hooks/D02_print_environment finished
+  -rwxr-xr-x 1 root root 767656 Apr 18  2019 bash
+  -rwxr-xr-x 3 root root  26052 Jul 11  2019 bunzip2
+  -rwxr-xr-x 3 root root  26052 Jul 11  2019 bzcat
+  lrwxrwxrwx 1 root root      6 Jul 11  2019 bzcmp -> bzdiff
+  -rwxr-xr-x 1 root root   2227 Jul 11  2019 bzdiff
+  lrwxrwxrwx 1 root root      6 Jul 11  2019 bzegrep -> bzgrep
+  -rwxr-xr-x 1 root root   4877 Jun 25  2019 bzexe
+  lrwxrwxrwx 1 root root      6 Jul 11  2019 bzfgrep -> bzgrep
+  -rwxr-xr-x 1 root root   3641 Jul 11  2019 bzgrep
+  -rwxr-xr-x 3 root root  26052 Jul 11  2019 bzip2
+  -rwxr-xr-x 1 root root   9636 Jul 11  2019 bzip2recover
+  lrwxrwxrwx 1 root root      6 Jul 11  2019 bzless -> bzmore
+  -rwxr-xr-x 1 root root   1297 Jul 11  2019 bzmore
+  -rwxr-xr-x 1 root root  22432 Mar  1  2019 cat
+  -rwxr-xr-x 1 root root  38868 Mar  1  2019 chgrp
+  -rwxr-xr-x 1 root root  38836 Mar  1  2019 chmod
+  -rwxr-xr-x 1 root root  42972 Mar  1  2019 chown
+  -rwxr-xr-x 1 root root  88376 Mar  1  2019 cp
+  -rwxr-xr-x 1 root root  75516 Jan 18  2019 dash
+  -rwxr-xr-x 1 root root  71648 Mar  1  2019 date
+  -rwxr-xr-x 1 root root  51212 Mar  1  2019 dd
+  -rwxr-xr-x 1 root root  55672 Mar  1  2019 df
+  -rwxr-xr-x 1 root root  88444 Mar  1  2019 dir
+  -rwxr-xr-x 1 root root  54872 Jan 10  2019 dmesg
+  lrwxrwxrwx 1 root root      8 Sep 27  2018 dnsdomainname -> hostname
+  lrwxrwxrwx 1 root root      8 Sep 27  2018 domainname -> hostname
+  -rwxr-xr-x 1 root root  22364 Mar  1  2019 echo
+  -rwxr-xr-x 1 root root     28 Jan  8  2019 egrep
+  -rwxr-xr-x 1 root root  18260 Mar  1  2019 false
+  -rwxr-xr-x 1 root root     28 Jan  8  2019 fgrep
+  -rwxr-xr-x 1 root root  47356 Jan 10  2019 findmnt
+  -rwsr-xr-x 1 root root  21980 Apr 23 09:38 fusermount
+  -rwxr-xr-x 1 root root 124508 Jan  8  2019 grep
+  -rwxr-xr-x 2 root root   2345 Jan  6  2019 gunzip
+  -rwxr-xr-x 1 root root   6375 Jan  6  2019 gzexe
+  -rwxr-xr-x 1 root root  64232 Jan  6  2019 gzip
+  -rwxr-xr-x 1 root root  13784 Sep 27  2018 hostname
+  -rwxr-xr-x 1 root root  43044 Mar  1  2019 ln
+  -rwxr-xr-x 1 root root  34932 Jul 27  2018 login
+  -rwxr-xr-x 1 root root  88444 Mar  1  2019 ls
+  -rwxr-xr-x 1 root root  67036 Jan 10  2019 lsblk
+  -rwxr-xr-x 1 root root  47168 Mar  1  2019 mkdir
+  -rwxr-xr-x 1 root root  43040 Mar  1  2019 mknod
+  -rwxr-xr-x 1 root root  26552 Mar  1  2019 mktemp
+  -rwxr-xr-x 1 root root  26024 Jan 10  2019 more
+  -rwsr-xr-x 1 root root  34268 Jan 10  2019 mount
+  -rwxr-xr-x 1 root root   9688 Jan 10  2019 mountpoint
+  -rwxr-xr-x 1 root root  84284 Mar  1  2019 mv
+  lrwxrwxrwx 1 root root      8 Sep 27  2018 nisdomainname -> hostname
+  lrwxrwxrwx 1 root root     14 Feb 15  2019 pidof -> /sbin/killall5
+  -rwxr-xr-x 1 root root  22416 Mar  1  2019 pwd
+  lrwxrwxrwx 1 root root      4 Apr 18  2019 rbash -> bash
+  -rwxr-xr-x 1 root root  26504 Mar  1  2019 readlink
+  -rwxr-xr-x 1 root root  42968 Mar  1  2019 rm
+  -rwxr-xr-x 1 root root  26496 Mar  1  2019 rmdir
+  -rwxr-xr-x 1 root root  14136 Jan 22  2019 run-parts
+  -rwxr-xr-x 1 root root  76012 Dec 23  2018 sed
+  lrwxrwxrwx 1 root root      4 Aug 13 02:30 sh -> bash
+  lrwxrwxrwx 1 root root      4 Aug 11 22:26 sh.distrib -> dash
+  -rwxr-xr-x 1 root root  22384 Mar  1  2019 sleep
+  -rwxr-xr-x 1 root root  51124 Mar  1  2019 stty
+  -rwsr-xr-x 1 root root  42472 Jan 10  2019 su
+  -rwxr-xr-x 1 root root  22392 Mar  1  2019 sync
+  -rwxr-xr-x 1 root root 283324 Apr 24  2019 tar
+  -rwxr-xr-x 1 root root   9808 Jan 22  2019 tempfile
+  -rwxr-xr-x 1 root root  63464 Mar  1  2019 touch
+  -rwxr-xr-x 1 root root  18260 Mar  1  2019 true
+  -rwxr-xr-x 1 root root   9636 Apr 23 09:38 ulockmgr_server
+  -rwsr-xr-x 1 root root  21976 Jan 10  2019 umount
+  -rwxr-xr-x 1 root root  22380 Mar  1  2019 uname
+  -rwxr-xr-x 2 root root   2345 Jan  6  2019 uncompress
+  -rwxr-xr-x 1 root root  88444 Mar  1  2019 vdir
+  -rwxr-xr-x 1 root root  21980 Jan 10  2019 wdctl
+  -rwxr-xr-x 1 root root    946 Jan 22  2019 which
+  lrwxrwxrwx 1 root root      8 Sep 27  2018 ypdomainname -> hostname
+  -rwxr-xr-x 1 root root   1983 Jan  6  2019 zcat
+  -rwxr-xr-x 1 root root   1677 Jan  6  2019 zcmp
+  -rwxr-xr-x 1 root root   5879 Jan  6  2019 zdiff
+  -rwxr-xr-x 1 root root     29 Jan  6  2019 zegrep
+  -rwxr-xr-x 1 root root     29 Jan  6  2019 zfgrep
+  -rwxr-xr-x 1 root root   2080 Jan  6  2019 zforce
+  -rwxr-xr-x 1 root root   7584 Jan  6  2019 zgrep
+  -rwxr-xr-x 1 root root   2205 Jan  6  2019 zless
+  -rwxr-xr-x 1 root root   1841 Jan  6  2019 zmore
+  -rwxr-xr-x 1 root root   4552 Jan  6  2019 znew
+I: user script /srv/workspace/pbuilder/7978/tmp/hooks/D02_print_environment finished
  -> Attempting to satisfy build-dependencies
  -> Creating pbuilder-satisfydepends-dummy package
 Package: pbuilder-satisfydepends-dummy
@@ -544,7 +578,7 @@
 Get: 349 http://deb.debian.org/debian buster/main armhf fonts-freefont-otf all 20120503-9 [3169 kB]
 Get: 350 http://deb.debian.org/debian buster/main armhf dh-octave-autopkgtest all 0.6.2 [8308 B]
 Get: 351 http://deb.debian.org/debian buster/main armhf dh-octave all 0.6.2 [19.9 kB]
-Fetched 111 MB in 32s (3489 kB/s)
+Fetched 111 MB in 26s (4257 kB/s)
 debconf: delaying package configuration, since apt-utils is not installed
 Selecting previously unselected package libbsd0:armhf.
 (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 ... 18932 files and directories currently installed.)
@@ -1995,7 +2029,7 @@
 Building tag database...
  -> Finished parsing the build-deps
 I: Building the package
-I: Running cd /build/octave-nurbs-1.3.13/ && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games" HOME="/nonexistent/first-build" dpkg-buildpackage -us -uc -b
+I: Running cd /build/octave-nurbs-1.3.13/ && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path" HOME="/nonexistent/second-build" dpkg-buildpackage -us -uc -b
 dpkg-buildpackage: info: source package octave-nurbs
 dpkg-buildpackage: info: source version 1.3.13-5
 dpkg-buildpackage: info: source distribution unstable
@@ -2033,9 +2067,9 @@
 /usr/bin/mkoctfile --verbose basisfun.cc low_level_functions.o
 /usr/bin/mkoctfile --verbose basisfunder.cc low_level_functions.o
 /usr/bin/mkoctfile --verbose bspderiv.cc low_level_functions.o
+/usr/bin/mkoctfile --verbose bspeval.cc low_level_functions.o
 g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security    basisfun.cc -o basisfun.o
 g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security    basisfunder.cc -o basisfunder.o
-/usr/bin/mkoctfile --verbose bspeval.cc low_level_functions.o
 g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security    bspderiv.cc -o bspderiv.o
 g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security    bspeval.cc -o bspeval.o
 bspeval.cc: In function 'bool bspeval_bad_arguments(const octave_value_list&)':
@@ -2052,17 +2086,17 @@
 /usr/include/octave-4.4.1/octave/../octave/ov.h:707:8: note: declared here
    bool is_real_type (void) const
         ^~~~~~~~~~~~
+g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o bspeval.oct  low_level_functions.o bspeval.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
 g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o basisfunder.oct  low_level_functions.o basisfunder.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
-g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o basisfun.oct  low_level_functions.o basisfun.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
 g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o bspderiv.oct  low_level_functions.o bspderiv.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
-g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o bspeval.oct  low_level_functions.o bspeval.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
+g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o basisfun.oct  low_level_functions.o basisfun.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
 /usr/bin/mkoctfile --verbose curvederivcpts.cc low_level_functions.o
 g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security    curvederivcpts.cc -o curvederivcpts.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-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security    nrb_srf_basisfun__.cc -o nrb_srf_basisfun__.o
 /usr/bin/mkoctfile --verbose nrb_srf_basisfun_der__.cc low_level_functions.o
-/usr/bin/mkoctfile --verbose nrbsurfderiveval.cc low_level_functions.o
 g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security    nrb_srf_basisfun_der__.cc -o nrb_srf_basisfun_der__.o
+/usr/bin/mkoctfile --verbose nrbsurfderiveval.cc low_level_functions.o
 g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security    nrbsurfderiveval.cc -o nrbsurfderiveval.o
 nrb_srf_basisfun__.cc: In function 'octave_value_list Fnrb_srf_basisfun__(const octave_value_list&, int)':
 nrb_srf_basisfun__.cc:69:64: warning: 'octave_value_list feval(const string&, const octave_value_list&, int)' is deprecated: [4.4]: use 'octave::feval' instead [-Wdeprecated-declarations]
@@ -2161,6 +2195,7 @@
 g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security    tbasisfun.cc -o tbasisfun.o
 g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o nrb_srf_basisfun_der__.oct  low_level_functions.o nrb_srf_basisfun_der__.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
 g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o surfderivcpts.oct  low_level_functions.o surfderivcpts.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
+g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o surfderiveval.oct  low_level_functions.o surfderiveval.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
 tbasisfun.cc: In function 'octave_value_list Ftbasisfun(const octave_value_list&, int)':
 tbasisfun.cc:159:26: warning: 'bool octave_value::is_cell() const' is deprecated: [4.4]: use 'iscell' instead [-Wdeprecated-declarations]
    if (! args(2).is_cell ())
@@ -2175,7 +2210,6 @@
 /usr/include/octave-4.4.1/octave/../octave/ov.h:540:8: note: declared here
    bool is_cell (void) const
         ^~~~~~~
-g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o surfderiveval.oct  low_level_functions.o surfderiveval.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
 g++ -I/usr/include/octave-4.4.1/octave/.. -I/usr/include/octave-4.4.1/octave  -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-nurbs-1.3.13=. -fstack-protector-strong -Wformat -Werror=format-security -shared -Wl,-Bsymbolic  -o tbasisfun.oct  low_level_functions.o tbasisfun.o   -L/usr/lib/arm-linux-gnueabihf/octave/4.4.1 -L/usr/lib/arm-linux-gnueabihf -loctinterp -loctave -Wl,-z,relro 
 make[1]: Leaving directory '/build/octave-nurbs-1.3.13/src'
 copyfile /build/octave-nurbs-1.3.13/./src/basisfun.oct /build/octave-nurbs-1.3.13/./src/basisfunder.oct /build/octave-nurbs-1.3.13/./src/bspderiv.oct /build/octave-nurbs-1.3.13/./src/bspeval.oct /build/octave-nurbs-1.3.13/./src/curvederivcpts.oct /build/octave-nurbs-1.3.13/./src/nrb_srf_basisfun__.oct /build/octave-nurbs-1.3.13/./src/nrb_srf_basisfun_der__.oct /build/octave-nurbs-1.3.13/./src/nrbsurfderiveval.oct /build/octave-nurbs-1.3.13/./src/surfderivcpts.oct /build/octave-nurbs-1.3.13/./src/surfderiveval.oct /build/octave-nurbs-1.3.13/./src/tbasisfun.oct /build/octave-nurbs-1.3.13/./inst/arm-unknown-linux-gnueabihf-api-v52
@@ -2186,425 +2220,7 @@
 Checking m files ...
 warning: function /build/octave-nurbs-1.3.13/debian/octave-nurbs/usr/share/octave/packages/nurbs-1.3.13/vecnorm.m shadows a core library function
 warning: called from
-    /tmp/fileJWalk0 at line 12 column 1
-[inst/kntrefine.m]
->>>>> /build/octave-nurbs-1.3.13/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/nrbdeval.m]
->>>>> /build/octave-nurbs-1.3.13/inst/nrbdeval.m
-***** demo
- crv = nrbtestcrv;
- nrbplot(crv,48);
- title('First derivatives along a test curve.');
- 
- tt = linspace(0.0,1.0,9);
- 
- dcrv = nrbderiv(crv);
- 
- [p1, dp] = nrbdeval(crv,dcrv,tt);
- 
- p2 = vecnorm(dp);
- 
- hold on;
- plot(p1(1,:),p1(2,:),'ro');
- h = quiver(p1(1,:),p1(2,:),p2(1,:),p2(2,:),0);
- set(h,'Color','black');
- hold off;
-***** 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,:,:)));
- set(h,'FaceColor','blue','EdgeColor','blue');
- title('First derivatives over a test surface.');
-
- npts = 5;
- tt = linspace(0.0,1.0,npts);
- dsrf = nrbderiv(srf);
- 
- [p1, dp] = nrbdeval(srf, dsrf, {tt, tt});
- 
- up2 = vecnorm(dp{1});
- vp2 = vecnorm(dp{2});
- 
- hold on;
- plot3(p1(1,:),p1(2,:),p1(3,:),'ro');
- h1 = quiver3(p1(1,:),p1(2,:),p1(3,:),up2(1,:),up2(2,:),up2(3,:));
- h2 = quiver3(p1(1,:),p1(2,:),p1(3,:),vp2(1,:),vp2(2,:),vp2(3,:));
- set(h1,'Color','black');
- set(h2,'Color','black');
- 
- 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]';
- ders = nrbderiv(nurbs);
- [points,jac] = nrbdeval(nurbs,ders,tt);
- assert(points,tt,1e-10)
- assert(jac{1}(1,:,:),ones(size(jac{1}(1,:,:))),1e-12)
- assert(jac{2}(2,:,:),ones(size(jac{2}(2,:,:))),1e-12)
- assert(jac{3}(3,:,:),ones(size(jac{3}(3,:,:))),1e-12)
- 
-***** test
- knots{1} = [0 0 0 1 1 1];
- knots{2} = [0 0 0 0 1 1 1 1];
- knots{3} = [0 0 0 1 1 1];
- cx = [0 0 1]; nx = length(cx);
- cy = [0 0 0 1]; ny = length(cy);
- cz = [0 0.5 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]';
- dnurbs = nrbderiv(nurbs);
- [points, jac] = nrbdeval(nurbs,dnurbs,tt);
- assert(points,[y.^3 x.^2 z]',1e-10);
- assert(jac{2}(1,:,:),3*y'.^2,1e-12)
- 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/basisfun.m]
->>>>> /build/octave-nurbs-1.3.13/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/nrbeval_der_p.m]
->>>>> /build/octave-nurbs-1.3.13/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:nrb.number(1) * nrb.number(2)
-   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:nrb.number(1) * nrb.number(2) * nrb.number(3)
-   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);
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/nrbmeasure.m]
->>>>> /build/octave-nurbs-1.3.13/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/nrbbasisfunder.m]
->>>>> /build/octave-nurbs-1.3.13/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)
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/nrbreverse.m]
->>>>> /build/octave-nurbs-1.3.13/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/nrbbasisfun.m]
->>>>> /build/octave-nurbs-1.3.13/inst/nrbbasisfun.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);
- B = nrbbasisfun (u, nrb);
- xplot = sum(bsxfun(@(x,y) x.*y, B, x),2);
- plot(xplot, B)
- title('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);
- B = nrbbasisfun (u, nrb);
- xplot = sum(bsxfun(@(x,y) x.*y, B, x),2);
-
- yy = y; yy(1) = 1;
- nrb2 = nrbmak ([x.*w;yy;y;w], U); 
- aux = nrbeval(nrb2,u);
- %figure, plot(xplot, B(:,1), aux(1,:).', w(1)*aux(2,:).')
- assert(B(:,1), w(1)*aux(2,:).', 1e-6)
- 
- yy = y; yy(2) = 1;
- nrb2 = nrbmak ([x.*w;yy;y;w], U);
- aux = nrbeval(nrb2, u);
- %figure, plot(xplot, B(:,2), aux(1,:).', w(2)*aux(2,:).')
- assert(B(:,2), w(2)*aux(2,:).', 1e-6)
-
- yy = y; yy(3) = 1;
- nrb2 = nrbmak ([x.*w;yy;y;w], U);
- aux = nrbeval(nrb2,u);
- %figure, plot(xplot, B(:,3), aux(1,:).', w(3)*aux(2,:).')
- assert(B(:,3), w(3)*aux(2,:).', 1e-6)
-
- yy = y; yy(4) = 1;
- nrb2 = nrbmak ([x.*w;yy;y;w], U);
- aux = nrbeval(nrb2,u);
- %figure, plot(xplot, B(:,4), aux(1,:).', w(4)*aux(2,:).')
- assert(B(:,4), w(4)*aux(2,:).', 1e-6)
-***** 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)});
- u = rand (1, 30); v = rand (1, 10);
- u = u - min (u); u = u / max (u);
- v = v - min (v); v = v / max (v);
- [B, N] = nrbbasisfun ({u, v}, nrb);
- assert (sum(B, 2), ones(300, 1), 1e-6)
- assert (all (all (B<=1)), true)
- assert (all (all (B>=0)), true)
- 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)
-***** test
- p1 = 2;   p2 = 3;   p3 = 2;
- n1 = 4;   n2 = 5;   n3 = 4;
- Lx = 1;  Ly = 1;   Lz = 1; 
- crv = nrbline([1 0], [2 0]);
- nrb = nrbtransp (nrbrevolve (crv, [], [0 0 1], pi/2));
- nrb = nrbextrude (nrb, [0 0 1]);
- nrb = nrbdegelev (nrb, [p1-1, p2-2, p3-1]);
- aux1 = linspace(0,1,n1); aux2 = linspace(0,1,n2); aux3 = linspace(0,1,n3);
- nrb = nrbkntins  (nrb, {aux1(2:end-1), aux2(2:end-1), aux3(2:end-1)});
- 
- u = rand (1, 12); v = rand (1, 10); w = rand (1, 15);
- u = u - min (u); u = u / max (u);
- v = v - min (v); v = v / max (v);
- w = w - min (w); w = w / max (w);
- [B, N] = nrbbasisfun ({u, v, w}, nrb);
- assert (all(sum(B, 2) - ones(numel(u)*numel(v)*numel(w),1) < 1e-6))
- assert (all (all (B <= 1)) == true)
- assert (all (all (B >= 0)) == true)
- 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)
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/crvkntremove.m]
->>>>> /build/octave-nurbs-1.3.13/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
+    /tmp/filelcx9fe at line 12 column 1
 [inst/nrbderiv.m]
 >>>>> /build/octave-nurbs-1.3.13/inst/nrbderiv.m
 ***** demo
@@ -2882,44 +2498,6 @@
  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;
 4 tests, 4 passed, 0 known failure, 0 skipped
-[inst/nrbtransp.m]
->>>>> /build/octave-nurbs-1.3.13/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/nrbnumbasisfun.m]
->>>>> /build/octave-nurbs-1.3.13/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/numbasisfun.m]
 >>>>> /build/octave-nurbs-1.3.13/inst/numbasisfun.m
 ***** test
@@ -2934,108 +2512,6 @@
   B = numbasisfun (s, u, p, U);
   assert (B, Bref)
 1 test, 1 passed, 0 known failure, 0 skipped
-[inst/nrbunclamp.m]
->>>>> /build/octave-nurbs-1.3.13/inst/nrbunclamp.m
-***** demo
- crv = nrbcirc (1,[],0,2*pi/3);
- crv = nrbdegelev (crv, 2);
- figure
- nrbctrlplot (crv); hold on
- nrbctrlplot (nrbtform (nrbunclamp (crv, 1), vectrans([-0.4, -0.4])));
- nrbctrlplot (nrbtform (nrbunclamp (crv, 2), vectrans([-0.8, -0.8])));
- nrbctrlplot (nrbtform (nrbunclamp (crv, 3), vectrans([-1.6, -1.6])));
- title ('Original curve and unclamped versions')
-***** test
- crv = nrbdegelev (nrbtestcrv,2);
- x = linspace (0, 1, 100);
- F = nrbeval (crv, x);
- ucrv = nrbunclamp (crv, 0);
- assert (F, nrbeval(ucrv, x));
- ucrv = nrbunclamp (crv, 1);
- assert (F, nrbeval(ucrv, x), 1e-14);
- ucrv = nrbunclamp (crv, 2);
- assert (F, nrbeval(ucrv, x), 1e-14);
- ucrv = nrbunclamp (crv, 3);
- assert (F, nrbeval(ucrv, x), 1e-14);
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/surfderivcpts.m]
->>>>> /build/octave-nurbs-1.3.13/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/nrbinverse.m]
->>>>> /build/octave-nurbs-1.3.13/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/kntbrkdegmult.m]
->>>>> /build/octave-nurbs-1.3.13/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/nrbcrvderiveval.m]
->>>>> /build/octave-nurbs-1.3.13/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/nrbsquare.m]
 >>>>> /build/octave-nurbs-1.3.13/inst/nrbsquare.m
 ***** test
@@ -3049,177 +2525,61 @@
  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/kntbrkdegreg.m]
->>>>> /build/octave-nurbs-1.3.13/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])
+[inst/aveknt.m]
+>>>>> /build/octave-nurbs-1.3.13/inst/aveknt.m
 ***** 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])
+ knt = [0 0 0 0.5 1 1 1];
+ pts = aveknt (knt, 3);
+ assert (pts - [0 1/4 3/4 1] < 1e-14)
+
 ***** 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])
+ 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
- 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/nrbeval_der_w.m]
->>>>> /build/octave-nurbs-1.3.13/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: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: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)
+ 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/octave-nurbs-1.3.13/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/nrbpermute.m]
->>>>> /build/octave-nurbs-1.3.13/inst/nrbpermute.m
-***** demo
- vol = nrbrevolve (nrb4surf ([1 0], [2 0], [1 1], [2 1]), [0 0 0], [0 1 0], pi/8);
- nrbplot(vol,[5 10 20]);
- title('NURBS volume and the same after reordering the directions')
- hold on
- vol.coefs(1,:,:) = vol.coefs(1,:,:) + 2;
- vol = nrbpermute(vol,[2 3 1]);
- nrbplot(vol,[5 10 20]);
- hold off
+[inst/nrbmeasure.m]
+>>>>> /build/octave-nurbs-1.3.13/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/surfderivcpts.m]
+>>>>> /build/octave-nurbs-1.3.13/inst/surfderivcpts.m
 ***** test
- vol = nrbrevolve (nrb4surf ([1 0], [2 0], [1 1], [2 1]), [0 0 0], [0 1 0], pi/8);
- perm1 = [1 3 2];
- perm2 = [2 1 3];
- vol2 = nrbpermute (vol, perm1);
- vol3 = nrbpermute (vol, perm2);
- assert (vol.number(perm1), vol2.number)
- assert (vol.order(perm1), vol2.order)
- assert ({vol.knots{perm1}}, vol2.knots)
- assert (permute(vol.coefs, [1, perm1+1]), vol2.coefs)
- assert (vol.number(perm2), vol3.number)
- assert (vol.order(perm2), vol3.order)
- assert ({vol.knots{perm2}}, vol3.knots)
- assert (permute(vol.coefs, [1, perm2+1]), vol3.coefs)
+ 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/basisfunder.m]
->>>>> /build/octave-nurbs-1.3.13/inst/basisfunder.m
-***** test
- k    = [0 0 0 0 1 1 1 1];
- p    = 3;
- u    = rand (1);
- i    = findspan (numel(k)-p-2, p, u, k);
- ders = basisfunder (i, p, u, k, 1);
- sumders = sum (squeeze(ders), 2);
- assert (sumders(1), 1, 1e-15);
- assert (sumders(2:end), 0, 1e-15);
-***** test
- k    = [0 0 0 0 1/3 2/3 1 1 1 1];
- p    = 3;
- u    = rand (1);
- i    = findspan (numel(k)-p-2, p, u, k);
- ders = basisfunder (i, p, u, k, 7); 
- sumders = sum (squeeze(ders), 2);
- assert (sumders(1), 1, 1e-15);
- assert (sumders(2:end), zeros(rows(squeeze(ders))-1, 1), 1e-13);
-***** test
- k    = [0 0 0 0 1/3 2/3 1 1 1 1];
- p    = 3;
- u    = rand (100, 1);
- i    = findspan (numel(k)-p-2, p, u, k);
- ders = basisfunder (i, p, u, k, 7);
- for ii=1:10
-   sumders = sum (squeeze(ders(ii,:,:)), 2);
-   assert (sumders(1), 1, 1e-15);
-   assert (sumders(2:end), zeros(rows(squeeze(ders(ii,:,:)))-1, 1), 1e-13);
- end
- 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/findspan.m]
 >>>>> /build/octave-nurbs-1.3.13/inst/findspan.m
 ***** test
@@ -3236,69 +2596,50 @@
  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/nrbclamp.m]
->>>>> /build/octave-nurbs-1.3.13/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/surfderiveval.m]
->>>>> /build/octave-nurbs-1.3.13/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/nrbeval.m]
->>>>> /build/octave-nurbs-1.3.13/inst/nrbeval.m
+[inst/nrbdeval.m]
+>>>>> /build/octave-nurbs-1.3.13/inst/nrbdeval.m
+***** demo
+ crv = nrbtestcrv;
+ nrbplot(crv,48);
+ title('First derivatives along a test curve.');
+ 
+ tt = linspace(0.0,1.0,9);
+ 
+ dcrv = nrbderiv(crv);
+ 
+ [p1, dp] = nrbdeval(crv,dcrv,tt);
+ 
+ p2 = vecnorm(dp);
+ 
+ hold on;
+ plot(p1(1,:),p1(2,:),'ro');
+ h = quiver(p1(1,:),p1(2,:),p2(1,:),p2(2,:),0);
+ set(h,'Color','black');
+ hold off;
 ***** 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
+ set(h,'FaceColor','blue','EdgeColor','blue');
+ title('First derivatives over a test surface.');
+
+ npts = 5;
+ tt = linspace(0.0,1.0,npts);
+ dsrf = nrbderiv(srf);
+ 
+ [p1, dp] = nrbdeval(srf, dsrf, {tt, tt});
+ 
+ up2 = vecnorm(dp{1});
+ vp2 = vecnorm(dp{2});
+ 
+ hold on;
+ plot3(p1(1,:),p1(2,:),p1(3,:),'ro');
+ h1 = quiver3(p1(1,:),p1(2,:),p1(3,:),up2(1,:),up2(2,:),up2(3,:));
+ h2 = quiver3(p1(1,:),p1(2,:),p1(3,:),vp2(1,:),vp2(2,:),vp2(3,:));
+ set(h1,'Color','black');
+ set(h2,'Color','black');
+ 
+ hold off;
 ***** test
  knots{1} = [0 0 0 1 1 1];
  knots{2} = [0 0 0 .5 1 1 1];
@@ -3313,34 +2654,20 @@
  nurbs = nrbmak(coefs, knots);
  x = rand(5,1); y = rand(5,1); z = rand(5,1);
  tt = [x y z]';
- points = nrbeval(nurbs,tt);
-
+ ders = nrbderiv(nurbs);
+ [points,jac] = nrbdeval(nurbs,ders,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);
-
+ assert(jac{1}(1,:,:),ones(size(jac{1}(1,:,:))),1e-12)
+ assert(jac{2}(2,:,:),ones(size(jac{2}(2,:,:))),1e-12)
+ assert(jac{3}(3,:,:),ones(size(jac{3}(3,:,:))),1e-12)
+ 
 ***** 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];
+ knots{3} = [0 0 0 1 1 1];
  cx = [0 0 1]; nx = length(cx);
  cy = [0 0 0 1]; ny = length(cy);
- cz = [0 1]; nz = length(cz);
+ cz = [0 0.5 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]);
@@ -3349,30 +2676,34 @@
  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]';
+ dnurbs = nrbderiv(nurbs);
+ [points, jac] = nrbdeval(nurbs,dnurbs,tt);
  assert(points,[y.^3 x.^2 z]',1e-10);
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/aveknt.m]
->>>>> /build/octave-nurbs-1.3.13/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);
-
+ assert(jac{2}(1,:,:),3*y'.^2,1e-12)
+ 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/basisfun.m]
+>>>>> /build/octave-nurbs-1.3.13/inst/basisfun.m
 ***** 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
+  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/nrbsurfderiveval.m]
 >>>>> /build/octave-nurbs-1.3.13/inst/nrbsurfderiveval.m
 ***** test
@@ -3582,14 +2913,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);
 21 tests, 21 passed, 0 known failure, 0 skipped
-[inst/curvederivcpts.m]
->>>>> /build/octave-nurbs-1.3.13/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/tbasisfun.m]
 >>>>> /build/octave-nurbs-1.3.13/inst/tbasisfun.m
 ***** demo
@@ -3623,6 +2946,717 @@
  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/nrbbasisfun.m]
+>>>>> /build/octave-nurbs-1.3.13/inst/nrbbasisfun.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);
+ B = nrbbasisfun (u, nrb);
+ xplot = sum(bsxfun(@(x,y) x.*y, B, x),2);
+ plot(xplot, B)
+ title('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);
+ B = nrbbasisfun (u, nrb);
+ xplot = sum(bsxfun(@(x,y) x.*y, B, x),2);
+
+ yy = y; yy(1) = 1;
+ nrb2 = nrbmak ([x.*w;yy;y;w], U); 
+ aux = nrbeval(nrb2,u);
+ %figure, plot(xplot, B(:,1), aux(1,:).', w(1)*aux(2,:).')
+ assert(B(:,1), w(1)*aux(2,:).', 1e-6)
+ 
+ yy = y; yy(2) = 1;
+ nrb2 = nrbmak ([x.*w;yy;y;w], U);
+ aux = nrbeval(nrb2, u);
+ %figure, plot(xplot, B(:,2), aux(1,:).', w(2)*aux(2,:).')
+ assert(B(:,2), w(2)*aux(2,:).', 1e-6)
+
+ yy = y; yy(3) = 1;
+ nrb2 = nrbmak ([x.*w;yy;y;w], U);
+ aux = nrbeval(nrb2,u);
+ %figure, plot(xplot, B(:,3), aux(1,:).', w(3)*aux(2,:).')
+ assert(B(:,3), w(3)*aux(2,:).', 1e-6)
+
+ yy = y; yy(4) = 1;
+ nrb2 = nrbmak ([x.*w;yy;y;w], U);
+ aux = nrbeval(nrb2,u);
+ %figure, plot(xplot, B(:,4), aux(1,:).', w(4)*aux(2,:).')
+ assert(B(:,4), w(4)*aux(2,:).', 1e-6)
+***** 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)});
+ u = rand (1, 30); v = rand (1, 10);
+ u = u - min (u); u = u / max (u);
+ v = v - min (v); v = v / max (v);
+ [B, N] = nrbbasisfun ({u, v}, nrb);
+ assert (sum(B, 2), ones(300, 1), 1e-6)
+ assert (all (all (B<=1)), true)
+ assert (all (all (B>=0)), true)
+ 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)
+***** test
+ p1 = 2;   p2 = 3;   p3 = 2;
+ n1 = 4;   n2 = 5;   n3 = 4;
+ Lx = 1;  Ly = 1;   Lz = 1; 
+ crv = nrbline([1 0], [2 0]);
+ nrb = nrbtransp (nrbrevolve (crv, [], [0 0 1], pi/2));
+ nrb = nrbextrude (nrb, [0 0 1]);
+ nrb = nrbdegelev (nrb, [p1-1, p2-2, p3-1]);
+ aux1 = linspace(0,1,n1); aux2 = linspace(0,1,n2); aux3 = linspace(0,1,n3);
+ nrb = nrbkntins  (nrb, {aux1(2:end-1), aux2(2:end-1), aux3(2:end-1)});
+ 
+ u = rand (1, 12); v = rand (1, 10); w = rand (1, 15);
+ u = u - min (u); u = u / max (u);
+ v = v - min (v); v = v / max (v);
+ w = w - min (w); w = w / max (w);
+ [B, N] = nrbbasisfun ({u, v, w}, nrb);
+ assert (all(sum(B, 2) - ones(numel(u)*numel(v)*numel(w),1) < 1e-6))
+ assert (all (all (B <= 1)) == true)
+ assert (all (all (B >= 0)) == true)
+ 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)
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/nrbreverse.m]
+>>>>> /build/octave-nurbs-1.3.13/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/kntbrkdegmult.m]
+>>>>> /build/octave-nurbs-1.3.13/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.m]
+>>>>> /build/octave-nurbs-1.3.13/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/nrbeval_der_w.m]
+>>>>> /build/octave-nurbs-1.3.13/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: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: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/nrbbasisfunder.m]
+>>>>> /build/octave-nurbs-1.3.13/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)
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/curvederiveval.m]
+>>>>> /build/octave-nurbs-1.3.13/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/kntrefine.m]
+>>>>> /build/octave-nurbs-1.3.13/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/nrbtransp.m]
+>>>>> /build/octave-nurbs-1.3.13/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/basisfunder.m]
+>>>>> /build/octave-nurbs-1.3.13/inst/basisfunder.m
+***** test
+ k    = [0 0 0 0 1 1 1 1];
+ p    = 3;
+ u    = rand (1);
+ i    = findspan (numel(k)-p-2, p, u, k);
+ ders = basisfunder (i, p, u, k, 1);
+ sumders = sum (squeeze(ders), 2);
+ assert (sumders(1), 1, 1e-15);
+ assert (sumders(2:end), 0, 1e-15);
+***** test
+ k    = [0 0 0 0 1/3 2/3 1 1 1 1];
+ p    = 3;
+ u    = rand (1);
+ i    = findspan (numel(k)-p-2, p, u, k);
+ ders = basisfunder (i, p, u, k, 7); 
+ sumders = sum (squeeze(ders), 2);
+ assert (sumders(1), 1, 1e-15);
+ assert (sumders(2:end), zeros(rows(squeeze(ders))-1, 1), 1e-13);
+***** test
+ k    = [0 0 0 0 1/3 2/3 1 1 1 1];
+ p    = 3;
+ u    = rand (100, 1);
+ i    = findspan (numel(k)-p-2, p, u, k);
+ ders = basisfunder (i, p, u, k, 7);
+ for ii=1:10
+   sumders = sum (squeeze(ders(ii,:,:)), 2);
+   assert (sumders(1), 1, 1e-15);
+   assert (sumders(2:end), zeros(rows(squeeze(ders(ii,:,:)))-1, 1), 1e-13);
+ end
+ 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/nrbcrvderiveval.m]
+>>>>> /build/octave-nurbs-1.3.13/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/nrbinverse.m]
+>>>>> /build/octave-nurbs-1.3.13/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/nrbeval_der_p.m]
+>>>>> /build/octave-nurbs-1.3.13/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:nrb.number(1) * nrb.number(2)
+   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:nrb.number(1) * nrb.number(2) * nrb.number(3)
+   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);
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/surfderiveval.m]
+>>>>> /build/octave-nurbs-1.3.13/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/curvederivcpts.m]
+>>>>> /build/octave-nurbs-1.3.13/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/crvkntremove.m]
+>>>>> /build/octave-nurbs-1.3.13/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/nrbclamp.m]
+>>>>> /build/octave-nurbs-1.3.13/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/nrbpermute.m]
+>>>>> /build/octave-nurbs-1.3.13/inst/nrbpermute.m
+***** demo
+ vol = nrbrevolve (nrb4surf ([1 0], [2 0], [1 1], [2 1]), [0 0 0], [0 1 0], pi/8);
+ nrbplot(vol,[5 10 20]);
+ title('NURBS volume and the same after reordering the directions')
+ hold on
+ vol.coefs(1,:,:) = vol.coefs(1,:,:) + 2;
+ vol = nrbpermute(vol,[2 3 1]);
+ nrbplot(vol,[5 10 20]);
+ hold off
+***** test
+ vol = nrbrevolve (nrb4surf ([1 0], [2 0], [1 1], [2 1]), [0 0 0], [0 1 0], pi/8);
+ perm1 = [1 3 2];
+ perm2 = [2 1 3];
+ vol2 = nrbpermute (vol, perm1);
+ vol3 = nrbpermute (vol, perm2);
+ assert (vol.number(perm1), vol2.number)
+ assert (vol.order(perm1), vol2.order)
+ assert ({vol.knots{perm1}}, vol2.knots)
+ assert (permute(vol.coefs, [1, perm1+1]), vol2.coefs)
+ assert (vol.number(perm2), vol3.number)
+ assert (vol.order(perm2), vol3.order)
+ 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/nrbnumbasisfun.m]
+>>>>> /build/octave-nurbs-1.3.13/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/nrbunclamp.m]
+>>>>> /build/octave-nurbs-1.3.13/inst/nrbunclamp.m
+***** demo
+ crv = nrbcirc (1,[],0,2*pi/3);
+ crv = nrbdegelev (crv, 2);
+ figure
+ nrbctrlplot (crv); hold on
+ nrbctrlplot (nrbtform (nrbunclamp (crv, 1), vectrans([-0.4, -0.4])));
+ nrbctrlplot (nrbtform (nrbunclamp (crv, 2), vectrans([-0.8, -0.8])));
+ nrbctrlplot (nrbtform (nrbunclamp (crv, 3), vectrans([-1.6, -1.6])));
+ title ('Original curve and unclamped versions')
+***** test
+ crv = nrbdegelev (nrbtestcrv,2);
+ x = linspace (0, 1, 100);
+ F = nrbeval (crv, x);
+ ucrv = nrbunclamp (crv, 0);
+ assert (F, nrbeval(ucrv, x));
+ ucrv = nrbunclamp (crv, 1);
+ assert (F, nrbeval(ucrv, x), 1e-14);
+ ucrv = nrbunclamp (crv, 2);
+ assert (F, nrbeval(ucrv, x), 1e-14);
+ ucrv = nrbunclamp (crv, 3);
+ assert (F, nrbeval(ucrv, x), 1e-14);
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/kntbrkdegreg.m]
+>>>>> /build/octave-nurbs-1.3.13/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
 Checking C++ files ...
 [src/basisfun.cc]
 >>>>> /build/octave-nurbs-1.3.13/src/basisfun.cc
@@ -3648,34 +3682,39 @@
   B = basisfun(s, u, p, U);
   assert (B, Bref, 1e-5);
 2 tests, 2 passed, 0 known failure, 0 skipped
-[src/surfderivcpts.cc]
->>>>> /build/octave-nurbs-1.3.13/src/surfderivcpts.cc
+[src/tbasisfun.cc]
+>>>>> /build/octave-nurbs-1.3.13/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
- plane = nrbdegelev(nrb4surf([0 0], [0 1], [1 0], [1 1]), [1, 1]);
- 
- 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,:,:)), 2);
- 
- 
- pkl2 = [  0   0   0   1   0   0   0   0   0   0   0   0   1   0 ...
- 	0   0   0   0   0   0   0   1   0   0   0   0   0 0.5   0 ...
- 	0   1   0   0   0   0   0 0.5   0   0   1   0   0   0   0 ...
- 	0 0.5  0   0   1   0   0   0   0   0   1   0   0   0   0 ...
- 	0   0   0   0   1   0   0   0   0   0   0   0   0   1   0 ...
- 	0   0   0   0   0   0   0]';
- 
- assert (pkl(:),pkl2);
-1 test, 1 passed, 0 known failure, 0 skipped
-[src/curvederivcpts.cc]
->>>>> /build/octave-nurbs-1.3.13/src/curvederivcpts.cc
+ 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
- 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
+ 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/surfderiveval.cc]
 >>>>> /build/octave-nurbs-1.3.13/src/surfderiveval.cc
 ***** shared srf
@@ -3958,39 +3997,34 @@
  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/octave-nurbs-1.3.13/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);
+[src/surfderivcpts.cc]
+>>>>> /build/octave-nurbs-1.3.13/src/surfderivcpts.cc
 ***** 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);
+ plane = nrbdegelev(nrb4surf([0 0], [0 1], [1 0], [1 1]), [1, 1]);
+ 
+ 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,:,:)), 2);
+ 
+ 
+ pkl2 = [  0   0   0   1   0   0   0   0   0   0   0   0   1   0 ...
+ 	0   0   0   0   0   0   0   1   0   0   0   0   0 0.5   0 ...
+ 	0   1   0   0   0   0   0 0.5   0   0   1   0   0   0   0 ...
+ 	0 0.5  0   0   1   0   0   0   0   0   1   0   0   0   0 ...
+ 	0   0   0   0   1   0   0   0   0   0   0   0   0   1   0 ...
+ 	0   0   0   0   0   0   0]';
+ 
+ assert (pkl(:),pkl2);
+1 test, 1 passed, 0 known failure, 0 skipped
+[src/curvederivcpts.cc]
+>>>>> /build/octave-nurbs-1.3.13/src/curvederivcpts.cc
 ***** 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
+ 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
 Summary: 123 tests, 123 passed, 0 known failures, 0 skipped
    dh_octave_substvar -O--buildsystem=octave
    dh_installdocs -O--buildsystem=octave
@@ -4009,8 +4043,8 @@
    dh_strip -O--buildsystem=octave
    dh_makeshlibs -O--buildsystem=octave
    dh_shlibdeps -O--buildsystem=octave
-dpkg-shlibdeps: warning: package could avoid a useless dependency if debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/surfderiveval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/tbasisfun.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrb_srf_basisfun__.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/basisfun.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/surfderivcpts.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/bspeval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/basisfunder.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrb_srf_basisfun_der__.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/bspderiv.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrbsurfderiveval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/curvederivcpts.oct were not linked against libpthread.so.0 (they use none of the library's symbols)
-dpkg-shlibdeps: warning: package could avoid a useless dependency if debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/surfderiveval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/tbasisfun.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrb_srf_basisfun__.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/basisfun.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/surfderivcpts.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/bspeval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/basisfunder.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrb_srf_basisfun_der__.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/bspderiv.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrbsurfderiveval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/curvederivcpts.oct were not linked against libgomp.so.1 (they use none of the library's symbols)
+dpkg-shlibdeps: warning: package could avoid a useless dependency if debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrbsurfderiveval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrb_srf_basisfun__.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/surfderiveval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/curvederivcpts.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrb_srf_basisfun_der__.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/surfderivcpts.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/basisfunder.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/tbasisfun.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/bspeval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/bspderiv.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/basisfun.oct were not linked against libgomp.so.1 (they use none of the library's symbols)
+dpkg-shlibdeps: warning: package could avoid a useless dependency if debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrbsurfderiveval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrb_srf_basisfun__.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/surfderiveval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/curvederivcpts.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/nrb_srf_basisfun_der__.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/surfderivcpts.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/basisfunder.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/tbasisfun.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/bspeval.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/bspderiv.oct debian/octave-nurbs/usr/lib/arm-linux-gnueabihf/octave/packages/nurbs-1.3.13/arm-unknown-linux-gnueabihf-api-v52/basisfun.oct were not linked against libpthread.so.0 (they use none of the library's symbols)
    dh_installdeb -O--buildsystem=octave
    dh_gencontrol -O--buildsystem=octave
 dpkg-gencontrol: warning: package octave-nurbs: substitution variable ${octave:Upstream-Description} unused, but is defined
@@ -4025,12 +4059,14 @@
  dpkg-source --after-build .
 dpkg-buildpackage: info: binary-only upload (no source included)
 I: copying local configuration
+I: user script /srv/workspace/pbuilder/7978/tmp/hooks/B01_cleanup starting
+I: user script /srv/workspace/pbuilder/7978/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/27930 and its subdirectories
-I: Current time: Wed Aug 12 00:29:32 -12 2020
-I: pbuilder-time-stamp: 1597235372
+I: removing directory /srv/workspace/pbuilder/7978 and its subdirectories
+I: Current time: Thu Aug 13 02:36:08 +14 2020
+I: pbuilder-time-stamp: 1597235768