Diff of the two buildlogs: -- --- b1/build.log 2023-06-08 17:23:25.565562942 +0000 +++ b2/build.log 2023-06-08 17:34:45.891779789 +0000 @@ -1,6 +1,6 @@ I: pbuilder: network access will be disabled during build -I: Current time: Thu Jun 8 05:12:42 -12 2023 -I: pbuilder-time-stamp: 1686244362 +I: Current time: Fri Jun 9 07:23:47 +14 2023 +I: pbuilder-time-stamp: 1686245027 I: Building the build Environment I: extracting base tarball [/var/cache/pbuilder/bookworm-reproducible-base.tgz] I: copying local configuration @@ -16,7 +16,7 @@ I: copying [./octave-stk_2.7.0.orig.tar.gz] I: copying [./octave-stk_2.7.0-3.debian.tar.xz] I: Extracting source -gpgv: Signature made Mon Dec 5 20:43:42 2022 -12 +gpgv: Signature made Tue Dec 6 22:43:42 2022 +14 gpgv: using RSA key 3F464391498FE874BDB5D98F2124AA1983785C90 gpgv: issuer "rafael@debian.org" gpgv: Can't check signature: No public key @@ -29,135 +29,167 @@ dpkg-source: info: applying octave7.patch I: Not using root during the build. I: Installing the build-deps -I: user script /srv/workspace/pbuilder/7156/tmp/hooks/D02_print_environment starting +I: user script /srv/workspace/pbuilder/9579/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 +'/bin/sh' -> '/bin/bash' +lrwxrwxrwx 1 root root 9 Jun 9 07:24 /bin/sh -> /bin/bash +I: Setting pbuilder2's login shell to /bin/bash +I: Setting pbuilder2's GECOS to second user,second room,second work-phone,second home-phone,second other +I: user script /srv/workspace/pbuilder/9579/tmp/hooks/D01_modify_environment finished +I: user script /srv/workspace/pbuilder/9579/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='bookworm' - HOME='/root' - HOST_ARCH='armhf' + BASH=/bin/sh + BASHOPTS=checkwinsize:cmdhist:complete_fullquote:extquote:force_fignore:globasciiranges:globskipdots:hostcomplete:interactive_comments:patsub_replacement:progcomp:promptvars:sourcepath + BASH_ALIASES=() + BASH_ARGC=() + BASH_ARGV=() + BASH_CMDS=() + BASH_LINENO=([0]="12" [1]="0") + BASH_LOADABLES_PATH=/usr/local/lib/bash:/usr/lib/bash:/opt/local/lib/bash:/usr/pkg/lib/bash:/opt/pkg/lib/bash:. + BASH_SOURCE=([0]="/tmp/hooks/D02_print_environment" [1]="/tmp/hooks/D02_print_environment") + BASH_VERSINFO=([0]="5" [1]="2" [2]="15" [3]="1" [4]="release" [5]="arm-unknown-linux-gnueabihf") + BASH_VERSION='5.2.15(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=bookworm + EUID=0 + FUNCNAME=([0]="Echo" [1]="main") + GROUPS=() + HOME=/root + HOSTNAME=i-capture-the-hostname + HOSTTYPE=arm + HOST_ARCH=armhf IFS=' ' - INVOCATION_ID='d5907f64c0ac4f748e59a6bab639a18a' - 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='7156' - PS1='# ' - PS2='> ' + INVOCATION_ID=334e0c390bea42adb9e18b59837e44ca + 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=9579 PS4='+ ' - PWD='/' - SHELL='/bin/bash' - SHLVL='2' - SUDO_COMMAND='/usr/bin/timeout -k 18.1h 18h /usr/bin/ionice -c 3 /usr/bin/nice /usr/sbin/pbuilder --build --configfile /srv/reproducible-results/rbuild-debian/r-b-build.E7tm2Pb6/pbuilderrc_jsl9 --distribution bookworm --hookdir /etc/pbuilder/first-build-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/bookworm-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.E7tm2Pb6/b1 --logfile b1/build.log octave-stk_2.7.0-3.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:3142/' + PWD=/ + SHELL=/bin/bash + SHELLOPTS=braceexpand:errexit:hashall:interactive-comments:posix + SHLVL=3 + SUDO_COMMAND='/usr/bin/timeout -k 24.1h 24h /usr/bin/ionice -c 3 /usr/bin/nice -n 11 /usr/bin/unshare --uts -- /usr/sbin/pbuilder --build --configfile /srv/reproducible-results/rbuild-debian/r-b-build.E7tm2Pb6/pbuilderrc_WhiK --distribution bookworm --hookdir /etc/pbuilder/rebuild-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/bookworm-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/r-b-build.E7tm2Pb6/b2 --logfile b2/build.log --extrapackages usrmerge octave-stk_2.7.0-3.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:3142/ I: uname -a - Linux virt64b 5.10.0-23-arm64 #1 SMP Debian 5.10.179-1 (2023-05-12) aarch64 GNU/Linux + Linux i-capture-the-hostname 5.10.0-23-armmp-lpae #1 SMP Debian 5.10.179-1 (2023-05-12) armv7l GNU/Linux I: ls -l /bin total 5072 - -rwxr-xr-x 1 root root 838488 Apr 23 09:24 bash - -rwxr-xr-x 3 root root 67144 Sep 18 2022 bunzip2 - -rwxr-xr-x 3 root root 67144 Sep 18 2022 bzcat - lrwxrwxrwx 1 root root 6 Sep 18 2022 bzcmp -> bzdiff - -rwxr-xr-x 1 root root 2225 Sep 18 2022 bzdiff - lrwxrwxrwx 1 root root 6 Sep 18 2022 bzegrep -> bzgrep - -rwxr-xr-x 1 root root 4893 Nov 27 2021 bzexe - lrwxrwxrwx 1 root root 6 Sep 18 2022 bzfgrep -> bzgrep - -rwxr-xr-x 1 root root 3775 Sep 18 2022 bzgrep - -rwxr-xr-x 3 root root 67144 Sep 18 2022 bzip2 - -rwxr-xr-x 1 root root 67112 Sep 18 2022 bzip2recover - lrwxrwxrwx 1 root root 6 Sep 18 2022 bzless -> bzmore - -rwxr-xr-x 1 root root 1297 Sep 18 2022 bzmore - -rwxr-xr-x 1 root root 67632 Sep 20 2022 cat - -rwxr-xr-x 1 root root 67676 Sep 20 2022 chgrp - -rwxr-xr-x 1 root root 67644 Sep 20 2022 chmod - -rwxr-xr-x 1 root root 67684 Sep 20 2022 chown - -rwxr-xr-x 1 root root 133532 Sep 20 2022 cp - -rwxr-xr-x 1 root root 132868 Jan 5 01:20 dash - -rwxr-xr-x 1 root root 133220 Sep 20 2022 date - -rwxr-xr-x 1 root root 67732 Sep 20 2022 dd - -rwxr-xr-x 1 root root 68104 Sep 20 2022 df - -rwxr-xr-x 1 root root 133632 Sep 20 2022 dir - -rwxr-xr-x 1 root root 59128 Mar 22 21:02 dmesg - lrwxrwxrwx 1 root root 8 Dec 19 01:33 dnsdomainname -> hostname - lrwxrwxrwx 1 root root 8 Dec 19 01:33 domainname -> hostname - -rwxr-xr-x 1 root root 67560 Sep 20 2022 echo - -rwxr-xr-x 1 root root 41 Jan 24 02:43 egrep - -rwxr-xr-x 1 root root 67548 Sep 20 2022 false - -rwxr-xr-x 1 root root 41 Jan 24 02:43 fgrep - -rwxr-xr-x 1 root root 55748 Mar 22 21:02 findmnt - -rwsr-xr-x 1 root root 26208 Mar 22 20:15 fusermount - -rwxr-xr-x 1 root root 128608 Jan 24 02:43 grep - -rwxr-xr-x 2 root root 2346 Apr 9 2022 gunzip - -rwxr-xr-x 1 root root 6447 Apr 9 2022 gzexe - -rwxr-xr-x 1 root root 64220 Apr 9 2022 gzip - -rwxr-xr-x 1 root root 67032 Dec 19 01:33 hostname - -rwxr-xr-x 1 root root 67720 Sep 20 2022 ln - -rwxr-xr-x 1 root root 35132 Mar 22 21:51 login - -rwxr-xr-x 1 root root 133632 Sep 20 2022 ls - -rwxr-xr-x 1 root root 136808 Mar 22 21:02 lsblk - -rwxr-xr-x 1 root root 67800 Sep 20 2022 mkdir - -rwxr-xr-x 1 root root 67764 Sep 20 2022 mknod - -rwxr-xr-x 1 root root 67596 Sep 20 2022 mktemp - -rwxr-xr-x 1 root root 38504 Mar 22 21:02 more - -rwsr-xr-x 1 root root 38496 Mar 22 21:02 mount - -rwxr-xr-x 1 root root 9824 Mar 22 21:02 mountpoint - -rwxr-xr-x 1 root root 133532 Sep 20 2022 mv - lrwxrwxrwx 1 root root 8 Dec 19 01:33 nisdomainname -> hostname - lrwxrwxrwx 1 root root 14 Apr 2 18:25 pidof -> /sbin/killall5 - -rwxr-xr-x 1 root root 67608 Sep 20 2022 pwd - lrwxrwxrwx 1 root root 4 Apr 23 09:24 rbash -> bash - -rwxr-xr-x 1 root root 67600 Sep 20 2022 readlink - -rwxr-xr-x 1 root root 67672 Sep 20 2022 rm - -rwxr-xr-x 1 root root 67600 Sep 20 2022 rmdir - -rwxr-xr-x 1 root root 67400 Nov 2 2022 run-parts - -rwxr-xr-x 1 root root 133372 Jan 5 07:55 sed - lrwxrwxrwx 1 root root 4 Jan 5 01:20 sh -> dash - -rwxr-xr-x 1 root root 67584 Sep 20 2022 sleep - -rwxr-xr-x 1 root root 67644 Sep 20 2022 stty - -rwsr-xr-x 1 root root 50800 Mar 22 21:02 su - -rwxr-xr-x 1 root root 67584 Sep 20 2022 sync - -rwxr-xr-x 1 root root 336764 Apr 6 02:25 tar - -rwxr-xr-x 1 root root 67144 Nov 2 2022 tempfile - -rwxr-xr-x 1 root root 133224 Sep 20 2022 touch - -rwxr-xr-x 1 root root 67548 Sep 20 2022 true - -rwxr-xr-x 1 root root 9768 Mar 22 20:15 ulockmgr_server - -rwsr-xr-x 1 root root 22108 Mar 22 21:02 umount - -rwxr-xr-x 1 root root 67572 Sep 20 2022 uname - -rwxr-xr-x 2 root root 2346 Apr 9 2022 uncompress - -rwxr-xr-x 1 root root 133632 Sep 20 2022 vdir - -rwxr-xr-x 1 root root 42608 Mar 22 21:02 wdctl - lrwxrwxrwx 1 root root 8 Dec 19 01:33 ypdomainname -> hostname - -rwxr-xr-x 1 root root 1984 Apr 9 2022 zcat - -rwxr-xr-x 1 root root 1678 Apr 9 2022 zcmp - -rwxr-xr-x 1 root root 6460 Apr 9 2022 zdiff - -rwxr-xr-x 1 root root 29 Apr 9 2022 zegrep - -rwxr-xr-x 1 root root 29 Apr 9 2022 zfgrep - -rwxr-xr-x 1 root root 2081 Apr 9 2022 zforce - -rwxr-xr-x 1 root root 8103 Apr 9 2022 zgrep - -rwxr-xr-x 1 root root 2206 Apr 9 2022 zless - -rwxr-xr-x 1 root root 1842 Apr 9 2022 zmore - -rwxr-xr-x 1 root root 4577 Apr 9 2022 znew -I: user script /srv/workspace/pbuilder/7156/tmp/hooks/D02_print_environment finished + -rwxr-xr-x 1 root root 838488 Apr 24 11:24 bash + -rwxr-xr-x 3 root root 67144 Sep 19 2022 bunzip2 + -rwxr-xr-x 3 root root 67144 Sep 19 2022 bzcat + lrwxrwxrwx 1 root root 6 Sep 19 2022 bzcmp -> bzdiff + -rwxr-xr-x 1 root root 2225 Sep 19 2022 bzdiff + lrwxrwxrwx 1 root root 6 Sep 19 2022 bzegrep -> bzgrep + -rwxr-xr-x 1 root root 4893 Nov 28 2021 bzexe + lrwxrwxrwx 1 root root 6 Sep 19 2022 bzfgrep -> bzgrep + -rwxr-xr-x 1 root root 3775 Sep 19 2022 bzgrep + -rwxr-xr-x 3 root root 67144 Sep 19 2022 bzip2 + -rwxr-xr-x 1 root root 67112 Sep 19 2022 bzip2recover + lrwxrwxrwx 1 root root 6 Sep 19 2022 bzless -> bzmore + -rwxr-xr-x 1 root root 1297 Sep 19 2022 bzmore + -rwxr-xr-x 1 root root 67632 Sep 21 2022 cat + -rwxr-xr-x 1 root root 67676 Sep 21 2022 chgrp + -rwxr-xr-x 1 root root 67644 Sep 21 2022 chmod + -rwxr-xr-x 1 root root 67684 Sep 21 2022 chown + -rwxr-xr-x 1 root root 133532 Sep 21 2022 cp + -rwxr-xr-x 1 root root 132868 Jan 6 03:20 dash + -rwxr-xr-x 1 root root 133220 Sep 21 2022 date + -rwxr-xr-x 1 root root 67732 Sep 21 2022 dd + -rwxr-xr-x 1 root root 68104 Sep 21 2022 df + -rwxr-xr-x 1 root root 133632 Sep 21 2022 dir + -rwxr-xr-x 1 root root 59128 Mar 23 23:02 dmesg + lrwxrwxrwx 1 root root 8 Dec 20 03:33 dnsdomainname -> hostname + lrwxrwxrwx 1 root root 8 Dec 20 03:33 domainname -> hostname + -rwxr-xr-x 1 root root 67560 Sep 21 2022 echo + -rwxr-xr-x 1 root root 41 Jan 25 04:43 egrep + -rwxr-xr-x 1 root root 67548 Sep 21 2022 false + -rwxr-xr-x 1 root root 41 Jan 25 04:43 fgrep + -rwxr-xr-x 1 root root 55748 Mar 23 23:02 findmnt + -rwsr-xr-x 1 root root 26208 Mar 23 22:15 fusermount + -rwxr-xr-x 1 root root 128608 Jan 25 04:43 grep + -rwxr-xr-x 2 root root 2346 Apr 10 2022 gunzip + -rwxr-xr-x 1 root root 6447 Apr 10 2022 gzexe + -rwxr-xr-x 1 root root 64220 Apr 10 2022 gzip + -rwxr-xr-x 1 root root 67032 Dec 20 03:33 hostname + -rwxr-xr-x 1 root root 67720 Sep 21 2022 ln + -rwxr-xr-x 1 root root 35132 Mar 23 23:51 login + -rwxr-xr-x 1 root root 133632 Sep 21 2022 ls + -rwxr-xr-x 1 root root 136808 Mar 23 23:02 lsblk + -rwxr-xr-x 1 root root 67800 Sep 21 2022 mkdir + -rwxr-xr-x 1 root root 67764 Sep 21 2022 mknod + -rwxr-xr-x 1 root root 67596 Sep 21 2022 mktemp + -rwxr-xr-x 1 root root 38504 Mar 23 23:02 more + -rwsr-xr-x 1 root root 38496 Mar 23 23:02 mount + -rwxr-xr-x 1 root root 9824 Mar 23 23:02 mountpoint + -rwxr-xr-x 1 root root 133532 Sep 21 2022 mv + lrwxrwxrwx 1 root root 8 Dec 20 03:33 nisdomainname -> hostname + lrwxrwxrwx 1 root root 14 Apr 3 20:25 pidof -> /sbin/killall5 + -rwxr-xr-x 1 root root 67608 Sep 21 2022 pwd + lrwxrwxrwx 1 root root 4 Apr 24 11:24 rbash -> bash + -rwxr-xr-x 1 root root 67600 Sep 21 2022 readlink + -rwxr-xr-x 1 root root 67672 Sep 21 2022 rm + -rwxr-xr-x 1 root root 67600 Sep 21 2022 rmdir + -rwxr-xr-x 1 root root 67400 Nov 3 2022 run-parts + -rwxr-xr-x 1 root root 133372 Jan 6 09:55 sed + lrwxrwxrwx 1 root root 9 Jun 9 07:24 sh -> /bin/bash + -rwxr-xr-x 1 root root 67584 Sep 21 2022 sleep + -rwxr-xr-x 1 root root 67644 Sep 21 2022 stty + -rwsr-xr-x 1 root root 50800 Mar 23 23:02 su + -rwxr-xr-x 1 root root 67584 Sep 21 2022 sync + -rwxr-xr-x 1 root root 336764 Apr 7 04:25 tar + -rwxr-xr-x 1 root root 67144 Nov 3 2022 tempfile + -rwxr-xr-x 1 root root 133224 Sep 21 2022 touch + -rwxr-xr-x 1 root root 67548 Sep 21 2022 true + -rwxr-xr-x 1 root root 9768 Mar 23 22:15 ulockmgr_server + -rwsr-xr-x 1 root root 22108 Mar 23 23:02 umount + -rwxr-xr-x 1 root root 67572 Sep 21 2022 uname + -rwxr-xr-x 2 root root 2346 Apr 10 2022 uncompress + -rwxr-xr-x 1 root root 133632 Sep 21 2022 vdir + -rwxr-xr-x 1 root root 42608 Mar 23 23:02 wdctl + lrwxrwxrwx 1 root root 8 Dec 20 03:33 ypdomainname -> hostname + -rwxr-xr-x 1 root root 1984 Apr 10 2022 zcat + -rwxr-xr-x 1 root root 1678 Apr 10 2022 zcmp + -rwxr-xr-x 1 root root 6460 Apr 10 2022 zdiff + -rwxr-xr-x 1 root root 29 Apr 10 2022 zegrep + -rwxr-xr-x 1 root root 29 Apr 10 2022 zfgrep + -rwxr-xr-x 1 root root 2081 Apr 10 2022 zforce + -rwxr-xr-x 1 root root 8103 Apr 10 2022 zgrep + -rwxr-xr-x 1 root root 2206 Apr 10 2022 zless + -rwxr-xr-x 1 root root 1842 Apr 10 2022 zmore + -rwxr-xr-x 1 root root 4577 Apr 10 2022 znew +I: user script /srv/workspace/pbuilder/9579/tmp/hooks/D02_print_environment finished -> Attempting to satisfy build-dependencies -> Creating pbuilder-satisfydepends-dummy package Package: pbuilder-satisfydepends-dummy @@ -700,7 +732,7 @@ Get: 503 http://deb.debian.org/debian bookworm/main armhf gnuplot-nox armhf 5.4.4+dfsg1-2+b2 [776 kB] Get: 504 http://deb.debian.org/debian bookworm/main armhf dh-octave-autopkgtest all 1.2.7 [9364 B] Get: 505 http://deb.debian.org/debian bookworm/main armhf dh-octave all 1.2.7 [22.8 kB] -Fetched 170 MB in 23s (7488 kB/s) +Fetched 170 MB in 27s (6292 kB/s) debconf: delaying package configuration, since apt-utils is not installed Selecting previously unselected package netbase. (Reading database ... (Reading database ... 5% (Reading database ... 10% (Reading database ... 15% (Reading database ... 20% (Reading database ... 25% (Reading database ... 30% (Reading database ... 35% (Reading database ... 40% (Reading database ... 45% (Reading database ... 50% (Reading database ... 55% (Reading database ... 60% (Reading database ... 65% (Reading database ... 70% (Reading database ... 75% (Reading database ... 80% (Reading database ... 85% (Reading database ... 90% (Reading database ... 95% (Reading database ... 100% (Reading database ... 19324 files and directories currently installed.) @@ -2771,8 +2803,19 @@ Writing extended state information... Building tag database... -> Finished parsing the build-deps +Reading package lists... +Building dependency tree... +Reading state information... +usrmerge is already the newest version (35). +0 upgraded, 0 newly installed, 0 to remove and 0 not upgraded. I: Building the package -I: Running cd /build/octave-stk-2.7.0/ && 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-stk_2.7.0-3_source.changes +I: user script /srv/workspace/pbuilder/9579/tmp/hooks/A99_set_merged_usr starting +Re-configuring usrmerge... +removed '/etc/unsupported-skip-usrmerge-conversion' +The system has been successfully converted. +I: user script /srv/workspace/pbuilder/9579/tmp/hooks/A99_set_merged_usr finished +hostname: Name or service not known +I: Running cd /build/octave-stk-2.7.0/ && 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-stk_2.7.0-3_source.changes dpkg-buildpackage: info: source package octave-stk dpkg-buildpackage: info: source version 2.7.0-3 dpkg-buildpackage: info: source distribution unstable @@ -2806,473 +2849,2138 @@ make[1]: Entering directory '/build/octave-stk-2.7.0/src' /usr/bin/mkoctfile --verbose --mex --output __stk_dist_pairwise__.mex stk_dist_pairwise.c /usr/bin/mkoctfile --verbose --mex --output __get_column_number__.mex get_column_number.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_dist_pairwise.c -o /tmp/oct-RkI216.o -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG get_column_number.c -o /tmp/oct-VYn4hI.o /usr/bin/mkoctfile --verbose --mex --output __stk_dist_matrixy__.mex stk_dist_matrixy.c +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_dist_pairwise.c -o /tmp/oct-mbETSv.o +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG get_column_number.c -o /tmp/oct-LgjSQo.o +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_dist_matrixy.c -o /tmp/oct-biOssP.o /usr/bin/mkoctfile --verbose --mex --output __stk_filldist_discr_mex__.mex stk_filldist_discr_mex.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_dist_matrixy.c -o /tmp/oct-J6MYJU.o -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_filldist_discr_mex.c -o /tmp/oct-7lvjEU.o -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_filldist_discr_mex__.mex /tmp/oct-7lvjEU.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_dist_matrixy__.mex /tmp/oct-J6MYJU.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __get_column_number__.mex /tmp/oct-VYn4hI.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_dist_pairwise__.mex /tmp/oct-RkI216.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_filldist_discr_mex.c -o /tmp/oct-5MzdVA.o +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_filldist_discr_mex__.mex /tmp/oct-5MzdVA.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_dist_matrixy__.mex /tmp/oct-biOssP.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_dist_pairwise__.mex /tmp/oct-mbETSv.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __get_column_number__.mex /tmp/oct-LgjSQo.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro /usr/bin/mkoctfile --verbose --mex --output __stk_gpquadform_matrixy__.mex stk_gpquadform_matrixy.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_gpquadform_matrixy.c -o /tmp/oct-lnz5J9.o /usr/bin/mkoctfile --verbose --mex --output __stk_isdominated_mex__.mex stk_isdominated_mex.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_isdominated_mex.c -o /tmp/oct-zPjAWE.o +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_gpquadform_matrixy.c -o /tmp/oct-74r4br.o +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_isdominated_mex.c -o /tmp/oct-MLFIWO.o /usr/bin/mkoctfile --verbose --mex --output __stk_mindist_mex__.mex stk_mindist_mex.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_mindist_mex.c -o /tmp/oct-r6XKAd.o +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_mindist_mex.c -o /tmp/oct-UDuBRu.o /usr/bin/mkoctfile --verbose --mex --output __stk_dist_matrixx__.mex stk_dist_matrixx.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_dist_matrixx.c -o /tmp/oct-hB0SkK.o -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_mindist_mex__.mex /tmp/oct-r6XKAd.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_gpquadform_matrixy__.mex /tmp/oct-lnz5J9.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_dist_matrixx.c -o /tmp/oct-lIXnNG.o +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_mindist_mex__.mex /tmp/oct-UDuBRu.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_gpquadform_matrixy__.mex /tmp/oct-74r4br.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro /usr/bin/mkoctfile --verbose --mex --output __stk_gpquadform_matrixx__.mex stk_gpquadform_matrixx.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_gpquadform_matrixx.c -o /tmp/oct-D4wPen.o -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_dist_matrixx__.mex /tmp/oct-hB0SkK.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro /usr/bin/mkoctfile --verbose --mex --output __stk_gpquadform_pairwise__.mex stk_gpquadform_pairwise.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_gpquadform_pairwise.c -o /tmp/oct-cstgwK.o +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_gpquadform_matrixx.c -o /tmp/oct-g40z1Q.o +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_gpquadform_pairwise.c -o /tmp/oct-TEQUNa.o +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_dist_matrixx__.mex /tmp/oct-lIXnNG.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_isdominated_mex__.mex /tmp/oct-MLFIWO.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro /usr/bin/mkoctfile --verbose --mex --output __stk_paretofind_mex__.mex stk_paretofind_mex.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_paretofind_mex.c -o /tmp/oct-5LgEfJ.o -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_gpquadform_matrixx__.mex /tmp/oct-D4wPen.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_gpquadform_pairwise__.mex /tmp/oct-cstgwK.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_isdominated_mex__.mex /tmp/oct-zPjAWE.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_paretofind_mex.c -o /tmp/oct-Kw8XQu.o /usr/bin/mkoctfile --verbose --mex --output __stk_sampling_vdc_rr2__.mex stk_sampling_vdc_rr2.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_sampling_vdc_rr2.c -o /tmp/oct-7Sf94T.o +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_sampling_vdc_rr2.c -o /tmp/oct-Qrutps.o +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_gpquadform_pairwise__.mex /tmp/oct-TEQUNa.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_gpquadform_matrixx__.mex /tmp/oct-g40z1Q.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro /usr/bin/mkoctfile --verbose --mex --output __stk_sampling_sobol_mex__.mex stk_sampling_sobol_mex.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_sampling_sobol_mex.c -o /tmp/oct-mnfpSW.o +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_sampling_sobol_mex.c -o /tmp/oct-7EjgSk.o /usr/bin/mkoctfile --verbose --mex --output __stk_distrib_bivnorm0_cdf__.mex stk_distrib_bivnorm0_cdf.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_distrib_bivnorm0_cdf.c -o /tmp/oct-r81x7K.o -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_paretofind_mex__.mex /tmp/oct-5LgEfJ.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_distrib_bivnorm0_cdf.c -o /tmp/oct-iqdYvA.o +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_paretofind_mex__.mex /tmp/oct-Kw8XQu.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro /usr/bin/mkoctfile --verbose --mex --output __stk_dominatedhv_mex__.mex stk_dominatedhv_mex.c wfg.c -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_dominatedhv_mex.c -o /tmp/oct-PnywwR.o -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_sampling_vdc_rr2__.mex /tmp/oct-7Sf94T.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_sampling_sobol_mex__.mex /tmp/oct-mnfpSW.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_distrib_bivnorm0_cdf__.mex /tmp/oct-r81x7K.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro -gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG wfg.c -o /tmp/oct-jkY6oB.o -g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_dominatedhv_mex__.mex /tmp/oct-PnywwR.o /tmp/oct-jkY6oB.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG stk_dominatedhv_mex.c -o /tmp/oct-jc9Pxr.o +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_sampling_vdc_rr2__.mex /tmp/oct-Qrutps.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_distrib_bivnorm0_cdf__.mex /tmp/oct-iqdYvA.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_sampling_sobol_mex__.mex /tmp/oct-7EjgSk.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro +gcc -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -I. -DMEX_DEBUG wfg.c -o /tmp/oct-u0cycN.o +g++ -I/usr/include/octave-7.3.0/octave/.. -I/usr/include/octave-7.3.0/octave -pthread -fopenmp -g -O2 -ffile-prefix-map=/build/octave-stk-2.7.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __stk_dominatedhv_mex__.mex /tmp/oct-jc9Pxr.o /tmp/oct-u0cycN.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/arm-linux-gnueabihf -shared -Wl,-Bsymbolic -Wl,-z,relro make[1]: Leaving directory '/build/octave-stk-2.7.0/src' copyfile /build/octave-stk-2.7.0/./src/__get_column_number__.mex /build/octave-stk-2.7.0/./src/__stk_dist_matrixx__.mex /build/octave-stk-2.7.0/./src/__stk_dist_matrixy__.mex /build/octave-stk-2.7.0/./src/__stk_dist_pairwise__.mex /build/octave-stk-2.7.0/./src/__stk_distrib_bivnorm0_cdf__.mex /build/octave-stk-2.7.0/./src/__stk_dominatedhv_mex__.mex /build/octave-stk-2.7.0/./src/__stk_filldist_discr_mex__.mex /build/octave-stk-2.7.0/./src/__stk_gpquadform_matrixx__.mex /build/octave-stk-2.7.0/./src/__stk_gpquadform_matrixy__.mex /build/octave-stk-2.7.0/./src/__stk_gpquadform_pairwise__.mex /build/octave-stk-2.7.0/./src/__stk_isdominated_mex__.mex /build/octave-stk-2.7.0/./src/__stk_mindist_mex__.mex /build/octave-stk-2.7.0/./src/__stk_paretofind_mex__.mex /build/octave-stk-2.7.0/./src/__stk_sampling_sobol_mex__.mex /build/octave-stk-2.7.0/./src/__stk_sampling_vdc_rr2__.mex /build/octave-stk-2.7.0/./inst/arm-unknown-linux-gnueabihf-api-v57 For information about changes from previous versions of the stk package, run 'news stk'. dh_octave_check -O--buildsystem=octave Checking package... Checking m files ... -[inst/lm/@stk_lm_cubic/stk_lm_cubic.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_cubic/stk_lm_cubic.m -***** test stk_test_class ('stk_lm_cubic') -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lm/@stk_lm_cubic/feval.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_cubic/feval.m -***** test - n = 15; d = 4; - x = stk_sampling_randunif (n, d); - P = feval (stk_lm_cubic (), x); - assert (isequal (size (P), [n, 1 + d * (11 + d * (6 + d)) / 6])) -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lm/@stk_lm_quadratic/stk_lm_quadratic.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_quadratic/stk_lm_quadratic.m -***** test stk_test_class ('stk_lm_quadratic') -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lm/@stk_lm_quadratic/feval.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_quadratic/feval.m +[inst/param/estim/stk_param_init_lnv.m] +>>>>> /build/octave-stk-2.7.0/inst/param/estim/stk_param_init_lnv.m ***** test - n = 15; d = 4; - x = stk_sampling_randunif (n, d); - P = feval (stk_lm_quadratic (), x); - assert (isequal (size (P), [n, 1 + d * (d + 3) / 2])); + f = @(x)(- (0.8 * x + sin (5 * x + 1) + 0.1 * sin (10 * x))); + ni = 20; + xi = (linspace (-1, 1, ni))' + 0.2 * (randn (ni, 1)); + zi = stk_feval (f, xi); + + model = stk_model (@stk_materncov_iso); + model.param = log ([1; 5/2; 1/0.4]); + model.lognoisevariance = nan; + lnv = stk_param_init_lnv (model, xi, zi); + + assert ((isscalar (lnv)) && (lnv > -30) && (lnv < 30)); 1 test, 1 passed, 0 known failure, 0 skipped -[inst/lm/@stk_lm_null/feval.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_null/feval.m +[inst/param/estim/stk_param_init.m] +>>>>> /build/octave-stk-2.7.0/inst/param/estim/stk_param_init.m ***** test - n = 15; d = 4; - x = stk_sampling_randunif (n, d); - P = feval (stk_lm_null (), x); - assert (isequal (size (P), [n, 0])); -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lm/@stk_lm_null/stk_lm_null.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_null/stk_lm_null.m -***** test stk_test_class ('stk_lm_null') -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lm/stk_lm_polynomial.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/stk_lm_polynomial.m -***** error lm = stk_lm_polynomial (); + xi = (1:10)'; zi = sin (xi); + model = stk_model (@stk_materncov52_iso); + model.param = stk_param_init (model, xi, zi, [1; 10], false); + xt = (1:9)' + 0.5; zt = sin (xt); + zp = stk_predict (model, xi, zi, xt); + assert (sum ((zt - zp.mean) .^ 2) < 1e-3); +***** test % check equivariance of parameter estimates + f = @(x) sin (x); + xi = stk_sampling_regulargrid (10, 1); zi = stk_feval (f, xi); + shift = 1000; scale = 0.01; + model = stk_model (@stk_materncov32_iso); + p1 = stk_param_init (model, xi, zi); + p2 = stk_param_init (model, xi, shift + scale .* zi); + assert (stk_isequal_tolabs (p2(1), p1(1) + log (scale^2), 1e-10)) + assert (stk_isequal_tolabs (p2(2), p1(2), eps)) +***** shared xi, zi, BOX, xt, zt + + f = @(x)(- (0.8 * x + sin (5 * x + 1) + 0.1 * sin (10 * x))); + DIM = 1; % Dimension of the factor space + BOX = [-1.0; 1.0]; % Factor space + + xi = stk_sampling_regulargrid (20, DIM, BOX); % Evaluation points + zi = stk_feval (f, xi); % Evaluation results + + NT = 400; % Number of points in the grid + xt = stk_sampling_regulargrid (NT, DIM, BOX); % Generate a regular grid + zt = stk_feval (f, xt); % Values of f on the grid +***** xtest % This test is known to fail in Debian buster i386 + % because of a known problem in Octave's sqp/qp functions + % http://savannah.gnu.org/bugs/?38875 + model = stk_model (@stk_materncov_iso); + [param0, model.lognoisevariance] = stk_param_init (model, xi, zi, BOX); + model.param = stk_param_estim (model, xi, zi, param0); + zp = stk_predict (model, xi, zi, xt); + assert (max ((zp.mean - zt) .^ 2) < 1e-3) ***** test - lm = stk_lm_polynomial (-1); - assert (isa (lm, 'stk_lm_null')); + model = stk_model (@stk_materncov_aniso); + [param0, model.lognoisevariance] = stk_param_init (model, xi, zi, BOX); + model.param = stk_param_estim (model, xi, zi, param0); + zp = stk_predict (model, xi, zi, xt); + assert (max ((zp.mean - zt) .^ 2) < 1e-3) ***** test - lm = stk_lm_polynomial (0); - assert (isa (lm, 'stk_lm_constant')); + model = stk_model (@stk_materncov32_iso); + [param0, model.lognoisevariance] = stk_param_init (model, xi, zi, BOX); + model.param = stk_param_estim (model, xi, zi, param0); + zp = stk_predict (model, xi, zi, xt); + assert (max ((zp.mean - zt) .^ 2) < 1e-3) ***** test - lm = stk_lm_polynomial (1); - assert (isa (lm, 'stk_lm_affine')); + model = stk_model (@stk_materncov32_aniso); + [param0, model.lognoisevariance] = stk_param_init (model, xi, zi, BOX); + model.param = stk_param_estim (model, xi, zi, param0); + zp = stk_predict (model, xi, zi, xt); + assert (max ((zp.mean - zt) .^ 2) < 1e-3) ***** test - lm = stk_lm_polynomial (2); - assert (isa (lm, 'stk_lm_quadratic')); + model = stk_model (@stk_materncov52_iso); + [param0, model.lognoisevariance] = stk_param_init (model, xi, zi, BOX); + model.param = stk_param_estim (model, xi, zi, param0); + zp = stk_predict (model, xi, zi, xt); + assert (max ((zp.mean - zt) .^ 2) < 1e-3) ***** test - lm = stk_lm_polynomial (3); - assert (isa (lm, 'stk_lm_cubic')); -6 tests, 6 passed, 0 known failure, 0 skipped -[inst/lm/@stk_lm_constant/feval.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_constant/feval.m + model = stk_model (@stk_materncov52_aniso); + [param0, model.lognoisevariance] = stk_param_init (model, xi, zi, BOX); + model.param = stk_param_estim (model, xi, zi, param0); + zp = stk_predict (model, xi, zi, xt); + assert (max ((zp.mean - zt) .^ 2) < 1e-3) ***** test - n = 15; d = 4; - x = stk_sampling_randunif (n, d); - P = feval (stk_lm_constant (), x); - assert (isequal (size (P), [n, 1])); -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lm/@stk_lm_constant/stk_lm_constant.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_constant/stk_lm_constant.m -***** test stk_test_class ('stk_lm_constant') -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lm/@stk_lm_affine/feval.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_affine/feval.m + model = stk_model (@stk_gausscov_iso); + [param0, model.lognoisevariance] = stk_param_init (model, xi, zi, BOX); + model.param = stk_param_estim (model, xi, zi, param0); + zp = stk_predict (model, xi, zi, xt); + assert (max ((zp.mean - zt) .^ 2) < 1e-3) +warning: matrix singular to machine precision, rcond = 1.21281e-17 +warning: called from + linsolve at line 113 column 7 + linsolve at line 42 column 3 + stk_set_righthandside at line 44 column 16 + stk_make_kreq at line 35 column 6 + stk_predict_ at line 112 column 10 + stk_predict at line 104 column 16 + __test__ at line 6 column 5 + test at line 683 column 11 + /tmp/tmp.9XTBtzs9Ty at line 22 column 31 + ***** test - n = 15; d = 4; - x = stk_sampling_randunif (n, d); - P = feval (stk_lm_affine (), x); - assert (isequal (size (P), [n, d + 1])); -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lm/@stk_lm_affine/stk_lm_affine.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_affine/stk_lm_affine.m -***** test stk_test_class ('stk_lm_affine') -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lm/@stk_lm_matrix/stk_lm_matrix.m] ->>>>> /build/octave-stk-2.7.0/inst/lm/@stk_lm_matrix/stk_lm_matrix.m -***** test stk_test_class ('stk_lm_matrix') -***** test %%% dim 1 - data = rand (10, 1); idx = 3:7; - lm = stk_lm_matrix (data); - assert (isa (lm, 'stk_lm_matrix')); - assert (isequal (data(idx, :), feval (lm, idx))); -***** test %%% dim 3 - data = rand (10, 3); idx = 3:7; - lm = stk_lm_matrix (data); - assert (isa (lm, 'stk_lm_matrix')); - assert (isequal (data(idx, :), feval (lm, idx))); -3 tests, 3 passed, 0 known failure, 0 skipped -[inst/sampling/stk_sampling_nesteddesign.m] ->>>>> /build/octave-stk-2.7.0/inst/sampling/stk_sampling_nesteddesign.m -***** shared x, n, dim, box, niter, levels - n = [23; 14; 5; 2]; dim = 2; box = [0, 0; 4, 4]; niter = 10; - levels = [10.1; 15.2; -9.3; 2.4; 17.5]; -***** error x = stk_sampling_nesteddesign (); -***** test x = stk_sampling_nesteddesign (n); -***** test x = stk_sampling_nesteddesign (n, dim); -***** test x = stk_sampling_nesteddesign (n, dim, box); -***** test x = stk_sampling_nesteddesign (n, dim, box, niter); -***** test x = stk_sampling_nesteddesign (n, dim, box, niter, levels); -***** assert ( isequal(size(x), [sum(n), dim + 1]) ); -***** assert ( isa(x, 'stk_dataframe') ); - cn = [0; cumsum(n)]; - for lev = 1:length(n), - y = x( (cn(lev) + 1):(cn(lev + 1)), 1:dim ); - assert (isequal (size (y), [n(lev) dim])); - if lev > 1 - assert ( isequal(z((end - n(lev) + 1):end, :), y) ); - end - if lev == length(n) - assert (stk_is_lhs (y, n(lev), dim, box)); - end - z = y; - end -***** assert (isequal (x.colnames{dim + 1}, 'Level')); - levels = stk_dataframe(levels, {'t'}); - box = stk_hrect(box, {'x1', 'x2', 'x3', 'x4'}); -***** test x = stk_sampling_nesteddesign (n, [], box, [], levels); -***** assert (isequal(x.colnames, {'x1', 'x2', 'x3', 'x4', 't'}) ); + model = stk_model (@stk_gausscov_aniso); + [param0, model.lognoisevariance] = stk_param_init (model, xi, zi, BOX); + model.param = stk_param_estim (model, xi, zi, param0); + zp = stk_predict (model, xi, zi, xt); + assert (max ((zp.mean - zt) .^ 2) < 1e-3) +warning: matrix singular to machine precision, rcond = 8.36758e-18 +warning: called from + linsolve at line 113 column 7 + linsolve at line 42 column 3 + stk_set_righthandside at line 44 column 16 + stk_make_kreq at line 35 column 6 + stk_predict_ at line 112 column 10 + stk_predict at line 104 column 16 + __test__ at line 6 column 5 + test at line 683 column 11 + /tmp/tmp.9XTBtzs9Ty at line 22 column 31 + +***** test % Homoscedastic case / do_estim_lnv = true + model = stk_model (@stk_materncov32_iso); + [model.param, model.lognoisevariance] = ... + stk_param_init (model, xi, zi, BOX, true); + [model.param, model.lognoisevariance] = ... + stk_param_estim (model, xi, zi); + zp = stk_predict (model, xi, zi, xt); + assert (model.lognoisevariance > -inf); + assert (max ((zp.mean - zt) .^ 2) < 1e-3) +warning: This way of setting do_estim_lnv is deprecated +warning: called from + stk_param_init>stk_param_init_ at line 250 column 5 + stk_param_init at line 87 column 18 + __test__ at line 4 column 40 + test at line 683 column 11 + /tmp/tmp.9XTBtzs9Ty at line 22 column 31 + +***** error % Homoscedastic case / do_estim_lnv = false / model.lnv = nan + model = stk_model (@stk_materncov32_iso); + model.lognoisevariance = nan; % not compatible with do_estim_lnv == false + [model.param, model.lognoisevariance] = ... + stk_param_init (model, xi, zi, BOX, false); +***** test % Heteroscedastic case / known noise variance + model = stk_model (@stk_materncov32_iso); + lnv = log ((100 + rand (size (zi))) / 1e6); + model.lognoisevariance = lnv; % here we say that lnv is known + [param0, model.lognoisevariance] = stk_param_init (model, xi, zi, BOX); + model.param = stk_param_estim (model, xi, zi, param0); + zp = stk_predict (model, xi, zi, xt); + assert (isequal (model.lognoisevariance, lnv)); % should be untouched + assert (max ((zp.mean - zt) .^ 2) < 1e-3) +***** shared model, x, z + model = stk_model (@stk_materncov52_iso); + n = 10; x = stk_sampling_regulargrid (n, 1, [0; 1]); z = ones (size (x)); +***** test % Constant response, noiseless model + [param, lnv] = stk_param_init (model, x, z); + assert ((all (isfinite (param))) && (length (param) == 2)); + assert (isequal (lnv, -inf)); +warning: Parameter estimation is impossible with constant-response data. +warning: called from + stk_param_init>paraminit_ at line 351 column 5 + stk_param_init>stk_param_init_ at line 323 column 14 + stk_param_init at line 87 column 18 + __test__ at line 3 column 15 + test at line 683 column 11 + /tmp/tmp.9XTBtzs9Ty at line 22 column 31 + +***** test % Constant response, noisy model + model.lognoisevariance = nan; + [param, lnv] = stk_param_init (model, x, z); + assert ((all (isfinite (param))) && (length (param) == 2)); + assert (isscalar (lnv) && isfinite (lnv)); +warning: Parameter estimation is impossible with constant-response data. +warning: called from + stk_param_init>paraminit_ at line 351 column 5 + stk_param_init>stk_param_init_ at line 323 column 14 + stk_param_init at line 87 column 18 + __test__ at line 4 column 15 + test at line 683 column 11 + /tmp/tmp.9XTBtzs9Ty at line 22 column 31 + +15 tests, 15 passed, 0 known failure, 0 skipped +[inst/param/estim/stk_param_estim.m] +>>>>> /build/octave-stk-2.7.0/inst/param/estim/stk_param_estim.m +***** shared f, xi, zi, NI, param0, param1, model + + f = @(x)(- (0.8 * x + sin (5 * x + 1) + 0.1 * sin (10 * x)) ); + DIM = 1; NI = 20; box = [-1.0; 1.0]; + xi = stk_sampling_regulargrid (NI, DIM, box); + + SIGMA2 = 1.0; % variance parameter + NU = 4.0; % regularity parameter + RHO1 = 0.4; % scale (range) parameter + param0 = log ([SIGMA2; NU; 1/RHO1]); + + model = stk_model (@stk_materncov_iso); +***** test % noiseless + zi = stk_feval (f, xi); + param1 = stk_param_estim (model, xi, zi, param0); + assert (isequal (size (param1), size (param0))) +***** test % same thing, with empty lnv0 (ok) + param2 = stk_param_estim (model, xi, zi, param0, []); + assert (stk_isequal_tolrel (param2, param1, 1e-2)) +***** error % same thing, with lnv0 == NaN (not ok as a starting point) + param2 = stk_param_estim (model, xi, zi, param0, nan); +***** error % same thing, with lnv0 == -inf (not ok as a starting point) + param2 = stk_param_estim (model, xi, zi, param0, -inf); +***** test % same thing, with explicit value for 'criterion' + param2 = stk_param_estim (model, xi, zi, param0, [], @stk_param_relik); + assert (stk_isequal_tolrel (param1, param2, 1e-2)) +***** test % noiseless + zi = stk_feval (f, xi); + param1 = stk_param_estim (model, xi, zi, param0); + param2 = stk_param_estim (model, xi, zi, param0, [], @stk_param_relik); + % We cannot assume a DETERMINISTIC optimization algorithm + % (for some reason, Octave's sqp is not exactly deterministic) + assert (stk_isequal_tolrel (param1, param2, 1e-2)) +***** test % noisy + NOISE_STD_TRUE = 0.1; + NOISE_STD_INIT = 1e-5; + zi = zi + NOISE_STD_TRUE * randn(NI, 1); + model.lognoisevariance = 2 * log(NOISE_STD_INIT); + [param, lnv] = stk_param_estim ... + (model, xi, zi, param0, model.lognoisevariance); +***** error param = stk_param_estim () +***** error param = stk_param_estim (model); +***** error param = stk_param_estim (model, xi); +***** test % Constant response + model = stk_model (@stk_materncov52_iso); + n = 10; x = stk_sampling_regulargrid (n, 1, [0; 1]); z = ones (size (x)); + param = stk_param_estim (model, x, z); + assert ((all (isfinite (param))) && (length (param) == 2)); +warning: Constant-response data: the output of stk_param_estim is likely to be unreliable. +warning: called from + stk_param_estim at line 95 column 5 + __test__ at line 5 column 8 + test at line 683 column 11 + /tmp/tmp.9XTBtzs9Ty at line 30 column 31 + +warning: Parameter estimation is impossible with constant-response data. +warning: called from + stk_param_init>paraminit_ at line 351 column 5 + stk_param_init>stk_param_init_ at line 323 column 14 + stk_param_init at line 87 column 18 + stk_param_estim>provide_starting_point at line 205 column 20 + stk_param_estim at line 110 column 30 + __test__ at line 5 column 8 + test at line 683 column 11 + /tmp/tmp.9XTBtzs9Ty at line 30 column 31 + 11 tests, 11 passed, 0 known failure, 0 skipped -[inst/sampling/stk_sampling_randunif.m] ->>>>> /build/octave-stk-2.7.0/inst/sampling/stk_sampling_randunif.m -***** shared x, n, dim, box - n = 10; dim = 2; box = [0, 0; 2, 2]; -***** error x = stk_sampling_randunif (); -***** test x = stk_sampling_randunif (n); -***** test x = stk_sampling_randunif (n, dim); -***** test x = stk_sampling_randunif (n, dim, box); -***** assert (isa(x, 'stk_dataframe')); -***** assert (isequal (x.colnames, {})); -***** test - cn = {'W', 'H'}; box = stk_hrect (box, cn); - x = stk_sampling_randunif (n, dim, box); - assert (isequal (x.colnames, cn)); -***** test - for dim = 1:5, - x = stk_sampling_randunif(n, dim); - assert(isequal(size(x), [n dim])); - u = double(x); u = u(:); - assert(~any(isnan(u) | isinf(u))); - assert((min(u) >= 0) && (max(u) <= 1)); +[inst/param/estim/stk_param_gls.m] +>>>>> /build/octave-stk-2.7.0/inst/param/estim/stk_param_gls.m +***** shared xi, zi, model, beta, sigma2 + xi = (1:10)'; zi = sin (xi); + model = stk_model (@stk_materncov52_iso); + model.param = [0.0 0.0]; +***** test + model.lm = stk_lm_constant (); + [beta, sigma2] = stk_param_gls (model, xi, zi); +***** assert (stk_isequal_tolabs (beta, 0.1346064, 1e-6)) +***** assert (stk_isequal_tolabs (sigma2, 0.4295288, 1e-6)) +***** test + model.lm = stk_lm_affine (); + [beta, sigma2] = stk_param_gls (model, xi, zi); +***** assert (stk_isequal_tolabs (beta, [0.4728342; -0.0614960], 1e-6)) +***** assert (stk_isequal_tolabs (sigma2, 0.4559431, 1e-6)) +***** test + model.lm = stk_lm_null (); + [beta, sigma2] = stk_param_gls (model, xi, zi); +***** assert (isequal (beta, zeros (0, 1))) +***** assert (stk_isequal_tolabs (sigma2, 0.3977993, 1e-6)) +9 tests, 9 passed, 0 known failure, 0 skipped +[inst/param/estim/stk_param_relik.m] +>>>>> /build/octave-stk-2.7.0/inst/param/estim/stk_param_relik.m +***** shared f, xi, zi, NI, model, C, dC1, dC2 + + f = @(x)(- (0.8 * x(:, 1) + sin (5 * x(:, 2) + 1) ... + + 0.1 * sin (10 * x(:, 3)))); + DIM = 3; NI = 20; box = repmat ([-1.0; 1.0], 1, DIM); + xi = stk_sampling_halton_rr2 (NI, DIM, box); + zi = stk_feval (f, xi); + + SIGMA2 = 1.0; % variance parameter + NU = 4.0; % regularity parameter + RHO1 = 0.4; % scale (range) parameter + + model = stk_model (@stk_materncov_aniso); + model.param = log ([SIGMA2; NU; 1/RHO1 * ones(DIM, 1)]); +***** error [C, dC1, dC2] = stk_param_relik (); +***** error [C, dC1, dC2] = stk_param_relik (model); +***** error [C, dC1, dC2] = stk_param_relik (model, xi); +***** test [C, dC1, dC2] = stk_param_relik (model, xi, zi); +***** test + TOL_REL = 0.01; + assert (stk_isequal_tolrel (C, 21.6, TOL_REL)); + assert (stk_isequal_tolrel (dC1, [4.387 -0.1803 0.7917 0.1392 2.580]', TOL_REL)); + assert (isequal (dC2, [])); +***** shared xi, zi, model, TOL_REL + xi = [-1 -.6 -.2 .2 .6 1]'; + zi = [-0.11 1.30 0.23 -1.14 0.36 -0.37]'; + model = stk_model (@stk_materncov_iso); + model.param = log ([1.0 4.0 2.5]); + model.lognoisevariance = log (0.01); + TOL_REL = 0.01; +***** test % Another simple 1D check + [C, dC1, dC2] = stk_param_relik (model, xi, zi); + assert (stk_isequal_tolrel (C, 6.327, TOL_REL)); + assert (stk_isequal_tolrel (dC1, [0.268 0.0149 -0.636]', TOL_REL)); + assert (stk_isequal_tolrel (dC2, -1.581e-04, TOL_REL)); +***** test % Same 1D test with simple kriging + model.lm = stk_lm_null; + [C, dC1, dC2] = stk_param_relik (model, xi, zi); + assert (stk_isequal_tolrel (C, 7.475, TOL_REL)); + assert (stk_isequal_tolrel (dC1, [0.765 0.0238 -1.019]', TOL_REL)); + assert (stk_isequal_tolrel (dC2, 3.0517e-03, TOL_REL)); +***** test % Check the gradient on a 2D test case + + f = @stk_testfun_braninhoo; + DIM = 2; + BOX = [[-5; 10], [0; 15]]; + NI = 20; + TOL_REL = 1e-2; + DELTA = 1e-6; + + model = stk_model (@stk_materncov52_iso, DIM); + model.param = [1 1]; + + xi = stk_sampling_halton_rr2 (NI, DIM, BOX); + zi = stk_feval (f, xi); + + for range = [0.3 2 10] + model.param(2) = - log (range); + for diff = 1:2 + assert (stk_test_critgrad ... + (@stk_param_relik, model, xi, zi, diff, 1e-6)); + end end 8 tests, 8 passed, 0 known failure, 0 skipped -[inst/sampling/stk_sampling_maximinlhs.m] ->>>>> /build/octave-stk-2.7.0/inst/sampling/stk_sampling_maximinlhs.m -***** shared x, n, dim, box, niter - n = 20; dim = 2; box = [0, 0; 1, 1]; niter = 1; -***** error x = stk_sampling_maximinlhs (); -***** test x = stk_sampling_maximinlhs (n); -***** test x = stk_sampling_maximinlhs (n, dim); -***** test x = stk_sampling_maximinlhs (n, dim, box); -***** test x = stk_sampling_maximinlhs (n, dim, box, niter); -***** assert (isa (x, 'stk_dataframe')); -***** assert (isequal (x.colnames, {})); -***** test - cn = {'W', 'H'}; box = stk_hrect (box, cn); - x = stk_sampling_maximinlhs (n, dim, box); - assert (isequal (x.colnames, cn)); -***** test - for dim = 1:5, - x = stk_sampling_randomlhs (n, dim); - assert (isequal (size (x), [n dim])); - u = double (x); u = u(:); - assert (~ any (isnan (u) | isinf (u))); - assert ((min (u) >= 0) && (max (u) <= 1)); - assert (stk_is_lhs (x, n, dim)); +[inst/covfcs/stk_noisecov.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_noisecov.m +***** shared ni, lognoisevariance, diff + ni = 5; + lognoisevariance = 0.0; + diff = -1; +***** error K = stk_noisecov (); +***** error K = stk_noisecov (ni); +***** test K = stk_noisecov (ni, lognoisevariance); +***** test K = stk_noisecov (ni, lognoisevariance, diff); +***** test K = stk_noisecov (ni, lognoisevariance, diff, true); +5 tests, 5 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_sphcov_iso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_sphcov_iso.m +***** shared param, x, y + dim = 1; + param = log ([1.0; 2.5]); + x = stk_sampling_randunif (5, dim); + y = stk_sampling_randunif (5, dim); +***** error K = stk_sphcov_iso ([param; 1.234], x, y); +***** error stk_sphcov_iso (); +***** error stk_sphcov_iso (param); +***** error stk_sphcov_iso (param, x); +***** test stk_sphcov_iso (param, x, y); +***** test stk_sphcov_iso (param, x, y, -1); +***** test stk_sphcov_iso (param, x, y, -1, false); +***** error stk_sphcov_iso (param, x, y, -2); +***** test stk_sphcov_iso (param, x, y, -1); +***** error stk_sphcov_iso (param, x, y, 0); +***** test stk_sphcov_iso (param, x, y, 1); +***** test stk_sphcov_iso (param, x, y, 2); +***** error stk_sphcov_iso (param, x, y, 3); +***** error stk_sphcov_iso (param, x, y, nan); +***** error stk_sphcov_iso (param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 2.5]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_sphcov_iso (param, x, y); + K2 = stk_sphcov_iso (param, x, y, -1); + assert (isequal (size (K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:2, + dK = stk_sphcov_iso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif (n, dim); + y = stk_sampling_randunif (n, dim); + + K1 = stk_sphcov_iso (param, x, y); + K2 = stk_sphcov_iso (param, x, y, -1, true); + assert (isequal (size (K1), [n n])); + assert (stk_isequal_tolabs (K2, diag (K1))); + + for i = 1:2, + dK1 = stk_sphcov_iso (param, x, y, i); + dK2 = stk_sphcov_iso (param, x, y, i, true); + assert (isequal (size (dK1), [n n])); + assert (stk_isequal_tolabs (dK2, diag (dK1))); + end +18 tests, 18 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_expcov_aniso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_expcov_aniso.m +***** shared param, x, y, K1, K2, K3 + dim = 1; + param = log ([1.0; 2.5]); + x = stk_sampling_randunif (5, dim); + y = stk_sampling_randunif (6, dim); +***** error K0 = stk_expcov_aniso (); +***** error K0 = stk_expcov_aniso (param); +***** error K0 = stk_expcov_aniso (param, x); +***** test K1 = stk_expcov_aniso (param, x, y); +***** test K2 = stk_expcov_aniso (param, x, y, -1); +***** test K3 = stk_expcov_aniso (param, x, y, -1, false); +***** assert (isequal (K1, K2)); +***** assert (isequal (K1, K3)); +***** test % df versus ordinary array + u = double (x); v = double (y); + K1 = stk_expcov_aniso (param, u, v, -1); + K2 = stk_expcov_aniso (param, stk_dataframe (u), stk_dataframe (v), -1); +***** error stk_expcov_aniso (param, x, y, -2); +***** test stk_expcov_aniso (param, x, y, -1); +***** error stk_expcov_aniso (param, x, y, 0); +***** test stk_expcov_aniso (param, x, y, 1); +***** test stk_expcov_aniso (param, x, y, 2); +***** error stk_expcov_aniso (param, x, y, 3); +***** error stk_expcov_aniso (param, x, y, nan); +***** error stk_expcov_aniso (param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 2.5; 2.4; 2.6]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_expcov_aniso (param, x, y); + K2 = stk_expcov_aniso (param, x, y, -1); + assert (isequal (size(K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:(dim + 1), + dK = stk_expcov_aniso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif (n, dim); + y = stk_sampling_randunif (n, dim); + + K1 = stk_expcov_aniso (param, x, y); + K2 = stk_expcov_aniso (param, x, y, -1, true); + assert (isequal (size (K1), [n n])); + assert (stk_isequal_tolabs (K2, diag (K1))); + + for i = 1:(dim + 1), + dK1 = stk_expcov_aniso (param, x, y, i); + dK2 = stk_expcov_aniso (param, x, y, i, true); + assert (isequal (size (dK1), [n n])); + assert (stk_isequal_tolabs (dK2, diag (dK1))); + end +20 tests, 20 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_materncov32_iso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_materncov32_iso.m +***** shared param, x, y + dim = 1; + param = log ([1.0; 2.5]); + x = stk_sampling_randunif (5, dim); + y = stk_sampling_randunif (5, dim); +***** error K = stk_materncov32_iso ([param; 1.234], x, y); +***** error stk_materncov32_iso (); +***** error stk_materncov32_iso (param); +***** error stk_materncov32_iso (param, x); +***** test stk_materncov32_iso (param, x, y); +***** test stk_materncov32_iso (param, x, y, -1); +***** test stk_materncov32_iso (param, x, y, -1, false); +***** error stk_materncov32_iso (param, x, y, -2); +***** test stk_materncov32_iso (param, x, y, -1); +***** error stk_materncov32_iso (param, x, y, 0); +***** test stk_materncov32_iso (param, x, y, 1); +***** test stk_materncov32_iso (param, x, y, 2); +***** error stk_materncov32_iso (param, x, y, 3); +***** error stk_materncov32_iso (param, x, y, nan); +***** error stk_materncov32_iso (param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 2.5]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_materncov32_iso (param, x, y); + K2 = stk_materncov32_iso (param, x, y, -1); + assert (isequal (size (K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:2, + dK = stk_materncov32_iso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif (n, dim); + y = stk_sampling_randunif (n, dim); + + K1 = stk_materncov32_iso (param, x, y); + K2 = stk_materncov32_iso (param, x, y, -1, true); + assert (isequal (size (K1), [n n])); + assert (stk_isequal_tolabs (K2, diag (K1))); + + for i = 1:2, + dK1 = stk_materncov32_iso (param, x, y, i); + dK2 = stk_materncov32_iso (param, x, y, i, true); + assert (isequal (size (dK1), [n n])); + assert (stk_isequal_tolabs (dK2, diag (dK1))); + end +18 tests, 18 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_materncov_aniso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_materncov_aniso.m +***** shared param, x, y + dim = 1; + param = log ([1.0; 1.5; 2.8]); + x = stk_sampling_randunif(5, dim); + y = stk_sampling_randunif(5, dim); +***** error stk_materncov_aniso(); +***** error stk_materncov_aniso(param); +***** error stk_materncov_aniso(param, x); +***** test stk_materncov_aniso(param, x, y); +***** test stk_materncov_aniso(param, x, y, -1); +***** test stk_materncov_aniso(param, x, y, -1, false); +***** error stk_materncov_aniso(param, x, y, -2); +***** test stk_materncov_aniso(param, x, y, -1); +***** error stk_materncov_aniso(param, x, y, 0); +***** test stk_materncov_aniso(param, x, y, 1); +***** test stk_materncov_aniso(param, x, y, 2); +***** test stk_materncov_aniso(param, x, y, 3); +***** error stk_materncov_aniso(param, x, y, 4); +***** error stk_materncov_aniso(param, x, y, nan); +***** error stk_materncov_aniso(param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 1.5; 2.8; 2.7; 2.9]); + nx = 4; ny = 10; + x = stk_sampling_randunif(nx, dim); + y = stk_sampling_randunif(ny, dim); +***** test + K1 = stk_materncov_aniso(param, x, y); + K2 = stk_materncov_aniso(param, x, y, -1); + assert(isequal(size(K1), [nx ny])); + assert(stk_isequal_tolabs(K1, K2)); +***** test + for i = 1:(dim+2), + dK = stk_materncov_aniso(param, x, y, i); + assert(isequal(size(dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif(n, dim); + y = stk_sampling_randunif(n, dim); + + K1 = stk_materncov_aniso(param, x, y); + K2 = stk_materncov_aniso(param, x, y, -1, true); + assert(isequal(size(K1), [n n])); + assert(stk_isequal_tolabs(K2, diag(K1))); + + for i = 1:(dim+2), + dK1 = stk_materncov_aniso(param, x, y, i); + dK2 = stk_materncov_aniso(param, x, y, i, true); + assert(isequal(size(dK1), [n n])); + assert(stk_isequal_tolabs(dK2, diag(dK1))); + end +18 tests, 18 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_gausscov_aniso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_gausscov_aniso.m +***** shared param, x, y, K1, K2, K3 + dim = 1; + param = log ([1.0; 2.5]); + x = stk_sampling_randunif (5, dim); + y = stk_sampling_randunif (6, dim); +***** error K0 = stk_gausscov_aniso (); +***** error K0 = stk_gausscov_aniso (param); +***** error K0 = stk_gausscov_aniso (param, x); +***** test K1 = stk_gausscov_aniso (param, x, y); +***** test K2 = stk_gausscov_aniso (param, x, y, -1); +***** test K3 = stk_gausscov_aniso (param, x, y, -1, false); +***** assert (isequal (K1, K2)); +***** assert (isequal (K1, K3)); +***** test % df versus ordinary array + u = double (x); v = double (y); + K1 = stk_gausscov_aniso (param, u, v, -1); + K2 = stk_gausscov_aniso (param, stk_dataframe (u), stk_dataframe (v), -1); +***** error stk_gausscov_aniso (param, x, y, -2); +***** test stk_gausscov_aniso (param, x, y, -1); +***** error stk_gausscov_aniso (param, x, y, 0); +***** test stk_gausscov_aniso (param, x, y, 1); +***** test stk_gausscov_aniso (param, x, y, 2); +***** error stk_gausscov_aniso (param, x, y, 3); +***** error stk_gausscov_aniso (param, x, y, nan); +***** error stk_gausscov_aniso (param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 2.5; 2.4; 2.6]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_gausscov_aniso (param, x, y); + K2 = stk_gausscov_aniso (param, x, y, -1); + assert (isequal (size (K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:(dim + 1), + dK = stk_gausscov_aniso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif (n, dim); + y = stk_sampling_randunif (n, dim); + + K1 = stk_gausscov_aniso (param, x, y); + K2 = stk_gausscov_aniso (param, x, y, -1, true); + assert (isequal (size (K1), [n n])); + assert (stk_isequal_tolabs (K2, diag (K1))); + + for i = 1:(dim + 1), + dK1 = stk_gausscov_aniso (param, x, y, i); + dK2 = stk_gausscov_aniso (param, x, y, i, true); + assert (isequal (size (dK1), [n n])); + assert (stk_isequal_tolabs (dK2, diag (dK1))); + end +20 tests, 20 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_materncov52_iso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_materncov52_iso.m +***** shared param, x, y + dim = 1; + param = log ([1.0; 2.5]); + x = stk_sampling_randunif (5, dim); + y = stk_sampling_randunif (5, dim); +***** error K = stk_materncov52_iso ([param; 1.234], x, y); +***** error stk_materncov52_iso (); +***** error stk_materncov52_iso (param); +***** error stk_materncov52_iso (param, x); +***** test stk_materncov52_iso (param, x, y); +***** test stk_materncov52_iso (param, x, y, -1); +***** test stk_materncov52_iso (param, x, y, -1, false); +***** error stk_materncov52_iso (param, x, y, -2); +***** test stk_materncov52_iso (param, x, y, -1); +***** error stk_materncov52_iso (param, x, y, 0); +***** test stk_materncov52_iso (param, x, y, 1); +***** test stk_materncov52_iso (param, x, y, 2); +***** error stk_materncov52_iso (param, x, y, 3); +***** error stk_materncov52_iso (param, x, y, nan); +***** error stk_materncov52_iso (param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 2.5]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_materncov52_iso (param, x, y); + K2 = stk_materncov52_iso (param, x, y, -1); + assert (isequal (size (K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:2, + dK = stk_materncov52_iso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif (n, dim); + y = stk_sampling_randunif (n, dim); + + K1 = stk_materncov52_iso (param, x, y); + K2 = stk_materncov52_iso (param, x, y, -1, true); + assert (isequal (size (K1), [n n])); + assert (stk_isequal_tolabs (K2, diag (K1))); + + for i = 1:2, + dK1 = stk_materncov52_iso (param, x, y, i); + dK2 = stk_materncov52_iso (param, x, y, i, true); + assert (isequal (size (dK1), [n n])); + assert (stk_isequal_tolabs (dK2, diag (dK1))); + end +18 tests, 18 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_gausscov_iso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_gausscov_iso.m +***** shared param, x, y + dim = 1; + param = log ([1.0; 2.5]); + x = stk_sampling_randunif (5, dim); + y = stk_sampling_randunif (5, dim); +***** error K = stk_gausscov_iso ([param; 1.234], x, y); +***** error stk_gausscov_iso (); +***** error stk_gausscov_iso (param); +***** error stk_gausscov_iso (param, x); +***** test stk_gausscov_iso (param, x, y); +***** test stk_gausscov_iso (param, x, y, -1); +***** test stk_gausscov_iso (param, x, y, -1, false); +***** error stk_gausscov_iso (param, x, y, -2); +***** test stk_gausscov_iso (param, x, y, -1); +***** error stk_gausscov_iso (param, x, y, 0); +***** test stk_gausscov_iso (param, x, y, 1); +***** test stk_gausscov_iso (param, x, y, 2); +***** error stk_gausscov_iso (param, x, y, 3); +***** error stk_gausscov_iso (param, x, y, nan); +***** error stk_gausscov_iso (param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 2.5]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_gausscov_iso (param, x, y); + K2 = stk_gausscov_iso (param, x, y, -1); + assert (isequal (size (K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:2, + dK = stk_gausscov_iso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif (n, dim); + y = stk_sampling_randunif (n, dim); + + K1 = stk_gausscov_iso (param, x, y); + K2 = stk_gausscov_iso (param, x, y, -1, true); + assert (isequal (size (K1), [n n])); + assert (stk_isequal_tolabs (K2, diag (K1))); + + for i = 1:2, + dK1 = stk_gausscov_iso (param, x, y, i); + dK2 = stk_gausscov_iso (param, x, y, i, true); + assert (isequal (size (dK1), [n n])); + assert (stk_isequal_tolabs (dK2, diag (dK1))); + end +18 tests, 18 passed, 0 known failure, 0 skipped +[inst/covfcs/rbf/stk_rbf_spherical.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/rbf/stk_rbf_spherical.m +***** shared h, diff + h = 1.0; diff = -1; +***** error stk_rbf_spherical (); +***** test stk_rbf_spherical (h); +***** test stk_rbf_spherical (h, diff); +***** test %% h = 0.0 => correlation = 1.0 + x = stk_rbf_spherical (0.0); + assert (stk_isequal_tolrel (x, 1.0, 1e-8)); +***** test %% check derivative numerically + h = [-1 -0.5 -0.1 0.1 0.5 1]; delta = 1e-9; + d1 = (stk_rbf_spherical (h + delta) - stk_rbf_spherical (h)) / delta; + d2 = stk_rbf_spherical (h, 1); + assert (stk_isequal_tolabs (d1, d2, 1e-4)); +***** assert (stk_rbf_spherical (inf) == 0) +6 tests, 6 passed, 0 known failure, 0 skipped +[inst/covfcs/rbf/stk_rbf_gauss.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/rbf/stk_rbf_gauss.m +***** shared h, diff + h = 1.0; diff = -1; +***** error stk_rbf_gauss (); +***** test stk_rbf_gauss (h); +***** test stk_rbf_gauss (h, diff); +***** test % h = 0.0 => correlation = 1.0 + x = stk_rbf_gauss (0.0); + assert (stk_isequal_tolrel (x, 1.0, 1e-8)); +4 tests, 4 passed, 0 known failure, 0 skipped +[inst/covfcs/rbf/stk_rbf_matern52.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/rbf/stk_rbf_matern52.m +***** shared h, diff + h = 1.0; diff = -1; +***** error stk_rbf_matern52 (); +***** test stk_rbf_matern52 (h); +***** test stk_rbf_matern52 (h, diff); +***** test %% h = 0.0 => correlation = 1.0 + x = stk_rbf_matern52 (0.0); + assert (stk_isequal_tolrel (x, 1.0, 1e-8)); +***** test %% consistency with stk_rbf_matern: function values + for h = 0.1:0.1:2.0, + x = stk_rbf_matern (5/2, h); + y = stk_rbf_matern52 (h); + assert (stk_isequal_tolrel (x, y, 1e-8)); + end +***** test %% consistency with stk_rbf_matern: derivatives + for h = 0.1:0.1:2.0, + x = stk_rbf_matern (5/2, h, 2); + y = stk_rbf_matern52 (h, 1); + assert (stk_isequal_tolrel (x, y, 1e-8)); + end +***** assert (stk_rbf_matern52 (inf) == 0) +7 tests, 7 passed, 0 known failure, 0 skipped +[inst/covfcs/rbf/stk_rbf_matern32.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/rbf/stk_rbf_matern32.m +***** shared h, diff + h = 1.0; diff = -1; +***** error stk_rbf_matern32 (); +***** test stk_rbf_matern32 (h); +***** test stk_rbf_matern32 (h, diff); +***** test %% h = 0.0 => correlation = 1.0 + x = stk_rbf_matern32 (0.0); + assert (stk_isequal_tolrel (x, 1.0, 1e-8)); +***** test %% consistency with stk_rbf_matern: function values + for h = 0.1:0.1:2.0, + x = stk_rbf_matern (3/2, h); + y = stk_rbf_matern32 (h); + assert (stk_isequal_tolrel (x, y, 1e-8)); + end +***** test %% consistency with stk_rbf_matern: derivatives + for h = 0.1:0.1:2.0, + x = stk_rbf_matern (3/2, h, 2); + y = stk_rbf_matern32 (h, 1); + assert (stk_isequal_tolrel (x, y, 1e-8)); + end +***** assert (stk_rbf_matern32 (inf) == 0) +7 tests, 7 passed, 0 known failure, 0 skipped +[inst/covfcs/rbf/stk_rbf_matern.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/rbf/stk_rbf_matern.m +***** shared nu, h, diff + nu = 1.0; h = 1.0; diff = -1; +***** error stk_rbf_matern (); +***** error stk_rbf_matern (nu); +***** test stk_rbf_matern (nu, h); +***** test stk_rbf_matern (nu, h, diff); +***** test %% h = 0.0 => correlation = 1.0 + for nu = 0.1:0.2:5.0, + x = stk_rbf_matern (nu, 0.0); + assert (stk_isequal_tolrel (x, 1.0, 1e-8)); end -9 tests, 9 passed, 0 known failure, 0 skipped -[inst/sampling/stk_sampling_olhs.m] ->>>>> /build/octave-stk-2.7.0/inst/sampling/stk_sampling_olhs.m -***** shared x, n, d, box, permut - n = 5; d = 2; box = [0 0; 1, 1]; permut = 1:2; -***** error x = stk_sampling_olhs(); -***** test x = stk_sampling_olhs(n); -***** test x = stk_sampling_olhs(n, d); -***** test x = stk_sampling_olhs(n, d, box); -***** test x = stk_sampling_olhs(n, d, box, permut); -***** error x = stk_sampling_olhs(n, d, box, permut, pi); -***** assert (isa (x, 'stk_dataframe')); -***** assert (isequal (x.colnames, {})); -***** test - cn = {'W', 'H'}; box = stk_hrect (box, cn); - x = stk_sampling_olhs (n, d, box); - assert (isequal (x.colnames, cn)); -***** test - for r = 1:5 - - n = 2 ^ (r + 1) + 1; d = 2 * r; - x = stk_sampling_olhs (n, d); - - assert (isequal (size (x), [n d])); - - box = repmat ([-1; 1], 1, d); - assert (stk_is_lhs (x, n, d, box)); - - x = double (x); w = x' * x; - assert (stk_isequal_tolabs (w / w(1,1), eye (d))); - - end -10 tests, 10 passed, 0 known failure, 0 skipped -[inst/sampling/stk_sampling_halton_rr2.m] ->>>>> /build/octave-stk-2.7.0/inst/sampling/stk_sampling_halton_rr2.m -***** error stk_sampling_halton_rr2 () % nargin < 1 -***** test - n = 300; d = 25; - x = stk_sampling_halton_rr2 (n, d); - assert (isequal (size (x), [n d])) -***** test - x = stk_sampling_halton_rr2 (1000, 3); - y = double (x(end, :)); - yref = [0.9052734375 0.028349336991312 0.74848]; - assert (stk_isequal_tolrel (y, yref, 1e-13)); -***** test - dim = 2; box = stk_hrect (dim); - x = stk_sampling_halton_rr2 (5, dim, box); - assert (isequal (x.colnames, {})); -***** test - dim = 2; cn = {'W', 'H'}; box = stk_hrect (dim, cn); - x = stk_sampling_halton_rr2 (5, dim, box); - assert (isequal (x.colnames, cn)); 5 tests, 5 passed, 0 known failure, 0 skipped -[inst/sampling/@stk_sampcrit_eqi/stk_sampcrit_eqi.m] ->>>>> /build/octave-stk-2.7.0/inst/sampling/@stk_sampcrit_eqi/stk_sampcrit_eqi.m -***** shared F, M, EQI - M = stk_model_gpposterior (stk_model, [1 2 3]', [1.234 3 2]'); -warning: Something went wrong during the optimization -crit0 = 2.594852, crit_opt = 2.594852: crit0 < crit_opt -***** test F = stk_sampcrit_eqi () % ending ";" omitted on purpose, to test disp - -F = -| -| model: -- (not instantiated) -| quantile_order: 0.5 -| point_batch_size: 1 -| current_minimum: Inf -| - -***** assert (isempty (F.model)) -***** assert (isempty (get (F, 'model'))) -***** assert (F.quantile_order == 0.5) -***** assert (get (F, 'quantile_order') == 0.5) -***** assert (F.current_minimum == +inf) -***** assert (get (F, 'current_minimum') == +inf) -***** error F.toto -***** error get (F, 'toto') -***** error F.current_min = 1.234; % read-only -***** error F = set (F, 'current_min', 1.234); % read-only -***** error F.toto = 1.234; % field does not exist -***** error F = set (F, 'toto', 1.234); % field does not exist -***** test F.quantile_order = 0.9; assert (F.quantile_order == 0.9) -***** test F = set (F, 'quantile_order', 0.8); assert (F.quantile_order == 0.8) -***** error F.quantile_order = 1.1; -***** error F.quantile_order = -0.1; -***** error F.quantile_order = [1 2]; -***** error F.current_minimum = 3.333; % read-only -***** error F.quantile_value = 2.222; % read-only -***** error EQI = feval (F, 1.0); -***** test F = stk_sampcrit_eqi (); F.model = M; - assert (~ isempty (F.model)); -***** test F = stk_sampcrit_eqi (); F = set (F, 'model', M); - assert (~ isempty (F.model)); -***** test F.model = []; % remove model - assert (isempty (F.model)); - assert (F.current_minimum == +inf); -***** test F = stk_sampcrit_eqi (M) % ending ";" omitted on purpose, to test disp - -F = -| -| model: -| quantile_order: 0.5 -| point_batch_size: 1 -| current_minimum: 1.234 -| - -***** assert (isequal (F.model, M)) -***** assert (stk_isequal_tolrel (F.current_minimum, 1.234, 10 * eps)); -***** test EQI = feval (F, [1.0; 1.1; 1.2]); -***** assert (isequal (size (EQI), [3 1])) -***** assert (all (EQI >= 0)) -***** test F.quantile_order = 0.9; assert (F.quantile_order == 0.9) -***** shared F, M, EQI - prior_model = stk_model (); - prior_model.lognoisevariance = 0.678; - M = stk_model_gpposterior (prior_model, [1 2 3]', [1.234 3 2]'); -***** test F = stk_sampcrit_eqi (M); -***** assert (isequal (F.model, M)) -***** assert (stk_isequal_tolrel (F.current_minimum, 2.077997, 1e-5)); -***** test EQI = feval (F, [1.0; 1.1; 1.2]); -***** assert (isequal (size (EQI), [3 1])) -***** assert (all (EQI >= 0)) -***** test F.quantile_order = 0.9; assert (F.quantile_order == 0.9) -***** shared F -***** test F = stk_sampcrit_eqi (stk_model ()); -***** assert (F.current_minimum == +inf); -***** error feval (F, 1.0); -***** shared F, M - M = stk_model_gpposterior (stk_model (), [1 2 3]', [1.234 3 2]'); -warning: Something went wrong during the optimization -crit0 = 2.594852, crit_opt = 2.594852: crit0 < crit_opt -***** error F = stk_sampcrit_eqi (M, [], 0); -***** error F = stk_sampcrit_eqi (M, [], 1.5); -***** error F = stk_sampcrit_eqi (M, [], nan); -***** error F = stk_sampcrit_eqi (M, [], [10 20]); -***** test F = stk_sampcrit_eqi (M, [], 10); -***** assert (isequal (F.quantile_order, 0.5)); -***** assert (isequal (F.point_batch_size, 10)); -***** error F = stk_sampcrit_eqi (M, 0.8, 0); -***** error F = stk_sampcrit_eqi (M, 0.8, 1.5); -***** error F = stk_sampcrit_eqi (M, 0.8, nan); -***** error F = stk_sampcrit_eqi (M, 0.8, [10 20]); -***** test F = stk_sampcrit_eqi (M, 0.8, 5); -***** assert (isequal (F.quantile_order, 0.8)); -***** assert (isequal (F.point_batch_size, 5)); -***** shared F, M, EQI - prior_model = stk_model (); - prior_model.lognoisevariance = 0.678; - M = stk_model_gpposterior (prior_model, [1 2 3]', [1.234 3 2]'); - F = stk_sampcrit_eqi (M); -***** test F.point_batch_size = 10; % numeric -***** assert (isequal (F.point_batch_size, 10)) -***** test EQI = feval (F, [1.0; 1.1; 1.2]); -***** test F.point_batch_size = @(x, n) 100 - n; % function handle -***** assert (isa (F.point_batch_size, 'function_handle')) -***** test EQI = feval (F, [1.0; 1.1; 1.2]); -***** test F.point_batch_size = 'sin'; % char -***** assert (isa (F.point_batch_size, 'function_handle')) -63 tests, 63 passed, 0 known failure, 0 skipped -[inst/sampling/stk_sampling_randomlhs.m] ->>>>> /build/octave-stk-2.7.0/inst/sampling/stk_sampling_randomlhs.m -***** shared x, n, dim, box - n = 10; dim = 2; box = [0, 0; 1, 1]; -***** error x = stk_sampling_randomlhs (); -***** test x = stk_sampling_randomlhs (n); -***** test x = stk_sampling_randomlhs (n, dim); -***** test x = stk_sampling_randomlhs (n, dim, box); -***** assert (isa(x, 'stk_dataframe')); -***** assert (isequal (x.colnames, {})); -***** test - cn = {'W', 'H'}; box = stk_hrect (box, cn); - x = stk_sampling_randomlhs (n, dim, box); - assert (isequal (x.colnames, cn)); -***** test - for dim = 1:5, - x = stk_sampling_randomlhs(n, dim); - assert(isequal(size(x), [n dim])); - u = double(x); u = u(:); - assert(~any(isnan(u) | isinf(u))); - assert((min(u) >= 0) && (max(u) <= 1)); - assert(stk_is_lhs(x, n, dim)); +[inst/covfcs/rbf/stk_rbf_exponential.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/rbf/stk_rbf_exponential.m +***** shared h, diff + h = 1.0; diff = -1; +***** error stk_rbf_exponential (); +***** test stk_rbf_exponential (h); +***** test stk_rbf_exponential (h, diff); +***** test %% h = 0.0 => correlation = 1.0 + x = stk_rbf_exponential (0.0); + assert (stk_isequal_tolrel (x, 1.0, 1e-8)); +***** test %% check derivative numerically + h = [-1 -0.5 -0.1 0.1 0.5 1]; delta = 1e-9; + d1 = (stk_rbf_exponential (h + delta) - stk_rbf_exponential (h)) / delta; + d2 = stk_rbf_exponential (h, 1); + assert (stk_isequal_tolabs (d1, d2, 1e-4)); +***** test %% consistency with stk_rbf_matern: function values + for h = 0.1:0.1:2.0, + x = stk_rbf_matern (1/2, h); + y = stk_rbf_exponential (h); + assert (stk_isequal_tolrel (x, y, 1e-8)); + end +***** test %% consistency with stk_rbf_matern: derivatives + for h = 0.1:0.1:2.0, + x = stk_rbf_matern (1/2, h, 2); + y = stk_rbf_exponential (h, 1); + assert (stk_isequal_tolrel (x, y, 1e-8)); end +***** assert (stk_rbf_exponential (inf) == 0) 8 tests, 8 passed, 0 known failure, 0 skipped -[inst/sampling/stk_sampcrit_ehvi_eval.m] ->>>>> /build/octave-stk-2.7.0/inst/sampling/stk_sampcrit_ehvi_eval.m -***** shared zr, zi - zr = [1 1]; - zi = [0.25 0.75; 0.5 0.5; 0.75 0.25]; -***** test % no improvement (1 computation) - zp_mean = [0.6 0.6]; zp_std = [0 0]; - EHVI = stk_sampcrit_ehvi_eval (zp_mean, zp_std, zi, zr); - assert (stk_isequal_tolabs (EHVI, 0, 1e-12)); -***** test % guaranteed improvement (1 computation) - zp_mean = [0 0]; zp_std = [0 0]; - EHVI = stk_sampcrit_ehvi_eval (zp_mean, zp_std, zi, zr); - assert (stk_isequal_tolabs (EHVI, 10 * 0.25 ^ 2)); -***** test % no improvement again (2 computations) - zp_mean = [0.5 0.5; 0.6 0.6]; zp_std = [0 0; 0 0]; - EHVI = stk_sampcrit_ehvi_eval (zp_mean, zp_std, zi, zr); - assert (stk_isequal_tolabs (EHVI, [0; 0], 1e-12)); -***** test % no observation -> EHVI wrt zr - zp_mean = [0.6 0.6]; zp_std = 0.01 * [1 1]; zi = []; - EHVI = stk_sampcrit_ehvi_eval (zp_mean, zp_std, zi, zr); - assert (stk_isequal_tolabs (EHVI, (1 - 0.6)^2, 1e-12)); -***** test % no observation below zr -> EHVI wrt zr - zp_mean = [0.6 0.6]; zp_std = 0.01 * [1 1]; zi = [2 2]; - EHVI = stk_sampcrit_ehvi_eval (zp_mean, zp_std, zi, zr); - assert (stk_isequal_tolabs (EHVI, (1 - 0.6)^2, 1e-12)); -5 tests, 5 passed, 0 known failure, 0 skipped -[inst/sampling/stk_sampling_nestedlhs.m] ->>>>> /build/octave-stk-2.7.0/inst/sampling/stk_sampling_nestedlhs.m -***** shared x, n, dim, box, niter, levels - n = [48; 12; 4; 2]; dim = 2; box = [0, 0; 4, 4]; niter = 10; - levels = [10.1; 15.2; -9.3; 2.4; 17.5]; -***** error x = stk_sampling_nestedlhs (); -***** test x = stk_sampling_nestedlhs (n); -***** test x = stk_sampling_nestedlhs (n, dim); -***** test x = stk_sampling_nestedlhs (n, dim, box); -***** test x = stk_sampling_nestedlhs (n, dim, box, niter); -***** test x = stk_sampling_nestedlhs (n, dim, box, niter, levels); -***** assert ( isequal(size(x), [sum(n), dim + 1]) ); -***** assert ( isa(x, 'stk_dataframe') ); - cn = [0; cumsum(n)]; - for lev = 1:length(n), - y = x( (cn(lev) + 1):(cn(lev + 1)), 1:dim ); - assert (isequal (size (y), [n(lev) dim])); - assert (stk_is_lhs (y, n(lev), dim, box)); - if lev > 1 - assert ( isequal(z((end - n(lev) + 1):end, :), y) ); - end - z = y; - end -***** assert (isequal (x.colnames{dim + 1}, 'Level')); - levels = stk_dataframe(levels, {'t'}); - box = stk_hrect(box, {'x1', 'x2', 'x3', 'x4'}); -***** test x = stk_sampling_nestedlhs (n, [], box, [], levels); -***** assert (isequal(x.colnames, {'x1', 'x2', 'x3', 'x4', 't'}) ); -11 tests, 11 passed, 0 known failure, 0 skipped -[inst/sampling/stk_sampling_vdc_rr2.m] ->>>>> /build/octave-stk-2.7.0/inst/sampling/stk_sampling_vdc_rr2.m -***** error __stk_sampling_vdc_rr2__() % two inputs required +[inst/covfcs/stk_materncov52_aniso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_materncov52_aniso.m +***** shared param, x, y + dim = 1; + param = log ([1.0; 2.5]); + x = stk_sampling_randunif (5, dim); + y = stk_sampling_randunif (5, dim); +***** error stk_materncov52_aniso (); +***** error stk_materncov52_aniso (param); +***** error stk_materncov52_aniso (param, x); +***** test stk_materncov52_aniso (param, x, y); +***** test stk_materncov52_aniso (param, x, y, -1); +***** test stk_materncov52_aniso (param, x, y, -1, false); +***** error stk_materncov52_aniso (param, x, y, -2); +***** test stk_materncov52_aniso (param, x, y, -1); +***** error stk_materncov52_aniso (param, x, y, 0); +***** test stk_materncov52_aniso (param, x, y, 1); +***** test stk_materncov52_aniso (param, x, y, 2); +***** error stk_materncov52_aniso (param, x, y, 3); +***** error stk_materncov52_aniso (param, x, y, nan); +***** error stk_materncov52_aniso (param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 2.5; 2.4; 2.6]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_materncov52_aniso (param, x, y); + K2 = stk_materncov52_aniso (param, x, y, -1); + assert (isequal (size (K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:(dim + 1), + dK = stk_materncov52_aniso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif (n, dim); + y = stk_sampling_randunif (n, dim); + + K1 = stk_materncov52_aniso (param, x, y); + K2 = stk_materncov52_aniso(param, x, y, -1, true); + assert (isequal (size (K1), [n n])); + assert (stk_isequal_tolabs (K2, diag (K1))); + + for i = 1:(dim + 1), + dK1 = stk_materncov52_aniso (param, x, y, i); + dK2 = stk_materncov52_aniso (param, x, y, i, true); + assert (isequal (size (dK1), [n n])); + assert (stk_isequal_tolabs (dK2, diag (dK1))); + end +17 tests, 17 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_sphcov_aniso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_sphcov_aniso.m +***** shared param, x, y, K1, K2, K3 + dim = 1; + param = log ([1.0; 2.5]); + x = stk_sampling_randunif (5, dim); + y = stk_sampling_randunif (6, dim); +***** error K0 = stk_sphcov_aniso (); +***** error K0 = stk_sphcov_aniso (param); +***** error K0 = stk_sphcov_aniso (param, x); +***** test K1 = stk_sphcov_aniso (param, x, y); +***** test K2 = stk_sphcov_aniso (param, x, y, -1); +***** test K3 = stk_sphcov_aniso (param, x, y, -1, false); +***** assert (isequal (K1, K2)); +***** assert (isequal (K1, K3)); +***** test % df versus ordinary array + u = double (x); v = double (y); + K1 = stk_sphcov_aniso (param, u, v, -1); + K2 = stk_sphcov_aniso (param, stk_dataframe (u), stk_dataframe (v), -1); + assert (isequal (K1, K2)); +***** error stk_sphcov_aniso (param, x, y, -2); +***** test stk_sphcov_aniso (param, x, y, -1); +***** error stk_sphcov_aniso (param, x, y, 0); +***** test stk_sphcov_aniso (param, x, y, 1); +***** test stk_sphcov_aniso (param, x, y, 2); +***** error stk_sphcov_aniso (param, x, y, 3); +***** error stk_sphcov_aniso (param, x, y, nan); +***** error stk_sphcov_aniso (param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 2.5; 2.4; 2.6]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_sphcov_aniso (param, x, y); + K2 = stk_sphcov_aniso (param, x, y, -1); + assert (isequal (size(K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:(dim + 1), + dK = stk_sphcov_aniso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif (n, dim); + y = stk_sampling_randunif (n, dim); + + K1 = stk_sphcov_aniso (param, x, y); + K2 = stk_sphcov_aniso (param, x, y, -1, true); + assert (isequal (size (K1), [n n])); + assert (stk_isequal_tolabs (K2, diag (K1))); + + for i = 1:(dim + 1), + dK1 = stk_sphcov_aniso (param, x, y, i); + dK2 = stk_sphcov_aniso (param, x, y, i, true); + assert (isequal (size (dK1), [n n])); + assert (stk_isequal_tolabs (dK2, diag (dK1))); + end +20 tests, 20 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_materncov_iso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_materncov_iso.m +***** shared param, x, y + dim = 1; + param = log ([1.0; 1.5; 2.9]); + x = stk_sampling_randunif(5, dim); + y = stk_sampling_randunif(5, dim); +***** error K = stk_materncov_iso ([param; 1.234], x, y); +***** error stk_materncov_iso(); +***** error stk_materncov_iso(param); +***** error stk_materncov_iso(param, x); +***** test stk_materncov_iso(param, x, y); +***** test stk_materncov_iso(param, x, y, -1); +***** test stk_materncov_iso(param, x, y, -1, false); +***** error stk_materncov_iso(param, x, y, -2); +***** test stk_materncov_iso(param, x, y, -1); +***** error stk_materncov_iso(param, x, y, 0); +***** test stk_materncov_iso(param, x, y, 1); +***** test stk_materncov_iso(param, x, y, 2); +***** test stk_materncov_iso(param, x, y, 3); +***** error stk_materncov_iso(param, x, y, 4); +***** error stk_materncov_iso(param, x, y, nan); +***** error stk_materncov_iso(param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 1.5; 2.9]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_materncov_iso (param, x, y); + K2 = stk_materncov_iso (param, x, y, -1); + assert (isequal (size (K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:3, + dK = stk_materncov_iso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif(n, dim); + y = stk_sampling_randunif(n, dim); + + K1 = stk_materncov_iso(param, x, y); + K2 = stk_materncov_iso(param, x, y, -1, true); + assert(isequal(size(K1), [n n])); + assert(stk_isequal_tolabs(K2, diag(K1))); + + for i = 1:3, + dK1 = stk_materncov_iso(param, x, y, i); + dK2 = stk_materncov_iso(param, x, y, i, true); + assert(isequal(size(dK1), [n n])); + assert(stk_isequal_tolabs(dK2, diag(dK1))); + end +19 tests, 19 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_materncov32_aniso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_materncov32_aniso.m +***** shared param, x, y, K1, K2, K3 + dim = 1; + param = log ([1.0; 2.5]); + x = stk_sampling_randunif (5, dim); + y = stk_sampling_randunif (6, dim); +***** error K0 = stk_materncov32_aniso (); +***** error K0 = stk_materncov32_aniso (param); +***** error K0 = stk_materncov32_aniso (param, x); +***** test K1 = stk_materncov32_aniso (param, x, y); +***** test K2 = stk_materncov32_aniso (param, x, y, -1); +***** test K3 = stk_materncov32_aniso (param, x, y, -1, false); +***** assert (isequal (K1, K2)); +***** assert (isequal (K1, K3)); +***** test % df versus ordinary array + u = double (x); v = double (y); + K1 = stk_materncov32_aniso (param, u, v, -1); + K2 = stk_materncov32_aniso (param, stk_dataframe (u), stk_dataframe (v), -1); + assert (isequal (K1, K2)); +***** error stk_materncov32_aniso (param, x, y, -2); +***** test stk_materncov32_aniso (param, x, y, -1); +***** error stk_materncov32_aniso (param, x, y, 0); +***** test stk_materncov32_aniso (param, x, y, 1); +***** test stk_materncov32_aniso (param, x, y, 2); +***** error stk_materncov32_aniso (param, x, y, 3); +***** error stk_materncov32_aniso (param, x, y, nan); +***** error stk_materncov32_aniso (param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 2.5; 2.4; 2.6]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_materncov32_aniso (param, x, y); + K2 = stk_materncov32_aniso (param, x, y, -1); + assert (isequal (size(K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:(dim + 1), + dK = stk_materncov32_aniso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif (n, dim); + y = stk_sampling_randunif (n, dim); + + K1 = stk_materncov32_aniso (param, x, y); + K2 = stk_materncov32_aniso (param, x, y, -1, true); + assert (isequal (size (K1), [n n])); + assert (stk_isequal_tolabs (K2, diag (K1))); + + for i = 1:(dim + 1), + dK1 = stk_materncov32_aniso (param, x, y, i); + dK2 = stk_materncov32_aniso (param, x, y, i, true); + assert (isequal (size (dK1), [n n])); + assert (stk_isequal_tolabs (dK2, diag (dK1))); + end +20 tests, 20 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_expcov_iso.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_expcov_iso.m +***** shared param, x, y + dim = 1; + param = log ([1.0; 2.5]); + x = stk_sampling_randunif (5, dim); + y = stk_sampling_randunif (5, dim); +***** error K = stk_expcov_iso ([param; 1.234], x, y); +***** error stk_expcov_iso (); +***** error stk_expcov_iso (param); +***** error stk_expcov_iso (param, x); +***** test stk_expcov_iso (param, x, y); +***** test stk_expcov_iso (param, x, y, -1); +***** test stk_expcov_iso (param, x, y, -1, false); +***** error stk_expcov_iso (param, x, y, -2); +***** test stk_expcov_iso (param, x, y, -1); +***** error stk_expcov_iso (param, x, y, 0); +***** test stk_expcov_iso (param, x, y, 1); +***** test stk_expcov_iso (param, x, y, 2); +***** error stk_expcov_iso (param, x, y, 3); +***** error stk_expcov_iso (param, x, y, nan); +***** error stk_expcov_iso (param, x, y, inf); +***** shared dim, param, x, y, nx, ny + dim = 3; + param = log ([1.0; 2.5]); + nx = 4; ny = 10; + x = stk_sampling_randunif (nx, dim); + y = stk_sampling_randunif (ny, dim); +***** test + K1 = stk_expcov_iso (param, x, y); + K2 = stk_expcov_iso (param, x, y, -1); + assert (isequal (size (K1), [nx ny])); + assert (stk_isequal_tolabs (K1, K2)); +***** test + for i = 1:2, + dK = stk_expcov_iso (param, x, y, i); + assert (isequal (size (dK), [nx ny])); + end +***** test + n = 7; + x = stk_sampling_randunif (n, dim); + y = stk_sampling_randunif (n, dim); + + K1 = stk_expcov_iso (param, x, y); + K2 = stk_expcov_iso (param, x, y, -1, true); + assert (isequal (size (K1), [n n])); + assert (stk_isequal_tolabs (K2, diag (K1))); + + for i = 1:2, + dK1 = stk_expcov_iso (param, x, y, i); + dK2 = stk_expcov_iso (param, x, y, i, true); + assert (isequal (size (dK1), [n n])); + assert (stk_isequal_tolabs (dK2, diag (dK1))); + end +18 tests, 18 passed, 0 known failure, 0 skipped +[inst/covfcs/stk_discretecov.m] +>>>>> /build/octave-stk-2.7.0/inst/covfcs/stk_discretecov.m +***** shared model, model2, x0 + n0 = 20; n1 = 10; dim = 4; + x0 = stk_sampling_randunif (n0, dim); + x1 = stk_sampling_randunif (n1, dim); + model = stk_model (@stk_materncov52_aniso, dim); + model.lm = stk_lm_affine (); + model.param = log ([1.0; 2.1; 2.2; 2.3; 2.4]); +***** test % without noise, pairwise = false + model.lognoisevariance = - inf; + model2 = stk_model (@stk_discretecov, model, x0); + idx = [1 4 9]; + [K1, P1] = stk_make_matcov (model, x0(idx, :)); + [K2, P2] = stk_make_matcov (model2, idx'); + assert (stk_isequal_tolrel (K1, K2)); + assert (stk_isequal_tolrel (P1, P2)); +***** test % without noise, pairwise = true + K1 = stk_make_matcov (model, x0([2 5 6], :), [], true); + K2 = stk_make_matcov (model2, [2 5 6]', [], true); + assert (stk_isequal_tolrel (K1, K2)); +***** test % with noise, pairwise = false + model.lognoisevariance = log (0.01); + model2 = stk_model (@stk_discretecov, model, x0); + idx = [1 4 9]; + [K1, P1] = stk_make_matcov (model, x0(idx, :)); + [K2, P2] = stk_make_matcov (model2, idx'); + assert (stk_isequal_tolrel (K1, K2)); + assert (stk_isequal_tolrel (P1, P2)); +3 tests, 3 passed, 0 known failure, 0 skipped +[inst/examples/01_kriging_basics/stk_example_kb09.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/01_kriging_basics/stk_example_kb09.m +***** test stk_example_kb09; close all; + +#========================# +# stk_example_kb09 # +#========================# + +'stk_example_kb09' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/01_kriging_basics/stk_example_kb09.m + + STK_EXAMPLE_KB09 Generation of sample paths conditioned on noisy observations + + A Matern Gaussian process model is used, with constant but unknown mean + (ordinary kriging) and known covariance parameters. + + Given noisy observations from the unknown function, a batch of conditioned + sample paths is drawn using the "conditioning by kriging" technique + (stk_generate_samplepaths function). + + See also: stk_generate_samplepaths, stk_conditioning, stk_example_kb05 + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. +warning: using the gnuplot graphics toolkit is discouraged + +The gnuplot graphics toolkit is not actively maintained and has a number +of limitations that are unlikely to be fixed. Communication with gnuplot +uses a one-directional pipe and limited information is passed back to the +Octave interpreter so most changes made interactively in the plot window +will not be reflected in the graphics properties managed by Octave. For +example, if the plot window is closed with a mouse click, Octave will not +be notified and will not update its internal list of open figure windows. +The qt toolkit is recommended instead. +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/01_kriging_basics/stk_example_kb02n.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/01_kriging_basics/stk_example_kb02n.m +***** test stk_example_kb02n; close all; + +#=========================# +# stk_example_kb02n # +#=========================# + +'stk_example_kb02n' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/01_kriging_basics/stk_example_kb02n.m + + STK_EXAMPLE_KB02N Noisy ordinary kriging in 1D with parameter estimation + + This example shows how to estimate covariance parameters and compute + ordinary kriging predictions on a one-dimensional noisy dataset. + + The model and data are the same as in stk_example_kb02, but this time the + parameters of the covariance function and the variance of the noise are + jointly estimated using the Restricted Maximum Likelihood (ReML) method. + + See also: stk_example_kb01n, stk_example_kb02 + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. +model = + + scalar structure containing the fields: + + covariance_type = @stk_materncov_iso + lm = + + + + dim = 1 + param = + + -0.098610 + 2.302585 + 0.613435 + + lognoisevariance = -3.0915 + +True noise variance = 0.0400 +Estimated noise variance = 0.0454 + + line 0: warning: iconv failed to convert degree sign + + + | True function and noisy observed data + 1.5 +|+ + | G ***** +--------------------+ + | *** G G** |***?***True function| + 1 +|+ G** **G +---G---Observations-+ + | ** *** + | ** ** + 0.5 +|+ ** G* G *** + G| ** ** **** G**G* G + |*G ** *** *** + F1 0 *|+ ** ** *** + *| ** **G **G + | **G ** ** + -0.5 +|+ G ** ** + | ** *G + | ***G ** + -1 +|+ ** *G + | ***G** + +|------------------------------------------------------------------ + -1.5 +-+ + + + + + -1 -0.5 0 0.5 1 + +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/01_kriging_basics/stk_example_kb05.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/01_kriging_basics/stk_example_kb05.m +***** test stk_example_kb05; close all; + +#========================# +# stk_example_kb05 # +#========================# + +'stk_example_kb05' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/01_kriging_basics/stk_example_kb05.m + + STK_EXAMPLE_KB05 Generation of conditioned sample paths + + A Matern Gaussian process model is used, with constant but unknown mean + (ordinary kriging) and known covariance parameters. + + Given noiseless observations from the unknown function, a batch of conditioned + sample paths is drawn using the "conditioning by kriging" technique. In short, + this means that unconditioned sample path are simulated first (using + stk_generate_samplepaths), and then conditioned on the observations by kriging + (using stk_conditioning). + + Note: in this example, for pedagogical purposes, conditioned samplepaths are + simulated in two steps: first, unconditioned samplepaths are simulated; + second, conditioned samplepaths are obtained using conditioning by kriging. + In practice, these two steps can be carried out all at once using + stk_generate_samplepath (see, e.g., stk_example_kb09). + + See also: stk_generate_samplepaths, stk_conditioning, stk_example_kb09 + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/01_kriging_basics/stk_example_kb07.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/01_kriging_basics/stk_example_kb07.m +***** test stk_example_kb07; close all; + +#========================# +# stk_example_kb07 # +#========================# + +'stk_example_kb07' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/01_kriging_basics/stk_example_kb07.m + + STK_EXAMPLE_KB07 Simulation of sample paths from a Matern process + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. + line 0: warning: iconv failed to convert degree sign + + |--Matern, nu = 0.5------| |--Matern, nu = 1.5------| + 4 +|+ + + + +-+| 3 +|+ + + + +-+| + 3 +|+ ### + @@@@ @ +&+| 2 +|+ + + + ***-+| + 2 $|######*#** @ @@@@@&&&| | ****%%%==##$$$ ***$$*&| + 1 $|$$$&$$$$$$$$$@@$$$**##==| 1 $|***&%*%##$$$$##$$*=&$$$$| +response z$|$$$$$&#$$$$$$##$$####$$response z*|$%$$$$$$$$******$$=$$==$| + 0 *|**==***####&$$$$$$##$$$$| -1 #|$$$####*=@**@@&&@$$$%@@+| + -1 =|=**=***&&&& $$##$$$$$**| #|####*=== *@***@@%#### | + -2 +|------------------------| -2 +|------------------------| + -3 +-+ *+ + + +-+ -3 +-+ == **** + + +-+ + -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 + |--Matern, nu = 2.5------| |--Matern, nu = 10.0-----| + 3 +|+ + + + +-+| 3 +|+ + + + +-+| + 2 +|+ + +*** + +-@| +| + + + +| + |** **@** $$| 2 +|+ #### $$$$$$$$####| + 1 =|*=&###&##$$$$$$*****$$$*| 1 =|$$$$$$$$$$$$#######$$$**| +response z#|####*##$$$####$$$$$$$##response z*|$******##$======&&&##$$*| + -1 #|##%*=*$$*=&&&###*&==@%$&| 0 $|&&&& $##**&=********##$&| + $|# %% $$%&&====%&&&%=====| -1 #|#@@&###@&**** ###@==@=#$| + -2 $|------------------------| @|------------------------| + -3 +-+ + + + +-+ -2 +-+ + + + +-+ + -1 -0.5 0 0.5 1 -1 -0.5 0 0.5 1 + input variable x input variable x +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/01_kriging_basics/stk_example_kb08.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/01_kriging_basics/stk_example_kb08.m +***** test stk_example_kb08; close all; + +#========================# +# stk_example_kb08 # +#========================# + +'stk_example_kb08' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/01_kriging_basics/stk_example_kb08.m + + STK_EXAMPLE_KB08 Generation of conditioned sample paths made easy + + It has been demonstrated, in stk_example_kb05, how to generate conditioned + sample paths using unconditioned sample paths and conditioning by kriging. + + This example shows how to do the same in a more concise way, letting STK + take care of the details. + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/01_kriging_basics/stk_example_kb06.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/01_kriging_basics/stk_example_kb06.m +***** test stk_example_kb06; close all; + +#========================# +# stk_example_kb06 # +#========================# + +'stk_example_kb06' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/01_kriging_basics/stk_example_kb06.m + + STK_EXAMPLE_KB06 Ordinary kriging VS kriging with a linear trend + + The same dataset is analyzed using two variants of kriging. + + The left panel shows the result of ordinary kriging, in other words, Gaussian + process interpolation assuming a constant (but unknown) mean. The right panel + shows the result of adding a linear trend in the mean of the Gaussian process. + + The difference with the left plot is clear in extrapolation: the first predic- + tor exhibits a "mean reverting" behaviour, while the second one captures an + increasing trend in the data. + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/01_kriging_basics/stk_example_kb01.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/01_kriging_basics/stk_example_kb01.m +***** test stk_example_kb01; close all; + +#========================# +# stk_example_kb01 # +#========================# + +'stk_example_kb01' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/01_kriging_basics/stk_example_kb01.m + + STK_EXAMPLE_KB01 Ordinary kriging in 1D, with noiseless data + + This example shows how to compute ordinary kriging predictions on a + one-dimensional noiseless dataset. + + The word 'ordinary' indicates that the mean function of the GP prior is + assumed to be constant and unknown. + + A Matern covariance function is used for the Gaussian Process (GP) prior. + The parameters of this covariance function are assumed to be known (i.e., + no parameter estimation is performed here). + + Note that the kriging predictor, which is the posterior mean of the GP, + interpolates the data in this noiseless example. + + See also: stk_example_kb01n, stk_example_kb02 + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. +model = + + scalar structure containing the fields: + + covariance_type = @stk_materncov_iso + lm = + + + + dim = 1 + param = + + -0.6931 + 1.3863 + 0.9163 + + lognoisevariance = -Inf + + line 0: warning: iconv failed to convert degree sign + + + | True function and observed data + 1.5 +|+ + | ****G** +--------------------+ + | *** **** |***?***True function| + 1 +|+ *** *** +---G---Observations-+ + | ** ** + | ** ** + 0.5 +|+ ** ** *****G + |*G ** *** *** + 0 +|* ** G* *** + G| ** *** *** + | ** ** ** + -0.5 +|+ ** ** * + | ** ** + | *G ** + -1 +|+ *** ** + | ******* + +|------------------------------------------------------------------ + -1.5 +-+ + + + + + -1 -0.5 0 0.5 1 + input variable x + +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/01_kriging_basics/stk_example_kb03.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/01_kriging_basics/stk_example_kb03.m +***** test stk_example_kb03; close all; + +#========================# +# stk_example_kb03 # +#========================# + +'stk_example_kb03' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/01_kriging_basics/stk_example_kb03.m + + STK_EXAMPLE_KB03 Ordinary kriging in 2D + + An anisotropic Matern covariance function is used for the Gaussian Process + (GP) prior. The parameters of this covariance function (variance, regularity + and ranges) are estimated using the Restricted Maximum Likelihood (ReML) + method. + + The mean function of the GP prior is assumed to be constant and unknown. This + default choice can be overridden by means of the model.lm property. + + The function is sampled on a space-filling Latin Hypercube design, and the + data is assumed to be noiseless. + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. + +BOX = 2-dimensional hyper-rectangle (stk_hrect object): + + : x_1 x_2 + lower_bounds : -5 0 + upper_bounds : 10 15 + +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/01_kriging_basics/stk_example_kb02.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/01_kriging_basics/stk_example_kb02.m +***** test stk_example_kb02; close all; + +#========================# +# stk_example_kb02 # +#========================# + +'stk_example_kb02' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/01_kriging_basics/stk_example_kb02.m + + STK_EXAMPLE_KB02 Ordinary kriging in 1D with parameter estimation + + This example shows how to estimate covariance parameters and compute + ordinary kriging predictions on a one-dimensional noiseless dataset. + + The model and data are the same as in stk_example_kb01, but this time the + parameters of the covariance function are estimated using the Restricted + Maximum Likelihood (ReML) method. + + See also: stk_example_kb01, stk_example_kb02n + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. +model = + + scalar structure containing the fields: + + covariance_type = @stk_materncov_iso + lm = + + + + dim = 1 + param = + + -0.2108 + 2.3026 + 0.7562 + + lognoisevariance = -Inf + + line 0: warning: iconv failed to convert degree sign + + + | True function and observed data + 1.5 +|+ + | ***G* +--------------------+ + | *** *** |***?***True function| + 1 +|+ ** *** +---G---Observations-+ + | ** *** + | ** ** + 0.5 +|+ ** ** *** + | G* ** **** *G** + |** ** G** *** + F1 0 *|+ ** ** *** + G| ** ** ** + | ** ** ** + -0.5 +|+ ** ** + | ** * + | *G* ** + -1 +|+ ** ** + | ****** + +|------------------------------------------------------------------ + -1.5 +-+ + + + + + -1 -0.5 0 0.5 1 + +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/01_kriging_basics/stk_example_kb01n.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/01_kriging_basics/stk_example_kb01n.m +***** test stk_example_kb01n; close all; + +#=========================# +# stk_example_kb01n # +#=========================# + +'stk_example_kb01n' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/01_kriging_basics/stk_example_kb01n.m + + STK_EXAMPLE_KB01N Ordinary kriging in 1D, with noisy data + + This example shows how to compute ordinary kriging predictions on a + one-dimensional noisy dataset. + + The Gaussian Process (GP) prior is the same as in stk_example_kb01. + + The observation noise is Gaussian and homoscedastic (constant variance). + Its variance is assumed to be known. + + Note that the kriging predictor, which is the posterior mean of the GP, + does NOT interpolate the data in this noisy example. + + See also: stk_example_kb01, stk_example_kb02n + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. +model = + + scalar structure containing the fields: + + covariance_type = @stk_materncov_iso + lm = + + + + dim = 1 + param = + + -0.6931 + 1.3863 + 0.9163 + + lognoisevariance = -3.2189 + + line 0: warning: iconv failed to convert degree sign + + + | True function and observed data + 1.5 +|+ + | G******* +--------------------+ + | *** G G*** |***?***True function| + 1 +|+ G** *G* +---G---Observations-+ + | ** ** + | ** *G G + 0.5 G|+ ** ** ******G G + |*G ** *** G *** + 0 +|* ** ** *** + *| ** **G *** + | ** ** *G + -0.5 +|+ G *G ** * + | ** *G + | ** G ** + -1 +|+ *** *G + | ***G*** + +|------------------------------------------------------------------ + -1.5 +-+ + + + + + -1 -0.5 0 0.5 1 + input variable x + +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/02_design_of_experiments/stk_example_doe02.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/02_design_of_experiments/stk_example_doe02.m +***** test stk_example_doe02; close all; + +#=========================# +# stk_example_doe02 # +#=========================# + +'stk_example_doe02' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/02_design_of_experiments/stk_example_doe02.m + + STK_EXAMPLE_DOE02 "Sequential Maximin" design + + In this example, a two-dimensional space-filling design is progressively + enriched with new points using a "sequential maximin" approach. More + precisely, the k-th point X(k, :) is selected to maximize the distance to the + set of all previously selected points X(1, :), X(2, :), ..., X(k-1, :). + + NOTES: + + * The resulting design is NOT optimal with respect to the maximin criterion + (separation distance). + + * This procedure is not truly a *sequential* design procedure, since the + choice of the k-th point X(k, :) does NOT depend on the response at the + previously selected locations X(i, :), i < k. + + REFERENCE + + [1] Emmanuel Vazquez and Julien Bect, "Sequential search based on kriging: + convergence analysis of some algorithms", In: ISI - 58th World + Statistics Congress of the International Statistical Institute (ISI'11), + Dublin, Ireland, August 21-26, 2011. + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. + line 0: warning: iconv failed to convert degree sign + +gnuplot> " font ":Bold,10" textcolor rgb " + ^ + line 0: invalid command + + + + |----------------------n = 10, fd = 3.26e-01-------------------------| + 1 +|+ + + + + +-+| + +| + +F + + +| + | | + | F | + 0.8 +|+F +-+| + | | + | F | + 0.6 +|+ +-+| + | F | + | | + | | + 0.4 +|+ F F +-+| + | | + | | + 0.2 +|+ F +-+| + | F | + | | + +|---------------------------------------------------------------------| + 0 F-+ + + + + +-+ + 0 0.2 0.4 0.6 0.8 1 + + line 0: warning: iconv failed to convert degree sign + +gnuplot> " font ":Bold,10" textcolor rgb " + ^ + line 0: invalid command + + + + |----------------------n = 11, fd = 3.07e-01-------------------------| + 1 +|+ + + + + +-+| + +| + +F + + +| + | | + | F | + 0.8 +|+F +-+| + | | + | F | + 0.6 +|+ +-F| + | F | + | | + | | + 0.4 +|+ F F +-+| + | | + | | + 0.2 +|+ F +-+| + | F | + | | + +|---------------------------------------------------------------------| + 0 F-+ + + + + +-+ + 0 0.2 0.4 0.6 0.8 1 + + line 0: warning: iconv failed to convert degree sign + +gnuplot> " font ":Bold,10" textcolor rgb " + ^ + line 0: invalid command + + + + |----------------------n = 12, fd = 2.93e-01-------------------------| + 1 +|+ + + + + +-F| + +| + +F + + +| + | | + | F | + 0.8 +|+F +-+| + | | + | F | + 0.6 +|+ +-F| + | F | + | | + | | + 0.4 +|+ F F +-+| + | | + | | + 0.2 +|+ F +-+| + | F | + | | + +|---------------------------------------------------------------------| + 0 F-+ + + + + +-+ + 0 0.2 0.4 0.6 0.8 1 + + line 0: warning: iconv failed to convert degree sign + +gnuplot> " font ":Bold,10" textcolor rgb " + ^ + line 0: invalid command + + + + |----------------------n = 13, fd = 2.76e-01-------------------------| + 1 +|+ + + + + +-F| + +| + +F + + +| + | | + | F | + 0.8 +|+F +-+| + | | + | F | + 0.6 +|+ +-F| + | F | + | | + | | + 0.4 +|+ F F +-+| + | | + | | + 0.2 +|+ F +-+| + | F | + | | + +|---------------------------------------------------------------------| + 0 F-+ + F + + + +-+ + 0 0.2 0.4 0.6 0.8 1 + + line 0: warning: iconv failed to convert degree sign + +gnuplot> " font ":Bold,10" textcolor rgb " + ^ + line 0: invalid command + + + + |----------------------n = 14, fd = 2.53e-01-------------------------| + 1 +|+ + + + + +-F| + +| + +F + + +| + | | + | F | + 0.8 +|+F +-+| + | | + | F | + 0.6 +|+ +-F| + | F | + | | + | | + 0.4 +|+ F F +-+| + F| | + | | + 0.2 +|+ F +-+| + | F | + | | + +|---------------------------------------------------------------------| + 0 F-+ + F + + + +-+ + 0 0.2 0.4 0.6 0.8 1 + + line 0: warning: iconv failed to convert degree sign + +gnuplot> " font ":Bold,10" textcolor rgb " + ^ + line 0: invalid command + + + + |----------------------n = 15, fd = 2.50e-01-------------------------| + 1 +|+ F + + + + +-F| + +| + +F + + +| + | | + | F | + 0.8 +|+F +-+| + | | + | F | + 0.6 +|+ +-F| + | F | + | | + | | + 0.4 +|+ F F +-+| + F| | + | | + 0.2 +|+ F +-+| + | F | + | | + +|---------------------------------------------------------------------| + 0 F-+ + F + + + +-+ + 0 0.2 0.4 0.6 0.8 1 + + line 0: warning: iconv failed to convert degree sign + +gnuplot> " font ":Bold,10" textcolor rgb " + ^ + line 0: invalid command + + + + |----------------------n = 16, fd = 2.37e-01-------------------------| + 1 +|+ F + + + + +-F| + +| + +F + + +| + | | + | F | + 0.8 +|+F +-+| + | | + | F | + 0.6 +|+ F +-F| + | F | + | | + | | + 0.4 +|+ F F +-+| + F| | + | | + 0.2 +|+ F +-+| + | F | + | | + +|---------------------------------------------------------------------| + 0 F-+ + F + + + +-+ + 0 0.2 0.4 0.6 0.8 1 + + line 0: warning: iconv failed to convert degree sign + +gnuplot> " font ":Bold,10" textcolor rgb " + ^ + line 0: invalid command + + + + |----------------------n = 17, fd = 2.28e-01-------------------------| + 1 +|+ F + + + F + +-F| + +| + +F + + +| + | | + | F | + 0.8 +|+F +-+| + | | + | F | + 0.6 +|+ F +-F| + | F | + | | + | | + 0.4 +|+ F F +-+| + F| | + | | + 0.2 +|+ F +-+| + | F | + | | + +|---------------------------------------------------------------------| + 0 F-+ + F + + + +-+ + 0 0.2 0.4 0.6 0.8 1 + + line 0: warning: iconv failed to convert degree sign + +gnuplot> " font ":Bold,10" textcolor rgb " + ^ + line 0: invalid command + + + + |----------------------n = 18, fd = 2.21e-01-------------------------| + 1 +|+ F + + + F + +-F| + +| + +F + + +| + | | + | F | + 0.8 +|+F +-+| + | | + | F | + 0.6 +|+ F +-F| + | F | + | | + | | + 0.4 +|+ F F +-+| + F| | + | | + 0.2 +|+ F +-+| + | F | + | | + +|---------------------------------------------------------------------| + 0 F-+ + F + + F + +-+ + 0 0.2 0.4 0.6 0.8 1 + + line 0: warning: iconv failed to convert degree sign + +gnuplot> " font ":Bold,10" textcolor rgb " + ^ + line 0: invalid command + + + + |----------------------n = 19, fd = 2.19e-01-------------------------| + 1 +|+ F + + + F + +-F| + +| + +F + + +| + | | + | F | + 0.8 +|+F +-+| + | | + | F | + 0.6 +|+ F +-F| + | F | + | | + | F | + 0.4 +|+ F F +-+| + F| | + | | + 0.2 +|+ F +-+| + | F | + | | + +|---------------------------------------------------------------------| + 0 F-+ + F + + F + +-+ + 0 0.2 0.4 0.6 0.8 1 + +1 test, 1 passed, 0 known failure, 0 skipped +[inst/examples/02_design_of_experiments/stk_example_doe05.m] +>>>>> /build/octave-stk-2.7.0/inst/examples/02_design_of_experiments/stk_example_doe05.m +***** test stk_example_doe05; close all; + +#=========================# +# stk_example_doe05 # +#=========================# + +'stk_example_doe05' is a script from the file /build/octave-stk-2.7.0/debian/octave-stk/usr/share/octave/packages/stk-2.7.0/examples/02_design_of_experiments/stk_example_doe05.m + + STK_EXAMPLE_DOE05 A simple illustration of 1D Bayesian optimization + + Our goal here is to minimize the one-dimensional function + + x |--> x * sin (x) + + over the interval [0; 4 * pi], using noisy evaluations. + + Evaluations points are chosen sequentially using either AKG criterion + (default) or the EQI criterion (set SAMPCRIT_NAME to 'EQI'); + + +Additional help for built-in functions and operators is +available in the online version of the manual. Use the command +'doc ' to search the manual index. + +Help and information about Octave is also available on the WWW +at https://www.octave.org and via the help@octave.org +mailing list. terminate called after throwing an instance of 'octave::execution_exception' - what(): __stk_sampling_vdc_rr2__: Two inputs required. + what(): __get_column_number__: There is no variable named type. fatal: caught signal Aborted -- stopping myself... -Aborted +/usr/bin/dh_octave_check: line 243: 29626 Aborted octave-cli --no-history --silent --no-init-file --no-window-system /tmp/tmp.9XTBtzs9Ty Summary: 0 tests, 0 passed, 0 known failures, 0 skipped dh_installdocs -O--buildsystem=octave dh_installchangelogs -O--buildsystem=octave @@ -3296,8 +5004,8 @@ dpkg-gencontrol: warning: package octave-stk: substitution variable ${octave:Upstream-Description} unused, but is defined dh_md5sums -O--buildsystem=octave dh_builddeb -O--buildsystem=octave -dpkg-deb: building package 'octave-stk-dbgsym' in '../octave-stk-dbgsym_2.7.0-3_armhf.deb'. dpkg-deb: building package 'octave-stk' in '../octave-stk_2.7.0-3_armhf.deb'. +dpkg-deb: building package 'octave-stk-dbgsym' in '../octave-stk-dbgsym_2.7.0-3_armhf.deb'. dpkg-genbuildinfo --build=binary -O../octave-stk_2.7.0-3_armhf.buildinfo dpkg-genchanges --build=binary -O../octave-stk_2.7.0-3_armhf.changes dpkg-genchanges: info: binary-only upload (no source code included) @@ -3305,12 +5013,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/9579/tmp/hooks/B01_cleanup starting +I: user script /srv/workspace/pbuilder/9579/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/7156 and its subdirectories -I: Current time: Thu Jun 8 05:23:16 -12 2023 -I: pbuilder-time-stamp: 1686244996 +I: removing directory /srv/workspace/pbuilder/9579 and its subdirectories +I: Current time: Fri Jun 9 07:34:37 +14 2023 +I: pbuilder-time-stamp: 1686245677