Diff of the two buildlogs:

--
--- b1/build.log	2021-08-10 13:20:12.408347505 +0000
+++ b2/build.log	2021-08-10 13:26:25.534479532 +0000
@@ -1,6 +1,6 @@
 I: pbuilder: network access will be disabled during build
-I: Current time: Tue Aug 10 01:16:05 -12 2021
-I: pbuilder-time-stamp: 1628601365
+I: Current time: Wed Aug 11 03:20:22 +14 2021
+I: pbuilder-time-stamp: 1628601622
 I: Building the build Environment
 I: extracting base tarball [/var/cache/pbuilder/bullseye-reproducible-base.tgz]
 I: copying local configuration
@@ -16,8 +16,8 @@
 I: copying [./octave-bim_1.1.5-7.debian.tar.xz]
 I: Extracting source
 gpgv: unknown type of key resource 'trustedkeys.kbx'
-gpgv: keyblock resource '/tmp/dpkg-verify-sig.Lw92ehy7/trustedkeys.kbx': General error
-gpgv: Signature made Sat Jul 25 20:41:50 2020 -12
+gpgv: keyblock resource '/tmp/dpkg-verify-sig.uBeFTYFQ/trustedkeys.kbx': General error
+gpgv: Signature made Sun Jul 26 22:41:50 2020 +14
 gpgv:                using RSA key 3F464391498FE874BDB5D98F2124AA1983785C90
 gpgv: Can't check signature: No public key
 dpkg-source: warning: failed to verify signature on ./octave-bim_1.1.5-7.dsc
@@ -30,135 +30,169 @@
 dpkg-source: info: applying xtest-in-bim1c_norm.patch
 I: Not using root during the build.
 I: Installing the build-deps
-I: user script /srv/workspace/pbuilder/12574/tmp/hooks/D02_print_environment starting
+I: user script /srv/workspace/pbuilder/28350/tmp/hooks/D01_modify_environment starting
+debug: Running on virt32b.
+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/28350/tmp/hooks/D01_modify_environment finished
+I: user script /srv/workspace/pbuilder/28350/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,-fixfilepath 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]="1" [2]="4" [3]="1" [4]="release" [5]="arm-unknown-linux-gnueabihf")
+  BASH_VERSION='5.1.4(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,-fixfilepath 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='a438d16796f241f6be609265a6696a6f'
-  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='12574'
-  PS1='# '
-  PS2='> '
+  INVOCATION_ID=f1f8be207f904255bb3ce6e4ff172d54
+  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=28350
   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.QZCxaR9SSf/pbuilderrc_fnJp --hookdir /etc/pbuilder/first-build-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/bullseye-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/tmp.QZCxaR9SSf/b1 --logfile b1/build.log octave-bim_1.1.5-7.dsc'
-  SUDO_GID='113'
-  SUDO_UID='107'
-  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.QZCxaR9SSf/pbuilderrc_9uue --hookdir /etc/pbuilder/rebuild-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/bullseye-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/tmp.QZCxaR9SSf/b2 --logfile b2/build.log --extrapackages usrmerge octave-bim_1.1.5-7.dsc'
+  SUDO_GID=112
+  SUDO_UID=106
+  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 virt64b 5.10.0-8-arm64 #1 SMP Debian 5.10.46-4 (2021-08-03) aarch64 GNU/Linux
+  Linux i-capture-the-hostname 5.10.0-8-armmp-lpae #1 SMP Debian 5.10.46-4 (2021-08-03) armv7l GNU/Linux
 I: ls -l /bin
   total 3580
-  -rwxr-xr-x 1 root root 816764 Aug  4 08:25 bash
-  -rwxr-xr-x 3 root root  26052 Jul 20  2020 bunzip2
-  -rwxr-xr-x 3 root root  26052 Jul 20  2020 bzcat
-  lrwxrwxrwx 1 root root      6 Jul 20  2020 bzcmp -> bzdiff
-  -rwxr-xr-x 1 root root   2225 Jul 20  2020 bzdiff
-  lrwxrwxrwx 1 root root      6 Jul 20  2020 bzegrep -> bzgrep
-  -rwxr-xr-x 1 root root   4877 Sep  4  2019 bzexe
-  lrwxrwxrwx 1 root root      6 Jul 20  2020 bzfgrep -> bzgrep
-  -rwxr-xr-x 1 root root   3775 Jul 20  2020 bzgrep
-  -rwxr-xr-x 3 root root  26052 Jul 20  2020 bzip2
-  -rwxr-xr-x 1 root root   9636 Jul 20  2020 bzip2recover
-  lrwxrwxrwx 1 root root      6 Jul 20  2020 bzless -> bzmore
-  -rwxr-xr-x 1 root root   1297 Jul 20  2020 bzmore
-  -rwxr-xr-x 1 root root  26668 Sep 22  2020 cat
-  -rwxr-xr-x 1 root root  43104 Sep 22  2020 chgrp
-  -rwxr-xr-x 1 root root  38984 Sep 22  2020 chmod
-  -rwxr-xr-x 1 root root  43112 Sep 22  2020 chown
-  -rwxr-xr-x 1 root root  92616 Sep 22  2020 cp
-  -rwxr-xr-x 1 root root  75524 Dec 10  2020 dash
-  -rwxr-xr-x 1 root root  75880 Sep 22  2020 date
-  -rwxr-xr-x 1 root root  55436 Sep 22  2020 dd
-  -rwxr-xr-x 1 root root  59912 Sep 22  2020 df
-  -rwxr-xr-x 1 root root  96764 Sep 22  2020 dir
-  -rwxr-xr-x 1 root root  55012 Jul 28 07:09 dmesg
-  lrwxrwxrwx 1 root root      8 Nov  6  2019 dnsdomainname -> hostname
-  lrwxrwxrwx 1 root root      8 Nov  6  2019 domainname -> hostname
-  -rwxr-xr-x 1 root root  22508 Sep 22  2020 echo
-  -rwxr-xr-x 1 root root     28 Nov  9  2020 egrep
-  -rwxr-xr-x 1 root root  22496 Sep 22  2020 false
-  -rwxr-xr-x 1 root root     28 Nov  9  2020 fgrep
-  -rwxr-xr-x 1 root root  47492 Jul 28 07:09 findmnt
-  -rwsr-xr-x 1 root root  26076 Feb 26 04:12 fusermount
-  -rwxr-xr-x 1 root root 124508 Nov  9  2020 grep
-  -rwxr-xr-x 2 root root   2346 Mar  2 11:30 gunzip
-  -rwxr-xr-x 1 root root   6376 Mar  2 11:30 gzexe
-  -rwxr-xr-x 1 root root  64212 Mar  2 11:30 gzip
-  -rwxr-xr-x 1 root root  13784 Nov  6  2019 hostname
-  -rwxr-xr-x 1 root root  43180 Sep 22  2020 ln
-  -rwxr-xr-x 1 root root  35068 Feb  7  2020 login
-  -rwxr-xr-x 1 root root  96764 Sep 22  2020 ls
-  -rwxr-xr-x 1 root root  99940 Jul 28 07:09 lsblk
-  -rwxr-xr-x 1 root root  51408 Sep 22  2020 mkdir
-  -rwxr-xr-x 1 root root  43184 Sep 22  2020 mknod
-  -rwxr-xr-x 1 root root  30780 Sep 22  2020 mktemp
-  -rwxr-xr-x 1 root root  34408 Jul 28 07:09 more
-  -rwsr-xr-x 1 root root  34400 Jul 28 07:09 mount
-  -rwxr-xr-x 1 root root   9824 Jul 28 07:09 mountpoint
-  -rwxr-xr-x 1 root root  88524 Sep 22  2020 mv
-  lrwxrwxrwx 1 root root      8 Nov  6  2019 nisdomainname -> hostname
-  lrwxrwxrwx 1 root root     14 Apr 18 03:38 pidof -> /sbin/killall5
-  -rwxr-xr-x 1 root root  26652 Sep 22  2020 pwd
-  lrwxrwxrwx 1 root root      4 Aug  4 08:25 rbash -> bash
-  -rwxr-xr-x 1 root root  30740 Sep 22  2020 readlink
-  -rwxr-xr-x 1 root root  43104 Sep 22  2020 rm
-  -rwxr-xr-x 1 root root  30732 Sep 22  2020 rmdir
-  -rwxr-xr-x 1 root root  14144 Sep 27  2020 run-parts
-  -rwxr-xr-x 1 root root  76012 Dec 22  2018 sed
-  lrwxrwxrwx 1 root root      4 Aug  8 21:24 sh -> dash
-  -rwxr-xr-x 1 root root  22532 Sep 22  2020 sleep
-  -rwxr-xr-x 1 root root  55360 Sep 22  2020 stty
-  -rwsr-xr-x 1 root root  46704 Jul 28 07:09 su
-  -rwxr-xr-x 1 root root  22532 Sep 22  2020 sync
-  -rwxr-xr-x 1 root root 340872 Feb 16 21:55 tar
-  -rwxr-xr-x 1 root root   9808 Sep 27  2020 tempfile
-  -rwxr-xr-x 1 root root  67696 Sep 22  2020 touch
-  -rwxr-xr-x 1 root root  22496 Sep 22  2020 true
-  -rwxr-xr-x 1 root root   9636 Feb 26 04:12 ulockmgr_server
-  -rwsr-xr-x 1 root root  22108 Jul 28 07:09 umount
-  -rwxr-xr-x 1 root root  22520 Sep 22  2020 uname
-  -rwxr-xr-x 2 root root   2346 Mar  2 11:30 uncompress
-  -rwxr-xr-x 1 root root  96764 Sep 22  2020 vdir
-  -rwxr-xr-x 1 root root  38512 Jul 28 07:09 wdctl
-  lrwxrwxrwx 1 root root      8 Nov  6  2019 ypdomainname -> hostname
-  -rwxr-xr-x 1 root root   1984 Mar  2 11:30 zcat
-  -rwxr-xr-x 1 root root   1678 Mar  2 11:30 zcmp
-  -rwxr-xr-x 1 root root   5880 Mar  2 11:30 zdiff
-  -rwxr-xr-x 1 root root     29 Mar  2 11:30 zegrep
-  -rwxr-xr-x 1 root root     29 Mar  2 11:30 zfgrep
-  -rwxr-xr-x 1 root root   2081 Mar  2 11:30 zforce
-  -rwxr-xr-x 1 root root   7585 Mar  2 11:30 zgrep
-  -rwxr-xr-x 1 root root   2206 Mar  2 11:30 zless
-  -rwxr-xr-x 1 root root   1842 Mar  2 11:30 zmore
-  -rwxr-xr-x 1 root root   4553 Mar  2 11:30 znew
-I: user script /srv/workspace/pbuilder/12574/tmp/hooks/D02_print_environment finished
+  -rwxr-xr-x 1 root root 816764 Aug  5 10:25 bash
+  -rwxr-xr-x 3 root root  26052 Jul 21  2020 bunzip2
+  -rwxr-xr-x 3 root root  26052 Jul 21  2020 bzcat
+  lrwxrwxrwx 1 root root      6 Jul 21  2020 bzcmp -> bzdiff
+  -rwxr-xr-x 1 root root   2225 Jul 21  2020 bzdiff
+  lrwxrwxrwx 1 root root      6 Jul 21  2020 bzegrep -> bzgrep
+  -rwxr-xr-x 1 root root   4877 Sep  5  2019 bzexe
+  lrwxrwxrwx 1 root root      6 Jul 21  2020 bzfgrep -> bzgrep
+  -rwxr-xr-x 1 root root   3775 Jul 21  2020 bzgrep
+  -rwxr-xr-x 3 root root  26052 Jul 21  2020 bzip2
+  -rwxr-xr-x 1 root root   9636 Jul 21  2020 bzip2recover
+  lrwxrwxrwx 1 root root      6 Jul 21  2020 bzless -> bzmore
+  -rwxr-xr-x 1 root root   1297 Jul 21  2020 bzmore
+  -rwxr-xr-x 1 root root  26668 Sep 23  2020 cat
+  -rwxr-xr-x 1 root root  43104 Sep 23  2020 chgrp
+  -rwxr-xr-x 1 root root  38984 Sep 23  2020 chmod
+  -rwxr-xr-x 1 root root  43112 Sep 23  2020 chown
+  -rwxr-xr-x 1 root root  92616 Sep 23  2020 cp
+  -rwxr-xr-x 1 root root  75524 Dec 11  2020 dash
+  -rwxr-xr-x 1 root root  75880 Sep 23  2020 date
+  -rwxr-xr-x 1 root root  55436 Sep 23  2020 dd
+  -rwxr-xr-x 1 root root  59912 Sep 23  2020 df
+  -rwxr-xr-x 1 root root  96764 Sep 23  2020 dir
+  -rwxr-xr-x 1 root root  55012 Jul 29 09:09 dmesg
+  lrwxrwxrwx 1 root root      8 Nov  8  2019 dnsdomainname -> hostname
+  lrwxrwxrwx 1 root root      8 Nov  8  2019 domainname -> hostname
+  -rwxr-xr-x 1 root root  22508 Sep 23  2020 echo
+  -rwxr-xr-x 1 root root     28 Nov 10  2020 egrep
+  -rwxr-xr-x 1 root root  22496 Sep 23  2020 false
+  -rwxr-xr-x 1 root root     28 Nov 10  2020 fgrep
+  -rwxr-xr-x 1 root root  47492 Jul 29 09:09 findmnt
+  -rwsr-xr-x 1 root root  26076 Feb 27 06:12 fusermount
+  -rwxr-xr-x 1 root root 124508 Nov 10  2020 grep
+  -rwxr-xr-x 2 root root   2346 Mar  3 13:30 gunzip
+  -rwxr-xr-x 1 root root   6376 Mar  3 13:30 gzexe
+  -rwxr-xr-x 1 root root  64212 Mar  3 13:30 gzip
+  -rwxr-xr-x 1 root root  13784 Nov  8  2019 hostname
+  -rwxr-xr-x 1 root root  43180 Sep 23  2020 ln
+  -rwxr-xr-x 1 root root  35068 Feb  8  2020 login
+  -rwxr-xr-x 1 root root  96764 Sep 23  2020 ls
+  -rwxr-xr-x 1 root root  99940 Jul 29 09:09 lsblk
+  -rwxr-xr-x 1 root root  51408 Sep 23  2020 mkdir
+  -rwxr-xr-x 1 root root  43184 Sep 23  2020 mknod
+  -rwxr-xr-x 1 root root  30780 Sep 23  2020 mktemp
+  -rwxr-xr-x 1 root root  34408 Jul 29 09:09 more
+  -rwsr-xr-x 1 root root  34400 Jul 29 09:09 mount
+  -rwxr-xr-x 1 root root   9824 Jul 29 09:09 mountpoint
+  -rwxr-xr-x 1 root root  88524 Sep 23  2020 mv
+  lrwxrwxrwx 1 root root      8 Nov  8  2019 nisdomainname -> hostname
+  lrwxrwxrwx 1 root root     14 Apr 19 05:38 pidof -> /sbin/killall5
+  -rwxr-xr-x 1 root root  26652 Sep 23  2020 pwd
+  lrwxrwxrwx 1 root root      4 Aug  5 10:25 rbash -> bash
+  -rwxr-xr-x 1 root root  30740 Sep 23  2020 readlink
+  -rwxr-xr-x 1 root root  43104 Sep 23  2020 rm
+  -rwxr-xr-x 1 root root  30732 Sep 23  2020 rmdir
+  -rwxr-xr-x 1 root root  14144 Sep 28  2020 run-parts
+  -rwxr-xr-x 1 root root  76012 Dec 23  2018 sed
+  lrwxrwxrwx 1 root root      4 Aug 11 03:20 sh -> bash
+  lrwxrwxrwx 1 root root      4 Aug  9 23:25 sh.distrib -> dash
+  -rwxr-xr-x 1 root root  22532 Sep 23  2020 sleep
+  -rwxr-xr-x 1 root root  55360 Sep 23  2020 stty
+  -rwsr-xr-x 1 root root  46704 Jul 29 09:09 su
+  -rwxr-xr-x 1 root root  22532 Sep 23  2020 sync
+  -rwxr-xr-x 1 root root 340872 Feb 17 23:55 tar
+  -rwxr-xr-x 1 root root   9808 Sep 28  2020 tempfile
+  -rwxr-xr-x 1 root root  67696 Sep 23  2020 touch
+  -rwxr-xr-x 1 root root  22496 Sep 23  2020 true
+  -rwxr-xr-x 1 root root   9636 Feb 27 06:12 ulockmgr_server
+  -rwsr-xr-x 1 root root  22108 Jul 29 09:09 umount
+  -rwxr-xr-x 1 root root  22520 Sep 23  2020 uname
+  -rwxr-xr-x 2 root root   2346 Mar  3 13:30 uncompress
+  -rwxr-xr-x 1 root root  96764 Sep 23  2020 vdir
+  -rwxr-xr-x 1 root root  38512 Jul 29 09:09 wdctl
+  lrwxrwxrwx 1 root root      8 Nov  8  2019 ypdomainname -> hostname
+  -rwxr-xr-x 1 root root   1984 Mar  3 13:30 zcat
+  -rwxr-xr-x 1 root root   1678 Mar  3 13:30 zcmp
+  -rwxr-xr-x 1 root root   5880 Mar  3 13:30 zdiff
+  -rwxr-xr-x 1 root root     29 Mar  3 13:30 zegrep
+  -rwxr-xr-x 1 root root     29 Mar  3 13:30 zfgrep
+  -rwxr-xr-x 1 root root   2081 Mar  3 13:30 zforce
+  -rwxr-xr-x 1 root root   7585 Mar  3 13:30 zgrep
+  -rwxr-xr-x 1 root root   2206 Mar  3 13:30 zless
+  -rwxr-xr-x 1 root root   1842 Mar  3 13:30 zmore
+  -rwxr-xr-x 1 root root   4553 Mar  3 13:30 znew
+I: user script /srv/workspace/pbuilder/28350/tmp/hooks/D02_print_environment finished
  -> Attempting to satisfy build-dependencies
  -> Creating pbuilder-satisfydepends-dummy package
 Package: pbuilder-satisfydepends-dummy
@@ -702,7 +736,7 @@
 Get: 500 http://deb.debian.org/debian bullseye/main armhf gmsh armhf 4.7.1+ds1-5 [46.0 kB]
 Get: 501 http://deb.debian.org/debian bullseye/main armhf octave-splines all 1.3.3-2 [22.3 kB]
 Get: 502 http://deb.debian.org/debian bullseye/main armhf octave-msh armhf 1.0.10-11 [32.1 kB]
-Fetched 168 MB in 18s (9142 kB/s)
+Fetched 168 MB in 18s (9543 kB/s)
 debconf: delaying package configuration, since apt-utils is not installed
 Selecting previously unselected package bsdextrautils.
 (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 ... 19398 files and directories currently installed.)
@@ -2763,8 +2797,26 @@
 Writing extended state information...
 Building tag database...
  -> Finished parsing the build-deps
+Reading package lists...
+Building dependency tree...
+Reading state information...
+The following NEW packages will be installed:
+  usrmerge
+0 upgraded, 1 newly installed, 0 to remove and 0 not upgraded.
+Need to get 13.0 kB of archives.
+After this operation, 43.0 kB of additional disk space will be used.
+Get:1 http://deb.debian.org/debian bullseye/main armhf usrmerge all 25 [13.0 kB]
+debconf: delaying package configuration, since apt-utils is not installed
+Fetched 13.0 kB in 0s (485 kB/s)
+Selecting previously unselected package usrmerge.
+(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 ... 38873 files and directories currently installed.)
+Preparing to unpack .../archives/usrmerge_25_all.deb ...
+Unpacking usrmerge (25) ...
+Setting up usrmerge (25) ...
+The system has been successfully converted.
 I: Building the package
-I: Running cd /build/octave-bim-1.1.5/ && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games" HOME="/nonexistent/first-build" dpkg-buildpackage -us -uc -b && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games" HOME="/nonexistent/first-build" dpkg-genchanges -S  > ../octave-bim_1.1.5-7_source.changes
+hostname: Name or service not known
+I: Running cd /build/octave-bim-1.1.5/ && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path" HOME="/nonexistent/second-build" dpkg-buildpackage -us -uc -b && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path" HOME="/nonexistent/second-build" dpkg-genchanges -S  > ../octave-bim_1.1.5-7_source.changes
 dpkg-buildpackage: info: source package octave-bim
 dpkg-buildpackage: info: source version 1.1.5-7
 dpkg-buildpackage: info: source distribution unstable
@@ -2805,101 +2857,297 @@
    dh_octave_check -O--buildsystem=octave
 Checking package...
 Checking m files ...
-[inst/bim2a_axisymmetric_reaction.m]
->>>>> /build/octave-bim-1.1.5/inst/bim2a_axisymmetric_reaction.m
+[inst/bim2c_intrp.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim2c_intrp.m
+***** test
+ msh = bim2c_mesh_properties (msh2m_structured_mesh (linspace (0, 1, 11), linspace (0, 1, 13), 1, 1:4));
+ x = y = linspace (0, 1, 100).';
+ u = msh.p(1, :).';
+ ui = bim2c_intrp (msh, u, [], [x, y]); 
+ assert (ui, linspace (0, 1, 100), 10*eps);
+warning: 'matrix' object indexed with empty index list
+warning: called from
+    msh2m_geometrical_properties>computearea at line 406 column 7
+    msh2m_geometrical_properties at line 167 column 13
+    bim2c_mesh_properties at line 44 column 40
+    __test__ at line 3 column 6
+    test at line 677 column 11
+    /tmp/tmp.RoVIMMmIw0 at line 14 column 31
+
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/bim3a_reaction.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim3a_reaction.m
 ***** shared mesh,delta,zeta,nnodes,nelem
 ***** test
 ***** test
 ***** test
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/bim3c_intrp.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim3c_intrp.m
 ***** test
- n = 20;
- [mesh] = msh2m_structured_mesh(linspace(1,2,n+1),linspace(0,1,n+1),1,1:4);
+ msh = bim3c_mesh_properties (msh3m_structured_mesh (linspace (0, 1, 11), linspace (0, 1, 9), linspace (0, 1, 13), 1, 1:6));
+ x = y = z = linspace (0, 1, 100).';
+ u = msh.p(1, :).';
+ ui = bim3c_intrp (msh, u, [], [x, y, z]); 
+ assert (ui, linspace (0, 1, 100), 10*eps);
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/bim3c_norm.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim3c_norm.m
+***** test
+***** shared L, V, x, y, z, m
+ L = rand (1); V = rand (1); x = linspace (0,L,4);  y = x;  z = x;
+ m = msh3m_structured_mesh (x,y,z,1,1:6);
+ m.area = msh3m_geometrical_properties (m, 'area');
+ m.shg  = msh3m_geometrical_properties (m, 'shg');
+ u    = V * ones (columns(m.p),1);
+ uinf = bim3c_norm (m, u, 'inf');
+ uL2  = bim3c_norm (m, u, 'L2');
+ uH1  = bim3c_norm (m, u, 'H1');
+ assert ([uinf, uL2, uH1], [V, V*sqrt(L^3), V*sqrt(L^3)], 1e-12);
+***** test
+ u    = V * (m.p(1,:) + 2*m.p(2,:) + 3*m.p(3,:))';
+ uinf = bim3c_norm (m, u, 'inf');
+ uL2  = bim3c_norm (m, u, 'L2');
+ uH1  = bim3c_norm (m, u, 'H1');
+ assert ([uinf, uL2, uH1], 
+         [6*L*V, V*sqrt(61/6*L^5), V*sqrt(61/6*L^5 + 14*L^3)],
+          1e-12);
+***** test
+ u    = V * ones (columns(m.t),1);
+ uinf = bim3c_norm (m, u, 'inf');
+ uL2  = bim3c_norm (m, u, 'L2');
+ assert ([uinf, uL2], [V, V*sqrt(L^3)], 1e-12);
+***** test
+ u     = V * ones (columns(m.t),1);
+ uvect = [u, 2*u, 3*u];
+ uinf  = bim3c_norm (m, uvect, 'inf');
+ uL2   = bim3c_norm (m, uvect, 'L2');
+ assert ([uinf, uL2], [3*V, V*sqrt(14*L^3)], 1e-12);
+4 tests, 4 passed, 0 known failure, 0 skipped
+[inst/bim2a_advection_diffusion.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim2a_advection_diffusion.m
+***** test 
+ [mesh] = msh2m_structured_mesh([0:1/3:1],[0:1/3:1],1,1:4);
  mesh   = bim2c_mesh_properties(mesh);
- uex    = @(r,z) exp(r) .* exp(1-z);
- duexdr = @(r,z) uex(r,z);
- d2uexdr2 = @(r,z) uex(r,z);
- duexdz = @(r,z) -uex(r,z);
- d2uexdz2 = @(r,z) uex(r,z);
- Dnodes = bim2c_unknowns_on_side(mesh,[1,2,3,4]);
+ x      = mesh.p(1,:)';
+ Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
  Nnodes    = columns(mesh.p);
  Nelements = columns(mesh.t);
  Varnodes  = setdiff(1:Nnodes,Dnodes);
- D = 1; vr = 1; vz = 1; sigma = 1;
- alpha  = D*ones(Nelements,1);
+ alpha  = ones(Nelements,1);
+ eta    = .1*ones(Nnodes,1);
+ beta   = [ones(1,Nelements);zeros(1,Nelements)];
  gamma  = ones(Nnodes,1);
- eta    = ones(Nnodes,1);
- delta  = sigma*ones(columns(mesh.t),1);
- zeta   = ones(columns(mesh.p),1);
- beta   = 1/D*[vr*ones(1,Nelements); vz*ones(1,Nelements)];
- f = @(r,z) -D./r.*duexdr(r,z) - D.*d2uexdr2(r,z) ...
-           + vr./r .* uex(r,z) + vr * duexdr(r,z) ...
-           - D.*d2uexdz2(r,z) + vz * duexdz(r,z) ...
-           + sigma * uex(r,z);
- rhs    = bim2a_axisymmetric_rhs(mesh, ones(Nelements,1), f(mesh.p(1,:), mesh.p(2,:)));
- S = bim2a_axisymmetric_advection_diffusion(mesh,alpha,gamma,eta,beta);
- C = bim2a_axisymmetric_reaction(mesh,delta,zeta);
- S += C;
- u = zeros(Nnodes,1); u(Dnodes) = uex(mesh.p(1,Dnodes), mesh.p(2,Dnodes));
- u(Varnodes) = S(Varnodes,Varnodes)\(rhs(Varnodes) - S(Varnodes,Dnodes)*u(Dnodes));
- assert(u,uex(mesh.p(1,:), mesh.p(2,:))',1e-3)
+ f      = bim2a_rhs(mesh,ones(Nelements,1),ones(Nnodes,1));
+ S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ u = zeros(Nnodes,1);
+ u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
+ uex = x - (exp(10*x)-1)/(exp(10)-1);
+ assert(u,uex,1e-7)
 warning: 'matrix' object indexed with empty index list
 warning: called from
     msh2m_geometrical_properties>computearea at line 406 column 7
     msh2m_geometrical_properties at line 167 column 13
     bim2c_mesh_properties at line 44 column 40
-    __test__ at line 5 column 9
+    __test__ at line 4 column 9
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 14 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 46 column 31
+
+***** test 
+ [mesh] = msh2m_structured_mesh([0:1/3:1],[0:1/3:1],1,1:4);
+ mesh   = bim2c_mesh_properties(mesh);
+ x      = mesh.p(1,:)';
+ Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
+ Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
+ Varnodes = setdiff(1:Nnodes,Dnodes);
+ alpha  = ones(Nelements,1);
+ eta    = .1*ones(Nnodes,1);
+ beta   = x;
+ gamma  = ones(Nnodes,1);
+ f      = bim2a_rhs(mesh,ones(Nelements,1),ones(Nnodes,1));
+ S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ u = zeros(Nnodes,1);
+ u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
+ uex = x - (exp(10*x)-1)/(exp(10)-1);
+ assert(u,uex,1e-7)
+warning: 'matrix' object indexed with empty index list
+warning: called from
+    msh2m_geometrical_properties>computearea at line 406 column 7
+    msh2m_geometrical_properties at line 167 column 13
+    bim2c_mesh_properties at line 44 column 40
+    __test__ at line 4 column 9
+    test at line 677 column 11
+    /tmp/tmp.RoVIMMmIw0 at line 46 column 31
 
-4 tests, 4 passed, 0 known failure, 0 skipped
-[inst/bim1c_norm.m]
->>>>> /build/octave-bim-1.1.5/inst/bim1c_norm.m
-***** xtest
-***** shared L, V, m
- L = rand (1); V = rand (1); m = linspace (0,1,5).^2;  m *= L;
- u    = V * ones (size (m))';
- uinf = bim1c_norm (m, u, 'inf');
- uL2  = bim1c_norm (m, u, 'L2');
- uH1  = bim1c_norm (m, u, 'H1');
- assert ([uinf, uL2, uH1], [V, V*sqrt(L), V*sqrt(L)], 1e-12);
 ***** test
- u    = V * m';
- uinf = bim1c_norm (m, u, 'inf');
- uL2  = bim1c_norm (m, u, 'L2');
- uH1  = bim1c_norm (m, u, 'H1');
- assert ([uinf, uL2, uH1], 
-         [L*V, V*sqrt(L^3/3), V*sqrt(L^3/3 + L)],
-          1e-12);
+ [mesh] = msh2m_structured_mesh([0:1/3:1],[0:1/3:1],1,1:4);
+ mesh   = bim2c_mesh_properties(mesh);
+ x      = mesh.p(1,:)';
+ Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
+ Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
+ Varnodes = setdiff(1:Nnodes,Dnodes);
+ alpha  = 10*ones(Nelements,1);
+ eta    = .01*ones(Nnodes,1);
+ beta   = x/10;
+ gamma  = ones(Nnodes,1);
+ f      = bim2a_rhs(mesh,ones(Nelements,1),ones(Nnodes,1));
+ S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ u = zeros(Nnodes,1);
+ u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
+ uex = x - (exp(10*x)-1)/(exp(10)-1);
+ assert(u,uex,1e-7)
+warning: 'matrix' object indexed with empty index list
+warning: called from
+    msh2m_geometrical_properties>computearea at line 406 column 7
+    msh2m_geometrical_properties at line 167 column 13
+    bim2c_mesh_properties at line 44 column 40
+    __test__ at line 4 column 9
+    test at line 677 column 11
+    /tmp/tmp.RoVIMMmIw0 at line 46 column 31
+
 ***** test
- u    = V * ones (size (diff (m)))';
- uinf = bim1c_norm (m, u, 'inf');
- uL2  = bim1c_norm (m, u, 'L2');
- assert ([uinf, uL2], [V, V*sqrt(L)], 1e-12);
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/bim1a_advection_upwind.m]
->>>>> /build/octave-bim-1.1.5/inst/bim1a_advection_upwind.m
+ [mesh] = msh2m_structured_mesh([0:1/3:1],[0:1/3:1],1,1:4);
+ mesh   = bim2c_mesh_properties(mesh);
+ x      = mesh.p(1,:)';
+ Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
+ Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
+ Varnodes = setdiff(1:Nnodes,Dnodes);
+ alpha  = 10*ones(Nelements,1); eta = .001*ones(Nnodes,1);
+ beta   = x/100;
+ gamma  = 10*ones(Nnodes,1);
+ f      = bim2a_rhs(mesh,ones(Nelements,1),ones(Nnodes,1));
+ S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ u = zeros(Nnodes,1);
+ u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
+ uex = x - (exp(10*x)-1)/(exp(10)-1);
+ assert(u,uex,1e-7)
+warning: 'matrix' object indexed with empty index list
+warning: called from
+    msh2m_geometrical_properties>computearea at line 406 column 7
+    msh2m_geometrical_properties at line 167 column 13
+    bim2c_mesh_properties at line 44 column 40
+    __test__ at line 4 column 9
+    test at line 677 column 11
+    /tmp/tmp.RoVIMMmIw0 at line 46 column 31
+
 ***** test
- n = 200;
- mesh      = linspace(0,1,n+1)';
- uex       = @(r) - r.^2 + 1;
- Nnodes    = numel(mesh);
- Nelements = Nnodes-1;
- D = 1;  v = 1;  sigma = 0;
- alpha  = D*ones(Nelements,1);
+ [mesh] = msh2m_structured_mesh([0:1/1e3:1],[0:1/2:1],1,1:4);
+ mesh   = bim2c_mesh_properties(mesh);
+ x      = mesh.p(1,:)';
+ Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
+ Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
+ Varnodes = setdiff(1:Nnodes,Dnodes);
+ alpha  = 3*ones(Nelements,1); eta = x+1;
+ beta   = [ones(1,Nelements);zeros(1,Nelements)];
+ gamma  = 2*x;
+ ff     = 2*(6*x.^2+6*x) - (6*x+6).*(1-2*x)+6*(x-x.^2);
+ f      = bim2a_rhs(mesh,ones(Nelements,1),ff);
+ S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ u = zeros(Nnodes,1);
+ u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
+ uex = x - x.^2;
+ assert(u,uex,5e-3)
+warning: 'matrix' object indexed with empty index list
+warning: called from
+    msh2m_geometrical_properties>computearea at line 406 column 7
+    msh2m_geometrical_properties at line 167 column 13
+    bim2c_mesh_properties at line 44 column 40
+    __test__ at line 4 column 9
+    test at line 677 column 11
+    /tmp/tmp.RoVIMMmIw0 at line 46 column 31
+
+***** test
+ [mesh] = msh2m_structured_mesh([0:1/1e3:1],[0:1/2:1],1,1:4);
+ mesh = bim2c_mesh_properties(mesh);
+ x = mesh.p(1,:)';
+ Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
+ Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
+ Varnodes = setdiff(1:Nnodes,Dnodes);
+ alpha  = ones(Nelements,1); eta = ones(Nnodes,1);
+ beta   = 0;
+ gamma  = x+1;
+ ff     = 4*x+1;
+ f      = bim2a_rhs(mesh,ones(Nelements,1),ff);
+ S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ u = zeros(Nnodes,1);
+ u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
+ uex = x - x.^2; 
+ assert(u,uex,1e-7)
+warning: 'matrix' object indexed with empty index list
+warning: called from
+    msh2m_geometrical_properties>computearea at line 406 column 7
+    msh2m_geometrical_properties at line 167 column 13
+    bim2c_mesh_properties at line 44 column 40
+    __test__ at line 4 column 7
+    test at line 677 column 11
+    /tmp/tmp.RoVIMMmIw0 at line 46 column 31
+
+***** test
+ [mesh] = msh2m_structured_mesh([0:.1:1],[0:.1:1],1,1:4);
+ mesh = bim2c_mesh_properties(mesh);
+ x = mesh.p(1,:)';y = mesh.p(2,:)';
+ Dnodes = bim2c_unknowns_on_side(mesh,[1:4]);
+ Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
+ Varnodes = setdiff(1:Nnodes,Dnodes);
+ alpha  = ones(Nelements,1); diff = 1e-2; eta=diff*ones(Nnodes,1);
+ beta   =[ones(1,Nelements);ones(1,Nelements)];
+ gamma  = x*0+1;
+ ux  = y.*(1-exp((y-1)/diff)) .* (1-exp((x-1)/diff)-x.*exp((x-1)/diff)/diff);
+ uy  = x.*(1-exp((x-1)/diff)) .* (1-exp((y-1)/diff)-y.*exp((y-1)/diff)/diff);
+ uxx = y.*(1-exp((y-1)/diff)) .* (-2*exp((x-1)/diff)/diff-x.*exp((x-1)/diff)/(diff^2));
+ uyy = x.*(1-exp((x-1)/diff)) .* (-2*exp((y-1)/diff)/diff-y.*exp((y-1)/diff)/(diff^2));
+ ff  = -diff*(uxx+uyy)+ux+uy;
+ f   = bim2a_rhs(mesh,ones(Nelements,1),ff);
+ S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ u = zeros(Nnodes,1);
+ u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
+ uex = x.*y.*(1-exp((x-1)/diff)).*(1-exp((y-1)/diff)); 
+ assert(u,uex,1e-7)
+warning: 'matrix' object indexed with empty index list
+warning: called from
+    msh2m_geometrical_properties>computearea at line 406 column 7
+    msh2m_geometrical_properties at line 167 column 13
+    bim2c_mesh_properties at line 44 column 40
+    __test__ at line 4 column 7
+    test at line 677 column 11
+    /tmp/tmp.RoVIMMmIw0 at line 46 column 31
+
+***** test
+ [mesh] = msh2m_structured_mesh([0:.1:1],[0:.1:1],1,1:4);
+ mesh   = bim2c_mesh_properties(mesh);
+ x = mesh.p(1,:)'; y = mesh.p(2,:)';
+ Dnodes = bim2c_unknowns_on_side(mesh,[1:4]);
+ Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
+ alpha  = ones(Nelements,1); eta=ones(Nnodes,1);
+ beta   = 0;
  gamma  = ones(Nnodes,1);
- eta    = ones(Nnodes,1);
- beta   = 1/D*v*ones(Nelements,1);
- delta  = ones(Nelements,1);
- zeta   = sigma*ones(Nnodes,1);
- f      = @(r) 2*D - 2*v.*r + sigma*uex(r);
- rhs    = bim1a_rhs(mesh, ones(Nelements,1), f(mesh));
- S = bim1a_laplacian(mesh,alpha,gamma);
- A = bim1a_advection_upwind(mesh, beta);
- R = bim1a_reaction(mesh, delta, zeta);
- S += (A+R);
- u = zeros(Nnodes,1); u([1 end]) = uex(mesh([1 end]));
- u(2:end-1) = S(2:end-1,2:end-1)\(rhs(2:end-1) - S(2:end-1,[1 end])*u([1 end]));
- assert(u,uex(mesh),1e-3)
+ A = bim2a_advection_diffusion(mesh,1,1,1,0);
+ B = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ assert(A,B)
+warning: 'matrix' object indexed with empty index list
+warning: called from
+    msh2m_geometrical_properties>computearea at line 406 column 7
+    msh2m_geometrical_properties at line 167 column 13
+    bim2c_mesh_properties at line 44 column 40
+    __test__ at line 4 column 9
+    test at line 677 column 11
+    /tmp/tmp.RoVIMMmIw0 at line 46 column 31
+
+8 tests, 8 passed, 0 known failure, 0 skipped
+[inst/bim1c_elem_to_nodes.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim1c_elem_to_nodes.m
+***** test
+ n = 10; msh = linspace (0, 1, n+1);
+ nel  = n;
+ nnod = n+1;
+ u_el = randn (nel, 1);
+ un1 = bim1c_elem_to_nodes (msh, u_el);
+ [un2, m] = bim1c_elem_to_nodes (msh, u_el);
+ un3 = bim1c_elem_to_nodes (m, u_el);
+ [un4, m] = bim1c_elem_to_nodes (m, u_el);
+ assert (un1, un2, 1e-10)
+ assert (un1, un3, 1e-10)
+ assert (un1, un4, 1e-10)
 1 test, 1 passed, 0 known failure, 0 skipped
 [inst/bim2a_axisymmetric_boundary_mass.m]
 >>>>> /build/octave-bim-1.1.5/inst/bim2a_axisymmetric_boundary_mass.m
@@ -2942,7 +3190,7 @@
     bim2c_mesh_properties at line 44 column 40
     __test__ at line 5 column 9
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 38 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 62 column 31
 
 ***** test
  n = 10;
@@ -2992,7 +3240,7 @@
     bim2c_mesh_properties at line 44 column 40
     __test__ at line 5 column 9
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 38 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 62 column 31
 
 2 tests, 2 passed, 0 known failure, 0 skipped
 [inst/bim1a_axisymmetric_advection_diffusion.m]
@@ -3062,160 +3310,6 @@
  B = bim1a_axisymmetric_advection_diffusion(x,alpha,gamma,eta,0);
  assert(A,B)
 4 tests, 4 passed, 0 known failure, 0 skipped
-[inst/bim3c_unknowns_on_faces.m]
->>>>> /build/octave-bim-1.1.5/inst/bim3c_unknowns_on_faces.m
-***** shared mesh
-***** test
-***** test
-***** test
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/bim3a_reaction.m]
->>>>> /build/octave-bim-1.1.5/inst/bim3a_reaction.m
-***** shared mesh,delta,zeta,nnodes,nelem
-***** test
-***** test
-***** test
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/bim1a_reaction.m]
->>>>> /build/octave-bim-1.1.5/inst/bim1a_reaction.m
-***** test
- x = linspace(0,1,101);
- A = bim1a_reaction(x,1,1);
- delta = ones(100,1);
- zeta  = ones(101,1);
- B = bim1a_reaction(x,delta,zeta);
- assert(A,B)
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/bim1a_rhs.m]
->>>>> /build/octave-bim-1.1.5/inst/bim1a_rhs.m
-***** test
- x = linspace(0,1,101);
- A = bim1a_rhs(x,1,1);
- delta = ones(100,1);
- zeta  = ones(101,1);
- B = bim1a_rhs(x,delta,zeta);
- assert(A,B)
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/bim1a_advection_diffusion.m]
->>>>> /build/octave-bim-1.1.5/inst/bim1a_advection_diffusion.m
-***** test
- x = linspace(0,1,101);
- A = bim1a_advection_diffusion(x,1,1,1,0);
- alpha = ones(100,1);
- gamma = ones(101,1);
- eta   = gamma;
- B = bim1a_advection_diffusion(x,alpha,gamma,eta,0);
- assert(A,B)
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/bim3a_osc_laplacian.m]
->>>>> /build/octave-bim-1.1.5/inst/bim3a_osc_laplacian.m
-***** shared msh, epsilon, M, nnodes, nelem, x, y, z
-***** test
- msh = bim3c_mesh_properties (msh3m_structured_mesh (0:5, 0:5, 0:5, 1, 1:6));
- x = msh.p (1, :).';
- y = msh.p (2, :).';
- z = msh.p (3, :).';
- u = ones (size (x));
- M = bim3a_osc_laplacian (msh, 1);
- assert (M * u, zeros (size (u)), eps * 100)
-***** test
- u = x;
- bnd = bim3c_unknowns_on_faces (msh, [1, 2]);
- int = setdiff (1:columns (msh.p), bnd);
- assert (M(int, int) * u(int), -M(int, bnd) * u(bnd), 100 * eps)
-***** test
- u = y;
- bnd = bim3c_unknowns_on_faces (msh, [3, 4]);
- int = setdiff (1:columns (msh.p), bnd);
- assert (M(int, int) * u(int), -M(int, bnd) * u(bnd), 100 * eps)
-***** test
- u = z;
- bnd = bim3c_unknowns_on_faces (msh, [5, 6]);
- int = setdiff (1:columns (msh.p), bnd);
- assert (M(int, int) * u(int), -M(int, bnd) * u(bnd), 100 * eps)
-***** test
- u = z;
- bnd = bim3c_unknowns_on_faces (msh, [5, 6]);
- int = setdiff (1:columns (msh.p), bnd);
- M = bim3a_osc_laplacian (msh, pi);
- assert (M(int, int) * u(int), -M(int, bnd) * u(bnd), 100 * eps)
-5 tests, 5 passed, 0 known failure, 0 skipped
-[inst/bim3c_norm.m]
->>>>> /build/octave-bim-1.1.5/inst/bim3c_norm.m
-***** test
-***** shared L, V, x, y, z, m
- L = rand (1); V = rand (1); x = linspace (0,L,4);  y = x;  z = x;
- m = msh3m_structured_mesh (x,y,z,1,1:6);
- m.area = msh3m_geometrical_properties (m, 'area');
- m.shg  = msh3m_geometrical_properties (m, 'shg');
- u    = V * ones (columns(m.p),1);
- uinf = bim3c_norm (m, u, 'inf');
- uL2  = bim3c_norm (m, u, 'L2');
- uH1  = bim3c_norm (m, u, 'H1');
- assert ([uinf, uL2, uH1], [V, V*sqrt(L^3), V*sqrt(L^3)], 1e-12);
-***** test
- u    = V * (m.p(1,:) + 2*m.p(2,:) + 3*m.p(3,:))';
- uinf = bim3c_norm (m, u, 'inf');
- uL2  = bim3c_norm (m, u, 'L2');
- uH1  = bim3c_norm (m, u, 'H1');
- assert ([uinf, uL2, uH1], 
-         [6*L*V, V*sqrt(61/6*L^5), V*sqrt(61/6*L^5 + 14*L^3)],
-          1e-12);
-***** test
- u    = V * ones (columns(m.t),1);
- uinf = bim3c_norm (m, u, 'inf');
- uL2  = bim3c_norm (m, u, 'L2');
- assert ([uinf, uL2], [V, V*sqrt(L^3)], 1e-12);
-***** test
- u     = V * ones (columns(m.t),1);
- uvect = [u, 2*u, 3*u];
- uinf  = bim3c_norm (m, uvect, 'inf');
- uL2   = bim3c_norm (m, uvect, 'L2');
- assert ([uinf, uL2], [3*V, V*sqrt(14*L^3)], 1e-12);
-4 tests, 4 passed, 0 known failure, 0 skipped
-[inst/bim3c_tri_to_nodes.m]
->>>>> /build/octave-bim-1.1.5/inst/bim3c_tri_to_nodes.m
-***** test
- msh = bim3c_mesh_properties (msh3m_structured_mesh (linspace (0, 1, 31), linspace (0, 1, 13), linspace (0, 1, 13), 1, 1:6));
- nel  = columns (msh.t);
- nnod = columns (msh.p);
- u_tri = randn (nel, 1);
- un1 = bim3c_tri_to_nodes (msh, u_tri);
- [un2, m] = bim3c_tri_to_nodes (msh, u_tri);
- assert (un1, un2, 1e-10)
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/bim3a_laplacian.m]
->>>>> /build/octave-bim-1.1.5/inst/bim3a_laplacian.m
-***** shared mesh,epsilon,kappa,nnodes,nelem
-***** test
-***** test
-***** test
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/bim2c_intrp.m]
->>>>> /build/octave-bim-1.1.5/inst/bim2c_intrp.m
-***** test
- msh = bim2c_mesh_properties (msh2m_structured_mesh (linspace (0, 1, 11), linspace (0, 1, 13), 1, 1:4));
- x = y = linspace (0, 1, 100).';
- u = msh.p(1, :).';
- ui = bim2c_intrp (msh, u, [], [x, y]); 
- assert (ui, linspace (0, 1, 100), 10*eps);
-warning: 'matrix' object indexed with empty index list
-warning: called from
-    msh2m_geometrical_properties>computearea at line 406 column 7
-    msh2m_geometrical_properties at line 167 column 13
-    bim2c_mesh_properties at line 44 column 40
-    __test__ at line 3 column 6
-    test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 126 column 31
-
-1 test, 1 passed, 0 known failure, 0 skipped
-[inst/bim3a_rhs.m]
->>>>> /build/octave-bim-1.1.5/inst/bim3a_rhs.m
-***** shared mesh,f,g,nnodes,nelem
-***** test
-***** test
-***** test
-3 tests, 3 passed, 0 known failure, 0 skipped
 [inst/bim3a_osc_advection_diffusion.m]
 >>>>> /build/octave-bim-1.1.5/inst/bim3a_osc_advection_diffusion.m
 ***** shared msh, epsilon, M, nnodes, nelem, x, y, z
@@ -3373,35 +3467,17 @@
  fpl_vtk_raw_write_field (fname_out, msh, {uosc, "u_osc"; ugal, "u_galerkin"}, {});
  unlink (fname);
 8 tests, 8 passed, 0 known failure, 0 skipped
-[inst/bim1c_elem_to_nodes.m]
->>>>> /build/octave-bim-1.1.5/inst/bim1c_elem_to_nodes.m
+[inst/bim1a_advection_diffusion.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim1a_advection_diffusion.m
 ***** test
- n = 10; msh = linspace (0, 1, n+1);
- nel  = n;
- nnod = n+1;
- u_el = randn (nel, 1);
- un1 = bim1c_elem_to_nodes (msh, u_el);
- [un2, m] = bim1c_elem_to_nodes (msh, u_el);
- un3 = bim1c_elem_to_nodes (m, u_el);
- [un4, m] = bim1c_elem_to_nodes (m, u_el);
- assert (un1, un2, 1e-10)
- assert (un1, un3, 1e-10)
- assert (un1, un4, 1e-10)
+ x = linspace(0,1,101);
+ A = bim1a_advection_diffusion(x,1,1,1,0);
+ alpha = ones(100,1);
+ gamma = ones(101,1);
+ eta   = gamma;
+ B = bim1a_advection_diffusion(x,alpha,gamma,eta,0);
+ assert(A,B)
 1 test, 1 passed, 0 known failure, 0 skipped
-[inst/bim2a_rhs.m]
->>>>> /build/octave-bim-1.1.5/inst/bim2a_rhs.m
-***** shared mesh,f,g,nnodes,nelem
-***** test
-***** test
-***** test
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/bim2a_reaction.m]
->>>>> /build/octave-bim-1.1.5/inst/bim2a_reaction.m
-***** shared mesh,delta,zeta,nnodes,nelem
-***** test
-***** test
-***** test
-3 tests, 3 passed, 0 known failure, 0 skipped
 [inst/bim2c_tri_to_nodes.m]
 >>>>> /build/octave-bim-1.1.5/inst/bim2c_tri_to_nodes.m
 ***** test
@@ -3419,9 +3495,180 @@
     bim2c_mesh_properties at line 44 column 40
     __test__ at line 3 column 6
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 174 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 94 column 31
 
 1 test, 1 passed, 0 known failure, 0 skipped
+[inst/bim3c_tri_to_nodes.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim3c_tri_to_nodes.m
+***** test
+ msh = bim3c_mesh_properties (msh3m_structured_mesh (linspace (0, 1, 31), linspace (0, 1, 13), linspace (0, 1, 13), 1, 1:6));
+ nel  = columns (msh.t);
+ nnod = columns (msh.p);
+ u_tri = randn (nel, 1);
+ un1 = bim3c_tri_to_nodes (msh, u_tri);
+ [un2, m] = bim3c_tri_to_nodes (msh, u_tri);
+ assert (un1, un2, 1e-10)
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/bim1a_axisymmetric_reaction.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim1a_axisymmetric_reaction.m
+***** test
+ n = 100;
+ mesh      = linspace(0,1,n+1)';
+ cm        = (mesh(1:end-1) + mesh(2:end))/2; 
+ uex       = @(r) - r.^2 + 1;
+ Nnodes    = numel(mesh);
+ Nelements = Nnodes-1;
+ D = 1;  v = cm;  sigma = 1;
+ alpha  = D*ones(Nelements,1);
+ gamma  = ones(Nnodes,1);
+ eta    = ones(Nnodes,1);
+ beta   = 0;
+ delta  = ones(Nelements,1);
+ zeta   = sigma*ones(Nnodes,1);
+ f      = @(r) 4*D + sigma*uex(r);
+ rhs    = bim1a_axisymmetric_rhs(mesh, ones(Nelements,1), f(mesh));
+ S = bim1a_axisymmetric_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ R = bim1a_axisymmetric_reaction(mesh, delta, zeta);
+ S += R;
+ u = zeros(Nnodes,1); u(end) = uex(mesh(end));
+ u(1:end-1) = S(1:end-1,1:end-1)\(rhs(1:end-1) - S(1:end-1,end)*u(end));
+ assert(u,uex(mesh),1e-3)
+***** test
+ n = 100;
+ mesh      = linspace(0,1,n+1)';
+ cm        = (mesh(1:end-1) + mesh(2:end))/2; 
+ uex       = @(r) - r.^2 + 1;
+ Nnodes    = numel(mesh);
+ Nelements = Nnodes-1;
+ D = 1;  v = cm;  sigma = 1;
+ alpha  = D*ones(Nelements,1);
+ gamma  = ones(Nnodes,1);
+ eta    = ones(Nnodes,1);
+ beta   = 1/D*v;
+ delta  = ones(Nelements,1);
+ zeta   = sigma*ones(Nnodes,1);
+ f      = @(r) 4*D + 2 - 4*r.^2 + sigma*uex(r);
+ rhs    = bim1a_axisymmetric_rhs(mesh, ones(Nelements,1), f(mesh));
+ S = bim1a_axisymmetric_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ R = bim1a_axisymmetric_reaction(mesh, delta, zeta);
+ S += R;
+ u = zeros(Nnodes,1); u(end) = uex(mesh(end));
+ u(1:end-1) = S(1:end-1,1:end-1)\(rhs(1:end-1) - S(1:end-1,end)*u(end));
+ assert(u,uex(mesh),1e-3)
+***** test
+ x = linspace(0,1,101);
+ A = bim1a_axisymmetric_reaction(x,1,1);
+ delta = ones(100,1);
+ zeta  = ones(101,1);
+ B = bim1a_axisymmetric_reaction(x,delta,zeta);
+ assert(A,B)
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/bim3c_unknowns_on_faces.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim3c_unknowns_on_faces.m
+***** shared mesh
+***** test
+***** test
+***** test
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/bim1a_reaction.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim1a_reaction.m
+***** test
+ x = linspace(0,1,101);
+ A = bim1a_reaction(x,1,1);
+ delta = ones(100,1);
+ zeta  = ones(101,1);
+ B = bim1a_reaction(x,delta,zeta);
+ assert(A,B)
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/bim1c_norm.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim1c_norm.m
+***** xtest
+***** shared L, V, m
+ L = rand (1); V = rand (1); m = linspace (0,1,5).^2;  m *= L;
+ u    = V * ones (size (m))';
+ uinf = bim1c_norm (m, u, 'inf');
+ uL2  = bim1c_norm (m, u, 'L2');
+ uH1  = bim1c_norm (m, u, 'H1');
+ assert ([uinf, uL2, uH1], [V, V*sqrt(L), V*sqrt(L)], 1e-12);
+***** test
+ u    = V * m';
+ uinf = bim1c_norm (m, u, 'inf');
+ uL2  = bim1c_norm (m, u, 'L2');
+ uH1  = bim1c_norm (m, u, 'H1');
+ assert ([uinf, uL2, uH1], 
+         [L*V, V*sqrt(L^3/3), V*sqrt(L^3/3 + L)],
+          1e-12);
+***** test
+ u    = V * ones (size (diff (m)))';
+ uinf = bim1c_norm (m, u, 'inf');
+ uL2  = bim1c_norm (m, u, 'L2');
+ assert ([uinf, uL2], [V, V*sqrt(L)], 1e-12);
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/bim1a_advection_upwind.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim1a_advection_upwind.m
+***** test
+ n = 200;
+ mesh      = linspace(0,1,n+1)';
+ uex       = @(r) - r.^2 + 1;
+ Nnodes    = numel(mesh);
+ Nelements = Nnodes-1;
+ D = 1;  v = 1;  sigma = 0;
+ alpha  = D*ones(Nelements,1);
+ gamma  = ones(Nnodes,1);
+ eta    = ones(Nnodes,1);
+ beta   = 1/D*v*ones(Nelements,1);
+ delta  = ones(Nelements,1);
+ zeta   = sigma*ones(Nnodes,1);
+ f      = @(r) 2*D - 2*v.*r + sigma*uex(r);
+ rhs    = bim1a_rhs(mesh, ones(Nelements,1), f(mesh));
+ S = bim1a_laplacian(mesh,alpha,gamma);
+ A = bim1a_advection_upwind(mesh, beta);
+ R = bim1a_reaction(mesh, delta, zeta);
+ S += (A+R);
+ u = zeros(Nnodes,1); u([1 end]) = uex(mesh([1 end]));
+ u(2:end-1) = S(2:end-1,2:end-1)\(rhs(2:end-1) - S(2:end-1,[1 end])*u([1 end]));
+ assert(u,uex(mesh),1e-3)
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/bim3a_osc_laplacian.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim3a_osc_laplacian.m
+***** shared msh, epsilon, M, nnodes, nelem, x, y, z
+***** test
+ msh = bim3c_mesh_properties (msh3m_structured_mesh (0:5, 0:5, 0:5, 1, 1:6));
+ x = msh.p (1, :).';
+ y = msh.p (2, :).';
+ z = msh.p (3, :).';
+ u = ones (size (x));
+ M = bim3a_osc_laplacian (msh, 1);
+ assert (M * u, zeros (size (u)), eps * 100)
+***** test
+ u = x;
+ bnd = bim3c_unknowns_on_faces (msh, [1, 2]);
+ int = setdiff (1:columns (msh.p), bnd);
+ assert (M(int, int) * u(int), -M(int, bnd) * u(bnd), 100 * eps)
+***** test
+ u = y;
+ bnd = bim3c_unknowns_on_faces (msh, [3, 4]);
+ int = setdiff (1:columns (msh.p), bnd);
+ assert (M(int, int) * u(int), -M(int, bnd) * u(bnd), 100 * eps)
+***** test
+ u = z;
+ bnd = bim3c_unknowns_on_faces (msh, [5, 6]);
+ int = setdiff (1:columns (msh.p), bnd);
+ assert (M(int, int) * u(int), -M(int, bnd) * u(bnd), 100 * eps)
+***** test
+ u = z;
+ bnd = bim3c_unknowns_on_faces (msh, [5, 6]);
+ int = setdiff (1:columns (msh.p), bnd);
+ M = bim3a_osc_laplacian (msh, pi);
+ assert (M(int, int) * u(int), -M(int, bnd) * u(bnd), 100 * eps)
+5 tests, 5 passed, 0 known failure, 0 skipped
+[inst/bim3a_laplacian.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim3a_laplacian.m
+***** shared mesh,epsilon,kappa,nnodes,nelem
+***** test
+***** test
+***** test
+3 tests, 3 passed, 0 known failure, 0 skipped
 [inst/bim1a_axisymmetric_advection_upwind.m]
 >>>>> /build/octave-bim-1.1.5/inst/bim1a_axisymmetric_advection_upwind.m
 ***** test
@@ -3446,22 +3693,53 @@
  uh = M \ rhs;  
  assert(uh, uex(mesh), 1e-3);
 1 test, 1 passed, 0 known failure, 0 skipped
-[inst/bim2a_axisymmetric_rhs.m]
->>>>> /build/octave-bim-1.1.5/inst/bim2a_axisymmetric_rhs.m
-***** shared mesh,f,g,nnodes,nelem
+[inst/bim2a_axisymmetric_reaction.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim2a_axisymmetric_reaction.m
+***** shared mesh,delta,zeta,nnodes,nelem
 ***** test
 ***** test
 ***** test
-3 tests, 3 passed, 0 known failure, 0 skipped
-[inst/bim3c_intrp.m]
->>>>> /build/octave-bim-1.1.5/inst/bim3c_intrp.m
 ***** test
- msh = bim3c_mesh_properties (msh3m_structured_mesh (linspace (0, 1, 11), linspace (0, 1, 9), linspace (0, 1, 13), 1, 1:6));
- x = y = z = linspace (0, 1, 100).';
- u = msh.p(1, :).';
- ui = bim3c_intrp (msh, u, [], [x, y, z]); 
- assert (ui, linspace (0, 1, 100), 10*eps);
-1 test, 1 passed, 0 known failure, 0 skipped
+ n = 20;
+ [mesh] = msh2m_structured_mesh(linspace(1,2,n+1),linspace(0,1,n+1),1,1:4);
+ mesh   = bim2c_mesh_properties(mesh);
+ uex    = @(r,z) exp(r) .* exp(1-z);
+ duexdr = @(r,z) uex(r,z);
+ d2uexdr2 = @(r,z) uex(r,z);
+ duexdz = @(r,z) -uex(r,z);
+ d2uexdz2 = @(r,z) uex(r,z);
+ Dnodes = bim2c_unknowns_on_side(mesh,[1,2,3,4]);
+ Nnodes    = columns(mesh.p);
+ Nelements = columns(mesh.t);
+ Varnodes  = setdiff(1:Nnodes,Dnodes);
+ D = 1; vr = 1; vz = 1; sigma = 1;
+ alpha  = D*ones(Nelements,1);
+ gamma  = ones(Nnodes,1);
+ eta    = ones(Nnodes,1);
+ delta  = sigma*ones(columns(mesh.t),1);
+ zeta   = ones(columns(mesh.p),1);
+ beta   = 1/D*[vr*ones(1,Nelements); vz*ones(1,Nelements)];
+ f = @(r,z) -D./r.*duexdr(r,z) - D.*d2uexdr2(r,z) ...
+           + vr./r .* uex(r,z) + vr * duexdr(r,z) ...
+           - D.*d2uexdz2(r,z) + vz * duexdz(r,z) ...
+           + sigma * uex(r,z);
+ rhs    = bim2a_axisymmetric_rhs(mesh, ones(Nelements,1), f(mesh.p(1,:), mesh.p(2,:)));
+ S = bim2a_axisymmetric_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ C = bim2a_axisymmetric_reaction(mesh,delta,zeta);
+ S += C;
+ u = zeros(Nnodes,1); u(Dnodes) = uex(mesh.p(1,Dnodes), mesh.p(2,Dnodes));
+ u(Varnodes) = S(Varnodes,Varnodes)\(rhs(Varnodes) - S(Varnodes,Dnodes)*u(Dnodes));
+ assert(u,uex(mesh.p(1,:), mesh.p(2,:))',1e-3)
+warning: 'matrix' object indexed with empty index list
+warning: called from
+    msh2m_geometrical_properties>computearea at line 406 column 7
+    msh2m_geometrical_properties at line 167 column 13
+    bim2c_mesh_properties at line 44 column 40
+    __test__ at line 5 column 9
+    test at line 677 column 11
+    /tmp/tmp.RoVIMMmIw0 at line 174 column 31
+
+4 tests, 4 passed, 0 known failure, 0 skipped
 [inst/bim2c_norm.m]
 >>>>> /build/octave-bim-1.1.5/inst/bim2c_norm.m
 ***** test
@@ -3495,225 +3773,15 @@
  uL2   = bim2c_norm (m, uvect, 'L2');
  assert ([uinf, uL2], [2*V, V*L*sqrt(5)], 1e-12);
 4 tests, 4 passed, 0 known failure, 0 skipped
-[inst/bim2a_advection_diffusion.m]
->>>>> /build/octave-bim-1.1.5/inst/bim2a_advection_diffusion.m
-***** test 
- [mesh] = msh2m_structured_mesh([0:1/3:1],[0:1/3:1],1,1:4);
- mesh   = bim2c_mesh_properties(mesh);
- x      = mesh.p(1,:)';
- Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
- Nnodes    = columns(mesh.p);
- Nelements = columns(mesh.t);
- Varnodes  = setdiff(1:Nnodes,Dnodes);
- alpha  = ones(Nelements,1);
- eta    = .1*ones(Nnodes,1);
- beta   = [ones(1,Nelements);zeros(1,Nelements)];
- gamma  = ones(Nnodes,1);
- f      = bim2a_rhs(mesh,ones(Nelements,1),ones(Nnodes,1));
- S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
- u = zeros(Nnodes,1);
- u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
- uex = x - (exp(10*x)-1)/(exp(10)-1);
- assert(u,uex,1e-7)
-warning: 'matrix' object indexed with empty index list
-warning: called from
-    msh2m_geometrical_properties>computearea at line 406 column 7
-    msh2m_geometrical_properties at line 167 column 13
-    bim2c_mesh_properties at line 44 column 40
-    __test__ at line 4 column 9
-    test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 214 column 31
-
-***** test 
- [mesh] = msh2m_structured_mesh([0:1/3:1],[0:1/3:1],1,1:4);
- mesh   = bim2c_mesh_properties(mesh);
- x      = mesh.p(1,:)';
- Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
- Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
- Varnodes = setdiff(1:Nnodes,Dnodes);
- alpha  = ones(Nelements,1);
- eta    = .1*ones(Nnodes,1);
- beta   = x;
- gamma  = ones(Nnodes,1);
- f      = bim2a_rhs(mesh,ones(Nelements,1),ones(Nnodes,1));
- S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
- u = zeros(Nnodes,1);
- u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
- uex = x - (exp(10*x)-1)/(exp(10)-1);
- assert(u,uex,1e-7)
-warning: 'matrix' object indexed with empty index list
-warning: called from
-    msh2m_geometrical_properties>computearea at line 406 column 7
-    msh2m_geometrical_properties at line 167 column 13
-    bim2c_mesh_properties at line 44 column 40
-    __test__ at line 4 column 9
-    test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 214 column 31
-
-***** test
- [mesh] = msh2m_structured_mesh([0:1/3:1],[0:1/3:1],1,1:4);
- mesh   = bim2c_mesh_properties(mesh);
- x      = mesh.p(1,:)';
- Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
- Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
- Varnodes = setdiff(1:Nnodes,Dnodes);
- alpha  = 10*ones(Nelements,1);
- eta    = .01*ones(Nnodes,1);
- beta   = x/10;
- gamma  = ones(Nnodes,1);
- f      = bim2a_rhs(mesh,ones(Nelements,1),ones(Nnodes,1));
- S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
- u = zeros(Nnodes,1);
- u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
- uex = x - (exp(10*x)-1)/(exp(10)-1);
- assert(u,uex,1e-7)
-warning: 'matrix' object indexed with empty index list
-warning: called from
-    msh2m_geometrical_properties>computearea at line 406 column 7
-    msh2m_geometrical_properties at line 167 column 13
-    bim2c_mesh_properties at line 44 column 40
-    __test__ at line 4 column 9
-    test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 214 column 31
-
-***** test
- [mesh] = msh2m_structured_mesh([0:1/3:1],[0:1/3:1],1,1:4);
- mesh   = bim2c_mesh_properties(mesh);
- x      = mesh.p(1,:)';
- Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
- Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
- Varnodes = setdiff(1:Nnodes,Dnodes);
- alpha  = 10*ones(Nelements,1); eta = .001*ones(Nnodes,1);
- beta   = x/100;
- gamma  = 10*ones(Nnodes,1);
- f      = bim2a_rhs(mesh,ones(Nelements,1),ones(Nnodes,1));
- S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
- u = zeros(Nnodes,1);
- u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
- uex = x - (exp(10*x)-1)/(exp(10)-1);
- assert(u,uex,1e-7)
-warning: 'matrix' object indexed with empty index list
-warning: called from
-    msh2m_geometrical_properties>computearea at line 406 column 7
-    msh2m_geometrical_properties at line 167 column 13
-    bim2c_mesh_properties at line 44 column 40
-    __test__ at line 4 column 9
-    test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 214 column 31
-
-***** test
- [mesh] = msh2m_structured_mesh([0:1/1e3:1],[0:1/2:1],1,1:4);
- mesh   = bim2c_mesh_properties(mesh);
- x      = mesh.p(1,:)';
- Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
- Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
- Varnodes = setdiff(1:Nnodes,Dnodes);
- alpha  = 3*ones(Nelements,1); eta = x+1;
- beta   = [ones(1,Nelements);zeros(1,Nelements)];
- gamma  = 2*x;
- ff     = 2*(6*x.^2+6*x) - (6*x+6).*(1-2*x)+6*(x-x.^2);
- f      = bim2a_rhs(mesh,ones(Nelements,1),ff);
- S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
- u = zeros(Nnodes,1);
- u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
- uex = x - x.^2;
- assert(u,uex,5e-3)
-warning: 'matrix' object indexed with empty index list
-warning: called from
-    msh2m_geometrical_properties>computearea at line 406 column 7
-    msh2m_geometrical_properties at line 167 column 13
-    bim2c_mesh_properties at line 44 column 40
-    __test__ at line 4 column 9
-    test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 214 column 31
-
-***** test
- [mesh] = msh2m_structured_mesh([0:1/1e3:1],[0:1/2:1],1,1:4);
- mesh = bim2c_mesh_properties(mesh);
- x = mesh.p(1,:)';
- Dnodes = bim2c_unknowns_on_side(mesh,[2,4]);
- Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
- Varnodes = setdiff(1:Nnodes,Dnodes);
- alpha  = ones(Nelements,1); eta = ones(Nnodes,1);
- beta   = 0;
- gamma  = x+1;
- ff     = 4*x+1;
- f      = bim2a_rhs(mesh,ones(Nelements,1),ff);
- S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
- u = zeros(Nnodes,1);
- u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
- uex = x - x.^2; 
- assert(u,uex,1e-7)
-warning: 'matrix' object indexed with empty index list
-warning: called from
-    msh2m_geometrical_properties>computearea at line 406 column 7
-    msh2m_geometrical_properties at line 167 column 13
-    bim2c_mesh_properties at line 44 column 40
-    __test__ at line 4 column 7
-    test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 214 column 31
-
-***** test
- [mesh] = msh2m_structured_mesh([0:.1:1],[0:.1:1],1,1:4);
- mesh = bim2c_mesh_properties(mesh);
- x = mesh.p(1,:)';y = mesh.p(2,:)';
- Dnodes = bim2c_unknowns_on_side(mesh,[1:4]);
- Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
- Varnodes = setdiff(1:Nnodes,Dnodes);
- alpha  = ones(Nelements,1); diff = 1e-2; eta=diff*ones(Nnodes,1);
- beta   =[ones(1,Nelements);ones(1,Nelements)];
- gamma  = x*0+1;
- ux  = y.*(1-exp((y-1)/diff)) .* (1-exp((x-1)/diff)-x.*exp((x-1)/diff)/diff);
- uy  = x.*(1-exp((x-1)/diff)) .* (1-exp((y-1)/diff)-y.*exp((y-1)/diff)/diff);
- uxx = y.*(1-exp((y-1)/diff)) .* (-2*exp((x-1)/diff)/diff-x.*exp((x-1)/diff)/(diff^2));
- uyy = x.*(1-exp((x-1)/diff)) .* (-2*exp((y-1)/diff)/diff-y.*exp((y-1)/diff)/(diff^2));
- ff  = -diff*(uxx+uyy)+ux+uy;
- f   = bim2a_rhs(mesh,ones(Nelements,1),ff);
- S = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
- u = zeros(Nnodes,1);
- u(Varnodes) = S(Varnodes,Varnodes)\f(Varnodes);
- uex = x.*y.*(1-exp((x-1)/diff)).*(1-exp((y-1)/diff)); 
- assert(u,uex,1e-7)
-warning: 'matrix' object indexed with empty index list
-warning: called from
-    msh2m_geometrical_properties>computearea at line 406 column 7
-    msh2m_geometrical_properties at line 167 column 13
-    bim2c_mesh_properties at line 44 column 40
-    __test__ at line 4 column 7
-    test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 214 column 31
-
+[inst/bim1a_rhs.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim1a_rhs.m
 ***** test
- [mesh] = msh2m_structured_mesh([0:.1:1],[0:.1:1],1,1:4);
- mesh   = bim2c_mesh_properties(mesh);
- x = mesh.p(1,:)'; y = mesh.p(2,:)';
- Dnodes = bim2c_unknowns_on_side(mesh,[1:4]);
- Nnodes = columns(mesh.p); Nelements = columns(mesh.t);
- alpha  = ones(Nelements,1); eta=ones(Nnodes,1);
- beta   = 0;
- gamma  = ones(Nnodes,1);
- A = bim2a_advection_diffusion(mesh,1,1,1,0);
- B = bim2a_advection_diffusion(mesh,alpha,gamma,eta,beta);
+ x = linspace(0,1,101);
+ A = bim1a_rhs(x,1,1);
+ delta = ones(100,1);
+ zeta  = ones(101,1);
+ B = bim1a_rhs(x,delta,zeta);
  assert(A,B)
-warning: 'matrix' object indexed with empty index list
-warning: called from
-    msh2m_geometrical_properties>computearea at line 406 column 7
-    msh2m_geometrical_properties at line 167 column 13
-    bim2c_mesh_properties at line 44 column 40
-    __test__ at line 4 column 9
-    test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 214 column 31
-
-8 tests, 8 passed, 0 known failure, 0 skipped
-[inst/bim3c_global_flux.m]
->>>>> /build/octave-bim-1.1.5/inst/bim3c_global_flux.m
-***** test
- N = 10; pp = linspace (0, 1, N); msh = bim3c_mesh_properties (msh3m_structured_mesh (pp, pp, pp, 1, 1:6));
- u = ones (N^3, 1);
- v = ones (N^3, 1);
- alpha = ones (columns (msh.t), 1);
- F =  bim3c_global_flux (msh, u, alpha, v);
- assert (norm (F(:), inf), 0, 100*eps);
 1 test, 1 passed, 0 known failure, 0 skipped
 [inst/bim3c_mesh_properties.m]
 >>>>> /build/octave-bim-1.1.5/inst/bim3c_mesh_properties.m
@@ -3733,6 +3801,23 @@
  B = bim1a_axisymmetric_rhs(x,delta,zeta);
  assert(A,B)
 1 test, 1 passed, 0 known failure, 0 skipped
+[inst/bim3c_global_flux.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim3c_global_flux.m
+***** test
+ N = 10; pp = linspace (0, 1, N); msh = bim3c_mesh_properties (msh3m_structured_mesh (pp, pp, pp, 1, 1:6));
+ u = ones (N^3, 1);
+ v = ones (N^3, 1);
+ alpha = ones (columns (msh.t), 1);
+ F =  bim3c_global_flux (msh, u, alpha, v);
+ assert (norm (F(:), inf), 0, 100*eps);
+1 test, 1 passed, 0 known failure, 0 skipped
+[inst/bim3a_rhs.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim3a_rhs.m
+***** shared mesh,f,g,nnodes,nelem
+***** test
+***** test
+***** test
+3 tests, 3 passed, 0 known failure, 0 skipped
 [inst/bim2a_axisymmetric_advection_diffusion.m]
 >>>>> /build/octave-bim-1.1.5/inst/bim2a_axisymmetric_advection_diffusion.m
 ***** test
@@ -3768,7 +3853,7 @@
     bim2c_mesh_properties at line 44 column 40
     __test__ at line 5 column 12
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 246 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 230 column 31
 
 ***** test
  n         = 20;
@@ -3803,7 +3888,7 @@
     bim2c_mesh_properties at line 44 column 40
     __test__ at line 5 column 12
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 246 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 230 column 31
 
 ***** test
  n         = 10;
@@ -3838,7 +3923,7 @@
     bim2c_mesh_properties at line 44 column 40
     __test__ at line 5 column 12
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 246 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 230 column 31
 
 ***** test
  n         = 10;
@@ -3871,7 +3956,7 @@
     bim2c_mesh_properties at line 44 column 40
     __test__ at line 5 column 12
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 246 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 230 column 31
 
 ***** test
  n         = 3;
@@ -3906,7 +3991,7 @@
     bim2c_mesh_properties at line 44 column 40
     __test__ at line 5 column 12
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 246 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 230 column 31
 
 ***** test
  n         = 10;
@@ -3939,7 +4024,7 @@
     bim2c_mesh_properties at line 44 column 40
     __test__ at line 5 column 12
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 246 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 230 column 31
 
 ***** test
  [mesh] = msh2m_structured_mesh([0:.1:1],[0:.1:1],1,1:4);
@@ -3960,62 +4045,15 @@
     bim2c_mesh_properties at line 44 column 40
     __test__ at line 4 column 9
     test at line 677 column 11
-    /tmp/tmp.uNQIAI2I9a at line 246 column 31
+    /tmp/tmp.RoVIMMmIw0 at line 230 column 31
 
 7 tests, 7 passed, 0 known failure, 0 skipped
-[inst/bim1a_axisymmetric_reaction.m]
->>>>> /build/octave-bim-1.1.5/inst/bim1a_axisymmetric_reaction.m
+[inst/bim2a_rhs.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim2a_rhs.m
+***** shared mesh,f,g,nnodes,nelem
 ***** test
- n = 100;
- mesh      = linspace(0,1,n+1)';
- cm        = (mesh(1:end-1) + mesh(2:end))/2; 
- uex       = @(r) - r.^2 + 1;
- Nnodes    = numel(mesh);
- Nelements = Nnodes-1;
- D = 1;  v = cm;  sigma = 1;
- alpha  = D*ones(Nelements,1);
- gamma  = ones(Nnodes,1);
- eta    = ones(Nnodes,1);
- beta   = 0;
- delta  = ones(Nelements,1);
- zeta   = sigma*ones(Nnodes,1);
- f      = @(r) 4*D + sigma*uex(r);
- rhs    = bim1a_axisymmetric_rhs(mesh, ones(Nelements,1), f(mesh));
- S = bim1a_axisymmetric_advection_diffusion(mesh,alpha,gamma,eta,beta);
- R = bim1a_axisymmetric_reaction(mesh, delta, zeta);
- S += R;
- u = zeros(Nnodes,1); u(end) = uex(mesh(end));
- u(1:end-1) = S(1:end-1,1:end-1)\(rhs(1:end-1) - S(1:end-1,end)*u(end));
- assert(u,uex(mesh),1e-3)
 ***** test
- n = 100;
- mesh      = linspace(0,1,n+1)';
- cm        = (mesh(1:end-1) + mesh(2:end))/2; 
- uex       = @(r) - r.^2 + 1;
- Nnodes    = numel(mesh);
- Nelements = Nnodes-1;
- D = 1;  v = cm;  sigma = 1;
- alpha  = D*ones(Nelements,1);
- gamma  = ones(Nnodes,1);
- eta    = ones(Nnodes,1);
- beta   = 1/D*v;
- delta  = ones(Nelements,1);
- zeta   = sigma*ones(Nnodes,1);
- f      = @(r) 4*D + 2 - 4*r.^2 + sigma*uex(r);
- rhs    = bim1a_axisymmetric_rhs(mesh, ones(Nelements,1), f(mesh));
- S = bim1a_axisymmetric_advection_diffusion(mesh,alpha,gamma,eta,beta);
- R = bim1a_axisymmetric_reaction(mesh, delta, zeta);
- S += R;
- u = zeros(Nnodes,1); u(end) = uex(mesh(end));
- u(1:end-1) = S(1:end-1,1:end-1)\(rhs(1:end-1) - S(1:end-1,end)*u(end));
- assert(u,uex(mesh),1e-3)
 ***** test
- x = linspace(0,1,101);
- A = bim1a_axisymmetric_reaction(x,1,1);
- delta = ones(100,1);
- zeta  = ones(101,1);
- B = bim1a_axisymmetric_reaction(x,delta,zeta);
- assert(A,B)
 3 tests, 3 passed, 0 known failure, 0 skipped
 [inst/bim2c_mesh_properties.m]
 >>>>> /build/octave-bim-1.1.5/inst/bim2c_mesh_properties.m
@@ -4024,6 +4062,20 @@
 ***** test
 ***** test
 3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/bim2a_axisymmetric_rhs.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim2a_axisymmetric_rhs.m
+***** shared mesh,f,g,nnodes,nelem
+***** test
+***** test
+***** test
+3 tests, 3 passed, 0 known failure, 0 skipped
+[inst/bim2a_reaction.m]
+>>>>> /build/octave-bim-1.1.5/inst/bim2a_reaction.m
+***** shared mesh,delta,zeta,nnodes,nelem
+***** test
+***** test
+***** test
+3 tests, 3 passed, 0 known failure, 0 skipped
 Checking C++ files ...
 Summary: 92 tests, 92 passed, 0 known failures, 0 skipped
    dh_installdocs -O--buildsystem=octave
@@ -4051,12 +4103,14 @@
 dpkg-buildpackage: info: binary-only upload (no source included)
 dpkg-genchanges: info: not including original source code in upload
 I: copying local configuration
+I: user script /srv/workspace/pbuilder/28350/tmp/hooks/B01_cleanup starting
+I: user script /srv/workspace/pbuilder/28350/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/12574 and its subdirectories
-I: Current time: Tue Aug 10 01:20:06 -12 2021
-I: pbuilder-time-stamp: 1628601606
+I: removing directory /srv/workspace/pbuilder/28350 and its subdirectories
+I: Current time: Wed Aug 11 03:26:21 +14 2021
+I: pbuilder-time-stamp: 1628601981