Diff of the two buildlogs: -- --- b1/build.log 2021-07-09 11:16:58.409860977 +0000 +++ b2/build.log 2021-07-09 11:19:44.595661546 +0000 @@ -1,6 +1,6 @@ I: pbuilder: network access will be disabled during build -I: Current time: Thu Jul 8 23:05:13 -12 2021 -I: pbuilder-time-stamp: 1625828713 +I: Current time: Fri Aug 12 07:40:04 +14 2022 +I: pbuilder-time-stamp: 1660239604 I: Building the build Environment I: extracting base tarball [/var/cache/pbuilder/bullseye-reproducible-base.tgz] I: copying local configuration @@ -17,8 +17,8 @@ I: copying [./octave-optim_1.6.0-6.debian.tar.xz] I: Extracting source gpgv: unknown type of key resource 'trustedkeys.kbx' -gpgv: keyblock resource '/tmp/dpkg-verify-sig.Qrq77Qve/trustedkeys.kbx': General error -gpgv: Signature made Wed Dec 9 21:30:51 2020 -12 +gpgv: keyblock resource '/tmp/dpkg-verify-sig.5H4TNco0/trustedkeys.kbx': General error +gpgv: Signature made Thu Dec 10 23:30:51 2020 +14 gpgv: using RSA key 53951D95272E0C5B82BE8C4A2CECE9350ECEBE4A gpgv: Can't check signature: No public key dpkg-source: warning: failed to verify signature on ./octave-optim_1.6.0-6.dsc @@ -31,137 +31,171 @@ dpkg-source: info: applying octave6.patch I: Not using root during the build. I: Installing the build-deps -I: user script /srv/workspace/pbuilder/30876/tmp/hooks/D02_print_environment starting +I: user script /srv/workspace/pbuilder/55977/tmp/hooks/D01_modify_environment starting +debug: Running on ionos6-i386. +I: Changing host+domainname to test build reproducibility +I: Adding a custom variable just for the fun of it... +I: Changing /bin/sh to bash +Removing 'diversion of /bin/sh to /bin/sh.distrib by dash' +Adding 'diversion of /bin/sh to /bin/sh.distrib by bash' +Removing 'diversion of /usr/share/man/man1/sh.1.gz to /usr/share/man/man1/sh.distrib.1.gz by dash' +Adding 'diversion of /usr/share/man/man1/sh.1.gz to /usr/share/man/man1/sh.distrib.1.gz by bash' +I: Setting pbuilder2's login shell to /bin/bash +I: Setting pbuilder2's GECOS to second user,second room,second work-phone,second home-phone,second other +I: user script /srv/workspace/pbuilder/55977/tmp/hooks/D01_modify_environment finished +I: user script /srv/workspace/pbuilder/55977/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='i386' - DEBIAN_FRONTEND='noninteractive' - DEB_BUILD_OPTIONS='buildinfo=+all reproducible=+all,-fixfilepath parallel=10' - DISTRIBUTION='' - HOME='/root' - HOST_ARCH='i386' + BASH=/bin/sh + BASHOPTS=checkwinsize:cmdhist:complete_fullquote:extquote:force_fignore:globasciiranges:hostcomplete:interactive_comments:progcomp:promptvars:sourcepath + BASH_ALIASES=() + BASH_ARGC=() + BASH_ARGV=() + BASH_CMDS=() + BASH_LINENO=([0]="12" [1]="0") + BASH_SOURCE=([0]="/tmp/hooks/D02_print_environment" [1]="/tmp/hooks/D02_print_environment") + BASH_VERSINFO=([0]="5" [1]="1" [2]="4" [3]="1" [4]="release" [5]="i686-pc-linux-gnu") + BASH_VERSION='5.1.4(1)-release' + BUILDDIR=/build + BUILDUSERGECOS='second user,second room,second work-phone,second home-phone,second other' + BUILDUSERNAME=pbuilder2 + BUILD_ARCH=i386 + DEBIAN_FRONTEND=noninteractive + DEB_BUILD_OPTIONS='buildinfo=+all reproducible=+all,-fixfilepath parallel=18' + DIRSTACK=() + DISTRIBUTION= + EUID=0 + FUNCNAME=([0]="Echo" [1]="main") + GROUPS=() + HOME=/root + HOSTNAME=i-capture-the-hostname + HOSTTYPE=i686 + HOST_ARCH=i386 IFS=' ' - INVOCATION_ID='95f27d8465e04f52bdd8fc1b1c59f07c' - LANG='C' - LANGUAGE='en_US:en' - LC_ALL='C' - LD_LIBRARY_PATH='/usr/lib/libeatmydata' - LD_PRELOAD='libeatmydata.so' - MAIL='/var/mail/root' - OPTIND='1' - PATH='/usr/sbin:/usr/bin:/sbin:/bin:/usr/games' - PBCURRENTCOMMANDLINEOPERATION='build' - PBUILDER_OPERATION='build' - PBUILDER_PKGDATADIR='/usr/share/pbuilder' - PBUILDER_PKGLIBDIR='/usr/lib/pbuilder' - PBUILDER_SYSCONFDIR='/etc' - PPID='30876' - PS1='# ' - PS2='> ' + INVOCATION_ID=00d7dfd9f01f490089ca8166d4b31b2a + LANG=C + LANGUAGE=de_CH:de + LC_ALL=C + LD_LIBRARY_PATH=/usr/lib/libeatmydata + LD_PRELOAD=libeatmydata.so + MACHTYPE=i686-pc-linux-gnu + MAIL=/var/mail/root + OPTERR=1 + OPTIND=1 + OSTYPE=linux-gnu + PATH=/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path + PBCURRENTCOMMANDLINEOPERATION=build + PBUILDER_OPERATION=build + PBUILDER_PKGDATADIR=/usr/share/pbuilder + PBUILDER_PKGLIBDIR=/usr/lib/pbuilder + PBUILDER_SYSCONFDIR=/etc + PIPESTATUS=([0]="0") + POSIXLY_CORRECT=y + PPID=55977 PS4='+ ' - PWD='/' - SHELL='/bin/bash' - SHLVL='2' - SUDO_COMMAND='/usr/bin/timeout -k 18.1h 18h /usr/bin/ionice -c 3 /usr/bin/nice /usr/sbin/pbuilder --build --configfile /srv/reproducible-results/rbuild-debian/tmp.aLtbgrlQtJ/pbuilderrc_srLF --hookdir /etc/pbuilder/first-build-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/bullseye-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/tmp.aLtbgrlQtJ/b1 --logfile b1/build.log octave-optim_1.6.0-6.dsc' - SUDO_GID='112' - SUDO_UID='107' - SUDO_USER='jenkins' - TERM='unknown' - TZ='/usr/share/zoneinfo/Etc/GMT+12' - USER='root' - _='/usr/bin/systemd-run' - http_proxy='http://78.137.99.97:3128' + PWD=/ + SHELL=/bin/bash + SHELLOPTS=braceexpand:errexit:hashall:interactive-comments:posix + SHLVL=3 + SUDO_COMMAND='/usr/bin/timeout -k 24.1h 24h /usr/bin/ionice -c 3 /usr/bin/nice -n 11 /usr/bin/unshare --uts -- /usr/sbin/pbuilder --build --configfile /srv/reproducible-results/rbuild-debian/tmp.aLtbgrlQtJ/pbuilderrc_N2r5 --hookdir /etc/pbuilder/rebuild-hooks --debbuildopts -b --basetgz /var/cache/pbuilder/bullseye-reproducible-base.tgz --buildresult /srv/reproducible-results/rbuild-debian/tmp.aLtbgrlQtJ/b2 --logfile b2/build.log --extrapackages usrmerge octave-optim_1.6.0-6.dsc' + SUDO_GID=112 + SUDO_UID=107 + SUDO_USER=jenkins + TERM=unknown + TZ=/usr/share/zoneinfo/Etc/GMT-14 + UID=0 + USER=root + _='I: set' + http_proxy=http://85.184.249.68:3128 I: uname -a - Linux ionos12-i386 4.19.0-17-686-pae #1 SMP Debian 4.19.194-2 (2021-06-21) i686 GNU/Linux + Linux i-capture-the-hostname 4.19.0-17-amd64 #1 SMP Debian 4.19.194-2 (2021-06-21) x86_64 GNU/Linux I: ls -l /bin - total 5804 - -rwxr-xr-x 1 root root 1367848 Jun 21 14:25 bash - -rwxr-xr-x 3 root root 38280 Jul 20 2020 bunzip2 - -rwxr-xr-x 3 root root 38280 Jul 20 2020 bzcat - lrwxrwxrwx 1 root root 6 Jul 20 2020 bzcmp -> bzdiff - -rwxr-xr-x 1 root root 2225 Jul 20 2020 bzdiff - lrwxrwxrwx 1 root root 6 Jul 20 2020 bzegrep -> bzgrep - -rwxr-xr-x 1 root root 4877 Sep 4 2019 bzexe - lrwxrwxrwx 1 root root 6 Jul 20 2020 bzfgrep -> bzgrep - -rwxr-xr-x 1 root root 3775 Jul 20 2020 bzgrep - -rwxr-xr-x 3 root root 38280 Jul 20 2020 bzip2 - -rwxr-xr-x 1 root root 17768 Jul 20 2020 bzip2recover - lrwxrwxrwx 1 root root 6 Jul 20 2020 bzless -> bzmore - -rwxr-xr-x 1 root root 1297 Jul 20 2020 bzmore - -rwxr-xr-x 1 root root 38824 Sep 22 2020 cat - -rwxr-xr-x 1 root root 71624 Sep 22 2020 chgrp - -rwxr-xr-x 1 root root 67528 Sep 22 2020 chmod - -rwxr-xr-x 1 root root 75752 Sep 22 2020 chown - -rwxr-xr-x 1 root root 157960 Sep 22 2020 cp - -rwxr-xr-x 1 root root 128724 Dec 10 2020 dash - -rwxr-xr-x 1 root root 124904 Sep 22 2020 date - -rwxr-xr-x 1 root root 92172 Sep 22 2020 dd - -rwxr-xr-x 1 root root 100752 Sep 22 2020 df - -rwxr-xr-x 1 root root 153964 Sep 22 2020 dir - -rwxr-xr-x 1 root root 83644 Feb 7 02:38 dmesg - lrwxrwxrwx 1 root root 8 Nov 6 2019 dnsdomainname -> hostname - lrwxrwxrwx 1 root root 8 Nov 6 2019 domainname -> hostname - -rwxr-xr-x 1 root root 34664 Sep 22 2020 echo - -rwxr-xr-x 1 root root 28 Nov 9 2020 egrep - -rwxr-xr-x 1 root root 34664 Sep 22 2020 false - -rwxr-xr-x 1 root root 28 Nov 9 2020 fgrep - -rwxr-xr-x 1 root root 71928 Feb 7 02:38 findmnt - -rwsr-xr-x 1 root root 30112 Feb 26 04:12 fusermount - -rwxr-xr-x 1 root root 210488 Nov 9 2020 grep - -rwxr-xr-x 2 root root 2346 Mar 2 11:30 gunzip - -rwxr-xr-x 1 root root 6376 Mar 2 11:30 gzexe - -rwxr-xr-x 1 root root 100952 Mar 2 11:30 gzip - -rwxr-xr-x 1 root root 21916 Nov 6 2019 hostname - -rwxr-xr-x 1 root root 83980 Sep 22 2020 ln - -rwxr-xr-x 1 root root 55572 Feb 7 2020 login - -rwxr-xr-x 1 root root 153964 Sep 22 2020 ls - -rwxr-xr-x 1 root root 153124 Feb 7 02:38 lsblk - -rwxr-xr-x 1 root root 96328 Sep 22 2020 mkdir - -rwxr-xr-x 1 root root 79912 Sep 22 2020 mknod - -rwxr-xr-x 1 root root 47048 Sep 22 2020 mktemp - -rwxr-xr-x 1 root root 58920 Feb 7 02:38 more - -rwsr-xr-x 1 root root 50720 Feb 7 02:38 mount - -rwxr-xr-x 1 root root 13856 Feb 7 02:38 mountpoint - -rwxr-xr-x 1 root root 157996 Sep 22 2020 mv - lrwxrwxrwx 1 root root 8 Nov 6 2019 nisdomainname -> hostname - lrwxrwxrwx 1 root root 14 Apr 18 03:38 pidof -> /sbin/killall5 - -rwxr-xr-x 1 root root 38824 Sep 22 2020 pwd - lrwxrwxrwx 1 root root 4 Jun 21 14:25 rbash -> bash - -rwxr-xr-x 1 root root 46984 Sep 22 2020 readlink - -rwxr-xr-x 1 root root 75720 Sep 22 2020 rm - -rwxr-xr-x 1 root root 46984 Sep 22 2020 rmdir - -rwxr-xr-x 1 root root 22292 Sep 27 2020 run-parts - -rwxr-xr-x 1 root root 125036 Dec 22 2018 sed - lrwxrwxrwx 1 root root 4 Jul 8 20:43 sh -> dash - -rwxr-xr-x 1 root root 34696 Sep 22 2020 sleep - -rwxr-xr-x 1 root root 83880 Sep 22 2020 stty - -rwsr-xr-x 1 root root 79396 Feb 7 02:38 su - -rwxr-xr-x 1 root root 34696 Sep 22 2020 sync - -rwxr-xr-x 1 root root 602584 Feb 16 21:55 tar - -rwxr-xr-x 1 root root 13860 Sep 27 2020 tempfile - -rwxr-xr-x 1 root root 108520 Sep 22 2020 touch - -rwxr-xr-x 1 root root 34664 Sep 22 2020 true - -rwxr-xr-x 1 root root 17768 Feb 26 04:12 ulockmgr_server - -rwsr-xr-x 1 root root 30236 Feb 7 02:38 umount - -rwxr-xr-x 1 root root 34664 Sep 22 2020 uname - -rwxr-xr-x 2 root root 2346 Mar 2 11:30 uncompress - -rwxr-xr-x 1 root root 153964 Sep 22 2020 vdir - -rwxr-xr-x 1 root root 63024 Feb 7 02:38 wdctl - lrwxrwxrwx 1 root root 8 Nov 6 2019 ypdomainname -> hostname - -rwxr-xr-x 1 root root 1984 Mar 2 11:30 zcat - -rwxr-xr-x 1 root root 1678 Mar 2 11:30 zcmp - -rwxr-xr-x 1 root root 5880 Mar 2 11:30 zdiff - -rwxr-xr-x 1 root root 29 Mar 2 11:30 zegrep - -rwxr-xr-x 1 root root 29 Mar 2 11:30 zfgrep - -rwxr-xr-x 1 root root 2081 Mar 2 11:30 zforce - -rwxr-xr-x 1 root root 7585 Mar 2 11:30 zgrep - -rwxr-xr-x 1 root root 2206 Mar 2 11:30 zless - -rwxr-xr-x 1 root root 1842 Mar 2 11:30 zmore - -rwxr-xr-x 1 root root 4553 Mar 2 11:30 znew -I: user script /srv/workspace/pbuilder/30876/tmp/hooks/D02_print_environment finished + total 5776 + -rwxr-xr-x 1 root root 1367848 Jun 22 2021 bash + -rwxr-xr-x 3 root root 38280 Jul 21 2020 bunzip2 + -rwxr-xr-x 3 root root 38280 Jul 21 2020 bzcat + lrwxrwxrwx 1 root root 6 Jul 21 2020 bzcmp -> bzdiff + -rwxr-xr-x 1 root root 2225 Jul 21 2020 bzdiff + lrwxrwxrwx 1 root root 6 Jul 21 2020 bzegrep -> bzgrep + -rwxr-xr-x 1 root root 4877 Sep 5 2019 bzexe + lrwxrwxrwx 1 root root 6 Jul 21 2020 bzfgrep -> bzgrep + -rwxr-xr-x 1 root root 3775 Jul 21 2020 bzgrep + -rwxr-xr-x 3 root root 38280 Jul 21 2020 bzip2 + -rwxr-xr-x 1 root root 17768 Jul 21 2020 bzip2recover + lrwxrwxrwx 1 root root 6 Jul 21 2020 bzless -> bzmore + -rwxr-xr-x 1 root root 1297 Jul 21 2020 bzmore + -rwxr-xr-x 1 root root 38824 Sep 23 2020 cat + -rwxr-xr-x 1 root root 71624 Sep 23 2020 chgrp + -rwxr-xr-x 1 root root 67528 Sep 23 2020 chmod + -rwxr-xr-x 1 root root 75752 Sep 23 2020 chown + -rwxr-xr-x 1 root root 157960 Sep 23 2020 cp + -rwxr-xr-x 1 root root 128724 Dec 11 2020 dash + -rwxr-xr-x 1 root root 124904 Sep 23 2020 date + -rwxr-xr-x 1 root root 92172 Sep 23 2020 dd + -rwxr-xr-x 1 root root 100752 Sep 23 2020 df + -rwxr-xr-x 1 root root 153964 Sep 23 2020 dir + -rwxr-xr-x 1 root root 83644 Feb 8 2021 dmesg + lrwxrwxrwx 1 root root 8 Nov 8 2019 dnsdomainname -> hostname + lrwxrwxrwx 1 root root 8 Nov 8 2019 domainname -> hostname + -rwxr-xr-x 1 root root 34664 Sep 23 2020 echo + -rwxr-xr-x 1 root root 28 Nov 10 2020 egrep + -rwxr-xr-x 1 root root 34664 Sep 23 2020 false + -rwxr-xr-x 1 root root 28 Nov 10 2020 fgrep + -rwxr-xr-x 1 root root 71928 Feb 8 2021 findmnt + -rwsr-xr-x 1 root root 30112 Feb 27 2021 fusermount + -rwxr-xr-x 1 root root 210488 Nov 10 2020 grep + -rwxr-xr-x 2 root root 2346 Mar 3 2021 gunzip + -rwxr-xr-x 1 root root 6376 Mar 3 2021 gzexe + -rwxr-xr-x 1 root root 100952 Mar 3 2021 gzip + -rwxr-xr-x 1 root root 21916 Nov 8 2019 hostname + -rwxr-xr-x 1 root root 83980 Sep 23 2020 ln + -rwxr-xr-x 1 root root 55572 Feb 8 2020 login + -rwxr-xr-x 1 root root 153964 Sep 23 2020 ls + -rwxr-xr-x 1 root root 153124 Feb 8 2021 lsblk + -rwxr-xr-x 1 root root 96328 Sep 23 2020 mkdir + -rwxr-xr-x 1 root root 79912 Sep 23 2020 mknod + -rwxr-xr-x 1 root root 47048 Sep 23 2020 mktemp + -rwxr-xr-x 1 root root 58920 Feb 8 2021 more + -rwsr-xr-x 1 root root 50720 Feb 8 2021 mount + -rwxr-xr-x 1 root root 13856 Feb 8 2021 mountpoint + -rwxr-xr-x 1 root root 157996 Sep 23 2020 mv + lrwxrwxrwx 1 root root 8 Nov 8 2019 nisdomainname -> hostname + lrwxrwxrwx 1 root root 14 Apr 19 2021 pidof -> /sbin/killall5 + -rwxr-xr-x 1 root root 38824 Sep 23 2020 pwd + lrwxrwxrwx 1 root root 4 Jun 22 2021 rbash -> bash + -rwxr-xr-x 1 root root 46984 Sep 23 2020 readlink + -rwxr-xr-x 1 root root 75720 Sep 23 2020 rm + -rwxr-xr-x 1 root root 46984 Sep 23 2020 rmdir + -rwxr-xr-x 1 root root 22292 Sep 28 2020 run-parts + -rwxr-xr-x 1 root root 125036 Dec 23 2018 sed + lrwxrwxrwx 1 root root 4 Aug 12 07:40 sh -> bash + lrwxrwxrwx 1 root root 4 Aug 12 05:48 sh.distrib -> dash + -rwxr-xr-x 1 root root 34696 Sep 23 2020 sleep + -rwxr-xr-x 1 root root 83880 Sep 23 2020 stty + -rwsr-xr-x 1 root root 79396 Feb 8 2021 su + -rwxr-xr-x 1 root root 34696 Sep 23 2020 sync + -rwxr-xr-x 1 root root 602584 Feb 17 2021 tar + -rwxr-xr-x 1 root root 13860 Sep 28 2020 tempfile + -rwxr-xr-x 1 root root 108520 Sep 23 2020 touch + -rwxr-xr-x 1 root root 34664 Sep 23 2020 true + -rwxr-xr-x 1 root root 17768 Feb 27 2021 ulockmgr_server + -rwsr-xr-x 1 root root 30236 Feb 8 2021 umount + -rwxr-xr-x 1 root root 34664 Sep 23 2020 uname + -rwxr-xr-x 2 root root 2346 Mar 3 2021 uncompress + -rwxr-xr-x 1 root root 153964 Sep 23 2020 vdir + -rwxr-xr-x 1 root root 63024 Feb 8 2021 wdctl + lrwxrwxrwx 1 root root 8 Nov 8 2019 ypdomainname -> hostname + -rwxr-xr-x 1 root root 1984 Mar 3 2021 zcat + -rwxr-xr-x 1 root root 1678 Mar 3 2021 zcmp + -rwxr-xr-x 1 root root 5880 Mar 3 2021 zdiff + -rwxr-xr-x 1 root root 29 Mar 3 2021 zegrep + -rwxr-xr-x 1 root root 29 Mar 3 2021 zfgrep + -rwxr-xr-x 1 root root 2081 Mar 3 2021 zforce + -rwxr-xr-x 1 root root 7585 Mar 3 2021 zgrep + -rwxr-xr-x 1 root root 2206 Mar 3 2021 zless + -rwxr-xr-x 1 root root 1842 Mar 3 2021 zmore + -rwxr-xr-x 1 root root 4553 Mar 3 2021 znew +I: user script /srv/workspace/pbuilder/55977/tmp/hooks/D02_print_environment finished -> Attempting to satisfy build-dependencies -> Creating pbuilder-satisfydepends-dummy package Package: pbuilder-satisfydepends-dummy @@ -685,7 +719,7 @@ Get: 476 http://deb.debian.org/debian bullseye/main i386 texlive-latex-base all 2020.20210202-3 [1120 kB] Get: 477 http://deb.debian.org/debian bullseye/main i386 texlive-latex-recommended all 2020.20210202-3 [14.5 MB] Get: 478 http://deb.debian.org/debian bullseye/main i386 texlive all 2020.20210202-3 [34.2 kB] -Fetched 232 MB in 1min 7s (3470 kB/s) +Fetched 232 MB in 3s (79.1 MB/s) debconf: delaying package configuration, since apt-utils is not installed Selecting previously unselected package bsdextrautils. (Reading database ... (Reading database ... 5% (Reading database ... 10% (Reading database ... 15% (Reading database ... 20% (Reading database ... 25% (Reading database ... 30% (Reading database ... 35% (Reading database ... 40% (Reading database ... 45% (Reading database ... 50% (Reading database ... 55% (Reading database ... 60% (Reading database ... 65% (Reading database ... 70% (Reading database ... 75% (Reading database ... 80% (Reading database ... 85% (Reading database ... 90% (Reading database ... 95% (Reading database ... 100% (Reading database ... 19675 files and directories currently installed.) @@ -2658,8 +2692,25 @@ Writing extended state information... Building tag database... -> Finished parsing the build-deps +Reading package lists... +Building dependency tree... +Reading state information... +The following NEW packages will be installed: + usrmerge +0 upgraded, 1 newly installed, 0 to remove and 0 not upgraded. +Need to get 13.0 kB of archives. +After this operation, 43.0 kB of additional disk space will be used. +Get:1 http://deb.debian.org/debian bullseye/main i386 usrmerge all 25 [13.0 kB] +debconf: delaying package configuration, since apt-utils is not installed +Fetched 13.0 kB in 0s (0 B/s) +Selecting previously unselected package usrmerge. +(Reading database ... (Reading database ... 5% (Reading database ... 10% (Reading database ... 15% (Reading database ... 20% (Reading database ... 25% (Reading database ... 30% (Reading database ... 35% (Reading database ... 40% (Reading database ... 45% (Reading database ... 50% (Reading database ... 55% (Reading database ... 60% (Reading database ... 65% (Reading database ... 70% (Reading database ... 75% (Reading database ... 80% (Reading database ... 85% (Reading database ... 90% (Reading database ... 95% (Reading database ... 100% (Reading database ... 46984 files and directories currently installed.) +Preparing to unpack .../archives/usrmerge_25_all.deb ... +Unpacking usrmerge (25) ... +Setting up usrmerge (25) ... +The system has been successfully converted. I: Building the package -I: Running cd /build/octave-optim-1.6.0/ && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games" HOME="/nonexistent/first-build" dpkg-buildpackage -us -uc -b +I: Running cd /build/octave-optim-1.6.0/ && env PATH="/usr/sbin:/usr/bin:/sbin:/bin:/usr/games:/i/capture/the/path" HOME="/nonexistent/second-build" dpkg-buildpackage -us -uc -b dpkg-buildpackage: info: source package octave-optim dpkg-buildpackage: info: source version 1.6.0-6 dpkg-buildpackage: info: source distribution unstable @@ -2700,7 +2751,7 @@ checking for mkoctfile... /usr/bin/mkoctfile --verbose checking for octave-config... /usr/bin/octave-config -checking for a sed that does not truncate output... /bin/sed +checking for a sed that does not truncate output... /usr/bin/sed checking whether the C++ compiler works... yes checking for C++ compiler default output file name... a.out checking for suffix of executables... @@ -2714,8 +2765,8 @@ checking whether gcc accepts -g... yes checking for gcc option to accept ISO C89... none needed checking how to run the C preprocessor... gcc -E -checking for grep that handles long lines and -e... /bin/grep -checking for egrep... /bin/grep -E +checking for grep that handles long lines and -e... /usr/bin/grep +checking for egrep... /usr/bin/grep -E checking for ANSI C header files... yes checking for sys/types.h... yes checking for sys/stat.h... yes @@ -2749,19 +2800,19 @@ if grep -q DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS numhessian.cc; then echo "The string 'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS' must not be present in source code, but is in numhessian.cc."; exit 1; fi if grep -q DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS __disna_optim__.cc; then echo "The string 'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS' must not be present in source code, but is in __disna_optim__.cc."; exit 1; fi CXXFLAGS="-g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security" /usr/bin/mkoctfile --verbose -c error-helpers.cc -if grep -q DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS numhessian.cc; then echo "The string 'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS' must not be present in source code, but is in numhessian.cc."; exit 1; fi -if grep -q DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS numgradient.cc; then echo "The string 'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS' must not be present in source code, but is in numgradient.cc."; exit 1; fi if grep -q DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS __bfgsmin.cc; then echo "The string 'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS' must not be present in source code, but is in __bfgsmin.cc."; exit 1; fi +if grep -q DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS numgradient.cc; then echo "The string 'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS' must not be present in source code, but is in numgradient.cc."; exit 1; fi if grep -q DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS __disna_optim__.cc; then echo "The string 'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS' must not be present in source code, but is in __disna_optim__.cc."; exit 1; fi +if grep -q DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS numhessian.cc; then echo "The string 'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS' must not be present in source code, but is in numhessian.cc."; exit 1; fi g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security error-helpers.cc -o error-helpers.o -(echo "#include "; echo "int main () {"; sed -e s/DEFUN_DLD/DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g -e s/DEFUNX_DLD/DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g numhessian.cc | g++ -E `/usr/bin/mkoctfile --verbose -p INCFLAGS` -x c++ -iquote '.' -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' -D'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,fname,gname,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' - | sed -e '/.*DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/!D'; echo "}";) | g++ -x c++ -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)=printf("%c" #name "\n@c " #name " numhessian.cc\n" doc "\n\n", 0x1D);' -o numhessian.bin - +(echo "#include "; echo "int main () {"; sed -e s/DEFUN_DLD/DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g -e s/DEFUNX_DLD/DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g __bfgsmin.cc | g++ -E `/usr/bin/mkoctfile --verbose -p INCFLAGS` -x c++ -iquote '.' -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' -D'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,fname,gname,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' - | sed -e '/.*DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/!D'; echo "}";) | g++ -x c++ -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)=printf("%c" #name "\n@c " #name " __bfgsmin.cc\n" doc "\n\n", 0x1D);' -o __bfgsmin.bin - (echo "#include "; echo "int main () {"; sed -e s/DEFUN_DLD/DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g -e s/DEFUNX_DLD/DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g numgradient.cc | g++ -E `/usr/bin/mkoctfile --verbose -p INCFLAGS` -x c++ -iquote '.' -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' -D'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,fname,gname,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' - | sed -e '/.*DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/!D'; echo "}";) | g++ -x c++ -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)=printf("%c" #name "\n@c " #name " numgradient.cc\n" doc "\n\n", 0x1D);' -o numgradient.bin - (echo "#include "; echo "int main () {"; sed -e s/DEFUN_DLD/DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g -e s/DEFUNX_DLD/DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g __disna_optim__.cc | g++ -E `/usr/bin/mkoctfile --verbose -p INCFLAGS` -x c++ -iquote '.' -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' -D'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,fname,gname,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' - | sed -e '/.*DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/!D'; echo "}";) | g++ -x c++ -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)=printf("%c" #name "\n@c " #name " __disna_optim__.cc\n" doc "\n\n", 0x1D);' -o __disna_optim__.bin - -(echo "#include "; echo "int main () {"; sed -e s/DEFUN_DLD/DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g -e s/DEFUNX_DLD/DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g __bfgsmin.cc | g++ -E `/usr/bin/mkoctfile --verbose -p INCFLAGS` -x c++ -iquote '.' -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' -D'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,fname,gname,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' - | sed -e '/.*DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/!D'; echo "}";) | g++ -x c++ -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)=printf("%c" #name "\n@c " #name " __bfgsmin.cc\n" doc "\n\n", 0x1D);' -o __bfgsmin.bin - -(echo "### This file is generated automatically from the"; echo "### corresponding .cc file by a Makefile rule."; echo ""; ./__bfgsmin.bin) > __bfgsmin.cc.docstrings +(echo "#include "; echo "int main () {"; sed -e s/DEFUN_DLD/DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g -e s/DEFUNX_DLD/DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/g numhessian.cc | g++ -E `/usr/bin/mkoctfile --verbose -p INCFLAGS` -x c++ -iquote '.' -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' -D'DEFUNX_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,fname,gname,args,nargout,doc)=DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)' - | sed -e '/.*DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS/!D'; echo "}";) | g++ -x c++ -D'DEFUN_DLD_REPLACEMENT_FOR_PACKAGE_DOCS(name,doc)=printf("%c" #name "\n@c " #name " numhessian.cc\n" doc "\n\n", 0x1D);' -o numhessian.bin - +(echo "### This file is generated automatically from the"; echo "### corresponding .cc file by a Makefile rule."; echo ""; ./numhessian.bin) > numhessian.cc.docstrings (echo "### This file is generated automatically from the"; echo "### corresponding .cc file by a Makefile rule."; echo ""; ./__disna_optim__.bin) > __disna_optim__.cc.docstrings +(echo "### This file is generated automatically from the"; echo "### corresponding .cc file by a Makefile rule."; echo ""; ./__bfgsmin.bin) > __bfgsmin.cc.docstrings (echo "### This file is generated automatically from the"; echo "### corresponding .cc file by a Makefile rule."; echo ""; ./numgradient.bin) > numgradient.cc.docstrings -(echo "### This file is generated automatically from the"; echo "### corresponding .cc file by a Makefile rule."; echo ""; ./numhessian.bin) > numhessian.cc.docstrings ./munge-texi.pl MFDOCSTRINGS __bfgsmin.cc.docstrings numgradient.cc.docstrings numhessian.cc.docstrings __disna_optim__.cc.docstrings < ../doc/optim.txi > ../doc/optim.texi makeinfo ../doc/optim.texi -o ../doc/optim.info ../doc/macros.texi:51: warning: redefining Texinfo language command: @seealso @@ -2769,59 +2820,10 @@ CXXFLAGS="-g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security" /usr/bin/mkoctfile --verbose numgradient.cc error-helpers.o CXXFLAGS="-g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security" /usr/bin/mkoctfile --verbose numhessian.cc error-helpers.o CXXFLAGS="-g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security" /usr/bin/mkoctfile --verbose -llapack -lblas __disna_optim__.cc error-helpers.o -g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security __disna_optim__.cc -o /tmp/oct-kkUHkj.o -g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security numgradient.cc -o /tmp/oct-z27AJW.o -g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security __bfgsmin.cc -o /tmp/oct-hL8tH4.o -g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security numhessian.cc -o /tmp/oct-DRdsoc.o -In file included from __bfgsmin.cc:30: -__bfgsmin.cc: In function 'int __bfgsmin_obj(double&, std::string, octave_value_list, ColumnVector, int)': -error-helpers.h:119:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] - 119 | if (error_state) \ - | ^~~~~~~~~~~ -__bfgsmin.cc:39:3: note: in expansion of macro 'SET_ERR' - 39 | SET_ERR (f_return = OCTAVE__FEVAL (f, f_args_new), err); - | ^~~~~~~ -In file included from /usr/include/octave-6.2.0/octave/oct.h:36, - from __bfgsmin.cc:23: -/usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here - 495 | extern OCTINTERP_API int error_state; - | ^~~~~~~~~~~ -In file included from __bfgsmin.cc:30: -error-helpers.h:121:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] - 121 | error_state = 0; \ - | ^~~~~~~~~~~ -__bfgsmin.cc:39:3: note: in expansion of macro 'SET_ERR' - 39 | SET_ERR (f_return = OCTAVE__FEVAL (f, f_args_new), err); - | ^~~~~~~ -In file included from /usr/include/octave-6.2.0/octave/oct.h:36, - from __bfgsmin.cc:23: -/usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here - 495 | extern OCTINTERP_API int error_state; - | ^~~~~~~~~~~ -In file included from __bfgsmin.cc:30: -error-helpers.h:119:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] - 119 | if (error_state) \ - | ^~~~~~~~~~~ -__bfgsmin.cc:42:7: note: in expansion of macro 'SET_ERR' - 42 | SET_ERR (obj = f_return(0).double_value(), err); - | ^~~~~~~ -In file included from /usr/include/octave-6.2.0/octave/oct.h:36, - from __bfgsmin.cc:23: -/usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here - 495 | extern OCTINTERP_API int error_state; - | ^~~~~~~~~~~ -In file included from __bfgsmin.cc:30: -error-helpers.h:121:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] - 121 | error_state = 0; \ - | ^~~~~~~~~~~ -__bfgsmin.cc:42:7: note: in expansion of macro 'SET_ERR' - 42 | SET_ERR (obj = f_return(0).double_value(), err); - | ^~~~~~~ -In file included from /usr/include/octave-6.2.0/octave/oct.h:36, - from __bfgsmin.cc:23: -/usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here - 495 | extern OCTINTERP_API int error_state; - | ^~~~~~~~~~~ +g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security __bfgsmin.cc -o /tmp/oct-dvbGdm.o +g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security numgradient.cc -o /tmp/oct-Y8f4ds.o +g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security numhessian.cc -o /tmp/oct-J90epu.o +g++ -c -Wdate-time -D_FORTIFY_SOURCE=2 -fPIC -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security __disna_optim__.cc -o /tmp/oct-qtq2iD.o In file included from __disna_optim__.cc:21: __disna_optim__.cc: In function 'octave_value_list F__disna_optim__(const octave_value_list&, int)': error-helpers.h:48:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] @@ -2847,6 +2849,31 @@ /usr/include/octave-6.2.0/octave/../octave/error.h:495:26: note: declared here 495 | extern OCTINTERP_API int error_state; | ^~~~~~~~~~~ +In file included from numgradient.cc:24: +numgradient.cc: In function 'bool any_bad_argument(const octave_value_list&)': +error-helpers.h:119:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] + 119 | if (error_state) \ + | ^~~~~~~~~~~ +numgradient.cc:47:7: note: in expansion of macro 'SET_ERR' + 47 | SET_ERR (tmp = args(2).int_value(), err); + | ^~~~~~~ +In file included from /usr/include/octave-6.2.0/octave/oct.h:36, + from numgradient.cc:18: +/usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here + 495 | extern OCTINTERP_API int error_state; + | ^~~~~~~~~~~ +In file included from numgradient.cc:24: +error-helpers.h:121:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] + 121 | error_state = 0; \ + | ^~~~~~~~~~~ +numgradient.cc:47:7: note: in expansion of macro 'SET_ERR' + 47 | SET_ERR (tmp = args(2).int_value(), err); + | ^~~~~~~ +In file included from /usr/include/octave-6.2.0/octave/oct.h:36, + from numgradient.cc:18: +/usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here + 495 | extern OCTINTERP_API int error_state; + | ^~~~~~~~~~~ In file included from __disna_optim__.cc:21: error-helpers.h:48:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] 48 | if (error_state) \ @@ -2883,60 +2910,84 @@ /usr/include/octave-6.2.0/octave/../octave/error.h:495:26: note: declared here 495 | extern OCTINTERP_API int error_state; | ^~~~~~~~~~~ -In file included from numhessian.cc:24: -numhessian.cc: In function 'bool any_bad_argument(const octave_value_list&)': +In file included from __bfgsmin.cc:30: +__bfgsmin.cc: In function 'int __bfgsmin_obj(double&, std::string, octave_value_list, ColumnVector, int)': error-helpers.h:119:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] 119 | if (error_state) \ | ^~~~~~~~~~~ -numhessian.cc:47:7: note: in expansion of macro 'SET_ERR' - 47 | SET_ERR (tmp = args(2).int_value(), err); +__bfgsmin.cc:39:3: note: in expansion of macro 'SET_ERR' + 39 | SET_ERR (f_return = OCTAVE__FEVAL (f, f_args_new), err); + | ^~~~~~~ +In file included from /usr/include/octave-6.2.0/octave/oct.h:36, + from __bfgsmin.cc:23: +/usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here + 495 | extern OCTINTERP_API int error_state; + | ^~~~~~~~~~~ +In file included from __bfgsmin.cc:30: +error-helpers.h:121:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] + 121 | error_state = 0; \ + | ^~~~~~~~~~~ +__bfgsmin.cc:39:3: note: in expansion of macro 'SET_ERR' + 39 | SET_ERR (f_return = OCTAVE__FEVAL (f, f_args_new), err); + | ^~~~~~~ +In file included from /usr/include/octave-6.2.0/octave/oct.h:36, + from __bfgsmin.cc:23: +/usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here + 495 | extern OCTINTERP_API int error_state; + | ^~~~~~~~~~~ +In file included from __bfgsmin.cc:30: +error-helpers.h:119:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] + 119 | if (error_state) \ + | ^~~~~~~~~~~ +__bfgsmin.cc:42:7: note: in expansion of macro 'SET_ERR' + 42 | SET_ERR (obj = f_return(0).double_value(), err); | ^~~~~~~ In file included from /usr/include/octave-6.2.0/octave/oct.h:36, - from numhessian.cc:18: + from __bfgsmin.cc:23: /usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here 495 | extern OCTINTERP_API int error_state; | ^~~~~~~~~~~ -In file included from numhessian.cc:24: +In file included from __bfgsmin.cc:30: error-helpers.h:121:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] 121 | error_state = 0; \ | ^~~~~~~~~~~ -numhessian.cc:47:7: note: in expansion of macro 'SET_ERR' - 47 | SET_ERR (tmp = args(2).int_value(), err); +__bfgsmin.cc:42:7: note: in expansion of macro 'SET_ERR' + 42 | SET_ERR (obj = f_return(0).double_value(), err); | ^~~~~~~ In file included from /usr/include/octave-6.2.0/octave/oct.h:36, - from numhessian.cc:18: + from __bfgsmin.cc:23: /usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here 495 | extern OCTINTERP_API int error_state; | ^~~~~~~~~~~ -In file included from numgradient.cc:24: -numgradient.cc: In function 'bool any_bad_argument(const octave_value_list&)': +In file included from numhessian.cc:24: +numhessian.cc: In function 'bool any_bad_argument(const octave_value_list&)': error-helpers.h:119:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] 119 | if (error_state) \ | ^~~~~~~~~~~ -numgradient.cc:47:7: note: in expansion of macro 'SET_ERR' +numhessian.cc:47:7: note: in expansion of macro 'SET_ERR' 47 | SET_ERR (tmp = args(2).int_value(), err); | ^~~~~~~ In file included from /usr/include/octave-6.2.0/octave/oct.h:36, - from numgradient.cc:18: + from numhessian.cc:18: /usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here 495 | extern OCTINTERP_API int error_state; | ^~~~~~~~~~~ -In file included from numgradient.cc:24: +In file included from numhessian.cc:24: error-helpers.h:121:13: warning: 'error_state' is deprecated: [6]: this variable is obsolete and always has the value 0 [-Wdeprecated-declarations] 121 | error_state = 0; \ | ^~~~~~~~~~~ -numgradient.cc:47:7: note: in expansion of macro 'SET_ERR' +numhessian.cc:47:7: note: in expansion of macro 'SET_ERR' 47 | SET_ERR (tmp = args(2).int_value(), err); | ^~~~~~~ In file included from /usr/include/octave-6.2.0/octave/oct.h:36, - from numgradient.cc:18: + from numhessian.cc:18: /usr/include/octave-6.2.0/octave/error.h:495:26: note: declared here 495 | extern OCTINTERP_API int error_state; | ^~~~~~~~~~~ -g++ -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __disna_optim__.oct error-helpers.o /tmp/oct-kkUHkj.o -shared -Wl,-Bsymbolic -Wl,-z,relro -llapack -lblas -L/usr/lib/i386-linux-gnu -Wl,-z,relro -g++ -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security -o numhessian.oct error-helpers.o /tmp/oct-DRdsoc.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/i386-linux-gnu -Wl,-z,relro -g++ -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security -o numgradient.oct error-helpers.o /tmp/oct-z27AJW.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/i386-linux-gnu -Wl,-z,relro -g++ -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __bfgsmin.oct error-helpers.o /tmp/oct-hL8tH4.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/i386-linux-gnu -Wl,-z,relro +g++ -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __disna_optim__.oct error-helpers.o /tmp/oct-qtq2iD.o -shared -Wl,-Bsymbolic -Wl,-z,relro -llapack -lblas -L/usr/lib/i386-linux-gnu -Wl,-z,relro +g++ -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security -o numgradient.oct error-helpers.o /tmp/oct-Y8f4ds.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/i386-linux-gnu -Wl,-z,relro +g++ -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security -o numhessian.oct error-helpers.o /tmp/oct-J90epu.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/i386-linux-gnu -Wl,-z,relro +g++ -I/usr/include/octave-6.2.0/octave/.. -I/usr/include/octave-6.2.0/octave -pthread -fopenmp -mieee-fp -g -O2 -fdebug-prefix-map=/build/octave-optim-1.6.0=. -fstack-protector-strong -Wformat -Werror=format-security -o __bfgsmin.oct error-helpers.o /tmp/oct-dvbGdm.o -shared -Wl,-Bsymbolic -Wl,-z,relro -L/usr/lib/i386-linux-gnu -Wl,-z,relro rm numhessian.bin numgradient.cc.docstrings MFDOCSTRINGS numgradient.bin __disna_optim__.cc.docstrings __disna_optim__.bin numhessian.cc.docstrings __bfgsmin.bin __bfgsmin.cc.docstrings make[2]: Leaving directory '/build/octave-optim-1.6.0/src' copyfile /build/octave-optim-1.6.0/./src/__bfgsmin.oct /build/octave-optim-1.6.0/./src/__disna_optim__.oct /build/octave-optim-1.6.0/./src/numgradient.oct /build/octave-optim-1.6.0/./src/numhessian.oct /build/octave-optim-1.6.0/./inst/i686-pc-linux-gnu-api-v55 @@ -2945,7 +2996,7 @@ ( cd src ; ./configure ; make doc ) checking for mkoctfile... mkoctfile checking for octave-config... octave-config -checking for a sed that does not truncate output... /bin/sed +checking for a sed that does not truncate output... /usr/bin/sed checking whether the C++ compiler works... yes checking for C++ compiler default output file name... a.out checking for suffix of executables... @@ -2959,8 +3010,8 @@ checking whether gcc accepts -g... yes checking for gcc option to accept ISO C89... none needed checking how to run the C preprocessor... gcc -E -checking for grep that handles long lines and -e... /bin/grep -checking for egrep... /bin/grep -E +checking for grep that handles long lines and -e... /usr/bin/grep +checking for egrep... /usr/bin/grep -E checking for ANSI C header files... yes checking for sys/types.h... yes checking for sys/stat.h... yes @@ -2997,181 +3048,6 @@ dh_octave_check -O--buildsystem=octave Checking package... Checking m files ... -[inst/quadprog.m] ->>>>> /build/octave-optim-1.6.0/inst/quadprog.m -***** test - H= diag([1; 0]); - f = [3; 4]; - A= [-1 -3; 2 5; 3 4]; - b = [-15; 100; 80]; - l= zeros(2,1); - [x,fval,exitflag,output] = quadprog(H,f,A,b,[],[],l,[]); - assert(x,[0;5]) - assert(fval,20) - assert(exitflag,1) - assert(output.iterations,1) -***** demo - C = [0.9501 0.7620 0.6153 0.4057 - 0.2311 0.4564 0.7919 0.9354 - 0.6068 0.0185 0.9218 0.9169 - 0.4859 0.8214 0.7382 0.4102 - 0.8912 0.4447 0.1762 0.8936]; - %% Linear Inequality Constraints - d = [0.0578; 0.3528; 0.8131; 0.0098; 0.1388]; - A =[0.2027 0.2721 0.7467 0.4659 - 0.1987 0.1988 0.4450 0.4186 - 0.6037 0.0152 0.9318 0.8462]; - b =[0.5251; 0.2026; 0.6721]; - %% Linear Equality Constraints - Aeq = [3 5 7 9]; - beq = 4; - %% Bound constraints - lb = -0.1*ones(4,1); - ub = ones(4,1); - H = C' * C; - f = -C' * d; - [x, obj, flag, output, lambda]=quadprog (H, f, A, b, Aeq, beq, lb, ub) -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lsqcurvefit.m] ->>>>> /build/octave-optim-1.6.0/inst/lsqcurvefit.m -***** test - xdata = [0 .3 .8 1.1 1.6 2.3]; - ydata = [.82 .72 .63 .60 .55 .50]; - yhat = @(p,x) p(1) + p(2)*exp(-x); - opt = optimset('TolFun',1e-100); - [p, resnorm, residual] = lsqcurvefit(yhat,[1 1], xdata, ydata,[0 0],[],opt); - assert (p, [ 0.47595; 0.34132], 1e-5) - assert (resnorm, 3.2419e-004, 1e-8) - assert(residual, [-2.7283e-003, 8.8079e-003, -6.8307e-004, -1.0432e-002, -5.1366e-003, 1.0172e-002], 1e-5) -***** test - problem.solver = "lsqcurvefit"; - problem.objective = @(p,x) p(1) + p(2)*exp(-x); - problem.x0 = [1, 1]; - problem.xdata = [0 .3 .8 1.1 1.6 2.3]; - problem.ydata = [.82 .72 .63 .60 .55 .50]; - problem.lb = [0, 0]; - problem.ub = []; - problem.options = optimset('TolFun',1e-100); - [p, resnorm, residual] = lsqcurvefit (problem); - assert (p, [ 0.47595; 0.34132], 1e-5) - assert (resnorm, 3.2419e-004, 1e-8) - assert(residual, [-2.7283e-003, 8.8079e-003, -6.8307e-004, -1.0432e-002, -5.1366e-003, 1.0172e-002], 1e-5) -***** demo - %% Example for user specified Jacobian. - - %% independents - x = [1:10:100]'; - %% observed data - y =[9.2160e-001, 3.3170e-001, 8.9789e-002, 2.8480e-002, 2.6055e-002,... - 8.3641e-003, 4.2362e-003, 3.1693e-003, 1.4739e-004, 2.9406e-004]'; - %% initial values: - p0=[0.8; 0.05]; - %% bounds - lb=[0; 0]; ub=[]; - %% Jacobian setting - opts = optimset ("Jacobian", "on") - - %% model function: - function [F,J] = myfun (p, x) - F = p(1) * exp (-p(2) * x); - if nargout > 1 - J = [exp(- p(2) * x), - p(1) * x .* exp(- p(2) * x)]; - endif - endfunction - - [c, resnorm, residual, flag, output, lambda, jacob] = ... - lsqcurvefit (@ (varargin) myfun(varargin{:}), p0, x, y, lb, ub, opts) -2 tests, 2 passed, 0 known failure, 0 skipped -[inst/jacobs.m] ->>>>> /build/octave-optim-1.6.0/inst/jacobs.m -***** assert (jacobs (1, @(x) x), 1) -***** assert (jacobs (6, @(x) x^2), 12) -***** assert (jacobs ([1; 1], @(x) [x(1)^2; x(1)*x(2)]), [2, 0; 1, 1]) -***** assert (jacobs ([1; 2], @(x) [x(1)^2 + x(2); x(2)*exp(x(1))]), [2, 1; 2*exp(1), exp(1)]) -***** error jacobs () -***** error jacobs (1) -***** error jacobs (1, 2, 3, 4) -***** error jacobs (@sin, 1, [1, 1]) -***** error jacobs (@sin, 1, ones(2, 2)) -***** demo - # Relative error against several h-values - k = 3:20; h = 10 .^ (-k); x = 0.3*pi; - err = zeros (1, numel (k)); - for count = 1 : numel (k) - err(count) = abs (jacobs (x, @sin, struct ("h", h(count))) - cos (x)) / abs (cos (x)) + eps; - endfor - loglog (h, err); grid minor; - xlabel ("h"); ylabel ("|Df(x) - cos(x)| / |cos(x)|") - title ("f(x)=sin(x), f'(x)=cos(x) at x = 0.3pi") -9 tests, 9 passed, 0 known failure, 0 skipped -[inst/statget.m] ->>>>> /build/octave-optim-1.6.0/inst/statget.m -***** shared opts - opts = statset ("maxit", 100); -***** assert (statget (opts, "maxit"), 100) -***** assert (statget (opts, "MaxITer"), 100) -***** assert (statget (opts, "TolFun"), []) -***** assert (statget (opts, "TolFun", 1e-3), 1e-3) -***** error statget () -***** error statget (1) -***** error statget (1,2,3,4,5) -***** error statget (1, "name") -***** error statget (struct (), 2) -9 tests, 9 passed, 0 known failure, 0 skipped -[inst/de_min.m] ->>>>> /build/octave-optim-1.6.0/inst/de_min.m -***** function result = f(x); - result = 100 * (x(2) - x(1)^2)^2 + (1 - x(1))^2; -***** test - tol = 1.0e-4; - ctl.tol = 0.0; - ctl.VTR = 1.0e-6; - ctl.XVmin = [-2 -2]; - ctl.XVmax = [ 2 2]; - rand("state", 11) - [x, obj_value, nfeval, convergence] = de_min (@f, ctl); - assert (convergence == 1); - assert (f(x) == obj_value); - assert (obj_value < ctl.VTR); -***** demo - ## define a simple example function - f = @(x) peaks(x(1), x(2)); - ## plot the function to see where the minimum might be - peaks() - ## first we set the region where we expect the minimum - ctl.XVmin = [-3 -3]; - ctl.XVmax = [ 3 3]; - ## and solve it with de_min - [x, obj_value, nfeval, convergence] = de_min (f, ctl) -1 test, 1 passed, 0 known failure, 0 skipped -[inst/__all_stat_opts__.m] ->>>>> /build/octave-optim-1.6.0/inst/__all_stat_opts__.m -***** assert (1) -1 test, 1 passed, 0 known failure, 0 skipped -[inst/nonlin_curvefit.m] ->>>>> /build/octave-optim-1.6.0/inst/nonlin_curvefit.m -***** demo - ## Example for linear inequality constraints - ## (see also the same example in 'demo nonlin_residmin') - - ## independents and observations - indep = 1:5; - obs = [1, 2, 4, 7, 14]; - ## model function: - f = @ (p, x) p(1) * exp (p(2) * x); - ## initial values: - init = [.25; .25]; - ## linear constraints, A.' * parametervector + B >= 0 - A = [1; -1]; B = 0; # p(1) >= p(2); - settings = optimset ("inequc", {A, B}); - - ## start optimization - [p, model_values, cvg, outp] = nonlin_curvefit (f, init, indep, obs, settings) -***** test - t = optim_problems ().curve.schittkowski_327; - cp = nonlin_curvefit (@ (p, x) t.f (x, p), t.init_p, t.data.x, t.data.y, optimset ("lbound", t.non_strict_inequc.bounds(:, 1), "ubound", t.non_strict_inequc.bounds(:, 2), "inequc", {t.non_strict_inequc.general}, "dfdp", @ (p, x) t.dfdp (x, p))); - assert (cp, t.result.p, [.0001; .001]); -1 test, 1 passed, 0 known failure, 0 skipped [inst/vfzero.m] >>>>> /build/octave-optim-1.6.0/inst/vfzero.m ***** shared opt0 @@ -3179,76 +3055,6 @@ ***** assert(vfzero(@cos, [0, 3], opt0), pi/2, 10*eps) ***** assert(vfzero(@(x) x^(1/3) - 1e-8, [0,1], opt0), 1e-24, 1e-22*eps) 2 tests, 2 passed, 0 known failure, 0 skipped -[inst/fmincon.m] ->>>>> /build/octave-optim-1.6.0/inst/fmincon.m -***** demo - ## Example for default optimization (Levenberg/Marquardt with - ## BFGS), one non-linear equality constraint. Constrained optimum is - ## at p = [0; 1]. - objective_function = @ (p) p(1)^2 + p(2)^2; - pin = [-2; 5]; - constraint_function = @ (p) - (p(1)^2 + 1 - p(2)); - [p, objf, cvg, outp] = fmincon (objective_function, pin, [], [], [], [], [], [], @ (p) {[], constraint_function(p)}{:}) -***** test - ## equality constraint - objective_function = @ (p) p(1)^2 + p(2)^2; - pin = [-2; 5]; - constraint_function = @ (p) - (p(1)^2 + 1 - p(2)); - [p, objf, cvg, outp] = fmincon (objective_function, pin, [], [], [], [], [], [], @ (p) {[], constraint_function(p)}{:}, optimset ("Algorithm", "lm_feasible")); - assert (p, [0; 1], 1e-6) -***** test - ## inequality constraint - objective_function = @ (p) p(1)^2 + p(2)^2; - pin = [2; 6]; - constraint_function = @ (p) p(1)^2 + 1 - p(2); - [p, objf, cvg, outp] = fmincon (objective_function, pin, [], [], [], [], [], [], @ (p) {constraint_function(p), []}{:}, optimset ("Algorithm", "lm_feasible")); - assert (p, [0; 1], 1e-6) -***** test - ## independents - indep = 1:5; - ## objective function: - f = @ (p) sumsq (p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]); - ## initial values: - init = [.25; .25]; - ## linear constraints, A.' * parametervector + B >= 0 - A = [1; -1]; B = 0; # p(1) >= p(2); - - assert (fmincon (f, init, -A, B), [.6203; .6203], .0001); -***** test - ## problem structure - indep = 1:5; - problem = struct ("objective", - @ (p) sumsq (p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]), - "x0", [.25; .25], - "Aineq", [-1; 1], - "bineq", 0, - "solver", "fmincon"); - assert (fmincon (problem), [.6203; .6203], .0001); -***** test - ## Octave sqp solver with a lot of inequality constraints - objf = @ (p) sumsq (p(4:9)); - init = [300; -100; -.1997; -127; -151; 379; 421; 460; 426]; - lbound = [-Inf; -Inf; -Inf; 0; 0; 0; 0; 0; 0]; - inequc = @ (p) vertcat ( ... - p(1) + p(2) * exp (-5 * p(3)) + p(4) - 127, - p(1) + p(2) * exp (-3 * p(3)) + p(5) - 151, - p(1) + p(2) * exp (-p(3)) + p(6) - 379, - p(1) + p(2) * exp (p(3)) + p(7) - 421, - p(1) + p(2) * exp (3 * p(3)) + p(8) - 460, - p(1) + p(2) * exp (5 * p(3)) + p(9) - 426, - -p(1) - p(2) * exp (-5 * p(3)) + p(4) + 127, - -p(1) - p(2) * exp (-3 * p(3)) + p(5) + 151, - -p(1) - p(2) * exp (-p(3)) + p(6) + 379, - -p(1) - p(2) * exp (p(3)) + p(7) + 421, - -p(1) - p(2) * exp (3 * p(3)) + p(8) + 460, - -p(1) - p(2) * exp (5 * p(3)) + p(9) + 426); - [p, objf, cvg, outp] = fmincon (objf, init, [], [], [], [], lbound, [], - @ (p) {- inequc(p), []}{:}, - optimset ("Algorithm", "octave_sqp")); - assert (p, [5.2330e+02; -1.5694e+02; -1.9966e-01; 2.9607e+01; - 8.6615e+01; 4.7326e+01; 2.6235e+01; 2.2915e+01; - 3.9470e+01], .01); -5 tests, 5 passed, 0 known failure, 0 skipped [inst/lsqnonlin.m] >>>>> /build/octave-optim-1.6.0/inst/lsqnonlin.m ***** test @@ -3307,18 +3113,6 @@ ***** error statset ("%NOT_A_REAL_FUNCTION_NAME%") ***** warning statset ("foobar", 13); 5 tests, 5 passed, 0 known failure, 0 skipped -[inst/polyconf.m] ->>>>> /build/octave-optim-1.6.0/inst/polyconf.m -***** test - # data from Hocking, RR, "Methods and Applications of Linear Models" - temperature=[40;40;40;45;45;45;50;50;50;55;55;55;60;60;60;65;65;65]; - strength=[66.3;64.84;64.36;69.70;66.26;72.06;73.23;71.4;68.85;75.78;72.57;76.64;78.87;77.37;75.94;78.82;77.13;77.09]; - [p,s] = polyfit(temperature,strength,1); - [y,dy] = polyconf(p,40,s,0.05,'ci'); - assert([y,dy],[66.15396825396826,1.71702862681486],200*eps); - [y,dy] = polyconf(p,40,s,0.05,'pi'); - assert(dy,4.45345484470743,200*eps); -1 test, 1 passed, 0 known failure, 0 skipped [inst/linprog.m] >>>>> /build/octave-optim-1.6.0/inst/linprog.m ***** test @@ -3364,91 +3158,27 @@ beq = 0; assert (linprog (f, A, b, Aeq, beq, lb, ub), x_exp); 3 tests, 3 passed, 0 known failure, 0 skipped -[inst/nlinfit.m] ->>>>> /build/octave-optim-1.6.0/inst/nlinfit.m -***** test - modelfun = @(b, x) (b(1) + b(2) * exp (- b(3) * x)); - b = [1;3;2]; - v = rande ("state"); - rande ("seed", 1234); - xdata = exprnd (2,100,1); - rande ("state", v); - v = randn ("state"); - randn ("seed", 1234); - ydata = modelfun (b,xdata) + normrnd (0,0.1,100,1); - randn ("state", v); - beta0 = [2;2;3]; - beta = nlinfit(xdata,ydata,modelfun,beta0); - assert (beta, [1;3;2], 1e-1) -***** demo - modelfun = @(b, x) (b(1) + b(2) * exp (- b(3) * x)); - %% actual value - beta_without_noise = [1; 3; 2]; - x = [3.49622; 0.33751; 1.25675; 3.66981; 0.26237; 5.51095; ... - 2.11407; 1.48774; 6.22436; 2.04519]; - y_actual = modelfun (beta_without_noise, x); - noise = [0.176110; -0.066850; 0.231000; -0.047570; -0.108230; ... - 0.122790; 0.062940; 0.151510; 0.116010; -0.097460]; - y_noisy = y_actual + noise; - %% initial guess - beta0 = [2; 2; 2]; - %% weights vector - weights = [5; 16; 1; 20; 12; 11; 17; 8; 11; 13]; - [beta, R, J, covb, mse] = nlinfit (x, y_noisy, modelfun, beta0) - [beta_w, R_w, J_w, covb_w, mse_w] = nlinfit (x, y_noisy, modelfun, beta0, [], "weights", weights) -1 test, 1 passed, 0 known failure, 0 skipped -[inst/lsqlin.m] ->>>>> /build/octave-optim-1.6.0/inst/lsqlin.m -***** test -***** shared C,d,A,b - C = [0.9501,0.7620,0.6153,0.4057;... - 0.2311,0.4564,0.7919,0.9354;... - 0.6068,0.0185,0.9218,0.9169;... - 0.4859,0.8214,0.7382,0.4102;... - 0.8912,0.4447,0.1762,0.8936]; - d = [0.0578; 0.3528; 0.8131; 0.0098; 0.1388]; - A =[0.2027, 0.2721, 0.7467, 0.4659;... - 0.1987, 0.1988, 0.4450, 0.4186;... - 0.6037 , 0.0152, 0.9318, 0.8462]; - b =[0.5251;0.2026;0.6721]; - Aeq = [3, 5, 7, 9]; - beq = 4; - lb = -0.1*ones(4,1); - ub = 2*ones(4,1); - [x,resnorm,residual,exitflag] = lsqlin(C,d,A,b,Aeq,beq,lb,ub); - assert(x,[-0.10000; -0.10000; 0.15991; 0.40896],10e-5) - assert(resnorm,0.16951,10e-5) - assert(residual, [0.035297; 0.087623; -0.353251; 0.145270; 0.121232],10e-5) - assert(exitflag,1) -warning: colon arguments should be scalars -warning: called from - null at line 67 column 14 - quadprog at line 302 column 13 - lsqlin at line 123 column 21 - __test__ at line 17 column 32 - test at line 677 column 11 - /tmp/tmp.KPwpj3PXRx at line 126 column 31 - +[inst/quadprog.m] +>>>>> /build/octave-optim-1.6.0/inst/quadprog.m ***** test - Aeq = []; - beq = []; - lb = []; - ub = []; - x0 = 0.1*ones(4,1); - x = lsqlin(C,d,A,b,Aeq,beq,lb,ub,x0); - [x,resnorm,residual,exitflag] = lsqlin(C,d,A,b,Aeq,beq,lb,ub,x0); - assert(x,[ 0.12986; -0.57569 ; 0.42510; 0.24384],10e-5) - assert(resnorm,0.017585,10e-5) - assert(residual, [-0.0126033; -0.0208040; -0.1295084; -0.0057389; 0.01372462],10e-5) + H= diag([1; 0]); + f = [3; 4]; + A= [-1 -3; 2 5; 3 4]; + b = [-15; 100; 80]; + l= zeros(2,1); + [x,fval,exitflag,output] = quadprog(H,f,A,b,[],[],l,[]); + assert(x,[0;5]) + assert(fval,20) assert(exitflag,1) + assert(output.iterations,1) ***** demo C = [0.9501 0.7620 0.6153 0.4057 0.2311 0.4564 0.7919 0.9354 0.6068 0.0185 0.9218 0.9169 0.4859 0.8214 0.7382 0.4102 0.8912 0.4447 0.1762 0.8936]; - d = [0.0578; 0.3528; 0.8131; 0.0098; 0.1388]; %% Linear Inequality Constraints + d = [0.0578; 0.3528; 0.8131; 0.0098; 0.1388]; A =[0.2027 0.2721 0.7467 0.4659 0.1987 0.1988 0.4450 0.4186 0.6037 0.0152 0.9318 0.8462]; @@ -3459,8 +3189,10 @@ %% Bound constraints lb = -0.1*ones(4,1); ub = ones(4,1); - [x, resnorm, residual, flag, output, lambda] = lsqlin (C, d, A, b, Aeq, beq, lb, ub) -2 tests, 2 passed, 0 known failure, 0 skipped + H = C' * C; + f = -C' * d; + [x, obj, flag, output, lambda]=quadprog (H, f, A, b, Aeq, beq, lb, ub) +1 test, 1 passed, 0 known failure, 0 skipped [inst/residmin_stat.m] >>>>> /build/octave-optim-1.6.0/inst/residmin_stat.m ***** test @@ -3586,45 +3318,300 @@ assert (info.corp, corp, .0001); 4 tests, 4 passed, 0 known failure, 0 skipped -[inst/nonlin_min.m] ->>>>> /build/octave-optim-1.6.0/inst/nonlin_min.m +[inst/__all_stat_opts__.m] +>>>>> /build/octave-optim-1.6.0/inst/__all_stat_opts__.m +***** assert (1) +1 test, 1 passed, 0 known failure, 0 skipped +[inst/lsqlin.m] +>>>>> /build/octave-optim-1.6.0/inst/lsqlin.m +***** test +***** shared C,d,A,b + C = [0.9501,0.7620,0.6153,0.4057;... + 0.2311,0.4564,0.7919,0.9354;... + 0.6068,0.0185,0.9218,0.9169;... + 0.4859,0.8214,0.7382,0.4102;... + 0.8912,0.4447,0.1762,0.8936]; + d = [0.0578; 0.3528; 0.8131; 0.0098; 0.1388]; + A =[0.2027, 0.2721, 0.7467, 0.4659;... + 0.1987, 0.1988, 0.4450, 0.4186;... + 0.6037 , 0.0152, 0.9318, 0.8462]; + b =[0.5251;0.2026;0.6721]; + Aeq = [3, 5, 7, 9]; + beq = 4; + lb = -0.1*ones(4,1); + ub = 2*ones(4,1); + [x,resnorm,residual,exitflag] = lsqlin(C,d,A,b,Aeq,beq,lb,ub); + assert(x,[-0.10000; -0.10000; 0.15991; 0.40896],10e-5) + assert(resnorm,0.16951,10e-5) + assert(residual, [0.035297; 0.087623; -0.353251; 0.145270; 0.121232],10e-5) + assert(exitflag,1) +warning: colon arguments should be scalars +warning: called from + null at line 67 column 14 + quadprog at line 302 column 13 + lsqlin at line 123 column 21 + __test__ at line 17 column 32 + test at line 677 column 11 + /tmp/tmp.i07XArg8KN at line 70 column 31 + +***** test + Aeq = []; + beq = []; + lb = []; + ub = []; + x0 = 0.1*ones(4,1); + x = lsqlin(C,d,A,b,Aeq,beq,lb,ub,x0); + [x,resnorm,residual,exitflag] = lsqlin(C,d,A,b,Aeq,beq,lb,ub,x0); + assert(x,[ 0.12986; -0.57569 ; 0.42510; 0.24384],10e-5) + assert(resnorm,0.017585,10e-5) + assert(residual, [-0.0126033; -0.0208040; -0.1295084; -0.0057389; 0.01372462],10e-5) + assert(exitflag,1) ***** demo - ## Example for default optimization (Levenberg/Marquardt with - ## BFGS), one non-linear equality constraint. Constrained optimum is - ## at p = [0; 1]. - objective_function = @ (p) p(1)^2 + p(2)^2; - pin = [-2; 5]; - constraint_function = @ (p) p(1)^2 + 1 - p(2); - [p, objf, cvg, outp] = nonlin_min (objective_function, pin, optimset ("equc", {constraint_function})) + C = [0.9501 0.7620 0.6153 0.4057 + 0.2311 0.4564 0.7919 0.9354 + 0.6068 0.0185 0.9218 0.9169 + 0.4859 0.8214 0.7382 0.4102 + 0.8912 0.4447 0.1762 0.8936]; + d = [0.0578; 0.3528; 0.8131; 0.0098; 0.1388]; + %% Linear Inequality Constraints + A =[0.2027 0.2721 0.7467 0.4659 + 0.1987 0.1988 0.4450 0.4186 + 0.6037 0.0152 0.9318 0.8462]; + b =[0.5251; 0.2026; 0.6721]; + %% Linear Equality Constraints + Aeq = [3 5 7 9]; + beq = 4; + %% Bound constraints + lb = -0.1*ones(4,1); + ub = ones(4,1); + [x, resnorm, residual, flag, output, lambda] = lsqlin (C, d, A, b, Aeq, beq, lb, ub) +2 tests, 2 passed, 0 known failure, 0 skipped +[inst/polyconf.m] +>>>>> /build/octave-optim-1.6.0/inst/polyconf.m +***** test + # data from Hocking, RR, "Methods and Applications of Linear Models" + temperature=[40;40;40;45;45;45;50;50;50;55;55;55;60;60;60;65;65;65]; + strength=[66.3;64.84;64.36;69.70;66.26;72.06;73.23;71.4;68.85;75.78;72.57;76.64;78.87;77.37;75.94;78.82;77.13;77.09]; + [p,s] = polyfit(temperature,strength,1); + [y,dy] = polyconf(p,40,s,0.05,'ci'); + assert([y,dy],[66.15396825396826,1.71702862681486],200*eps); + [y,dy] = polyconf(p,40,s,0.05,'pi'); + assert(dy,4.45345484470743,200*eps); +1 test, 1 passed, 0 known failure, 0 skipped +[inst/statget.m] +>>>>> /build/octave-optim-1.6.0/inst/statget.m +***** shared opts + opts = statset ("maxit", 100); +***** assert (statget (opts, "maxit"), 100) +***** assert (statget (opts, "MaxITer"), 100) +***** assert (statget (opts, "TolFun"), []) +***** assert (statget (opts, "TolFun", 1e-3), 1e-3) +***** error statget () +***** error statget (1) +***** error statget (1,2,3,4,5) +***** error statget (1, "name") +***** error statget (struct (), 2) +9 tests, 9 passed, 0 known failure, 0 skipped +[inst/nlinfit.m] +>>>>> /build/octave-optim-1.6.0/inst/nlinfit.m +***** test + modelfun = @(b, x) (b(1) + b(2) * exp (- b(3) * x)); + b = [1;3;2]; + v = rande ("state"); + rande ("seed", 1234); + xdata = exprnd (2,100,1); + rande ("state", v); + v = randn ("state"); + randn ("seed", 1234); + ydata = modelfun (b,xdata) + normrnd (0,0.1,100,1); + randn ("state", v); + beta0 = [2;2;3]; + beta = nlinfit(xdata,ydata,modelfun,beta0); + assert (beta, [1;3;2], 1e-1) ***** demo - ## Example for simulated annealing, two parameters, "trace_steps" - ## is true; - t_init = .2; - t_min = .002; - mu_t = 1.002; - iters_fixed_t = 10; - init_p = [2; 2]; - max_rand_step = [.2; .2]; - [p, objf, cvg, outp] = nonlin_min (@ (p) (p(1)/10)^2 + (p(2)/10)^2 + .1 * (-cos(4*p(1)) - cos(4*p(2))), init_p, optimset ("algorithm", "siman", "max_rand_step", max_rand_step, "t_init", t_init, "T_min", t_min, "mu_t", mu_t, "iters_fixed_T", iters_fixed_t, "trace_steps", true)); - p - objf - x = (outp.trace(:, 1) - 1) * iters_fixed_t + outp.trace(:, 2); - x(1) = 0; - plot (x, cat (2, outp.trace(:, 3:end), t_init ./ (mu_t .^ outp.trace(:, 1)))) - legend ({"objective function value", "p(1)", "p(2)", "Temperature"}) - xlabel ("subiteration") + modelfun = @(b, x) (b(1) + b(2) * exp (- b(3) * x)); + %% actual value + beta_without_noise = [1; 3; 2]; + x = [3.49622; 0.33751; 1.25675; 3.66981; 0.26237; 5.51095; ... + 2.11407; 1.48774; 6.22436; 2.04519]; + y_actual = modelfun (beta_without_noise, x); + noise = [0.176110; -0.066850; 0.231000; -0.047570; -0.108230; ... + 0.122790; 0.062940; 0.151510; 0.116010; -0.097460]; + y_noisy = y_actual + noise; + %% initial guess + beta0 = [2; 2; 2]; + %% weights vector + weights = [5; 16; 1; 20; 12; 11; 17; 8; 11; 13]; + [beta, R, J, covb, mse] = nlinfit (x, y_noisy, modelfun, beta0) + [beta_w, R_w, J_w, covb_w, mse_w] = nlinfit (x, y_noisy, modelfun, beta0, [], "weights", weights) +1 test, 1 passed, 0 known failure, 0 skipped +[inst/lsqcurvefit.m] +>>>>> /build/octave-optim-1.6.0/inst/lsqcurvefit.m +***** test + xdata = [0 .3 .8 1.1 1.6 2.3]; + ydata = [.82 .72 .63 .60 .55 .50]; + yhat = @(p,x) p(1) + p(2)*exp(-x); + opt = optimset('TolFun',1e-100); + [p, resnorm, residual] = lsqcurvefit(yhat,[1 1], xdata, ydata,[0 0],[],opt); + assert (p, [ 0.47595; 0.34132], 1e-5) + assert (resnorm, 3.2419e-004, 1e-8) + assert(residual, [-2.7283e-003, 8.8079e-003, -6.8307e-004, -1.0432e-002, -5.1366e-003, 1.0172e-002], 1e-5) +***** test + problem.solver = "lsqcurvefit"; + problem.objective = @(p,x) p(1) + p(2)*exp(-x); + problem.x0 = [1, 1]; + problem.xdata = [0 .3 .8 1.1 1.6 2.3]; + problem.ydata = [.82 .72 .63 .60 .55 .50]; + problem.lb = [0, 0]; + problem.ub = []; + problem.options = optimset('TolFun',1e-100); + [p, resnorm, residual] = lsqcurvefit (problem); + assert (p, [ 0.47595; 0.34132], 1e-5) + assert (resnorm, 3.2419e-004, 1e-8) + assert(residual, [-2.7283e-003, 8.8079e-003, -6.8307e-004, -1.0432e-002, -5.1366e-003, 1.0172e-002], 1e-5) +***** demo + %% Example for user specified Jacobian. + + %% independents + x = [1:10:100]'; + %% observed data + y =[9.2160e-001, 3.3170e-001, 8.9789e-002, 2.8480e-002, 2.6055e-002,... + 8.3641e-003, 4.2362e-003, 3.1693e-003, 1.4739e-004, 2.9406e-004]'; + %% initial values: + p0=[0.8; 0.05]; + %% bounds + lb=[0; 0]; ub=[]; + %% Jacobian setting + opts = optimset ("Jacobian", "on") + + %% model function: + function [F,J] = myfun (p, x) + F = p(1) * exp (-p(2) * x); + if nargout > 1 + J = [exp(- p(2) * x), - p(1) * x .* exp(- p(2) * x)]; + endif + endfunction + + [c, resnorm, residual, flag, output, lambda, jacob] = ... + lsqcurvefit (@ (varargin) myfun(varargin{:}), p0, x, y, lb, ub, opts) +2 tests, 2 passed, 0 known failure, 0 skipped +[inst/jacobs.m] +>>>>> /build/octave-optim-1.6.0/inst/jacobs.m +***** assert (jacobs (1, @(x) x), 1) +***** assert (jacobs (6, @(x) x^2), 12) +***** assert (jacobs ([1; 1], @(x) [x(1)^2; x(1)*x(2)]), [2, 0; 1, 1]) +***** assert (jacobs ([1; 2], @(x) [x(1)^2 + x(2); x(2)*exp(x(1))]), [2, 1; 2*exp(1), exp(1)]) +***** error jacobs () +***** error jacobs (1) +***** error jacobs (1, 2, 3, 4) +***** error jacobs (@sin, 1, [1, 1]) +***** error jacobs (@sin, 1, ones(2, 2)) +***** demo + # Relative error against several h-values + k = 3:20; h = 10 .^ (-k); x = 0.3*pi; + err = zeros (1, numel (k)); + for count = 1 : numel (k) + err(count) = abs (jacobs (x, @sin, struct ("h", h(count))) - cos (x)) / abs (cos (x)) + eps; + endfor + loglog (h, err); grid minor; + xlabel ("h"); ylabel ("|Df(x) - cos(x)| / |cos(x)|") + title ("f(x)=sin(x), f'(x)=cos(x) at x = 0.3pi") +9 tests, 9 passed, 0 known failure, 0 skipped +[inst/de_min.m] +>>>>> /build/octave-optim-1.6.0/inst/de_min.m +***** function result = f(x); + result = 100 * (x(2) - x(1)^2)^2 + (1 - x(1))^2; +***** test + tol = 1.0e-4; + ctl.tol = 0.0; + ctl.VTR = 1.0e-6; + ctl.XVmin = [-2 -2]; + ctl.XVmax = [ 2 2]; + rand("state", 11) + [x, obj_value, nfeval, convergence] = de_min (@f, ctl); + assert (convergence == 1); + assert (f(x) == obj_value); + assert (obj_value < ctl.VTR); +***** demo + ## define a simple example function + f = @(x) peaks(x(1), x(2)); + ## plot the function to see where the minimum might be + peaks() + ## first we set the region where we expect the minimum + ctl.XVmin = [-3 -3]; + ctl.XVmax = [ 3 3]; + ## and solve it with de_min + [x, obj_value, nfeval, convergence] = de_min (f, ctl) +1 test, 1 passed, 0 known failure, 0 skipped +[inst/nonlin_curvefit.m] +>>>>> /build/octave-optim-1.6.0/inst/nonlin_curvefit.m +***** demo + ## Example for linear inequality constraints + ## (see also the same example in 'demo nonlin_residmin') + + ## independents and observations + indep = 1:5; + obs = [1, 2, 4, 7, 14]; + ## model function: + f = @ (p, x) p(1) * exp (p(2) * x); + ## initial values: + init = [.25; .25]; + ## linear constraints, A.' * parametervector + B >= 0 + A = [1; -1]; B = 0; # p(1) >= p(2); + settings = optimset ("inequc", {A, B}); + + ## start optimization + [p, model_values, cvg, outp] = nonlin_curvefit (f, init, indep, obs, settings) +***** test + t = optim_problems ().curve.schittkowski_327; + cp = nonlin_curvefit (@ (p, x) t.f (x, p), t.init_p, t.data.x, t.data.y, optimset ("lbound", t.non_strict_inequc.bounds(:, 1), "ubound", t.non_strict_inequc.bounds(:, 2), "inequc", {t.non_strict_inequc.general}, "dfdp", @ (p, x) t.dfdp (x, p))); + assert (cp, t.result.p, [.0001; .001]); +1 test, 1 passed, 0 known failure, 0 skipped +[inst/nonlin_residmin.m] +>>>>> /build/octave-optim-1.6.0/inst/nonlin_residmin.m +***** demo + ## Example for linear inequality constraints + ## (see also the same example in 'demo nonlin_curvefit') + + ## independents + indep = 1:5; + ## residual function: + f = @ (p) p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]; + ## initial values: + init = [.25; .25]; + ## linear constraints, A.' * parametervector + B >= 0 + A = [1; -1]; B = 0; # p(1) >= p(2); + settings = optimset ("inequc", {A, B}); + + ## start optimization + [p, residuals, cvg, outp] = nonlin_residmin (f, init, settings) ***** test ## independents indep = 1:5; - ## objective function: - f = @ (p) sumsq (p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]); + ## residual function: + f = @ (p) p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]; ## initial values: init = [.25; .25]; ## linear constraints, A.' * parametervector + B >= 0 A = [1; -1]; B = 0; # p(1) >= p(2); settings = optimset ("inequc", {A, B}); - assert (nonlin_min (f, init, settings), [.6203; .6203], .0001); + assert (nonlin_residmin (f, init, settings), [.6203; .6203], .0001); +***** test + ## independents + indep = 1:5; + ## residual function: + f = @ (p) p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]; + ## initial values: + init = single ([.25; .25]); + ## linear constraints, A.' * parametervector + B >= 0 + A = [1; -1]; B = 0; # p(1) >= p(2); + settings = optimset ("inequc", {A, B}, + "complex_step_derivative_f", true); + + result = nonlin_residmin (f, init, settings); + assert (result, [.6203; .6203], .0001); + assert (isa (result, "single")); ***** test ***** shared x, misc x = [.871, .643, .550; @@ -3677,20 +3664,29 @@ pconf.c.ubound = .3; pconf.d.fixed = true; - settings = optimset ("param_config", pconf, "objf_pstruct", true); + settings = optimset ("param_config", pconf, "f_pstruct", true); - f = @ (p) sumsq (( ... + f = @ (p) ... subsasgn (x, struct ("type", "()", "subs", {{9, 3}}), p.c) ... * horzcat (p.a, p.a([3, 1, 2]), p.a([3, 2, 1])) ... - - p.d ... - * subsasgn (misc, struct ("type", "()", "subs", {{3, 2}}), p.b))(:)); + - p.d * subsasgn (misc, struct ("type", "()", "subs", {{3, 2}}), p.b); - [p, ~, ~, outp] = nonlin_min (f, pin, settings); + [p, ~, ~, outp] = nonlin_residmin (f, pin, settings); - assert (p.a, [1.0590; 1.9266; 4.0456], .005); - assert (p.b, 2.7061, .005); + assert (p.a, [1.0590; 1.9266; 4.0456], .0001); + assert (p.b, 2.7061, .0001); assert (p.c, .3, .000001); assert (p.d, 1); + assert (isempty (outp.lambda.ineqlin)); + assert (isempty (outp.lambda.eqlin)); + assert (isempty (outp.lambda.ineqnonlin)); + assert (isempty (outp.lambda.eqnonlin)); + assert (! any (outp.lambda.lower.a)); + assert (! outp.lambda.lower.b); + assert (! outp.lambda.lower.c); + assert (! any (outp.lambda.upper.a)); + assert (! outp.lambda.upper.b); + assert (outp.lambda.upper.c > 0); ***** test pin = zeros (6, 1); pin(6) = 1; @@ -3700,61 +3696,124 @@ "diff_onesided", true, "fixed", [false; false; false; false; false; true]); - f = @ (p) sumsq (( ... + f = @ (p) ... subsasgn (x, struct ("type", "()", "subs", {{9, 3}}), p(5)) ... * horzcat (p([1, 2, 3]), p([3, 1, 2]), p([3, 2, 1])) ... - - p(6) ... - * subsasgn (misc, struct ("type", "()", "subs", {{3, 2}}), p(4)))(:)); + - p(6) * subsasgn (misc, struct ("type", "()", "subs", {{3, 2}}), p(4)); - p = nonlin_min (f, pin, settings); + p = nonlin_residmin (f, pin, settings); - assert (p, [1.0590; 1.9266; 4.0456; 2.7061; .3; 1], .005); -3 tests, 3 passed, 0 known failure, 0 skipped -[inst/nonlin_residmin.m] ->>>>> /build/octave-optim-1.6.0/inst/nonlin_residmin.m + assert (p, [1.0590; 1.9266; 4.0456; 2.7061; .3; 1], .0001); +4 tests, 4 passed, 0 known failure, 0 skipped +[inst/fmincon.m] +>>>>> /build/octave-optim-1.6.0/inst/fmincon.m ***** demo - ## Example for linear inequality constraints - ## (see also the same example in 'demo nonlin_curvefit') - - ## independents - indep = 1:5; - ## residual function: - f = @ (p) p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]; - ## initial values: - init = [.25; .25]; - ## linear constraints, A.' * parametervector + B >= 0 - A = [1; -1]; B = 0; # p(1) >= p(2); - settings = optimset ("inequc", {A, B}); - - ## start optimization - [p, residuals, cvg, outp] = nonlin_residmin (f, init, settings) + ## Example for default optimization (Levenberg/Marquardt with + ## BFGS), one non-linear equality constraint. Constrained optimum is + ## at p = [0; 1]. + objective_function = @ (p) p(1)^2 + p(2)^2; + pin = [-2; 5]; + constraint_function = @ (p) - (p(1)^2 + 1 - p(2)); + [p, objf, cvg, outp] = fmincon (objective_function, pin, [], [], [], [], [], [], @ (p) {[], constraint_function(p)}{:}) +***** test + ## equality constraint + objective_function = @ (p) p(1)^2 + p(2)^2; + pin = [-2; 5]; + constraint_function = @ (p) - (p(1)^2 + 1 - p(2)); + [p, objf, cvg, outp] = fmincon (objective_function, pin, [], [], [], [], [], [], @ (p) {[], constraint_function(p)}{:}, optimset ("Algorithm", "lm_feasible")); + assert (p, [0; 1], 1e-6) +***** test + ## inequality constraint + objective_function = @ (p) p(1)^2 + p(2)^2; + pin = [2; 6]; + constraint_function = @ (p) p(1)^2 + 1 - p(2); + [p, objf, cvg, outp] = fmincon (objective_function, pin, [], [], [], [], [], [], @ (p) {constraint_function(p), []}{:}, optimset ("Algorithm", "lm_feasible")); + assert (p, [0; 1], 1e-6) ***** test ## independents indep = 1:5; - ## residual function: - f = @ (p) p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]; + ## objective function: + f = @ (p) sumsq (p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]); ## initial values: init = [.25; .25]; ## linear constraints, A.' * parametervector + B >= 0 A = [1; -1]; B = 0; # p(1) >= p(2); - settings = optimset ("inequc", {A, B}); - assert (nonlin_residmin (f, init, settings), [.6203; .6203], .0001); + assert (fmincon (f, init, -A, B), [.6203; .6203], .0001); +***** test + ## problem structure + indep = 1:5; + problem = struct ("objective", + @ (p) sumsq (p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]), + "x0", [.25; .25], + "Aineq", [-1; 1], + "bineq", 0, + "solver", "fmincon"); + assert (fmincon (problem), [.6203; .6203], .0001); +***** test + ## Octave sqp solver with a lot of inequality constraints + objf = @ (p) sumsq (p(4:9)); + init = [300; -100; -.1997; -127; -151; 379; 421; 460; 426]; + lbound = [-Inf; -Inf; -Inf; 0; 0; 0; 0; 0; 0]; + inequc = @ (p) vertcat ( ... + p(1) + p(2) * exp (-5 * p(3)) + p(4) - 127, + p(1) + p(2) * exp (-3 * p(3)) + p(5) - 151, + p(1) + p(2) * exp (-p(3)) + p(6) - 379, + p(1) + p(2) * exp (p(3)) + p(7) - 421, + p(1) + p(2) * exp (3 * p(3)) + p(8) - 460, + p(1) + p(2) * exp (5 * p(3)) + p(9) - 426, + -p(1) - p(2) * exp (-5 * p(3)) + p(4) + 127, + -p(1) - p(2) * exp (-3 * p(3)) + p(5) + 151, + -p(1) - p(2) * exp (-p(3)) + p(6) + 379, + -p(1) - p(2) * exp (p(3)) + p(7) + 421, + -p(1) - p(2) * exp (3 * p(3)) + p(8) + 460, + -p(1) - p(2) * exp (5 * p(3)) + p(9) + 426); + [p, objf, cvg, outp] = fmincon (objf, init, [], [], [], [], lbound, [], + @ (p) {- inequc(p), []}{:}, + optimset ("Algorithm", "octave_sqp")); + assert (p, [5.2330e+02; -1.5694e+02; -1.9966e-01; 2.9607e+01; + 8.6615e+01; 4.7326e+01; 2.6235e+01; 2.2915e+01; + 3.9470e+01], .01); +5 tests, 5 passed, 0 known failure, 0 skipped +[inst/nonlin_min.m] +>>>>> /build/octave-optim-1.6.0/inst/nonlin_min.m +***** demo + ## Example for default optimization (Levenberg/Marquardt with + ## BFGS), one non-linear equality constraint. Constrained optimum is + ## at p = [0; 1]. + objective_function = @ (p) p(1)^2 + p(2)^2; + pin = [-2; 5]; + constraint_function = @ (p) p(1)^2 + 1 - p(2); + [p, objf, cvg, outp] = nonlin_min (objective_function, pin, optimset ("equc", {constraint_function})) +***** demo + ## Example for simulated annealing, two parameters, "trace_steps" + ## is true; + t_init = .2; + t_min = .002; + mu_t = 1.002; + iters_fixed_t = 10; + init_p = [2; 2]; + max_rand_step = [.2; .2]; + [p, objf, cvg, outp] = nonlin_min (@ (p) (p(1)/10)^2 + (p(2)/10)^2 + .1 * (-cos(4*p(1)) - cos(4*p(2))), init_p, optimset ("algorithm", "siman", "max_rand_step", max_rand_step, "t_init", t_init, "T_min", t_min, "mu_t", mu_t, "iters_fixed_T", iters_fixed_t, "trace_steps", true)); + p + objf + x = (outp.trace(:, 1) - 1) * iters_fixed_t + outp.trace(:, 2); + x(1) = 0; + plot (x, cat (2, outp.trace(:, 3:end), t_init ./ (mu_t .^ outp.trace(:, 1)))) + legend ({"objective function value", "p(1)", "p(2)", "Temperature"}) + xlabel ("subiteration") ***** test ## independents indep = 1:5; - ## residual function: - f = @ (p) p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]; + ## objective function: + f = @ (p) sumsq (p(1) * exp (p(2) * indep) - [1, 2, 4, 7, 14]); ## initial values: - init = single ([.25; .25]); + init = [.25; .25]; ## linear constraints, A.' * parametervector + B >= 0 A = [1; -1]; B = 0; # p(1) >= p(2); - settings = optimset ("inequc", {A, B}, - "complex_step_derivative_f", true); + settings = optimset ("inequc", {A, B}); - result = nonlin_residmin (f, init, settings); - assert (result, [.6203; .6203], .0001); - assert (isa (result, "single")); + assert (nonlin_min (f, init, settings), [.6203; .6203], .0001); ***** test ***** shared x, misc x = [.871, .643, .550; @@ -3807,29 +3866,20 @@ pconf.c.ubound = .3; pconf.d.fixed = true; - settings = optimset ("param_config", pconf, "f_pstruct", true); + settings = optimset ("param_config", pconf, "objf_pstruct", true); - f = @ (p) ... + f = @ (p) sumsq (( ... subsasgn (x, struct ("type", "()", "subs", {{9, 3}}), p.c) ... * horzcat (p.a, p.a([3, 1, 2]), p.a([3, 2, 1])) ... - - p.d * subsasgn (misc, struct ("type", "()", "subs", {{3, 2}}), p.b); + - p.d ... + * subsasgn (misc, struct ("type", "()", "subs", {{3, 2}}), p.b))(:)); - [p, ~, ~, outp] = nonlin_residmin (f, pin, settings); + [p, ~, ~, outp] = nonlin_min (f, pin, settings); - assert (p.a, [1.0590; 1.9266; 4.0456], .0001); - assert (p.b, 2.7061, .0001); + assert (p.a, [1.0590; 1.9266; 4.0456], .005); + assert (p.b, 2.7061, .005); assert (p.c, .3, .000001); assert (p.d, 1); - assert (isempty (outp.lambda.ineqlin)); - assert (isempty (outp.lambda.eqlin)); - assert (isempty (outp.lambda.ineqnonlin)); - assert (isempty (outp.lambda.eqnonlin)); - assert (! any (outp.lambda.lower.a)); - assert (! outp.lambda.lower.b); - assert (! outp.lambda.lower.c); - assert (! any (outp.lambda.upper.a)); - assert (! outp.lambda.upper.b); - assert (outp.lambda.upper.c > 0); ***** test pin = zeros (6, 1); pin(6) = 1; @@ -3839,15 +3889,16 @@ "diff_onesided", true, "fixed", [false; false; false; false; false; true]); - f = @ (p) ... + f = @ (p) sumsq (( ... subsasgn (x, struct ("type", "()", "subs", {{9, 3}}), p(5)) ... * horzcat (p([1, 2, 3]), p([3, 1, 2]), p([3, 2, 1])) ... - - p(6) * subsasgn (misc, struct ("type", "()", "subs", {{3, 2}}), p(4)); + - p(6) ... + * subsasgn (misc, struct ("type", "()", "subs", {{3, 2}}), p(4)))(:)); - p = nonlin_residmin (f, pin, settings); + p = nonlin_min (f, pin, settings); - assert (p, [1.0590; 1.9266; 4.0456; 2.7061; .3; 1], .0001); -4 tests, 4 passed, 0 known failure, 0 skipped + assert (p, [1.0590; 1.9266; 4.0456; 2.7061; .3; 1], .005); +3 tests, 3 passed, 0 known failure, 0 skipped Checking C++ files ... Summary: 56 tests, 56 passed, 0 known failures, 0 skipped dh_installdocs -O--buildsystem=octave @@ -3865,7 +3916,7 @@ dh_strip -a -O--buildsystem=octave dh_makeshlibs -a -O--buildsystem=octave dh_shlibdeps -a -O--buildsystem=octave -dpkg-shlibdeps: warning: package could avoid a useless dependency if debian/octave-optim/usr/lib/i386-linux-gnu/octave/packages/optim-1.6.0/i686-pc-linux-gnu-api-v55/__bfgsmin.oct debian/octave-optim/usr/lib/i386-linux-gnu/octave/packages/optim-1.6.0/i686-pc-linux-gnu-api-v55/__disna_optim__.oct debian/octave-optim/usr/lib/i386-linux-gnu/octave/packages/optim-1.6.0/i686-pc-linux-gnu-api-v55/numgradient.oct debian/octave-optim/usr/lib/i386-linux-gnu/octave/packages/optim-1.6.0/i686-pc-linux-gnu-api-v55/numhessian.oct were not linked against libpthread.so.0 (they use none of the library's symbols) +dpkg-shlibdeps: warning: package could avoid a useless dependency if debian/octave-optim/usr/lib/i386-linux-gnu/octave/packages/optim-1.6.0/i686-pc-linux-gnu-api-v55/__bfgsmin.oct debian/octave-optim/usr/lib/i386-linux-gnu/octave/packages/optim-1.6.0/i686-pc-linux-gnu-api-v55/numhessian.oct debian/octave-optim/usr/lib/i386-linux-gnu/octave/packages/optim-1.6.0/i686-pc-linux-gnu-api-v55/__disna_optim__.oct debian/octave-optim/usr/lib/i386-linux-gnu/octave/packages/optim-1.6.0/i686-pc-linux-gnu-api-v55/numgradient.oct were not linked against libpthread.so.0 (they use none of the library's symbols) dh_octave_substvar -O--buildsystem=octave dh_installdeb -O--buildsystem=octave dh_gencontrol -O--buildsystem=octave @@ -3873,8 +3924,8 @@ dpkg-gencontrol: warning: package octave-optim: substitution variable ${octave:Upstream-Description} unused, but is defined dh_md5sums -O--buildsystem=octave dh_builddeb -O--buildsystem=octave -dpkg-deb: building package 'octave-optim-dbgsym' in '../octave-optim-dbgsym_1.6.0-6_i386.deb'. dpkg-deb: building package 'octave-optim' in '../octave-optim_1.6.0-6_i386.deb'. +dpkg-deb: building package 'octave-optim-dbgsym' in '../octave-optim-dbgsym_1.6.0-6_i386.deb'. dpkg-genbuildinfo --build=binary dpkg-genchanges --build=binary >../octave-optim_1.6.0-6_i386.changes dpkg-genchanges: info: binary-only upload (no source code included) @@ -3882,12 +3933,14 @@ dpkg-source: info: using options from octave-optim-1.6.0/debian/source/options: --extend-diff-ignore=doc/optim\.info|doc/optim\.texi dpkg-buildpackage: info: binary-only upload (no source included) I: copying local configuration +I: user script /srv/workspace/pbuilder/55977/tmp/hooks/B01_cleanup starting +I: user script /srv/workspace/pbuilder/55977/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/30876 and its subdirectories -I: Current time: Thu Jul 8 23:16:57 -12 2021 -I: pbuilder-time-stamp: 1625829417 +I: removing directory /srv/workspace/pbuilder/55977 and its subdirectories +I: Current time: Fri Aug 12 07:42:43 +14 2022 +I: pbuilder-time-stamp: 1660239763