>>> mpb: Building testing/mpb 1.12.0-r0 (using abuild 3.18.0_rc1-r0) started Tue, 30 Jun 2026 03:59:26 +0000 >>> mpb: Validating /home/buildozer/aports/testing/mpb/APKBUILD... >>> mpb: Analyzing dependencies... >>> mpb: Installing for build: build-base fftw-dev gfortran guile-dev hdf5-dev libctl-dev openblas-dev ( 1/31) Installing fftw-double-libs (3.3.11-r0) ( 2/31) Installing fftw-long-double-libs (3.3.11-r0) ( 3/31) Installing fftw-single-libs (3.3.11-r0) ( 4/31) Installing fftw-dev (3.3.11-r0) ( 5/31) Installing libgfortran (15.2.0-r6) ( 6/31) Installing gfortran (15.2.0-r6) ( 7/31) Installing libffi (3.5.2-r1) ( 8/31) Installing gc (8.2.12-r0) ( 9/31) Installing guile-libs (3.0.9-r2) (10/31) Installing guile (3.0.9-r2) (11/31) Installing libgc++ (8.2.12-r0) (12/31) Installing gc-dev (8.2.12-r0) (13/31) Installing libgmpxx (6.3.0-r4) (14/31) Installing gmp-dev (6.3.0-r4) (15/31) Installing guile-dev (3.0.9-r2) (16/31) Installing zlib-dev (1.3.2-r0) (17/31) Installing libsz (1.1.6-r0) (18/31) Installing hdf5 (1.14.4.2-r1) (19/31) Installing hdf5-cpp (1.14.4.2-r1) (20/31) Installing hdf5-fortran (1.14.4.2-r1) (21/31) Installing hdf5-hl (1.14.4.2-r1) (22/31) Installing hdf5-hl-cpp (1.14.4.2-r1) (23/31) Installing hdf5-hl-fortran (1.14.4.2-r1) (24/31) Installing hdf5-dev (1.14.4.2-r1) (25/31) Installing libctl (4.5.1-r1) (26/31) Installing libctl-dev (4.5.1-r1) (27/31) Installing openblas (0.3.30-r2) (28/31) Installing liblapack (0.3.30-r2) (29/31) Installing liblapacke (0.3.30-r2) (30/31) Installing openblas-dev (0.3.30-r2) (31/31) Installing .makedepends-mpb (20260630.035927) Executing busybox-1.38.0-r1.trigger OK: 516.2 MiB in 138 packages >>> mpb: Cleaning up srcdir >>> mpb: Cleaning up pkgdir >>> mpb: Cleaning up tmpdir >>> mpb: Fetching https://distfiles.alpinelinux.org/distfiles/edge//mpb-1.12.0.tar.gz Connecting to distfiles.alpinelinux.org (172.105.82.32:443) wget: server returned error: HTTP/1.1 404 Not Found >>> mpb: Fetching https://github.com/NanoComp/mpb/releases/download/v1.12.0/mpb-1.12.0.tar.gz Connecting to github.com (20.26.156.215:443) Connecting to release-assets.githubusercontent.com (185.199.111.133:443) saving to '/var/cache/distfiles/edge/mpb-1.12.0.tar.gz.part' mpb-1.12.0.tar.gz.pa 100% |********************************| 651k 0:00:00 ETA '/var/cache/distfiles/edge/mpb-1.12.0.tar.gz.part' saved /var/cache/distfiles/edge/mpb-1.12.0.tar.gz: OK >>> mpb: Fetching https://distfiles.alpinelinux.org/distfiles/edge//mpb-1.12.0.tar.gz /var/cache/distfiles/edge/mpb-1.12.0.tar.gz: OK >>> mpb: Unpacking /var/cache/distfiles/edge/mpb-1.12.0.tar.gz... checking whether to enable maintainer-specific portions of Makefiles... no checking for a BSD-compatible install... /usr/bin/install -c checking whether sleep supports fractional seconds... yes checking filesystem timestamp resolution... 2 checking whether build environment is sane... yes checking for a race-free mkdir -p... /bin/mkdir -p checking for gawk... no checking for mawk... no checking for nawk... no checking for awk... awk checking whether make sets $(MAKE)... yes checking whether make supports nested variables... yes checking xargs -n works... yes checking build system type... s390x-alpine-linux-musl checking host system type... s390x-alpine-linux-musl checking how to print strings... printf checking whether make supports the include directive... yes (GNU style) checking for s390x-alpine-linux-musl-gcc... cc checking whether the C compiler works... yes checking for C compiler default output file name... a.out checking for suffix of executables... checking whether we are cross compiling... no checking for suffix of object files... o checking whether the compiler supports GNU C... yes checking whether cc accepts -g... yes checking for cc option to enable C11 features... none needed checking whether cc understands -c and -o together... yes checking dependency style of cc... gcc3 checking for a sed that does not truncate output... /bin/sed checking for grep that handles long lines and -e... /bin/grep checking for egrep... /bin/grep -E checking for fgrep... /bin/grep -F checking for ld used by cc... /usr/s390x-alpine-linux-musl/bin/ld checking if the linker (/usr/s390x-alpine-linux-musl/bin/ld) is GNU ld... yes checking for BSD- or MS-compatible name lister (nm)... /usr/bin/nm -B checking the name lister (/usr/bin/nm -B) interface... BSD nm checking whether ln -s works... yes checking the maximum length of command line arguments... 98304 checking how to convert s390x-alpine-linux-musl file names to s390x-alpine-linux-musl format... func_convert_file_noop checking how to convert s390x-alpine-linux-musl file names to toolchain format... func_convert_file_noop checking for /usr/s390x-alpine-linux-musl/bin/ld option to reload object files... -r checking for file... file checking for s390x-alpine-linux-musl-objdump... no checking for objdump... objdump checking how to recognize dependent libraries... pass_all checking for s390x-alpine-linux-musl-dlltool... no checking for dlltool... no checking how to associate runtime and link libraries... printf %s\n checking for s390x-alpine-linux-musl-ranlib... no checking for ranlib... ranlib checking for s390x-alpine-linux-musl-ar... no checking for ar... ar checking for archiver @FILE support... @ checking for s390x-alpine-linux-musl-strip... no checking for strip... strip checking command to parse /usr/bin/nm -B output from cc object... ok checking for sysroot... no checking for a working dd... /bin/dd checking how to truncate binary pipes... /bin/dd bs=4096 count=1 checking for s390x-alpine-linux-musl-mt... no checking for mt... no checking if : is a manifest tool... no checking for stdio.h... yes checking for stdlib.h... yes checking for string.h... yes checking for inttypes.h... yes checking for stdint.h... yes checking for strings.h... yes checking for sys/stat.h... yes checking for sys/types.h... yes checking for unistd.h... yes checking for dlfcn.h... yes checking for objdir... .libs checking if cc supports -fno-rtti -fno-exceptions... no checking for cc option to produce PIC... -fPIC -DPIC checking if cc PIC flag -fPIC -DPIC works... yes checking if cc static flag -static works... yes checking if cc supports -c -o file.o... yes checking if cc supports -c -o file.o... (cached) yes checking whether the cc linker (/usr/s390x-alpine-linux-musl/bin/ld -m elf64_s390) supports shared libraries... yes checking whether -lc should be explicitly linked in... no checking dynamic linker characteristics... GNU/Linux ld.so checking how to hardcode library paths into programs... immediate checking whether stripping libraries is possible... yes checking if libtool supports shared libraries... yes checking whether to build shared libraries... yes checking whether to build static libraries... yes checking for s390x-alpine-linux-musl-g77... no checking for s390x-alpine-linux-musl-xlf... no checking for s390x-alpine-linux-musl-f77... no checking for s390x-alpine-linux-musl-frt... no checking for s390x-alpine-linux-musl-pgf77... no checking for s390x-alpine-linux-musl-cf77... no checking for s390x-alpine-linux-musl-fort77... no checking for s390x-alpine-linux-musl-fl32... no checking for s390x-alpine-linux-musl-af77... no checking for s390x-alpine-linux-musl-xlf90... no checking for s390x-alpine-linux-musl-f90... no checking for s390x-alpine-linux-musl-pgf90... no checking for s390x-alpine-linux-musl-pghpf... no checking for s390x-alpine-linux-musl-epcf90... no checking for s390x-alpine-linux-musl-gfortran... s390x-alpine-linux-musl-gfortran checking whether the compiler supports GNU Fortran 77... yes checking whether s390x-alpine-linux-musl-gfortran accepts -g... yes checking if libtool supports shared libraries... yes checking whether to build shared libraries... yes checking whether to build static libraries... yes checking for s390x-alpine-linux-musl-gfortran option to produce PIC... -fPIC checking if s390x-alpine-linux-musl-gfortran PIC flag -fPIC works... yes checking if s390x-alpine-linux-musl-gfortran static flag -static works... yes checking if s390x-alpine-linux-musl-gfortran supports -c -o file.o... yes checking if s390x-alpine-linux-musl-gfortran supports -c -o file.o... (cached) yes checking whether the s390x-alpine-linux-musl-gfortran linker (/usr/s390x-alpine-linux-musl/bin/ld -m elf64_s390) supports shared libraries... yes checking dynamic linker characteristics... (cached) GNU/Linux ld.so checking how to hardcode library paths into programs... immediate checking for perl... ./missing checking for vendor cc to be used instead of gcc... checking for cc... (cached) cc checking for s390x-alpine-linux-musl-gcc... (cached) cc checking whether the compiler supports GNU C... (cached) yes checking whether cc accepts -g... (cached) yes checking for cc option to enable C11 features... (cached) none needed checking whether cc understands -c and -o together... (cached) yes checking dependency style of cc... (cached) gcc3 checking for sqrt in -lm... yes checking how to get verbose linking output from s390x-alpine-linux-musl-gfortran... -v checking for Fortran 77 libraries of s390x-alpine-linux-musl-gfortran... -L/usr/lib/gcc/s390x-alpine-linux-musl/15.2.0 -L/usr/lib/gcc/s390x-alpine-linux-musl/15.2.0/../../../../s390x-alpine-linux-musl/lib/../lib -L/usr/lib/gcc/s390x-alpine-linux-musl/15.2.0/../../../../lib -L/lib/../lib -L/usr/lib/../lib -L/usr/lib/gcc/s390x-alpine-linux-musl/15.2.0/../../../../s390x-alpine-linux-musl/lib -L/usr/lib/gcc/s390x-alpine-linux-musl/15.2.0/../../.. -L/lib -L/usr/lib -lgfortran -lm -lssp_nonshared checking for dummy main to link with Fortran 77 libraries... none checking for Fortran 77 name-mangling scheme... lower case, underscore, no extra underscore checking for fftw_execute in -lfftw3... yes checking for sgemm_ in -lopenblas... yes checking for cheev_ in -lopenblas... yes checking for deflate in -lz... yes checking for H5Pcreate in -lhdf5... yes checking for hdf5.h... yes checking for guile... yes checking for guile-config... guile-config checking if linking to guile works... yes checking for libguile.h... yes checking for guile/gh.h... no checking for scm_make_smob_type... yes checking for scm_new_smob... yes checking for scm_array_get_handle... yes checking for scm_is_array... yes checking for SCM_SMOB_PREDICATE... no checking for SCM_SMOB_DATA... no checking for SCM_NEWSMOB... no checking how to activate readline in Guile... cannot checking for libctl dir... /usr/share/libctl checking for gen-ctl-io... gen-ctl-io checking for ctl_get_vector3 in -lctl... yes checking for ctl.h... yes checking whether libctl version is at least 4.0.0... checking how to run the C preprocessor... cc -E checking for egrep -e... (cached) /bin/grep -E ok checking for libctl_quiet feature... yes checking for vector3_dot in -lctlgeom... yes checking for ctlgeom.h... yes checking for nlopt.h... no checking for nlopt_minimize in -lnlopt... no checking for unistd.h... (cached) yes checking for getopt.h... yes checking for nlopt.h... (cached) no checking for an ANSI C-conforming const... yes checking for inline... inline checking for getopt... yes checking for strncmp... yes checking whether calling BLAS zdotc works... yes checking for bash... sh configure: WARNING: Couldn't find bash or ksh shell on your system. The mpb-split command may not work if your sh doesn't support job control. checking for feenableexcept... no checking whether feenableexcept declaration is usable... no checking that generated files are newer than configure... done configure: creating ./config.status config.status: creating Makefile config.status: creating examples/Makefile config.status: creating tests/Makefile config.status: creating src/Makefile config.status: creating src/util/Makefile config.status: creating src/matrices/Makefile config.status: creating src/matrixio/Makefile config.status: creating src/maxwell/Makefile config.status: creating mpb/Makefile config.status: creating mpb/mpb.scm config.status: creating mpb/mpb-split-preinstall config.status: creating utils/Makefile config.status: creating config.h config.status: creating src/mpbconf.h config.status: executing depfiles commands config.status: executing libtool commands make all-recursive make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0' Making all in src make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' cp -f mpbconf.h mpb.h make all-recursive make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' Making all in util make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' make all-am make[5]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' CC sphere-quad.o CC debug_malloc.lo CC mpi_utils.lo CC verbosity.lo mpi_utils.c: In function 'mpi_begin_critical_section': mpi_utils.c:202:22: warning: unused variable 'status' [-Wunused-variable] 202 | MPI_Status status; | ^~~~~~ mpi_utils.c: In function 'mpi_end_critical_section': mpi_utils.c:210:35: warning: unused parameter 'tag' [-Wunused-parameter] 210 | void mpi_end_critical_section(int tag) | ~~~~^~~ CCLD libutil.la CCLD sphere_quad make[5]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' Making all in matrices make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrices' CC libmatrices_la-blasglue.lo CC libmatrices_la-eigensolver.lo CC libmatrices_la-eigensolver_davidson.lo CC libmatrices_la-eigensolver_utils.lo CC libmatrices_la-evectmatrix.lo CC libmatrices_la-linmin.lo CC libmatrices_la-matrices.lo CC libmatrices_la-minpack2-linmin.lo CC libmatrices_la-sqmatrix.lo CCLD libmatrices.la make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrices' Making all in matrixio make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrixio' CC libmatrixio_a-evectmatrixio.o CC libmatrixio_a-fieldio.o CC libmatrixio_a-matrixio.o fieldio.c:41:46: warning: argument 8 of type 'const real *' {aka 'const double *'} declared as a pointer [-Warray-parameter=] 41 | const real *kvector, | ~~~~~~~~~~~~^~~~~~~ In file included from fieldio.c:28: matrixio.h:95:52: note: previously declared as an array 'const real[3]' {aka 'const double[3]'} 95 | const real kvector[3], | ~~~~~~~~~~~^~~~~~~~~~ matrixio.c: In function 'matrixio_create_sub': matrixio.c:469:24: warning: statement with no effect [-Wunused-value] 469 | IF_EXCLUSIVE(0,if (id.parallel) MPI_Barrier(mpb_comm)); | ^ matrixio.c:58:32: note: in definition of macro 'IF_EXCLUSIVE' 58 | # define IF_EXCLUSIVE(yes,no) yes | ^~~ matrixio.c:472:24: warning: statement with no effect [-Wunused-value] 472 | IF_EXCLUSIVE(0,if (id.parallel) MPI_Barrier(mpb_comm)); | ^ matrixio.c:58:32: note: in definition of macro 'IF_EXCLUSIVE' 58 | # define IF_EXCLUSIVE(yes,no) yes | ^~~ In file included from matrixio.c:29: matrixio.c: In function 'matrixio_open_dataset': matrixio.c:516:30: warning: comparison of integer expressions of different signedness: 'hsize_t' {aka 'long unsigned int'} and 'int' [-Wsign-compare] 516 | CHECK(dims_copy[i] == dims[i], | ^~ ./../util/check.h:37:12: note: in definition of macro 'CHECK' 37 | if (!(condition)) { \ | ^~~~~~~~~ matrixio.c: In function 'matrixio_create_dataset': matrixio.c:552:24: warning: statement with no effect [-Wunused-value] 552 | IF_EXCLUSIVE(0,if (id.parallel) MPI_Barrier(mpb_comm)); | ^ matrixio.c:58:32: note: in definition of macro 'IF_EXCLUSIVE' 58 | # define IF_EXCLUSIVE(yes,no) yes | ^~~ matrixio.c: In function 'matrixio_read_real_data': matrixio.c:837:35: warning: comparison of integer expressions of different signedness: 'hsize_t' {aka 'long unsigned int'} and 'int' [-Wsign-compare] 837 | CHECK(dims_copy[i] == dims[i], | ^~ ./../util/check.h:37:12: note: in definition of macro 'CHECK' 37 | if (!(condition)) { \ | ^~~~~~~~~ AR libmatrixio.a make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrixio' Making all in maxwell make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/maxwell' CC libmaxwell_la-maxwell.lo CC libmaxwell_la-maxwell_constraints.lo CC libmaxwell_la-maxwell_eps.lo CC libmaxwell_la-maxwell_op.lo CC libmaxwell_la-maxwell_pre.lo maxwell_op.c: In function 'maxwell_compute_fft': maxwell_op.c:187:21: warning: variable 'np' set but not used [-Wunused-but-set-variable] 187 | ptrdiff_t np[3]; | ^~ maxwell_op.c:177:12: warning: unused variable 'rarray_out' [-Wunused-variable] 177 | real *rarray_out = (real *) array_out; | ^~~~~~~~~~ maxwell_op.c:175:12: warning: unused variable 'rarray_in' [-Wunused-variable] 175 | real *rarray_in = (real *) array_in; | ^~~~~~~~~ maxwell_eps.c:181:12: warning: 'sym_matrix_eq' defined but not used [-Wunused-function] 181 | static int sym_matrix_eq(symmetric_matrix V1, symmetric_matrix V2, double tol) | ^~~~~~~~~~~~~ maxwell_op.c: In function 'maxwell_vectorfield_otherhalf': maxwell_op.c:637:50: warning: unused parameter 'd' [-Wunused-parameter] 637 | void maxwell_vectorfield_otherhalf(maxwell_data *d, scalar_complex *field, | ~~~~~~~~~~~~~~^ maxwell_op.c:637:69: warning: unused parameter 'field' [-Wunused-parameter] 637 | void maxwell_vectorfield_otherhalf(maxwell_data *d, scalar_complex *field, | ~~~~~~~~~~~~~~~~^~~~~ maxwell_op.c:638:41: warning: unused parameter 'phasex' [-Wunused-parameter] 638 | real phasex, real phasey, real phasez) | ~~~~~^~~~~~ maxwell_op.c:638:54: warning: unused parameter 'phasey' [-Wunused-parameter] 638 | real phasex, real phasey, real phasez) | ~~~~~^~~~~~ maxwell_op.c:638:67: warning: unused parameter 'phasez' [-Wunused-parameter] 638 | real phasex, real phasey, real phasez) | ~~~~~^~~~~~ maxwell_op.c: In function 'maxwell_cscalarfield_otherhalf': maxwell_op.c:841:51: warning: unused parameter 'd' [-Wunused-parameter] 841 | void maxwell_cscalarfield_otherhalf(maxwell_data *d, scalar_complex *field, | ~~~~~~~~~~~~~~^ maxwell_op.c:841:70: warning: unused parameter 'field' [-Wunused-parameter] 841 | void maxwell_cscalarfield_otherhalf(maxwell_data *d, scalar_complex *field, | ~~~~~~~~~~~~~~~~^~~~~ maxwell_op.c:842:42: warning: unused parameter 'phasex' [-Wunused-parameter] 842 | real phasex, real phasey, real phasez) | ~~~~~^~~~~~ maxwell_op.c:842:55: warning: unused parameter 'phasey' [-Wunused-parameter] 842 | real phasex, real phasey, real phasez) | ~~~~~^~~~~~ maxwell_op.c:842:68: warning: unused parameter 'phasez' [-Wunused-parameter] 842 | real phasex, real phasey, real phasez) | ~~~~~^~~~~~ maxwell_op.c: In function 'maxwell_scalarfield_otherhalf': maxwell_op.c:1042:50: warning: unused parameter 'd' [-Wunused-parameter] 1042 | void maxwell_scalarfield_otherhalf(maxwell_data *d, real *field) | ~~~~~~~~~~~~~~^ maxwell_op.c:1042:59: warning: unused parameter 'field' [-Wunused-parameter] 1042 | void maxwell_scalarfield_otherhalf(maxwell_data *d, real *field) | ~~~~~~^~~~~ CCLD libmaxwell.la make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/maxwell' make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' CCLD libmpb.la make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' Making all in tests make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/tests' CC malloctest.o CC blastest.o CC eigs_test.o CC maxwell_test.o CC normal_vectors-normal_vectors.o eigs_test.c: In function 'Aop': eigs_test.c:293:51: warning: unused parameter 'data' [-Wunused-parameter] 293 | void Aop(evectmatrix Xin, evectmatrix Xout, void *data, | ~~~~~~^~~~ eigs_test.c:294:14: warning: unused parameter 'is_current_eigenvector' [-Wunused-parameter] 294 | int is_current_eigenvector, evectmatrix Work) | ~~~~^~~~~~~~~~~~~~~~~~~~~~ eigs_test.c:294:50: warning: unused parameter 'Work' [-Wunused-parameter] 294 | int is_current_eigenvector, evectmatrix Work) | ~~~~~~~~~~~~^~~~ eigs_test.c: In function 'Bop': eigs_test.c:303:51: warning: unused parameter 'data' [-Wunused-parameter] 303 | void Bop(evectmatrix Xin, evectmatrix Xout, void *data, | ~~~~~~^~~~ eigs_test.c:304:14: warning: unused parameter 'is_current_eigenvector' [-Wunused-parameter] 304 | int is_current_eigenvector, evectmatrix Work) | ~~~~^~~~~~~~~~~~~~~~~~~~~~ eigs_test.c:304:50: warning: unused parameter 'Work' [-Wunused-parameter] 304 | int is_current_eigenvector, evectmatrix Work) | ~~~~~~~~~~~~^~~~ eigs_test.c: In function 'Ainvop': eigs_test.c:313:54: warning: unused parameter 'data' [-Wunused-parameter] 313 | void Ainvop(evectmatrix Xin, evectmatrix Xout, void *data, | ~~~~~~^~~~ eigs_test.c:314:25: warning: unused parameter 'Y' [-Wunused-parameter] 314 | evectmatrix Y, real *eigenvals, sqmatrix YtY) | ~~~~~~~~~~~~^ eigs_test.c:314:34: warning: unused parameter 'eigenvals' [-Wunused-parameter] 314 | evectmatrix Y, real *eigenvals, sqmatrix YtY) | ~~~~~~^~~~~~~~~ eigs_test.c:314:54: warning: unused parameter 'YtY' [-Wunused-parameter] 314 | evectmatrix Y, real *eigenvals, sqmatrix YtY) | ~~~~~~~~~^~~ eigs_test.c: In function 'Cop_old': eigs_test.c:324:55: warning: unused parameter 'data' [-Wunused-parameter] 324 | void Cop_old(evectmatrix Xin, evectmatrix Xout, void *data, | ~~~~~~^~~~ eigs_test.c:325:22: warning: unused parameter 'Y' [-Wunused-parameter] 325 | evectmatrix Y, real *eigenvals, sqmatrix YtY) | ~~~~~~~~~~~~^ eigs_test.c:325:31: warning: unused parameter 'eigenvals' [-Wunused-parameter] 325 | evectmatrix Y, real *eigenvals, sqmatrix YtY) | ~~~~~~^~~~~~~~~ eigs_test.c: In function 'Cop': eigs_test.c:349:51: warning: unused parameter 'data' [-Wunused-parameter] 349 | void Cop(evectmatrix Xin, evectmatrix Xout, void *data, | ~~~~~~^~~~ eigs_test.c:350:23: warning: unused parameter 'Y' [-Wunused-parameter] 350 | evectmatrix Y, real *eigenvals, sqmatrix YtY) | ~~~~~~~~~~~~^ CCLD malloctest CCLD blastest CCLD normal_vectors CCLD maxwell_test CCLD eigs_test make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/tests' Making all in mpb make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' cp -f /usr/share/libctl/base/main.c main.c cp -f /usr/share/libctl/utils/geom.c geom.c gen-ctl-io --header -o ctl-io.h mpb.scm /usr/share/libctl gen-ctl-io --code -o ctl-io.c mpb.scm /usr/share/libctl ;;; note: auto-compilation is enabled, set GUILE_AUTO_COMPILE=0 ;;; or pass the --no-auto-compile argument to disable. ;;; compiling /usr/share/libctl/base/include.scm ;;; compiled /home/buildozer/.cache/guile/ccache/3.0-BE-8-4.6/usr/share/libctl/base/include.scm.go ;;; note: auto-compilation is enabled, set GUILE_AUTO_COMPILE=0 ;;; or pass the --no-auto-compile argument to disable. ;;; compiling /usr/share/libctl/base/include.scm ;;; compiled /home/buildozer/.cache/guile/ccache/3.0-BE-8-4.6/usr/share/libctl/base/include.scm.go make all-am make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' CC mpb-transform.o CC mpb-medium.o CC mpb-epsilon_file.o CC mpb-field-smob.o CC mpb-fields.o CC mpb-material_grid.o CC mpb-material_grid_opt.o CC mpb-matrix-smob.o In file included from /usr/include/guile/3.0/libguile/procs.h:25, from /usr/include/guile/3.0/libguile/threads.h:25, from /usr/include/guile/3.0/libguile/async.h:25, from /usr/include/guile/3.0/libguile.h:35, from /usr/include/ctl.h:28, from ctl-io.h:7, from field-smob.h:21, from field-smob.c:27: field-smob.c: In function 'field_p': my-smob.h:56:44: warning: comparison of integer expressions of different signedness: 'long unsigned int' and 'long int' [-Wsign-compare] 56 | && SCM_TYP16 (X) == (scm_tc16_smob_ ## T)) | ^~ field-smob.h:45:20: note: in expansion of macro 'T_SMOB_P' 45 | #define FIELD_P(X) T_SMOB_P(field_smob, X) | ^~~~~~~~ field-smob.c:42:40: note: in expansion of macro 'FIELD_P' 42 | return ctl_convert_boolean_to_scm(FIELD_P(obj)); | ^~~~~~~ field-smob.c: In function 'rscalar_field_p': my-smob.h:56:44: warning: comparison of integer expressions of different signedness: 'long unsigned int' and 'long int' [-Wsign-compare] 56 | && SCM_TYP16 (X) == (scm_tc16_smob_ ## T)) | ^~ field-smob.h:45:20: note: in expansion of macro 'T_SMOB_P' 45 | #define FIELD_P(X) T_SMOB_P(field_smob, X) | ^~~~~~~~ field-smob.h:49:29: note: in expansion of macro 'FIELD_P' 49 | #define RSCALAR_FIELD_P(X) (FIELD_P(X) && ((FIELD(X))->type == RSCALAR_FIELD_SMOB)) | ^~~~~~~ field-smob.c:47:40: note: in expansion of macro 'RSCALAR_FIELD_P' 47 | return ctl_convert_boolean_to_scm(RSCALAR_FIELD_P(obj)); | ^~~~~~~~~~~~~~~ field-smob.c: In function 'cscalar_field_p': my-smob.h:56:44: warning: comparison of integer expressions of different signedness: 'long unsigned int' and 'long int' [-Wsign-compare] 56 | && SCM_TYP16 (X) == (scm_tc16_smob_ ## T)) | ^~ field-smob.h:45:20: note: in expansion of macro 'T_SMOB_P' 45 | #define FIELD_P(X) T_SMOB_P(field_smob, X) | ^~~~~~~~ field-smob.h:50:29: note: in expansion of macro 'FIELD_P' 50 | #define CSCALAR_FIELD_P(X) (FIELD_P(X) && ((FIELD(X))->type == CSCALAR_FIELD_SMOB)) | ^~~~~~~ field-smob.c:52:40: note: in expansion of macro 'CSCALAR_FIELD_P' 52 | return ctl_convert_boolean_to_scm(CSCALAR_FIELD_P(obj)); | ^~~~~~~~~~~~~~~ material_grid_opt.c: In function 'material_grids_min_tetm_gap': material_grid_opt.c:90:43: warning: unused parameter 'func_tol' [-Wunused-parameter] 90 | number func_tol, number eps_tol, | ~~~~~~~^~~~~~~~ material_grid_opt.c:91:44: warning: unused parameter 'maxeval' [-Wunused-parameter] 91 | integer maxeval, number maxtime) | ~~~~~~~~^~~~~~~ material_grid_opt.c:91:60: warning: unused parameter 'maxtime' [-Wunused-parameter] 91 | integer maxeval, number maxtime) | ~~~~~~~^~~~~~~ field-smob.c: In function 'cvector_field_p': my-smob.h:56:44: warning: comparison of integer expressions of different signedness: 'long unsigned int' and 'long int' [-Wsign-compare] 56 | && SCM_TYP16 (X) == (scm_tc16_smob_ ## T)) | ^~ field-smob.h:45:20: note: in expansion of macro 'T_SMOB_P' 45 | #define FIELD_P(X) T_SMOB_P(field_smob, X) | ^~~~~~~~ field-smob.h:51:29: note: in expansion of macro 'FIELD_P' 51 | #define CVECTOR_FIELD_P(X) (FIELD_P(X) && ((FIELD(X))->type == CVECTOR_FIELD_SMOB)) | ^~~~~~~ field-smob.c:57:40: note: in expansion of macro 'CVECTOR_FIELD_P' 57 | return ctl_convert_boolean_to_scm(CVECTOR_FIELD_P(obj)); | ^~~~~~~~~~~~~~~ material_grid_opt.c: In function 'material_grids_maxmin_gap': material_grid_opt.c:258:48: warning: unused parameter 'func_tol' [-Wunused-parameter] 258 | number func_tol, number eps_tol, | ~~~~~~~^~~~~~~~ material_grid_opt.c:259:49: warning: unused parameter 'maxeval' [-Wunused-parameter] 259 | integer maxeval, number maxtime) | ~~~~~~~~^~~~~~~ material_grid_opt.c:259:65: warning: unused parameter 'maxtime' [-Wunused-parameter] 259 | integer maxeval, number maxtime) | ~~~~~~~^~~~~~~ In file included from field-smob.h:22: field-smob.c: In function 'assert_field_smob': my-smob.h:56:44: warning: comparison of integer expressions of different signedness: 'long unsigned int' and 'long int' [-Wsign-compare] 56 | && SCM_TYP16 (X) == (scm_tc16_smob_ ## T)) | ^~ my-smob.h:69:28: note: in expansion of macro 'T_SMOB_P' 69 | #define SAFE_T_SMOB(T, X) (T_SMOB_P(T,X) ? T_SMOB(T,X) : NULL) | ^~~~~~~~ field-smob.h:47:65: note: in expansion of macro 'SAFE_T_SMOB' 47 | #define SAFE_FIELD(X) (cur_fieldp(X) ? update_curfield_smob() : SAFE_T_SMOB(field_smob, X)) | ^~~~~~~~~~~ field-smob.c:187:22: note: in expansion of macro 'SAFE_FIELD' 187 | field_smob *f = SAFE_FIELD(fo); | ^~~~~~~~~~ In file included from medium.c:75: epsilon.c: In function 'get_epsilon': epsilon.c:415:51: warning: variable 'local_y_start' set but not used [-Wunused-but-set-variable] 415 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~~~~~~ epsilon.c:415:47: warning: variable 'nz' set but not used [-Wunused-but-set-variable] 415 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~ epsilon.c:415:43: warning: variable 'nx' set but not used [-Wunused-but-set-variable] 415 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~ epsilon.c:415:26: warning: variable 'last_dim_stored' set but not used [-Wunused-but-set-variable] 415 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~~~~~~~~ epsilon.c:415:16: warning: variable 'last_dim' set but not used [-Wunused-but-set-variable] 415 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~ field-smob.c: In function 'integrate_fieldL': field-smob.c:407:54: warning: variable 'last_dim' set but not used [-Wunused-but-set-variable] 407 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~~~ field-smob.c:407:48: warning: variable 'rank' set but not used [-Wunused-but-set-variable] 407 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~ field-smob.c:407:40: warning: variable 'n_last' set but not used [-Wunused-but-set-variable] 407 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~ field-smob.c:407:31: warning: variable 'n_other' set but not used [-Wunused-but-set-variable] 407 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~~ In file included from medium.c:76: mu.c: In function 'get_mu': mu.c:415:51: warning: variable 'local_y_start' set but not used [-Wunused-but-set-variable] 415 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~~~~~~ mu.c:415:47: warning: variable 'nz' set but not used [-Wunused-but-set-variable] 415 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~ mu.c:415:43: warning: variable 'nx' set but not used [-Wunused-but-set-variable] 415 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~ mu.c:415:26: warning: variable 'last_dim_stored' set but not used [-Wunused-but-set-variable] 415 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~~~~~~~~ mu.c:415:16: warning: variable 'last_dim' set but not used [-Wunused-but-set-variable] 415 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~ field-smob.c: At top level: field-smob.c:395:17: warning: 'cvector3_conj' defined but not used [-Wunused-function] 395 | static cvector3 cvector3_conj(cvector3 c) | ^~~~~~~~~~~~~ In file included from /usr/include/guile/3.0/libguile/procs.h:25, from /usr/include/guile/3.0/libguile/threads.h:25, from /usr/include/guile/3.0/libguile/async.h:25, from /usr/include/guile/3.0/libguile.h:35, from /usr/include/ctl.h:28, from ctl-io.h:7, from matrix-smob.h:21, from matrix-smob.c:29: matrix-smob.c: In function 'evectmatrix_p': my-smob.h:56:44: warning: comparison of integer expressions of different signedness: 'long unsigned int' and 'long int' [-Wsign-compare] 56 | && SCM_TYP16 (X) == (scm_tc16_smob_ ## T)) | ^~ matrix-smob.h:29:26: note: in expansion of macro 'T_SMOB_P' 29 | #define EVECTMATRIX_P(X) T_SMOB_P(evectmatrix, X) | ^~~~~~~~ matrix-smob.c:44:40: note: in expansion of macro 'EVECTMATRIX_P' 44 | return ctl_convert_boolean_to_scm(EVECTMATRIX_P(obj)); | ^~~~~~~~~~~~~ matrix-smob.c: In function 'sqmatrix_p': my-smob.h:56:44: warning: comparison of integer expressions of different signedness: 'long unsigned int' and 'long int' [-Wsign-compare] 56 | && SCM_TYP16 (X) == (scm_tc16_smob_ ## T)) | ^~ matrix-smob.h:33:23: note: in expansion of macro 'T_SMOB_P' 33 | #define SQMATRIX_P(X) T_SMOB_P(sqmatrix, X) | ^~~~~~~~ matrix-smob.c:85:40: note: in expansion of macro 'SQMATRIX_P' 85 | return ctl_convert_boolean_to_scm(SQMATRIX_P(obj)); | ^~~~~~~~~~ In file included from matrix-smob.h:22: matrix-smob.c: In function 'assert_sqmatrix_smob': my-smob.h:56:44: warning: comparison of integer expressions of different signedness: 'long unsigned int' and 'long int' [-Wsign-compare] 56 | && SCM_TYP16 (X) == (scm_tc16_smob_ ## T)) | ^~ my-smob.h:69:28: note: in expansion of macro 'T_SMOB_P' 69 | #define SAFE_T_SMOB(T, X) (T_SMOB_P(T,X) ? T_SMOB(T,X) : NULL) | ^~~~~~~~ matrix-smob.h:35:26: note: in expansion of macro 'SAFE_T_SMOB' 35 | #define SAFE_SQMATRIX(X) SAFE_T_SMOB(sqmatrix, X) | ^~~~~~~~~~~ matrix-smob.c:186:20: note: in expansion of macro 'SAFE_SQMATRIX' 186 | sqmatrix *m = SAFE_SQMATRIX(mo); | ^~~~~~~~~~~~~ matrix-smob.c: In function 'assert_evectmatrix_smob': my-smob.h:56:44: warning: comparison of integer expressions of different signedness: 'long unsigned int' and 'long int' [-Wsign-compare] 56 | && SCM_TYP16 (X) == (scm_tc16_smob_ ## T)) | ^~ my-smob.h:69:28: note: in expansion of macro 'T_SMOB_P' 69 | #define SAFE_T_SMOB(T, X) (T_SMOB_P(T,X) ? T_SMOB(T,X) : NULL) | ^~~~~~~~ matrix-smob.h:31:29: note: in expansion of macro 'SAFE_T_SMOB' 31 | #define SAFE_EVECTMATRIX(X) SAFE_T_SMOB(evectmatrix, X) | ^~~~~~~~~~~ matrix-smob.c:193:23: note: in expansion of macro 'SAFE_EVECTMATRIX' 193 | evectmatrix *m = SAFE_EVECTMATRIX(mo); | ^~~~~~~~~~~~~~~~ material_grid.c: In function 'material_grids_addgradient': material_grid.c:390:54: warning: variable 'last_dim' set but not used [-Wunused-but-set-variable] 390 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~~~ material_grid.c:390:48: warning: variable 'rank' set but not used [-Wunused-but-set-variable] 390 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~ material_grid.c:390:40: warning: variable 'n_last' set but not used [-Wunused-but-set-variable] 390 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~ material_grid.c:390:31: warning: variable 'n_other' set but not used [-Wunused-but-set-variable] 390 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~~ material_grid.c:390:16: warning: unused variable 'k' [-Wunused-variable] 390 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^ material_grid.c:390:13: warning: unused variable 'j' [-Wunused-variable] 390 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^ material_grid.c:390:10: warning: unused variable 'i' [-Wunused-variable] 390 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^ material_grid.c: In function 'print_material_grids_deps_du_numeric': material_grid.c:634:14: warning: unused variable 'foo' [-Wunused-variable] 634 | double *foo; | ^~~ In file included from material_grid.c:25: material_grid.c: In function 'synchronize_material_grid': ../src/util/mpiglue.h:58:40: warning: statement with no effect [-Wunused-value] 58 | #define MPI_Bcast(b, n, t, root, comm) 0 | ^ material_grid.c:693:6: note: in expansion of macro 'MPI_Bcast' 693 | MPI_Bcast(grid, n, MPI_DOUBLE, 0, mpb_comm); | ^~~~~~~~~ material_grid.c:691:10: warning: unused variable 'n' [-Wunused-variable] 691 | int n = ((int) g->size.x) * ((int) g->size.y) * ((int) g->size.z); | ^ material_grid.c:690:14: warning: variable 'grid' set but not used [-Wunused-but-set-variable] 690 | double *grid; | ^~~~ material_grid.c: In function 'match_eps_func': material_grid.c:816:54: warning: variable 'last_dim' set but not used [-Wunused-but-set-variable] 816 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~~~ material_grid.c:816:48: warning: variable 'rank' set but not used [-Wunused-but-set-variable] 816 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~ material_grid.c:816:40: warning: variable 'n_last' set but not used [-Wunused-but-set-variable] 816 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~ material_grid.c:816:31: warning: variable 'n_other' set but not used [-Wunused-but-set-variable] 816 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~~ material_grid.c: At top level: material_grid.c:807:15: warning: 'match_eps_func' defined but not used [-Wunused-function] 807 | static double match_eps_func(int n, const double *u, double *grad, void *data) | ^~~~~~~~~~~~~~ fields.c: In function 'compute_field_energy_internal': fields.c:229:50: warning: variable 'local_y_start' set but not used [-Wunused-but-set-variable] 229 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~~~~~~ fields.c:229:46: warning: variable 'nz' set but not used [-Wunused-but-set-variable] 229 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~ fields.c:229:42: warning: variable 'nx' set but not used [-Wunused-but-set-variable] 229 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~ fields.c:229:25: warning: variable 'last_dim_stored' set but not used [-Wunused-but-set-variable] 229 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~~~~~~~~ fields.c:229:15: warning: variable 'last_dim' set but not used [-Wunused-but-set-variable] 229 | int i, N, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~ In file included from fields.c:25: fields.c: In function 'fix_field_phase': ../src/util/mpiglue.h:54:17: warning: comparison between two arrays [-Warray-compare] 54 | CHECK((sb) != (rb), "MPI_Allreduce doesn't work for sendbuf == recvbuf");\ | ^~ ../src/util/check.h:37:12: note: in definition of macro 'CHECK' 37 | if (!(condition)) { \ | ^~~~~~~~~ fields.c:465:6: note: in expansion of macro 'mpi_allreduce' 465 | mpi_allreduce(sq_sum2, sq_sum, 2, real, SCALAR_MPI_TYPE, | ^~~~~~~~~~~~~ ../src/util/mpiglue.h:54:17: note: use '&sq_sum2[0] != &sq_sum[0]' to compare the addresses 54 | CHECK((sb) != (rb), "MPI_Allreduce doesn't work for sendbuf == recvbuf");\ | ^~ ../src/util/check.h:37:12: note: in definition of macro 'CHECK' 37 | if (!(condition)) { \ | ^~~~~~~~~ fields.c:465:6: note: in expansion of macro 'mpi_allreduce' 465 | mpi_allreduce(sq_sum2, sq_sum, 2, real, SCALAR_MPI_TYPE, | ^~~~~~~~~~~~~ fields.c: In function 'get_val': fields.c:557:25: warning: unused parameter 'nx' [-Wunused-parameter] 557 | int nx, int ny, int nz, int last_dim_size, | ~~~~^~ fields.c:557:49: warning: unused parameter 'last_dim_size' [-Wunused-parameter] 557 | int nx, int ny, int nz, int last_dim_size, | ~~~~^~~~~~~~~~~~~ fields.c: At top level: fields.c:709:44: warning: argument 1 of type 'scalar_complex[3]' with mismatched bound [-Warray-parameter=] 709 | void get_bloch_field_point_(scalar_complex field[3], vector3 p) | ~~~~~~~~~~~~~~~^~~~~~~~ In file included from fields.c:35: mpb.h:77:45: note: previously declared as 'scalar_complex *' 77 | void get_bloch_field_point_(scalar_complex *field, vector3 p); | ~~~~~~~~~~~~~~~~^~~~~ fields.c: In function 'get_bloch_field_point': fields.c:721:15: warning: unused variable 'F' [-Wunused-variable] 721 | cvector3 F; | ^ fields.c: In function 'compute_energy_in_dielectric': fields.c:870:51: warning: variable 'local_y_start' set but not used [-Wunused-but-set-variable] 870 | int N, i, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~~~~~~ fields.c:870:47: warning: variable 'nz' set but not used [-Wunused-but-set-variable] 870 | int N, i, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~ fields.c:870:43: warning: variable 'nx' set but not used [-Wunused-but-set-variable] 870 | int N, i, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~ fields.c:870:26: warning: variable 'last_dim_stored' set but not used [-Wunused-but-set-variable] 870 | int N, i, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~~~~~~~~ fields.c:870:16: warning: variable 'last_dim' set but not used [-Wunused-but-set-variable] 870 | int N, i, last_dim, last_dim_stored, nx, nz, local_y_start; | ^~~~~~~~ fields.c: In function 'output_scalarfield': fields.c:943:36: warning: unused parameter 'last_dim_index' [-Wunused-parameter] 943 | int last_dim_index, | ~~~~^~~~~~~~~~~~~~ fields.c:944:36: warning: unused parameter 'last_dim_start' [-Wunused-parameter] 944 | int last_dim_start, int last_dim_size, | ~~~~^~~~~~~~~~~~~~ fields.c:944:56: warning: unused parameter 'last_dim_size' [-Wunused-parameter] 944 | int last_dim_start, int last_dim_size, | ~~~~^~~~~~~~~~~~~ fields.c:945:36: warning: unused parameter 'first_dim_start' [-Wunused-parameter] 945 | int first_dim_start, int first_dim_size, | ~~~~^~~~~~~~~~~~~~~ fields.c:945:57: warning: unused parameter 'first_dim_size' [-Wunused-parameter] 945 | int first_dim_start, int first_dim_size, | ~~~~^~~~~~~~~~~~~~ fields.c:946:36: warning: unused parameter 'write_start0_special' [-Wunused-parameter] 946 | int write_start0_special) | ~~~~^~~~~~~~~~~~~~~~~~~~ fields.c: In function 'compute_energy_in_object_list': fields.c:1319:6: warning: 'geom_fix_objects0' is deprecated [-Wdeprecated-declarations] 1319 | geom_fix_objects0(objects); | ^~~~~~~~~~~~~~~~~ In file included from fields.c:33: /usr/include/ctlgeom.h:65:13: note: declared here 65 | extern void geom_fix_objects0(GEOMETRIC_OBJECT_LIST geometry) CTLGEOM_DEPRECATED; | ^~~~~~~~~~~~~~~~~ fields.c:1306:54: warning: variable 'last_dim' set but not used [-Wunused-but-set-variable] 1306 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~~~ fields.c:1306:48: warning: variable 'rank' set but not used [-Wunused-but-set-variable] 1306 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~ fields.c:1306:40: warning: variable 'n_last' set but not used [-Wunused-but-set-variable] 1306 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~ fields.c:1306:31: warning: variable 'n_other' set but not used [-Wunused-but-set-variable] 1306 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~~ fields.c:1306:16: warning: unused variable 'k' [-Wunused-variable] 1306 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^ fields.c:1306:13: warning: unused variable 'j' [-Wunused-variable] 1306 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^ fields.c:1306:10: warning: unused variable 'i' [-Wunused-variable] 1306 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^ fields.c: In function 'compute_field_integral': fields.c:1374:54: warning: variable 'last_dim' set but not used [-Wunused-but-set-variable] 1374 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~~~ fields.c:1374:48: warning: variable 'rank' set but not used [-Wunused-but-set-variable] 1374 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~ fields.c:1374:40: warning: variable 'n_last' set but not used [-Wunused-but-set-variable] 1374 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~ fields.c:1374:31: warning: variable 'n_other' set but not used [-Wunused-but-set-variable] 1374 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^~~~~~~ fields.c:1374:16: warning: unused variable 'k' [-Wunused-variable] 1374 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^ fields.c:1374:13: warning: unused variable 'j' [-Wunused-variable] 1374 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^ fields.c:1374:10: warning: unused variable 'i' [-Wunused-variable] 1374 | int i, j, k, n1, n2, n3, n_other, n_last, rank, last_dim; | ^ CC mpb-mpb.o CC mpb-main.o CC mpb-geom.o In file included from mpb.c:35: mpb.c: In function 'ctl_start_hook': ../src/util/mpiglue.h:42:29: warning: statement with no effect [-Wunused-value] 42 | #define MPI_Init(argc,argv) 0 | ^ mpb.c:99:6: note: in expansion of macro 'MPI_Init' 99 | MPI_Init(argc, argv); | ^~~~~~~~ mpb.c:97:26: warning: unused parameter 'argc' [-Wunused-parameter] 97 | void ctl_start_hook(int *argc, char ***argv) | ~~~~~^~~~ mpb.c:97:40: warning: unused parameter 'argv' [-Wunused-parameter] 97 | void ctl_start_hook(int *argc, char ***argv) | ~~~~~~~~^~~~ mpb.c: In function 'ctl_stop_hook': ../src/util/mpiglue.h:43:24: warning: statement with no effect [-Wunused-value] 43 | #define MPI_Finalize() 0 | ^ mpb.c:143:6: note: in expansion of macro 'MPI_Finalize' 143 | MPI_Finalize(); | ^~~~~~~~~~~~ In file included from ctl-io.h:7, from main.c:37: main.c: In function 'main_entry': main.c:194:31: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 194 | gh_new_procedure("subplex", (SCM(*)(void))subplex_scm, 7, 0, 0); | ^ main.c:202:44: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 202 | gh_new_procedure("adaptive-integration", (SCM(*)(void))adaptive_integration_scm, 6, 0, 0); | ^ main.c:205:45: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 205 | gh_new_procedure("cadaptive-integration", (SCM(*)(void))cadaptive_integration_scm, 6, 0, 0); | ^ geom.c: In function 'overlap_integrand': geom.c:1156:10: warning: unused variable 's' [-Wunused-variable] 1156 | double s[2]; | ^ CC mpb-ctl-io.o geom.c: In function 'point_in_or_on_prism': geom.c:2110:10: warning: unused variable 'height' [-Wunused-variable] 2110 | double height = prsm->height; | ^~~~~~ geom.c: In function 'dcmp': geom.c:2137:17: warning: cast discards 'const' qualifier from pointer target type [-Wcast-qual] 2137 | double d1 = *((double *)pd1), d2 = *((double *)pd2); | ^ geom.c:2137:40: warning: cast discards 'const' qualifier from pointer target type [-Wcast-qual] 2137 | double d1 = *((double *)pd1), d2 = *((double *)pd2); | ^ geom.c: In function 'normal_to_prism': geom.c:2357:11: warning: unused variable 'axisp' [-Wunused-variable] 2357 | vector3 axisp = vector3_scale(height, zhatp); | ^~~~~ geom.c: In function 'triangle_normal': geom.c:2548:10: warning: unused variable 'nvnorm' [-Wunused-variable] 2548 | double nvnorm = vector3_norm(nv); | ^~~~~~ sed "s%\./mpb%"/usr/bin/"mpb%g" \ mpb-split-preinstall > mpb-split chmod u+x mpb-split ctl-io.c: In function 'block_equal': ctl-io.c:639:6: warning: suggest braces around empty body in an 'else' statement [-Wempty-body] 639 | else ; | ^ ctl-io.c: In function 'cylinder_equal': ctl-io.c:679:6: warning: suggest braces around empty body in an 'else' statement [-Wempty-body] 679 | else ; | ^ ctl-io.c: In function 'geometric_object_equal': ctl-io.c:716:6: warning: suggest braces around empty body in an 'else' statement [-Wempty-body] 716 | else ; | ^ ctl-io.c: In function 'material_type_equal': ctl-io.c:776:6: warning: suggest braces around empty body in an 'else' statement [-Wempty-body] 776 | else ; | ^ ctl-io.c: In function 'lattice_destroy': ctl-io.c:782:30: warning: unused parameter 'o' [-Wunused-parameter] 782 | void lattice_destroy(lattice o) | ~~~~~~~~^ ctl-io.c: In function 'ellipsoid_destroy': ctl-io.c:786:34: warning: unused parameter 'o' [-Wunused-parameter] 786 | void ellipsoid_destroy(ellipsoid o) | ~~~~~~~~~~^ ctl-io.c: In function 'sphere_destroy': ctl-io.c:845:28: warning: unused parameter 'o' [-Wunused-parameter] 845 | void sphere_destroy(sphere o) | ~~~~~~~^ ctl-io.c: In function 'wedge_destroy': ctl-io.c:849:26: warning: unused parameter 'o' [-Wunused-parameter] 849 | void wedge_destroy(wedge o) | ~~~~~~^ ctl-io.c: In function 'cone_destroy': ctl-io.c:853:24: warning: unused parameter 'o' [-Wunused-parameter] 853 | void cone_destroy(cone o) | ~~~~~^ ctl-io.c: In function 'material_grid_destroy': ctl-io.c:907:42: warning: unused parameter 'o' [-Wunused-parameter] 907 | void material_grid_destroy(material_grid o) | ~~~~~~~~~~~~~~^ ctl-io.c: In function 'material_function_destroy': ctl-io.c:911:50: warning: unused parameter 'o' [-Wunused-parameter] 911 | void material_function_destroy(material_function o) | ~~~~~~~~~~~~~~~~~~^ ctl-io.c: In function 'medium_anisotropic_destroy': ctl-io.c:915:52: warning: unused parameter 'o' [-Wunused-parameter] 915 | void medium_anisotropic_destroy(medium_anisotropic o) | ~~~~~~~~~~~~~~~~~~~^ ctl-io.c: In function 'medium_destroy': ctl-io.c:919:28: warning: unused parameter 'o' [-Wunused-parameter] 919 | void medium_destroy(medium o) | ~~~~~~~^ In file included from ctl-io.h:7, from ctl-io.c:7: ctl-io.c: In function 'export_external_functions': ctl-io.c:2885:61: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2885 | gh_new_procedure("compute-1-group-velocity-reciprocal-aux", (SCM (*)()) compute_1_group_velocity_reciprocal_aux, 1, 0, 0); | ^ ctl-io.c:2886:50: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2886 | gh_new_procedure("compute-1-group-velocity-aux", (SCM (*)()) compute_1_group_velocity_aux, 1, 0, 0); | ^ ctl-io.c:2887:60: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2887 | gh_new_procedure("compute-1-group-velocity-component-aux", (SCM (*)()) compute_1_group_velocity_component_aux, 2, 0, 0); | ^ ctl-io.c:2888:58: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2888 | gh_new_procedure("compute-group-velocity-component-aux", (SCM (*)()) compute_group_velocity_component_aux, 1, 0, 0); | ^ ctl-io.c:2891:53: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2891 | gh_new_procedure("material-grids-min-tetm-gap-aux", (SCM (*)()) material_grids_min_tetm_gap_aux, 6, 0, 0); | ^ ctl-io.c:2892:47: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2892 | gh_new_procedure("material-grids-mingap-aux", (SCM (*)()) material_grids_mingap_aux, 7, 0, 0); | ^ ctl-io.c:2893:47: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2893 | gh_new_procedure("material-grids-maxgap-aux", (SCM (*)()) material_grids_maxgap_aux, 7, 0, 0); | ^ ctl-io.c:2894:56: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2894 | gh_new_procedure("material-grids-approx-gradient-aux", (SCM (*)()) material_grids_approx_gradient_aux, 4, 0, 0); | ^ ctl-io.c:2895:62: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2895 | gh_new_procedure("print-material-grids-deps-du-numeric-aux", (SCM (*)()) print_material_grids_deps_du_numeric_aux, 1, 0, 0); | ^ ctl-io.c:2897:55: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2897 | gh_new_procedure("print-material-grids-gradient-aux", (SCM (*)()) print_material_grids_gradient_aux, 1, 0, 0); | ^ ctl-io.c:2898:60: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2898 | gh_new_procedure("material-grids-match-epsilon-file!-aux", (SCM (*)()) material_grids_match_epsilon_fileB_aux, 2, 0, 0); | ^ ctl-io.c:2899:45: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2899 | gh_new_procedure("load-material-grid!-aux", (SCM (*)()) load_material_gridB_aux, 3, 0, 0); | ^ ctl-io.c:2900:44: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2900 | gh_new_procedure("save-material-grid-aux", (SCM (*)()) save_material_grid_aux, 2, 0, 0); | ^ ctl-io.c:2901:50: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2901 | gh_new_procedure("randomize-material-grid!-aux", (SCM (*)()) randomize_material_gridB_aux, 2, 0, 0); | ^ ctl-io.c:2902:55: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2902 | gh_new_procedure("cvector-field-get-point-bloch-aux", (SCM (*)()) cvector_field_get_point_bloch_aux, 2, 0, 0); | ^ ctl-io.c:2903:49: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2903 | gh_new_procedure("cvector-field-get-point-aux", (SCM (*)()) cvector_field_get_point_aux, 2, 0, 0); | ^ ctl-io.c:2904:49: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2904 | gh_new_procedure("cscalar-field-get-point-aux", (SCM (*)()) cscalar_field_get_point_aux, 2, 0, 0); | ^ ctl-io.c:2905:49: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2905 | gh_new_procedure("rscalar-field-get-point-aux", (SCM (*)()) rscalar_field_get_point_aux, 2, 0, 0); | ^ ctl-io.c:2906:42: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2906 | gh_new_procedure("integrate-fieldL-aux", (SCM (*)()) integrate_fieldL_aux, 2, 0, 0); | ^ ctl-io.c:2907:37: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2907 | gh_new_procedure("field-mapL!-aux", (SCM (*)()) field_mapLB_aux, 3, 0, 0); | ^ ctl-io.c:2908:36: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2908 | gh_new_procedure("field-load-aux", (SCM (*)()) field_load_aux, 1, 0, 0); | ^ ctl-io.c:2909:36: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2909 | gh_new_procedure("field-set!-aux", (SCM (*)()) field_setB_aux, 2, 0, 0); | ^ ctl-io.c:2910:41: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2910 | gh_new_procedure("fields-conform?-aux", (SCM (*)()) fields_conformp_aux, 2, 0, 0); | ^ ctl-io.c:2911:36: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2911 | gh_new_procedure("field-make-aux", (SCM (*)()) field_make_aux, 1, 0, 0); | ^ ctl-io.c:2912:49: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2912 | gh_new_procedure("cvector-field-nonbloch!-aux", (SCM (*)()) cvector_field_nonblochB_aux, 1, 0, 0); | ^ ctl-io.c:2913:44: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2913 | gh_new_procedure("cvector-field-make-aux", (SCM (*)()) cvector_field_make_aux, 1, 0, 0); | ^ ctl-io.c:2914:44: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2914 | gh_new_procedure("rscalar-field-make-aux", (SCM (*)()) rscalar_field_make_aux, 1, 0, 0); | ^ ctl-io.c:2915:36: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2915 | gh_new_procedure("cur-field?-aux", (SCM (*)()) cur_fieldp_aux, 1, 0, 0); | ^ ctl-io.c:2916:48: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2916 | gh_new_procedure("get-dominant-planewave-aux", (SCM (*)()) get_dominant_planewave_aux, 1, 0, 0); | ^ ctl-io.c:2917:43: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2917 | gh_new_procedure("load-eigenvectors-aux", (SCM (*)()) load_eigenvectors_aux, 1, 0, 0); | ^ ctl-io.c:2918:43: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2918 | gh_new_procedure("save-eigenvectors-aux", (SCM (*)()) save_eigenvectors_aux, 1, 0, 0); | ^ ctl-io.c:2919:44: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2919 | gh_new_procedure("input-eigenvectors-aux", (SCM (*)()) input_eigenvectors_aux, 2, 0, 0); | ^ ctl-io.c:2920:45: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2920 | gh_new_procedure("output-eigenvectors-aux", (SCM (*)()) output_eigenvectors_aux, 2, 0, 0); | ^ ctl-io.c:2921:43: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2921 | gh_new_procedure("scale-eigenvector-aux", (SCM (*)()) scale_eigenvector_aux, 2, 0, 0); | ^ ctl-io.c:2922:42: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2922 | gh_new_procedure("dot-eigenvectors-aux", (SCM (*)()) dot_eigenvectors_aux, 2, 0, 0); | ^ ctl-io.c:2923:42: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2923 | gh_new_procedure("set-eigenvectors-aux", (SCM (*)()) set_eigenvectors_aux, 2, 0, 0); | ^ ctl-io.c:2924:42: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2924 | gh_new_procedure("get-eigenvectors-aux", (SCM (*)()) get_eigenvectors_aux, 2, 0, 0); | ^ ctl-io.c:2925:42: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2925 | gh_new_procedure("sqmatrix-eigvals-aux", (SCM (*)()) sqmatrix_eigvals_aux, 1, 0, 0); | ^ ctl-io.c:2926:40: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2926 | gh_new_procedure("sqmatrix-diagm-aux", (SCM (*)()) sqmatrix_diagm_aux, 1, 0, 0); | ^ ctl-io.c:2927:39: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2927 | gh_new_procedure("sqmatrix-mult-aux", (SCM (*)()) sqmatrix_mult_aux, 2, 0, 0); | ^ ctl-io.c:2928:38: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2928 | gh_new_procedure("sqmatrix-set-aux", (SCM (*)()) sqmatrix_set_aux, 4, 0, 0); | ^ ctl-io.c:2929:38: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2929 | gh_new_procedure("sqmatrix-ref-aux", (SCM (*)()) sqmatrix_ref_aux, 3, 0, 0); | ^ ctl-io.c:2930:39: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2930 | gh_new_procedure("sqmatrix-size-aux", (SCM (*)()) sqmatrix_size_aux, 1, 0, 0); | ^ ctl-io.c:2931:42: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2931 | gh_new_procedure("set-kpoint-index-aux", (SCM (*)()) set_kpoint_index_aux, 1, 0, 0); | ^ ctl-io.c:2935:33: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2935 | gh_new_procedure("mpi-max-aux", (SCM (*)()) mpi_max_aux, 1, 0, 0); | ^ ctl-io.c:2940:46: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2940 | gh_new_procedure("output-field-to-file-aux", (SCM (*)()) output_field_to_file_aux, 2, 0, 0); | ^ ctl-io.c:2941:42: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2941 | gh_new_procedure("compute-symmetry-aux", (SCM (*)()) compute_symmetry_aux, 3, 0, 0); | ^ ctl-io.c:2942:44: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2942 | gh_new_procedure("compute-symmetries-aux", (SCM (*)()) compute_symmetries_aux, 2, 0, 0); | ^ ctl-io.c:2943:45: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2943 | gh_new_procedure("transformed-overlap-aux", (SCM (*)()) transformed_overlap_aux, 2, 0, 0); | ^ ctl-io.c:2944:55: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2944 | gh_new_procedure("compute-energy-in-object-list-aux", (SCM (*)()) compute_energy_in_object_list_aux, 1, 0, 0); | ^ ctl-io.c:2945:49: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2945 | gh_new_procedure("compute-energy-integral-aux", (SCM (*)()) compute_energy_integral_aux, 1, 0, 0); | ^ ctl-io.c:2946:48: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2946 | gh_new_procedure("compute-field-integral-aux", (SCM (*)()) compute_field_integral_aux, 1, 0, 0); | ^ ctl-io.c:2947:54: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2947 | gh_new_procedure("compute-energy-in-dielectric-aux", (SCM (*)()) compute_energy_in_dielectric_aux, 2, 0, 0); | ^ ctl-io.c:2948:43: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2948 | gh_new_procedure("get-cscalar-point-aux", (SCM (*)()) get_cscalar_point_aux, 1, 0, 0); | ^ ctl-io.c:2949:49: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2949 | gh_new_procedure("get-bloch-cscalar-point-aux", (SCM (*)()) get_bloch_cscalar_point_aux, 1, 0, 0); | ^ ctl-io.c:2950:41: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2950 | gh_new_procedure("get-field-point-aux", (SCM (*)()) get_field_point_aux, 1, 0, 0); | ^ ctl-io.c:2951:47: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2951 | gh_new_procedure("get-bloch-field-point-aux", (SCM (*)()) get_bloch_field_point_aux, 1, 0, 0); | ^ ctl-io.c:2952:42: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2952 | gh_new_procedure("get-energy-point-aux", (SCM (*)()) get_energy_point_aux, 1, 0, 0); | ^ ctl-io.c:2953:58: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2953 | gh_new_procedure("get-epsilon-inverse-tensor-point-aux", (SCM (*)()) get_epsilon_inverse_tensor_point_aux, 1, 0, 0); | ^ ctl-io.c:2954:43: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2954 | gh_new_procedure("get-epsilon-point-aux", (SCM (*)()) get_epsilon_point_aux, 1, 0, 0); | ^ ctl-io.c:2961:36: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2961 | gh_new_procedure("get-bfield-aux", (SCM (*)()) get_bfield_aux, 1, 0, 0); | ^ ctl-io.c:2962:36: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2962 | gh_new_procedure("get-hfield-aux", (SCM (*)()) get_hfield_aux, 1, 0, 0); | ^ ctl-io.c:2963:36: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2963 | gh_new_procedure("get-dfield-aux", (SCM (*)()) get_dfield_aux, 1, 0, 0); | ^ ctl-io.c:2964:38: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2964 | gh_new_procedure("solve-kpoint-aux", (SCM (*)()) solve_kpoint_aux, 1, 0, 0); | ^ ctl-io.c:2966:36: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2966 | gh_new_procedure("set-parity-aux", (SCM (*)()) set_parity_aux, 1, 0, 0); | ^ ctl-io.c:2968:37: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2968 | gh_new_procedure("init-params-aux", (SCM (*)()) init_params_aux, 2, 0, 0); | ^ ctl-io.c:2969:38: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2969 | gh_new_procedure("square-basis-aux", (SCM (*)()) square_basis_aux, 2, 0, 0); | ^ ctl-io.c:2970:51: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2970 | gh_new_procedure("range-overlap-with-object-aux", (SCM (*)()) range_overlap_with_object_aux, 5, 0, 0); | ^ ctl-io.c:2971:55: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2971 | gh_new_procedure("display-geometric-object-info-aux", (SCM (*)()) display_geometric_object_info_aux, 2, 0, 0); | ^ ctl-io.c:2972:51: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2972 | gh_new_procedure("point-in-periodic-object?-aux", (SCM (*)()) point_in_periodic_objectp_aux, 2, 0, 0); | ^ ctl-io.c:2973:42: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2973 | gh_new_procedure("normal-to-object-aux", (SCM (*)()) normal_to_object_aux, 2, 0, 0); | ^ ctl-io.c:2974:42: warning: cast between incompatible function types from 'struct scm_unused_struct * (*)(struct scm_unused_struct *, struct scm_unused_struct *)' to 'struct scm_unused_struct * (*)(void)' [-Wcast-function-type] 2974 | gh_new_procedure("point-in-object?-aux", (SCM (*)()) point_in_objectp_aux, 2, 0, 0); | ^ CCLD mpb make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' Making all in examples make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/examples' make[2]: Nothing to be done for 'all'. make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/examples' Making all in utils make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/utils' CC mpb_data-mpb-data.o CCLD mpb-data make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/utils' make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0' make[2]: Nothing to be done for 'all-am'. make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0' make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0' Making check in src make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make check-recursive make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' Making check in util make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' make check-am make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' make[4]: Nothing to be done for 'check-am'. make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' Making check in matrices make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrices' make[3]: Nothing to be done for 'check'. make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrices' Making check in matrixio make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrixio' make[3]: Nothing to be done for 'check'. make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrixio' Making check in maxwell make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/maxwell' make[3]: Nothing to be done for 'check'. make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/maxwell' make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make[3]: Nothing to be done for 'check-am'. make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' Making check in tests make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/tests' make check-local make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/tests' ./blastest | sed 's/\-0\.000\([ ,)]\)/ 0.000\1/g' | sed 's/\-0\.000$/ 0.000/g' > blastest.out ./maxwell_test -1 -c 1e-9 -x 256 -E 1e-3 > maxwell_test.out ./maxwell_test -1 -c 1e-9 -x 256 -E 1e-3 -e -k 0.4 -n 1 > maxwell_test_2.out diff ./blastest.complex.out blastest.out ********************************************************** PASSED tests. ********************************************************** make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/tests' make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/tests' Making check in mpb make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' make check-am make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' make check-local make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' ./mpb ../examples/check.ctl ************************************************************************** Test case: 1d quarter-wave stack. ************************************************************************** init-params: initializing eigensolver data Computing 8 bands with 1.000000e-13 tolerance. Working in 1 dimensions. Grid size is 32 x 1 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 1e+20, height 0.25, axis (1, 0, 0) epsilon = 9, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Allocating fields... 6 k-points: (0,0,0) (0.1,0,0) (0.2,0,0) (0.3,0,0) (0.4,0,0) (0.5,0,0) Solving for band polarization: tm. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-9, mean 3, harm. mean 1.32353, 28.125% > 1, 25% "fill" Outputting check-epsilon... Warning! Library version information error. The HDF5 library version information are not consistent in its source code. This is NOT a fatal error but should be corrected. Setting the environment variable 'HDF5_DISABLE_VERSION_CHECK' to a value of 1 will suppress this warning. Library version information are: H5_VERS_MAJOR=1, H5_VERS_MINOR=14, H5_VERS_RELEASE=4, H5_VERS_SUBRELEASE=2, H5_VERS_INFO=HDF5 library version: 1.14.4 solve_kpoint (0,0,0): tmfreqs:, k index, k1, k2, k3, kmag/2pi, tm band 1, tm band 2, tm band 3, tm band 4, tm band 5, tm band 6, tm band 7, tm band 8 Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 30 iterations. tmfreqs:, 1, 0, 0, 0, 0, 0, 0.666384, 0.666668, 1.33099, 1.33336, 1.99162, 2.00024, 2.64509 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 74 iterations. tmfreqs:, 2, 0.1, 0, 0, 0.1, 0.0574931, 0.608788, 0.724374, 1.2736, 1.39098, 1.93585, 2.05634, 2.59377 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tmfreqs:, 3, 0.2, 0, 0, 0.2, 0.113352, 0.552761, 0.780756, 1.21673, 1.44856, 1.87756, 2.11569, 2.53392 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 11 iterations. tmfreqs:, 4, 0.3, 0, 0, 0.3, 0.164802, 0.501202, 0.832987, 1.16418, 1.5025, 1.82314, 2.17225, 2.4769 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 11 iterations. tmfreqs:, 5, 0.4, 0, 0, 0.4, 0.205537, 0.460405, 0.874781, 1.1222, 1.54664, 1.77873, 2.22033, 2.42855 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 6, 0.5, 0, 0, 0.5, 0.222451, 0.443471, 0.892368, 1.10456, 1.56579, 1.75949, 2.24248, 2.40631 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.22245099319196915 at #(0.5 0.0 0.0) Band 2 range: 0.44347071815372147 at #(0.5 0.0 0.0) to 0.6663842825289141 at #(0.0 0.0 0.0) Band 3 range: 0.6666675180317079 at #(0.0 0.0 0.0) to 0.8923678144024804 at #(0.5 0.0 0.0) Band 4 range: 1.104558052488972 at #(0.5 0.0 0.0) to 1.3309933766567763 at #(0.0 0.0 0.0) Band 5 range: 1.3333608767284046 at #(0.0 0.0 0.0) to 1.565793736926709 at #(0.5 0.0 0.0) Band 6 range: 1.7594893273125134 at #(0.5 0.0 0.0) to 1.9916198017301796 at #(0.0 0.0 0.0) Band 7 range: 2.000243026425641 at #(0.0 0.0 0.0) to 2.2424804310085236 at #(0.5 0.0 0.0) Band 8 range: 2.406311900435369 at #(0.5 0.0 0.0) to 2.645093785208238 at #(0.0 0.0 0.0) Gap from band 1 (0.22245099319196915) to band 2 (0.44347071815372147), 66.38009279352583% Gap from band 2 (0.6663842825289141) to band 3 (0.6666675180317079), 0.0424942980722321% Gap from band 3 (0.8923678144024804) to band 4 (1.104558052488972), 21.25168906913916% Gap from band 4 (1.3309933766567763) to band 5 (1.3333608767284046), 0.17771661321840537% Gap from band 5 (1.565793736926709) to band 6 (1.7594893273125134), 11.649870801607634% Gap from band 6 (1.9916198017301796) to band 7 (2.000243026425641), 0.43204013097039173% Gap from band 7 (2.2424804310085236) to band 8 (2.406311900435369), 7.048345365685982% total elapsed time for run: 0 seconds. done. check-freqs: PASSED init-params: initializing eigensolver data Computing 8 bands with 1.000000e-13 tolerance. Working in 1 dimensions. Grid size is 32 x 1 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 1e+20, height 0.25, axis (1, 0, 0) epsilon = 9, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... 6 k-points: (0,0,0) (0.1,0,0) (0.2,0,0) (0.3,0,0) (0.4,0,0) (0.5,0,0) Solving for band polarization: te. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-9, mean 3, harm. mean 1.32353, 28.125% > 1, 25% "fill" Outputting check-epsilon... solve_kpoint (0,0,0): tefreqs:, k index, k1, k2, k3, kmag/2pi, te band 1, te band 2, te band 3, te band 4, te band 5, te band 6, te band 7, te band 8 Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 31 iterations. tefreqs:, 1, 0, 0, 0, 0, 0, 0.666384, 0.666668, 1.33099, 1.33336, 1.99162, 2.00024, 2.64509 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 74 iterations. tefreqs:, 2, 0.1, 0, 0, 0.1, 0.0574931, 0.608788, 0.724374, 1.2736, 1.39098, 1.93585, 2.05634, 2.59377 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tefreqs:, 3, 0.2, 0, 0, 0.2, 0.113352, 0.552761, 0.780756, 1.21673, 1.44856, 1.87756, 2.11569, 2.53392 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 11 iterations. tefreqs:, 4, 0.3, 0, 0, 0.3, 0.164802, 0.501202, 0.832987, 1.16418, 1.5025, 1.82314, 2.17225, 2.4769 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 11 iterations. tefreqs:, 5, 0.4, 0, 0, 0.4, 0.205537, 0.460405, 0.874781, 1.1222, 1.54664, 1.77873, 2.22033, 2.42855 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tefreqs:, 6, 0.5, 0, 0, 0.5, 0.222451, 0.443471, 0.892368, 1.10456, 1.56579, 1.75949, 2.24248, 2.40631 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.2224509931919693 at #(0.5 0.0 0.0) Band 2 range: 0.44347071815372185 at #(0.5 0.0 0.0) to 0.6663842825289168 at #(0.0 0.0 0.0) Band 3 range: 0.666667518031825 at #(0.0 0.0 0.0) to 0.892367814402481 at #(0.5 0.0 0.0) Band 4 range: 1.1045580524889718 at #(0.5 0.0 0.0) to 1.3309933766567714 at #(0.0 0.0 0.0) Band 5 range: 1.333360876728407 at #(0.0 0.0 0.0) to 1.5657937369267116 at #(0.5 0.0 0.0) Band 6 range: 1.7594893273125127 at #(0.5 0.0 0.0) to 1.991619801730292 at #(0.0 0.0 0.0) Band 7 range: 2.0002430264256734 at #(0.0 0.0 0.0) to 2.242480431008526 at #(0.5 0.0 0.0) Band 8 range: 2.406311900435368 at #(0.5 0.0 0.0) to 2.6450937852081897 at #(0.0 0.0 0.0) Gap from band 1 (0.2224509931919693) to band 2 (0.44347071815372185), 66.38009279352585% Gap from band 2 (0.6663842825289168) to band 3 (0.666667518031825), 0.04249429808940151% Gap from band 3 (0.892367814402481) to band 4 (1.1045580524889718), 21.25168906913907% Gap from band 4 (1.3309933766567714) to band 5 (1.333360876728407), 0.17771661321895554% Gap from band 5 (1.5657937369267116) to band 6 (1.7594893273125127), 11.649870801607427% Gap from band 6 (1.991619801730292) to band 7 (2.0002430264256734), 0.4320401309663711% Gap from band 7 (2.242480431008526) to band 8 (2.406311900435368), 7.048345365685846% total elapsed time for run: 0 seconds. done. check-freqs: PASSED ************************************************************************** Test case: 1d sinusoidal Bragg mirrors. ************************************************************************** init-params: initializing eigensolver data Computing 8 bands with 1.000000e-13 tolerance. Working in 1 dimensions. Grid size is 32 x 1 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 1e+20, height 0.25, axis (1, 0, 0) epsilon = 9, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... 11 k-points: (0,0,0) (0.05,0,0) (0.1,0,0) (0.15,0,0) (0.2,0,0) (0.25,0,0) (0.3,0,0) (0.35,0,0) (0.4,0,0) (0.45,0,0) (0.5,0,0) Solving for band polarization: tm. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1.00285-9, mean 4.66295, harm. mean 2.61927, 100% > 1, 45.7676% "fill" Outputting check-epsilon... solve_kpoint (0,0,0): tmfreqs:, k index, k1, k2, k3, kmag/2pi, tm band 1, tm band 2, tm band 3, tm band 4, tm band 5, tm band 6, tm band 7, tm band 8 Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 35 iterations. tmfreqs:, 1, 0, 0, 0, 0, 0, 0.460648, 0.542428, 0.968587, 1.01617, 1.48337, 1.48387, 1.96764 elapsed time for k point: 0 seconds. solve_kpoint (0.05,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 126 iterations. tmfreqs:, 2, 0.05, 0, 0, 0.05, 0.0231425, 0.454294, 0.548894, 0.95836, 1.02642, 1.45914, 1.50811, 1.94949 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tmfreqs:, 3, 0.1, 0, 0, 0.1, 0.0462091, 0.439084, 0.564453, 0.938267, 1.04658, 1.43467, 1.53261, 1.92548 elapsed time for k point: 0 seconds. solve_kpoint (0.15,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 4, 0.15, 0, 0, 0.15, 0.0691102, 0.420016, 0.584145, 0.91575, 1.06922, 1.41023, 1.55711, 1.9011 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 5, 0.2, 0, 0, 0.2, 0.0917239, 0.399481, 0.605653, 0.892484, 1.09267, 1.38581, 1.58161, 1.87665 elapsed time for k point: 0 seconds. solve_kpoint (0.25,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 6, 0.25, 0, 0, 0.25, 0.113863, 0.378568, 0.628026, 0.868958, 1.11645, 1.36143, 1.6061, 1.85216 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 7, 0.3, 0, 0, 0.3, 0.135213, 0.357979, 0.650805, 0.845399, 1.14036, 1.33713, 1.63058, 1.82769 elapsed time for k point: 0 seconds. solve_kpoint (0.35,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 8, 0.35, 0, 0, 0.35, 0.155194, 0.338479, 0.673671, 0.822003, 1.16428, 1.31295, 1.65502, 1.80325 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 9, 0.4, 0, 0, 0.4, 0.17268, 0.32129, 0.696193, 0.799138, 1.188, 1.28906, 1.67935, 1.77891 elapsed time for k point: 0 seconds. solve_kpoint (0.45,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 8 iterations. tmfreqs:, 10, 0.45, 0, 0, 0.45, 0.185503, 0.308628, 0.717034, 0.778102, 1.21076, 1.26621, 1.70326, 1.755 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 8 iterations. tmfreqs:, 11, 0.5, 0, 0, 0.5, 0.190416, 0.303766, 0.72859, 0.766483, 1.2253, 1.25164, 1.721, 1.73726 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.19041596916807002 at #(0.5 0.0 0.0) Band 2 range: 0.30376617480859147 at #(0.5 0.0 0.0) to 0.4606482752180713 at #(0.0 0.0 0.0) Band 3 range: 0.5424277685176458 at #(0.0 0.0 0.0) to 0.7285900510563766 at #(0.5 0.0 0.0) Band 4 range: 0.7664833049941233 at #(0.5 0.0 0.0) to 0.9685865863610064 at #(0.0 0.0 0.0) Band 5 range: 1.0161706920387945 at #(0.0 0.0 0.0) to 1.2253024188808161 at #(0.5 0.0 0.0) Band 6 range: 1.2516393338141771 at #(0.5 0.0 0.0) to 1.4833733379476344 at #(0.0 0.0 0.0) Band 7 range: 1.4838659375381698 at #(0.0 0.0 0.0) to 1.720998655926316 at #(0.5 0.0 0.0) Band 8 range: 1.737259239629598 at #(0.5 0.0 0.0) to 1.9676376863470306 at #(0.0 0.0 0.0) Gap from band 1 (0.19041596916807002) to band 2 (0.30376617480859147), 45.87385724963571% Gap from band 2 (0.4606482752180713) to band 3 (0.5424277685176458), 16.3057414859608% Gap from band 3 (0.7285900510563766) to band 4 (0.7664833049941233), 5.069082902774536% Gap from band 4 (0.9685865863610064) to band 5 (1.0161706920387945), 4.794954647164973% Gap from band 5 (1.2253024188808161) to band 6 (1.2516393338141771), 2.126567159256415% Gap from band 6 (1.4833733379476344) to band 7 (1.4838659375381698), 0.03320255259527266% Gap from band 7 (1.720998655926316) to band 8 (1.737259239629598), 0.940391618807724% total elapsed time for run: 0 seconds. done. check-freqs: PASSED ************************************************************************** Test case: get-dominant-planewave. ************************************************************************** init-params: initializing eigensolver data Computing 8 bands with 1.000000e-09 tolerance. Working in 1 dimensions. Grid size is 32 x 1 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: Geometric object tree has depth 1 and 0 object nodes (vs. 0 actual objects) Initializing epsilon function... 1 k-points: (0.4,0,0) Solving for band polarization: te. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-1, mean 1, harm. mean 1, 0% > 1, 100% "fill" Outputting check-epsilon... solve_kpoint (0.4,0,0): tefreqs:, k index, k1, k2, k3, kmag/2pi, te band 1, te band 2, te band 3, te band 4, te band 5, te band 6, te band 7, te band 8 Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tefreqs:, 1, 0.4, 0, 0, 0.4, 0.4, 0.6, 1.4, 1.6, 2.4, 2.6, 3.4, 3.6 elapsed time for k point: 0 seconds. total elapsed time for run: 0 seconds. done. check-almost-equal: PASSED ************************************************************************** Test case: Square lattice of dielectric rods in air. ************************************************************************** init-params: initializing eigensolver data Computing 8 bands with 1.000000e-09 tolerance. Working in 2 dimensions. Grid size is 32 x 32 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 11.56, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Allocating fields... 16 k-points: (0,0,0) (0.1,0,0) (0.2,0,0) (0.3,0,0) (0.4,0,0) (0.5,0,0) (0.5,0.1,0) (0.5,0.2,0) (0.5,0.3,0) (0.5,0.4,0) (0.5,0.5,0) (0.4,0.4,0) (0.3,0.3,0) (0.2,0.2,0) (0.1,0.1,0) (0,0,0) Solving for band polarization: te. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-11.56, mean 2.327, harm. mean 1.1441, 14.5508% > 1, 12.5663% "fill" Outputting check-epsilon... solve_kpoint (0,0,0): tefreqs:, k index, k1, k2, k3, kmag/2pi, te band 1, te band 2, te band 3, te band 4, te band 5, te band 6, te band 7, te band 8 Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 42 iterations. tefreqs:, 1, 0, 0, 0, 0, 0, 0.561945, 0.780843, 0.780846, 0.924372, 1.00804, 1.00804, 1.09859 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 19 iterations. tefreqs:, 2, 0.1, 0, 0, 0.1, 0.0897998, 0.56021, 0.767785, 0.782438, 0.912603, 1.00889, 1.00949, 1.12115 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 93 iterations. tefreqs:, 3, 0.2, 0, 0, 0.2, 0.178853, 0.553301, 0.732522, 0.786672, 0.890463, 1.01386, 1.02078, 1.11301 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 15 iterations. tefreqs:, 4, 0.3, 0, 0, 0.3, 0.266124, 0.534865, 0.689376, 0.79205, 0.872926, 1.0209, 1.04285, 1.10753 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 20 iterations. tefreqs:, 5, 0.4, 0, 0, 0.4, 0.349588, 0.49478, 0.658508, 0.796524, 0.862684, 1.02898, 1.07044, 1.10003 elapsed time for k point: 1 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 34 iterations. tefreqs:, 6, 0.5, 0, 0, 0.5, 0.413346, 0.444623, 0.648672, 0.798265, 0.859327, 1.03311, 1.0958, 1.09742 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.1,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 19 iterations. tefreqs:, 7, 0.5, 0.1, 0, 0.509902, 0.424298, 0.448983, 0.644549, 0.802551, 0.854853, 0.98984, 1.05831, 1.11679 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.2,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 15 iterations. tefreqs:, 8, 0.5, 0.2, 0, 0.538516, 0.455354, 0.461161, 0.63323, 0.814782, 0.834251, 0.934558, 1.01274, 1.12367 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.3,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 14 iterations. tefreqs:, 9, 0.5, 0.3, 0, 0.583095, 0.478427, 0.501912, 0.61745, 0.784489, 0.83371, 0.906541, 0.967734, 1.12636 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.4,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 14 iterations. tefreqs:, 10, 0.5, 0.4, 0, 0.640312, 0.495605, 0.556867, 0.601739, 0.720584, 0.858432, 0.897772, 0.926611, 1.12762 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.5,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tefreqs:, 11, 0.5, 0.5, 0, 0.707107, 0.503606, 0.594344, 0.594357, 0.679389, 0.883665, 0.895782, 0.895787, 1.12801 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0.4,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tefreqs:, 12, 0.4, 0.4, 0, 0.565685, 0.474671, 0.549917, 0.608032, 0.745035, 0.848355, 0.895854, 0.945891, 1.12708 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0.3,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 15 iterations. tefreqs:, 13, 0.3, 0.3, 0, 0.424264, 0.373254, 0.543225, 0.646251, 0.817389, 0.830457, 0.896516, 1.0176, 1.12198 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0.2,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 20 iterations. tefreqs:, 14, 0.2, 0.2, 0, 0.282843, 0.252369, 0.551073, 0.700804, 0.797552, 0.899586, 0.90339, 1.08791, 1.10099 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0.1,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 24 iterations. tefreqs:, 15, 0.1, 0.1, 0, 0.141421, 0.12694, 0.558854, 0.755651, 0.785047, 0.909884, 0.968813, 1.04814, 1.11332 elapsed time for k point: 0 seconds. solve_kpoint (0,0,0): Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 60 iterations. tefreqs:, 16, 0, 0, 0, 0, 0, 0.561945, 0.780843, 0.780846, 0.924372, 1.00804, 1.00804, 1.09859 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.5036058015317006 at #(0.5 0.5 0.0) Band 2 range: 0.4446229135155304 at #(0.5 0.0 0.0) to 0.5943440246130545 at #(0.5 0.5 0.0) Band 3 range: 0.5943566395317802 at #(0.5 0.5 0.0) to 0.7808428128477852 at #(0.0 0.0 0.0) Band 4 range: 0.6793887413186156 at #(0.5 0.5 0.0) to 0.8173893720817341 at #(0.30000000000000004 0.30000000000000004 0.0) Band 5 range: 0.8304573822717966 at #(0.30000000000000004 0.30000000000000004 0.0) to 0.9243716835680688 at #(0.0 0.0 0.0) Band 6 range: 0.8957817684697049 at #(0.5 0.5 0.0) to 1.0331104139500322 at #(0.5 0.0 0.0) Band 7 range: 0.8957868745561722 at #(0.5 0.5 0.0) to 1.0958021498670216 at #(0.5 0.0 0.0) Band 8 range: 1.097416809770438 at #(0.5 0.0 0.0) to 1.1280127657657506 at #(0.5 0.5 0.0) Gap from band 2 (0.5943440246130545) to band 3 (0.5943566395317802), 0.002122471889893968% Gap from band 4 (0.8173893720817341) to band 5 (0.8304573822717966), 1.586071053699306% Gap from band 7 (1.0958021498670216) to band 8 (1.097416809770438), 0.14724110388716175% total elapsed time for run: 1 seconds. done. check-freqs: PASSED init-params: initializing eigensolver data Computing 8 bands with 1.000000e-09 tolerance. Working in 2 dimensions. Grid size is 32 x 32 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 11.56, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... 16 k-points: (0,0,0) (0.1,0,0) (0.2,0,0) (0.3,0,0) (0.4,0,0) (0.5,0,0) (0.5,0.1,0) (0.5,0.2,0) (0.5,0.3,0) (0.5,0.4,0) (0.5,0.5,0) (0.4,0.4,0) (0.3,0.3,0) (0.2,0.2,0) (0.1,0.1,0) (0,0,0) Solving for band polarization: tm. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-11.56, mean 2.327, harm. mean 1.1441, 14.5508% > 1, 12.5663% "fill" Outputting check-epsilon... solve_kpoint (0,0,0): tmfreqs:, k index, k1, k2, k3, kmag/2pi, tm band 1, tm band 2, tm band 3, tm band 4, tm band 5, tm band 6, tm band 7, tm band 8 Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 25 iterations. tmfreqs:, 1, 0, 0, 0, 0, 0, 0.550336, 0.561338, 0.56134, 0.822948, 0.868842, 0.965325, 1.08938 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 21 iterations. tmfreqs:, 2, 0.1, 0, 0, 0.1, 0.0651416, 0.525004, 0.561885, 0.586307, 0.823536, 0.867346, 0.954632, 1.05688 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 3, 0.2, 0, 0, 0.2, 0.127665, 0.49365, 0.563323, 0.617312, 0.822737, 0.86353, 0.924306, 1.03882 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 4, 0.3, 0, 0, 0.3, 0.184047, 0.461592, 0.565122, 0.651397, 0.810839, 0.858984, 0.892949, 1.03589 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 5, 0.4, 0, 0, 0.4, 0.227778, 0.433361, 0.566596, 0.689045, 0.778792, 0.855425, 0.879398, 1.03874 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 6, 0.5, 0, 0, 0.5, 0.245809, 0.420657, 0.567163, 0.720092, 0.747203, 0.85409, 0.877012, 1.0408 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.1,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 8 iterations. tmfreqs:, 7, 0.5, 0.1, 0, 0.509902, 0.249299, 0.427308, 0.560221, 0.718025, 0.756122, 0.855013, 0.877106, 1.02917 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.2,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 8, 0.5, 0.2, 0, 0.538516, 0.258693, 0.445541, 0.543385, 0.711698, 0.77935, 0.858136, 0.877357, 1.00092 elapsed time for k point: 1 seconds. solve_kpoint (0.5,0.3,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 9, 0.5, 0.3, 0, 0.583095, 0.270963, 0.470322, 0.524053, 0.701396, 0.810963, 0.864394, 0.877683, 0.965108 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.4,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 10, 0.5, 0.4, 0, 0.640312, 0.281614, 0.493207, 0.508823, 0.689961, 0.846785, 0.874595, 0.878068, 0.926581 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.5,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 12 iterations. tmfreqs:, 11, 0.5, 0.5, 0, 0.707107, 0.285906, 0.502981, 0.502983, 0.684476, 0.874359, 0.883317, 0.883317, 0.892993 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0.4,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tmfreqs:, 12, 0.4, 0.4, 0, 0.565685, 0.276089, 0.491352, 0.508683, 0.69258, 0.839723, 0.856431, 0.907219, 0.907347 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0.3,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 13, 0.3, 0.3, 0, 0.424264, 0.240239, 0.47923, 0.523499, 0.685362, 0.829266, 0.840451, 0.910791, 0.941595 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0.2,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 14, 0.2, 0.2, 0, 0.282843, 0.17524, 0.488533, 0.541561, 0.64751, 0.83002, 0.850443, 0.922657, 0.983925 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0.1,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 15, 0.1, 0.1, 0, 0.141421, 0.0915259, 0.516393, 0.555924, 0.601211, 0.824561, 0.86529, 0.948412, 1.03526 elapsed time for k point: 0 seconds. solve_kpoint (0,0,0): Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 22 iterations. tmfreqs:, 16, 0, 0, 0, 0, 0, 0.550336, 0.561338, 0.56134, 0.822948, 0.868842, 0.965325, 1.08938 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.2859057791271609 at #(0.5 0.5 0.0) Band 2 range: 0.4206573384061859 at #(0.5 0.0 0.0) to 0.5503360754972212 at #(0.0 0.0 0.0) Band 3 range: 0.5029830978387367 at #(0.5 0.5 0.0) to 0.5671632878212804 at #(0.5 0.0 0.0) Band 4 range: 0.5613397939938073 at #(0.0 0.0 0.0) to 0.7200918204690926 at #(0.5 0.0 0.0) Band 5 range: 0.7472029910634781 at #(0.5 0.0 0.0) to 0.8743593805271221 at #(0.5 0.5 0.0) Band 6 range: 0.8404509697715068 at #(0.30000000000000004 0.30000000000000004 0.0) to 0.8833173725850519 at #(0.5 0.5 0.0) Band 7 range: 0.877011871859037 at #(0.5 0.0 0.0) to 0.965325380977264 at #(0.0 0.0 0.0) Band 8 range: 0.8929933495601443 at #(0.5 0.5 0.0) to 1.0893779022325405 at #(0.0 0.0 0.0) Gap from band 1 (0.2859057791271609) to band 2 (0.4206573384061859), 38.14282289442182% Gap from band 4 (0.7200918204690926) to band 5 (0.7472029910634781), 3.6953951423120253% total elapsed time for run: 1 seconds. done. check-freqs: PASSED **************************************************************************** Test case: square lattice of magneto-electric rods in air. **************************************************************************** init-params: initializing eigensolver data Computing 10 bands with 1.000000e-09 tolerance. Working in 2 dimensions. Grid size is 32 x 32 x 1. Solving for 10 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 0.11, height 1e+20, axis (0, 0, 1) Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... Allocating fields... 10 k-points: (0,0,0) (0.166667,0,0) (0.333333,0,0) (0.5,0,0) (0.5,0.166667,0) (0.5,0.333333,0) (0.5,0.5,0) (0.333333,0.333333,0) (0.166667,0.166667,0) (0,0,0) Solving for band polarization: tm. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-15, mean 1.53164, harm. mean 1.04189, 4.39453% > 1, 3.79742% "fill" Outputting check-epsilon... mu: 1-13.9136, mean 1.42353, harm. mean 1.03906, 4.39453% > 1, 3.27973% "fill" Outputting check-mu... solve_kpoint (0,0,0): tmfreqs:, k index, k1, k2, k3, kmag/2pi, tm band 1, tm band 2, tm band 3, tm band 4, tm band 5, tm band 6, tm band 7, tm band 8, tm band 9, tm band 10 Solving for bands 2 to 10... iteration 66: trace = 3.332654113437661 (0.373881% change) iteration 137: trace = 2.922869083991767 (0.077839% change) iteration 212: trace = 2.621670309258063 (0.271461% change) iteration 265: trace = 2.452275916605919 (0.111455% change) iteration 335: trace = 2.413533242590518 (0.00345526% change) iteration 400: trace = 2.411709252043468 (0.000386865% change) iteration 473: trace = 2.410405775002738 (0.00150795% change) iteration 539: trace = 2.410031462228026 (0.000344194% change) iteration 611: trace = 2.409798474370296 (6.74858e-05% change) iteration 672: trace = 2.409687426571257 (1.04426e-05% change) iteration 717: trace = 2.40958267847257 (0.000190448% change) iteration 776: trace = 2.409450048283543 (5.05514e-05% change) iteration 838: trace = 2.40925473080393 (9.56854e-05% change) iteration 895: trace = 2.408991357730128 (0.000462023% change) iteration 955: trace = 2.408611050961347 (0.000184701% change) iteration 1014: trace = 2.407035592955471 (0.00105454% change) iteration 1089: trace = 2.40543420112172 (0.000546933% change) iteration 1163: trace = 2.401918608364619 (0.000705386% change) iteration 1234: trace = 2.399633565733665 (0.000746489% change) iteration 1304: trace = 2.395658675540483 (0.000618623% change) iteration 1377: trace = 2.39464206184472 (0.000709406% change) iteration 1434: trace = 2.393865823149387 (0.000672246% change) iteration 1503: trace = 2.393436581058496 (0.00024049% change) iteration 1575: trace = 2.393204239272516 (0.000152796% change) iteration 1648: trace = 2.393117989281403 (6.23159e-05% change) iteration 1707: trace = 2.393081549246357 (2.23961e-05% change) iteration 1767: trace = 2.393062516420064 (1.42036e-05% change) iteration 1817: trace = 2.393053162681409 (5.93485e-06% change) iteration 1878: trace = 2.393047016634214 (2.96446e-06% change) iteration 1934: trace = 2.393042321934401 (1.30026e-06% change) iteration 1992: trace = 2.393041049462053 (1.96343e-06% change) Finished solving for bands 2 to 10 after 2043 iterations. tmfreqs:, 1, 0, 0, 0, 0, 0, 0.265885, 0.356852, 0.36899, 0.503897, 0.506552, 0.539911, 0.635681, 0.645852, 0.660017 elapsed time for k point: 16 seconds. solve_kpoint (0.166667,0,0): Solving for bands 1 to 10... iteration 63: trace = 2.413664630369193 (0.00193653% change) iteration 127: trace = 2.412764087072741 (0.000141728% change) iteration 186: trace = 2.412679012174157 (2.72984e-05% change) iteration 251: trace = 2.412659193531059 (1.28657e-05% change) iteration 303: trace = 2.412646622748995 (4.38531e-06% change) iteration 366: trace = 2.412640644654746 (4.19394e-06% change) iteration 431: trace = 2.412636295764047 (2.06358e-06% change) iteration 497: trace = 2.412633820852816 (1.56746e-06% change) iteration 566: trace = 2.412632152069356 (3.05584e-07% change) Finished solving for bands 1 to 10 after 581 iterations. tmfreqs:, 2, 0.166667, 0, 0, 0.166667, 0.123306, 0.278292, 0.357366, 0.369461, 0.503797, 0.506642, 0.538217, 0.63557, 0.645595, 0.659432 elapsed time for k point: 5 seconds. solve_kpoint (0.333333,0,0): Solving for bands 1 to 10... iteration 55: trace = 2.455292826574757 (0.0016247% change) iteration 112: trace = 2.453461871283896 (0.000942013% change) iteration 138: trace = 2.452850294093757 (0.000516431% change) iteration 157: trace = 2.452594543112414 (0.000190195% change) iteration 212: trace = 2.452261647114671 (0.000166614% change) iteration 226: trace = 2.45222218084879 (0.000271425% change) iteration 255: trace = 2.452111509759052 (0.00010316% change) iteration 274: trace = 2.452043731962384 (4.89594e-05% change) iteration 338: trace = 2.451920855778976 (1.40795e-05% change) iteration 368: trace = 2.451893675155902 (5.3292e-05% change) iteration 388: trace = 2.451877011963294 (2.27752e-05% change) iteration 417: trace = 2.451851890520292 (7.3478e-06% change) iteration 441: trace = 2.451837518829035 (9.8259e-06% change) iteration 488: trace = 2.451821380847629 (2.13056e-05% change) iteration 525: trace = 2.451813300926426 (1.29873e-05% change) iteration 561: trace = 2.451804588379279 (3.06408e-06% change) iteration 611: trace = 2.45179865320105 (2.22218e-06% change) iteration 642: trace = 2.451796997197523 (8.22103e-07% change) iteration 670: trace = 2.451795163210634 (4.48774e-06% change) iteration 737: trace = 2.451792191422107 (1.35098e-06% change) iteration 799: trace = 2.451791130767661 (1.20237e-07% change) iteration 860: trace = 2.451790293354878 (4.86966e-07% change) Finished solving for bands 1 to 10 after 864 iterations. tmfreqs:, 3, 0.333333, 0, 0, 0.333333, 0.194655, 0.331458, 0.361679, 0.377947, 0.502693, 0.506755, 0.527227, 0.633335, 0.644215, 0.653287 elapsed time for k point: 11 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 10... iteration 57: trace = 2.455408918283909 (0.00479795% change) iteration 82: trace = 2.452103611846625 (0.00276672% change) iteration 121: trace = 2.450566888260983 (0.000296549% change) iteration 163: trace = 2.449770134987228 (0.000114968% change) iteration 204: trace = 2.449533444934879 (0.000215786% change) iteration 220: trace = 2.449426206730259 (6.73248e-05% change) iteration 269: trace = 2.449367617211575 (2.93637e-05% change) iteration 311: trace = 2.449343425249795 (2.08503e-05% change) iteration 362: trace = 2.449316336941978 (9.9765e-06% change) iteration 426: trace = 2.449305638980764 (1.50805e-06% change) iteration 492: trace = 2.449298594858241 (7.61886e-06% change) iteration 561: trace = 2.449295688781586 (4.0645e-07% change) iteration 626: trace = 2.449294268379298 (1.38997e-06% change) Finished solving for bands 1 to 10 after 661 iterations. tmfreqs:, 4, 0.5, 0, 0, 0.5, 0.208799, 0.348978, 0.365278, 0.443856, 0.477901, 0.506845, 0.509076, 0.617527, 0.63971, 0.646891 elapsed time for k point: 8 seconds. solve_kpoint (0.5,0.166667,0): Solving for bands 1 to 10... iteration 44: trace = 2.49096335243478 (0.0018121% change) iteration 97: trace = 2.487536804185604 (0.00237053% change) iteration 159: trace = 2.486231868717033 (0.000391527% change) iteration 228: trace = 2.485995305196467 (0.000103246% change) iteration 300: trace = 2.485901827753271 (1.65081e-05% change) iteration 363: trace = 2.485880557678704 (3.04459e-06% change) iteration 415: trace = 2.485871623123662 (8.95371e-06% change) iteration 478: trace = 2.485867638277805 (3.84048e-07% change) iteration 537: trace = 2.485866239230515 (3.9297e-07% change) Finished solving for bands 1 to 10 after 596 iterations. tmfreqs:, 5, 0.5, 0.166667, 0, 0.527046, 0.210673, 0.351103, 0.364915, 0.453355, 0.500502, 0.507018, 0.510125, 0.616821, 0.642478, 0.646894 elapsed time for k point: 5 seconds. solve_kpoint (0.5,0.333333,0): Solving for bands 1 to 10... iteration 69: trace = 2.582909536690412 (0.00344747% change) iteration 136: trace = 2.579501830383394 (0.000797224% change) iteration 199: trace = 2.578998167676302 (0.00013433% change) iteration 253: trace = 2.578914819807851 (3.50211e-05% change) iteration 309: trace = 2.578887808757357 (9.83849e-06% change) iteration 364: trace = 2.578883446709924 (1.15118e-06% change) iteration 433: trace = 2.57888190545302 (8.33347e-08% change) Finished solving for bands 1 to 10 after 433 iterations. tmfreqs:, 6, 0.5, 0.333333, 0, 0.600925, 0.214081, 0.353217, 0.365001, 0.471977, 0.507165, 0.51079, 0.563894, 0.615052, 0.646178, 0.648357 elapsed time for k point: 3 seconds. solve_kpoint (0.5,0.5,0): Solving for bands 1 to 10... iteration 59: trace = 2.660613872519377 (0.00995321% change) iteration 128: trace = 2.654135912253834 (0.000379581% change) iteration 200: trace = 2.653873703305905 (2.19124e-05% change) iteration 271: trace = 2.65386168176228 (4.19101e-07% change) Finished solving for bands 1 to 10 after 308 iterations. tmfreqs:, 7, 0.5, 0.5, 0, 0.707107, 0.215649, 0.353602, 0.365269, 0.480748, 0.507335, 0.512534, 0.60792, 0.619784, 0.648908, 0.650323 elapsed time for k point: 3 seconds. solve_kpoint (0.333333,0.333333,0): Solving for bands 1 to 10... iteration 57: trace = 2.534381609709547 (0.00247526% change) iteration 128: trace = 2.532451791307959 (6.94251e-05% change) iteration 199: trace = 2.532257044589847 (2.27282e-05% change) iteration 271: trace = 2.532229995337595 (5.01968e-06% change) iteration 343: trace = 2.532227108793403 (1.14039e-06% change) Finished solving for bands 1 to 10 after 371 iterations. tmfreqs:, 8, 0.333333, 0.333333, 0, 0.471405, 0.210879, 0.35463, 0.361245, 0.439486, 0.507022, 0.508252, 0.543815, 0.620599, 0.646589, 0.651154 elapsed time for k point: 2 seconds. solve_kpoint (0.166667,0.166667,0): Solving for bands 1 to 10... iteration 69: trace = 2.430575292235021 (0.00249951% change) iteration 107: trace = 2.429540484813175 (0.00108424% change) iteration 178: trace = 2.42921900002811 (4.9784e-05% change) iteration 250: trace = 2.429176559079801 (6.62091e-06% change) iteration 321: trace = 2.429164711786565 (3.83483e-06% change) iteration 389: trace = 2.429161384946178 (1.69304e-06% change) Finished solving for bands 1 to 10 after 443 iterations. tmfreqs:, 9, 0.166667, 0.166667, 0, 0.235702, 0.162353, 0.295875, 0.356214, 0.372161, 0.504258, 0.506735, 0.537355, 0.633231, 0.645916, 0.65717 elapsed time for k point: 4 seconds. solve_kpoint (0,0,0): Solving for bands 2 to 10... iteration 75: trace = 2.393309825162555 (0.000587191% change) iteration 141: trace = 2.393048997522238 (4.72568e-05% change) iteration 161: trace = 2.3930455381777 (2.03921e-05% change) iteration 232: trace = 2.393040971215491 (7.58165e-07% change) Finished solving for bands 2 to 10 after 277 iterations. tmfreqs:, 10, 0, 0, 0, 0, 0, 0.265885, 0.356852, 0.36899, 0.503897, 0.506552, 0.539911, 0.635681, 0.645852, 0.660016 elapsed time for k point: 2 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.2156485348682738 at #(0.5 0.5 0.0) Band 2 range: 0.265884981507677 at #(0.0 0.0 0.0) to 0.354630241126045 at #(0.3333333333333333 0.3333333333333333 0.0) Band 3 range: 0.3562142390294759 at #(0.16666666666666666 0.16666666666666666 0.0) to 0.3652783377830165 at #(0.5 0.0 0.0) Band 4 range: 0.36899013867629754 at #(0.0 0.0 0.0) to 0.48074847340165067 at #(0.5 0.5 0.0) Band 5 range: 0.47790149106730545 at #(0.5 0.0 0.0) to 0.5073352932597683 at #(0.5 0.5 0.0) Band 6 range: 0.5065518607417917 at #(0.0 0.0 0.0) to 0.5125343708174586 at #(0.5 0.5 0.0) Band 7 range: 0.5090758750097967 at #(0.5 0.0 0.0) to 0.6079200973337775 at #(0.5 0.5 0.0) Band 8 range: 0.6150518914912332 at #(0.5 0.33333333333333337 0.0) to 0.6356807918573262 at #(0.0 0.0 0.0) Band 9 range: 0.6397097808807445 at #(0.5 0.0 0.0) to 0.6489079032967237 at #(0.5 0.5 0.0) Band 10 range: 0.6468905461913745 at #(0.5 0.0 0.0) to 0.6600167739023248 at #(0.0 0.0 0.0) Gap from band 1 (0.2156485348682738) to band 2 (0.265884981507677), 20.865192112684326% Gap from band 2 (0.354630241126045) to band 3 (0.3562142390294759), 0.4456665129014822% Gap from band 3 (0.3652783377830165) to band 4 (0.36899013867629754), 1.0110200865981802% Gap from band 7 (0.6079200973337775) to band 8 (0.6150518914912332), 1.1663053974453932% Gap from band 8 (0.6356807918573262) to band 9 (0.6397097808807445), 0.6318047364531911% total elapsed time for run: 59 seconds. done. check-freqs: PASSED ************************************************************************** Test case: 3x3 triangular lattice of rods in air, dipole defect states. ************************************************************************** init-params: initializing eigensolver data Computing 2 bands with 1.000000e-09 tolerance. Target frequency is 0.35 Working in 2 dimensions. Grid size is 80 x 80 x 1. Solving for 2 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (2.59808, 1.5, 0) (2.59808, -1.5, 0) (0, 0, 1) Cell volume = 7.79423 Reciprocal lattice vectors (/ 2 pi): (0.19245, 0.333333, -0) (0.19245, -0.333333, 0) (-0, 0, 1) Geometric objects: cylinder, center = (1,1,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 12, mu = 1 cylinder, center = (0,1,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 12, mu = 1 cylinder, center = (-1,1,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 12, mu = 1 cylinder, center = (1,0,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 12, mu = 1 cylinder, center = (0,0,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 12, mu = 1 cylinder, center = (-1,0,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 12, mu = 1 cylinder, center = (1,-1,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 12, mu = 1 cylinder, center = (0,-1,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 12, mu = 1 cylinder, center = (-1,-1,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 12, mu = 1 cylinder, center = (0,0,0) radius 0.33, height 1e+20, axis (0, 0, 1) epsilon = 12, mu = 1 Geometric object tree has depth 4 and 11 object nodes (vs. 10 actual objects) Initializing epsilon function... Allocating fields... 1 k-points: (0,0.5,0) Solving for band polarization: tm. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-12, mean 2.90162, harm. mean 1.21599, 20.9844% > 1, 17.2875% "fill" Outputting check-epsilon... solve_kpoint (0,0.5,0): tmfreqs:, k index, k1, k2, k3, kmag/2pi, tm band 1, tm band 2 Solving for bands 1 to 2... iteration 59: trace = 0.007854138053471868 (0.658942% change) iteration 117: trace = 0.002855014262991655 (3.47332% change) iteration 176: trace = 0.0001815753140648246 (1.02747% change) iteration 214: trace = 0.0001491753155659719 (0.0129766% change) iteration 268: trace = 0.0001480701578613175 (8.0246e-06% change) Finished solving for bands 1 to 2 after 307 iterations. tmfreqs:, 1, 0, 0.5, 0, 0.19245, 0.33627, 0.338821 elapsed time for k point: 3 seconds. total elapsed time for run: 3 seconds. done. check-freqs: PASSED ************************************************************************** Test case: fcc lattice of air spheres in dielectric. ************************************************************************** init-params: initializing eigensolver data Computing 10 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 10 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (0, 0.707107, 0.707107) (0.707107, 0, 0.707107) (0.707107, 0.707107, 0) Cell volume = 0.707107 Reciprocal lattice vectors (/ 2 pi): (-0.707107, 0.707107, 0.707107) (0.707107, -0.707107, 0.707107) (0.707107, 0.707107, -0.707107) Geometric objects: sphere, center = (0,0,0) radius 0.5 epsilon = 1, mu = 1 Geometric object tree has depth 4 and 64 object nodes (vs. 1 actual objects) Initializing epsilon function... Allocating fields... 13 k-points: (0,0.5,0.5) (0,0.5625,0.4375) (0,0.625,0.375) (0,0.5625,0.1875) (0,0.5,0) (0,0.25,0) (0,0,0) (0,0.25,0.25) (0,0.5,0.5) (0.125,0.625,0.5) (0.25,0.75,0.5) (0.3125,0.75,0.4375) (0.375,0.75,0.375) Solving for band polarization: . Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-11.56, mean 3.72661, harm. mean 1.46352, 45.5811% > 1, 25.8202% "fill" Outputting check-epsilon... solve_kpoint (0,0.5,0.5): freqs:, k index, k1, k2, k3, kmag/2pi, band 1, band 2, band 3, band 4, band 5, band 6, band 7, band 8, band 9, band 10 Solving for bands 1 to 10... iteration 33: trace = 2.417928316057217 (1.68003e-05% change) Finished solving for bands 1 to 10 after 47 iterations. freqs:, 1, 0, 0.5, 0.5, 0.707107, 0.370307, 0.372085, 0.38252, 0.3836, 0.493989, 0.512361, 0.522307, 0.523822, 0.595103, 0.664079 elapsed time for k point: 1 seconds. solve_kpoint (0,0.5625,0.4375): Solving for bands 1 to 10... Finished solving for bands 1 to 10 after 31 iterations. freqs:, 2, 0, 0.5625, 0.4375, 0.71807, 0.368359, 0.377487, 0.384992, 0.387438, 0.471895, 0.506546, 0.523725, 0.531091, 0.610293, 0.649714 elapsed time for k point: 0 seconds. solve_kpoint (0,0.625,0.375): Solving for bands 1 to 10... iteration 23: trace = 2.416823960556714 (8.65146e-08% change) Finished solving for bands 1 to 10 after 23 iterations. freqs:, 3, 0, 0.625, 0.375, 0.75, 0.357533, 0.38095, 0.392789, 0.401454, 0.43935, 0.495158, 0.527478, 0.540892, 0.635515, 0.644621 elapsed time for k point: 1 seconds. solve_kpoint (0,0.5625,0.1875): Solving for bands 1 to 10... iteration 25: trace = 2.408973249392727 (4.40312e-06% change) Finished solving for bands 1 to 10 after 31 iterations. freqs:, 4, 0, 0.5625, 0.1875, 0.649519, 0.323742, 0.331653, 0.39821, 0.401206, 0.462868, 0.514292, 0.533282, 0.546116, 0.631046, 0.644236 elapsed time for k point: 1 seconds. solve_kpoint (0,0.5,0): Solving for bands 1 to 10... Finished solving for bands 1 to 10 after 18 iterations. freqs:, 5, 0, 0.5, 0, 0.612372, 0.307529, 0.308891, 0.387977, 0.389913, 0.491606, 0.536329, 0.536613, 0.538985, 0.624681, 0.629124 elapsed time for k point: 0 seconds. solve_kpoint (0,0.25,0): Solving for bands 1 to 10... Finished solving for bands 1 to 10 after 27 iterations. freqs:, 6, 0, 0.25, 0, 0.306186, 0.179755, 0.180421, 0.473402, 0.475968, 0.50346, 0.53577, 0.538144, 0.540345, 0.624073, 0.627676 elapsed time for k point: 0 seconds. solve_kpoint (0,0,0): Solving for bands 3 to 10... iteration 32: trace = 2.454760663552489 (5.02919e-06% change) Finished solving for bands 3 to 10 after 39 iterations. freqs:, 7, 0, 0, 0, 0, 0, 0, 0.518283, 0.521124, 0.521132, 0.543942, 0.543952, 0.546505, 0.61287, 0.61397 elapsed time for k point: 1 seconds. solve_kpoint (0,0.25,0.25): Solving for bands 1 to 10... iteration 26: trace = 2.428538456531416 (0.00426264% change) iteration 51: trace = 2.428131957456567 (1.932e-06% change) Finished solving for bands 1 to 10 after 57 iterations. freqs:, 8, 0, 0.25, 0.25, 0.353553, 0.207258, 0.207841, 0.473224, 0.475268, 0.506873, 0.526874, 0.530366, 0.532023, 0.604689, 0.653955 elapsed time for k point: 1 seconds. solve_kpoint (0,0.5,0.5): Solving for bands 1 to 10... iteration 25: trace = 2.418948340062556 (0.00984979% change) iteration 54: trace = 2.417940224603838 (3.48511e-05% change) iteration 75: trace = 2.417934601164372 (6.92443e-07% change) Finished solving for bands 1 to 10 after 82 iterations. freqs:, 9, 0, 0.5, 0.5, 0.707107, 0.370307, 0.372085, 0.38252, 0.3836, 0.493989, 0.512361, 0.522307, 0.523822, 0.595103, 0.664085 elapsed time for k point: 2 seconds. solve_kpoint (0.125,0.625,0.5): Solving for bands 1 to 10... Finished solving for bands 1 to 10 after 27 iterations. freqs:, 10, 0.125, 0.625, 0.5, 0.728869, 0.372309, 0.377154, 0.385968, 0.392562, 0.463569, 0.502026, 0.505399, 0.549256, 0.621502, 0.639203 elapsed time for k point: 0 seconds. solve_kpoint (0.25,0.75,0.5): Solving for bands 1 to 10... iteration 24: trace = 2.422730091067978 (0.00066389% change) iteration 38: trace = 2.42269331624072 (4.36128e-06% change) Finished solving for bands 1 to 10 after 45 iterations. freqs:, 11, 0.25, 0.75, 0.5, 0.790569, 0.373709, 0.386404, 0.386499, 0.410208, 0.435229, 0.490809, 0.492518, 0.566583, 0.627971, 0.654078 elapsed time for k point: 2 seconds. solve_kpoint (0.3125,0.75,0.4375): Solving for bands 1 to 10... iteration 32: trace = 2.419550704426088 (6.43583e-06% change) Finished solving for bands 1 to 10 after 40 iterations. freqs:, 12, 0.3125, 0.75, 0.4375, 0.760345, 0.364261, 0.382754, 0.390011, 0.405967, 0.438258, 0.493602, 0.505736, 0.559219, 0.634896, 0.645093 elapsed time for k point: 0 seconds. solve_kpoint (0.375,0.75,0.375): Solving for bands 1 to 10... Finished solving for bands 1 to 10 after 25 iterations. freqs:, 13, 0.375, 0.75, 0.375, 0.75, 0.359122, 0.380258, 0.391276, 0.402681, 0.438005, 0.49516, 0.529122, 0.541679, 0.637583, 0.645624 elapsed time for k point: 1 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.37370938045583396 at #(0.25 0.75 0.5) Band 2 range: 0.0 at #(0.0 0.0 0.0) to 0.38640412149814385 at #(0.25 0.75 0.5) Band 3 range: 0.382520138319054 at #(0.0 0.5 0.5) to 0.5182834488052886 at #(0.0 0.0 0.0) Band 4 range: 0.3836003495578335 at #(0.0 0.5 0.5) to 0.5211243460581921 at #(0.0 0.0 0.0) Band 5 range: 0.4352292385396697 at #(0.25 0.75 0.5) to 0.5211324741392463 at #(0.0 0.0 0.0) Band 6 range: 0.4908086394112833 at #(0.25 0.75 0.5) to 0.5439419096018505 at #(0.0 0.0 0.0) Band 7 range: 0.492518348487665 at #(0.25 0.75 0.5) to 0.5439517867826291 at #(0.0 0.0 0.0) Band 8 range: 0.5238217225058702 at #(0.0 0.5 0.5) to 0.5665833313384165 at #(0.25 0.75 0.5) Band 9 range: 0.5951030360714451 at #(0.0 0.5 0.5) to 0.6375830839471329 at #(0.375 0.75 0.375) Band 10 range: 0.6139704161698623 at #(0.0 0.0 0.0) to 0.6640850263391456 at #(0.0 0.5 0.5) Gap from band 8 (0.5665833313384165) to band 9 (0.5951030360714451), 4.910052408829956% total elapsed time for run: 10 seconds. done. check-freqs: PASSED ************************************************************************** Test case: simple cubic lattice with anisotropic dielectric. ************************************************************************** init-params: initializing eigensolver data Computing 3 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 3 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: block, center = (0,0,0) size (0.313,0.313,1) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (0.313,1,0.313) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (1,0.313,0.313) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (0.313,0.313,0.313) axes (1,0,0), (0,1,0), (0,0,1) epsilon = 12, mu = 1 Geometric object tree has depth 7 and 16 object nodes (vs. 4 actual objects) Initializing epsilon function... Allocating fields... 4 k-points: (0,0,0) (0.5,0,0) (0.5,0.5,0) (0.5,0.5,0.5) Solving for band polarization: . Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-12, mean 1.51138, harm. mean 1.14094, 34.9121% > 1, 4.64894% "fill" Outputting check-epsilon... solve_kpoint (0,0,0): freqs:, k index, k1, k2, k3, kmag/2pi, band 1, band 2, band 3 Solving for bands 3 to 3... Finished solving for bands 3 to 3 after 41 iterations. freqs:, 1, 0, 0, 0, 0, 0, 0, 0.546635 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 3... Finished solving for bands 1 to 3 after 30 iterations. freqs:, 2, 0.5, 0, 0, 0.5, 0.259951, 0.259951, 0.444658 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.5,0): Solving for bands 1 to 3... Finished solving for bands 1 to 3 after 41 iterations. freqs:, 3, 0.5, 0.5, 0, 0.707107, 0.300692, 0.345674, 0.497693 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.5,0.5): Solving for bands 1 to 3... Finished solving for bands 1 to 3 after 36 iterations. freqs:, 4, 0.5, 0.5, 0.5, 0.866025, 0.362782, 0.362782, 0.502237 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.3627824325771372 at #(0.5 0.5 0.5) Band 2 range: 0.0 at #(0.0 0.0 0.0) to 0.3627824326347085 at #(0.5 0.5 0.5) Band 3 range: 0.4446580171512465 at #(0.5 0.0 0.0) to 0.5466349527111005 at #(0.0 0.0 0.0) Gap from band 2 (0.3627824326347085) to band 3 (0.4446580171512465), 20.280278140200295% total elapsed time for run: 0 seconds. done. check-freqs: PASSED ******************************************************************************* Test case: group velocity in simple cubic lattice with anisotropic dielectric. ******************************************************************************* init-params: initializing eigensolver data Computing 3 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 3 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: block, center = (0,0,0) size (0.313,0.313,1) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (0.313,1,0.313) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (1,0.313,0.313) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (0.313,0.313,0.313) axes (1,0,0), (0,1,0), (0,0,1) epsilon = 12, mu = 1 Geometric object tree has depth 7 and 16 object nodes (vs. 4 actual objects) Initializing epsilon function... 1 k-points: (0.12,0.34,0.41) Solving for band polarization: . Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-12, mean 1.51138, harm. mean 1.14094, 34.9121% > 1, 4.64894% "fill" Outputting check-epsilon... solve_kpoint (0.12,0.34,0.41): freqs:, k index, k1, k2, k3, kmag/2pi, band 1, band 2, band 3 Solving for bands 1 to 3... Finished solving for bands 1 to 3 after 51 iterations. freqs:, 1, 0.12, 0.34, 0.41, 0.545985, 0.291848, 0.315613, 0.491642 elapsed time for k point: 1 seconds. total elapsed time for run: 1 seconds. done. check-almost-equal: PASSED check-almost-equal: PASSED ************************************************************************** Test case: symmetry transformed overlaps & inversion/mirror eigenvalues. ************************************************************************** init-params: initializing eigensolver data Computing 6 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 6 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: sphere, center = (0,0,0) radius 0.25 epsilon = 13, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Allocating fields... 1 k-points: (0.5,0.5,0.5) Solving for band polarization: . Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-13, mean 1.78533, harm. mean 1.08703, 11.2549% > 1, 6.54442% "fill" Outputting check-epsilon... solve_kpoint (0.5,0.5,0.5): freqs:, k index, k1, k2, k3, kmag/2pi, band 1, band 2, band 3, band 4, band 5, band 6 Solving for bands 1 to 6... Finished solving for bands 1 to 6 after 22 iterations. freqs:, 1, 0.5, 0.5, 0.5, 0.866025, 0.529861, 0.529881, 0.529888, 0.638532, 0.638544, 0.638573 elapsed time for k point: 0 seconds. total elapsed time for run: 0 seconds. done. check-almost-equal: PASSED check-almost-equal: PASSED init-params: initializing eigensolver data Computing 6 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 6 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: sphere, center = (0,0,0) radius 0.25 epsilon = 13, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... 1 k-points: (0,0,0) Solving for band polarization: zeven. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-13, mean 1.78533, harm. mean 1.08703, 11.2549% > 1, 6.54442% "fill" Outputting check-epsilon... solve_kpoint (0,0,0): zevenfreqs:, k index, k1, k2, k3, kmag/2pi, zeven band 1, zeven band 2, zeven band 3, zeven band 4, zeven band 5, zeven band 6 Solving for bands 2 to 6... Finished solving for bands 2 to 6 after 35 iterations. zevenfreqs:, 1, 0, 0, 0, 0, 0, 0.552074, 0.713314, 0.713341, 0.772993, 0.773109 elapsed time for k point: 0 seconds. total elapsed time for run: 0 seconds. done. check-almost-equal: PASSED check-almost-equal: PASSED init-params: initializing eigensolver data Computing 6 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 6 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: sphere, center = (0,0,0) radius 0.25 epsilon = 13, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... 1 k-points: (0,0,0) Solving for band polarization: zodd. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-13, mean 1.78533, harm. mean 1.08703, 11.2549% > 1, 6.54442% "fill" Outputting check-epsilon... solve_kpoint (0,0,0): zoddfreqs:, k index, k1, k2, k3, kmag/2pi, zodd band 1, zodd band 2, zodd band 3, zodd band 4, zodd band 5, zodd band 6 Solving for bands 2 to 6... iteration 30: trace = 2.328723746322124 (4.33192e-07% change) Finished solving for bands 2 to 6 after 31 iterations. zoddfreqs:, 1, 0, 0, 0, 0, 0, 0.552078, 0.552093, 0.71276, 0.770963, 0.785313 elapsed time for k point: 1 seconds. total elapsed time for run: 1 seconds. done. check-almost-equal: PASSED check-almost-equal: PASSED eigensolver iterations for 88 k-points: 8-18387/10, mean = 97.01410984848485, median = 20.0 mean flops per iteration = 4756300.0629984895 mean time per iteration = 0.008902172735338991 s Relative error ranged from 0.0 to 1.968577710495782e-4, with a mean of 3.7346024026232894e-6 Absolute error ranged to 1.9322975859081826e-4, with a mean of 2.7087067285005226e-6 PASSED all tests. ./mpb force-mu?=true ../examples/check.ctl command-line param: force-mu?=true ************************************************************************** Test case: 1d quarter-wave stack. ************************************************************************** init-params: initializing eigensolver data Computing 8 bands with 1.000000e-13 tolerance. Working in 1 dimensions. Grid size is 32 x 1 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 1e+20, height 0.25, axis (1, 0, 0) epsilon = 9, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... Allocating fields... 6 k-points: (0,0,0) (0.1,0,0) (0.2,0,0) (0.3,0,0) (0.4,0,0) (0.5,0,0) Solving for band polarization: tm. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-9, mean 3, harm. mean 1.32353, 28.125% > 1, 25% "fill" Outputting check-epsilon... Warning! Library version information error. The HDF5 library version information are not consistent in its source code. This is NOT a fatal error but should be corrected. Setting the environment variable 'HDF5_DISABLE_VERSION_CHECK' to a value of 1 will suppress this warning. Library version information are: H5_VERS_MAJOR=1, H5_VERS_MINOR=14, H5_VERS_RELEASE=4, H5_VERS_SUBRELEASE=2, H5_VERS_INFO=HDF5 library version: 1.14.4 mu: 1-1, mean 1, harm. mean 1, 0% > 1, 100% "fill" Outputting check-mu... solve_kpoint (0,0,0): tmfreqs:, k index, k1, k2, k3, kmag/2pi, tm band 1, tm band 2, tm band 3, tm band 4, tm band 5, tm band 6, tm band 7, tm band 8 Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 25 iterations. tmfreqs:, 1, 0, 0, 0, 0, 0, 0.666384, 0.666668, 1.33099, 1.33336, 1.99162, 2.00024, 2.64509 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 74 iterations. tmfreqs:, 2, 0.1, 0, 0, 0.1, 0.0574931, 0.608788, 0.724374, 1.2736, 1.39098, 1.93585, 2.05634, 2.59377 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tmfreqs:, 3, 0.2, 0, 0, 0.2, 0.113352, 0.552761, 0.780756, 1.21673, 1.44856, 1.87756, 2.11569, 2.53392 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 11 iterations. tmfreqs:, 4, 0.3, 0, 0, 0.3, 0.164802, 0.501202, 0.832987, 1.16418, 1.5025, 1.82314, 2.17225, 2.4769 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 11 iterations. tmfreqs:, 5, 0.4, 0, 0, 0.4, 0.205537, 0.460405, 0.874781, 1.1222, 1.54664, 1.77873, 2.22033, 2.42855 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 6, 0.5, 0, 0, 0.5, 0.222451, 0.443471, 0.892368, 1.10456, 1.56579, 1.75949, 2.24248, 2.40631 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.22245099319196907 at #(0.5 0.0 0.0) Band 2 range: 0.44347071815372 at #(0.5 0.0 0.0) to 0.6663842825289061 at #(0.0 0.0 0.0) Band 3 range: 0.666667518031661 at #(0.0 0.0 0.0) to 0.8923678144024793 at #(0.5 0.0 0.0) Band 4 range: 1.1045580524889709 at #(0.5 0.0 0.0) to 1.3309933766567734 at #(0.0 0.0 0.0) Band 5 range: 1.333360876728432 at #(0.0 0.0 0.0) to 1.5657937369267099 at #(0.5 0.0 0.0) Band 6 range: 1.7594893273125085 at #(0.5 0.0 0.0) to 1.9916198017299604 at #(0.0 0.0 0.0) Band 7 range: 2.000243026425645 at #(0.0 0.0 0.0) to 2.2424804310085227 at #(0.5 0.0 0.0) Band 8 range: 2.4063119004353655 at #(0.5 0.0 0.0) to 2.645093785208385 at #(0.0 0.0 0.0) Gap from band 1 (0.22245099319196907) to band 2 (0.44347071815372), 66.38009279352558% Gap from band 2 (0.6663842825289061) to band 3 (0.666667518031661), 0.04249429806640395% Gap from band 3 (0.8923678144024793) to band 4 (1.1045580524889709), 21.251689069139186% Gap from band 4 (1.3309933766567734) to band 5 (1.333360876728432), 0.17771661322067053% Gap from band 5 (1.5657937369267099) to band 6 (1.7594893273125085), 11.649870801607303% Gap from band 6 (1.9916198017299604) to band 7 (2.000243026425645), 0.4320401309815955% Gap from band 7 (2.2424804310085227) to band 8 (2.4063119004353655), 7.048345365685874% total elapsed time for run: 0 seconds. done. check-freqs: PASSED init-params: initializing eigensolver data Computing 8 bands with 1.000000e-13 tolerance. Working in 1 dimensions. Grid size is 32 x 1 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 1e+20, height 0.25, axis (1, 0, 0) epsilon = 9, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... 6 k-points: (0,0,0) (0.1,0,0) (0.2,0,0) (0.3,0,0) (0.4,0,0) (0.5,0,0) Solving for band polarization: te. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-9, mean 3, harm. mean 1.32353, 28.125% > 1, 25% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 100% "fill" Outputting check-mu... solve_kpoint (0,0,0): tefreqs:, k index, k1, k2, k3, kmag/2pi, te band 1, te band 2, te band 3, te band 4, te band 5, te band 6, te band 7, te band 8 Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 25 iterations. tefreqs:, 1, 0, 0, 0, 0, 0, 0.666384, 0.666668, 1.33099, 1.33336, 1.99162, 2.00024, 2.64509 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 74 iterations. tefreqs:, 2, 0.1, 0, 0, 0.1, 0.0574931, 0.608788, 0.724374, 1.2736, 1.39098, 1.93585, 2.05634, 2.59377 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tefreqs:, 3, 0.2, 0, 0, 0.2, 0.113352, 0.552761, 0.780756, 1.21673, 1.44856, 1.87756, 2.11569, 2.53392 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 11 iterations. tefreqs:, 4, 0.3, 0, 0, 0.3, 0.164802, 0.501202, 0.832987, 1.16418, 1.5025, 1.82314, 2.17225, 2.4769 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 11 iterations. tefreqs:, 5, 0.4, 0, 0, 0.4, 0.205537, 0.460405, 0.874781, 1.1222, 1.54664, 1.77873, 2.22033, 2.42855 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tefreqs:, 6, 0.5, 0, 0, 0.5, 0.222451, 0.443471, 0.892368, 1.10456, 1.56579, 1.75949, 2.24248, 2.40631 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.22245099319196965 at #(0.5 0.0 0.0) Band 2 range: 0.44347071815372274 at #(0.5 0.0 0.0) to 0.6663842825289686 at #(0.0 0.0 0.0) Band 3 range: 0.6666675180318095 at #(0.0 0.0 0.0) to 0.8923678144024818 at #(0.5 0.0 0.0) Band 4 range: 1.1045580524889715 at #(0.5 0.0 0.0) to 1.3309933766567872 at #(0.0 0.0 0.0) Band 5 range: 1.3333608767284217 at #(0.0 0.0 0.0) to 1.5657937369267096 at #(0.5 0.0 0.0) Band 6 range: 1.7594893273125132 at #(0.5 0.0 0.0) to 1.991619801730251 at #(0.0 0.0 0.0) Band 7 range: 2.0002430264256823 at #(0.0 0.0 0.0) to 2.2424804310085253 at #(0.5 0.0 0.0) Band 8 range: 2.4063119004353677 at #(0.5 0.0 0.0) to 2.6450937852113277 at #(0.0 0.0 0.0) Gap from band 1 (0.22245099319196965) to band 2 (0.44347071815372274), 66.38009279352589% Gap from band 2 (0.6663842825289686) to band 3 (0.6666675180318095), 0.04249429807928964% Gap from band 3 (0.8923678144024818) to band 4 (1.1045580524889715), 21.251689069138962% Gap from band 4 (1.3309933766567872) to band 5 (1.3333608767284217), 0.17771661321887017% Gap from band 5 (1.5657937369267096) to band 6 (1.7594893273125132), 11.64987080160758% Gap from band 6 (1.991619801730251) to band 7 (2.0002430264256823), 0.43204013096886656% Gap from band 7 (2.2424804310085253) to band 8 (2.4063119004353677), 7.048345365685847% total elapsed time for run: 0 seconds. done. check-freqs: PASSED ************************************************************************** Test case: 1d sinusoidal Bragg mirrors. ************************************************************************** init-params: initializing eigensolver data Computing 8 bands with 1.000000e-13 tolerance. Working in 1 dimensions. Grid size is 32 x 1 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 1e+20, height 0.25, axis (1, 0, 0) epsilon = 9, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... 11 k-points: (0,0,0) (0.05,0,0) (0.1,0,0) (0.15,0,0) (0.2,0,0) (0.25,0,0) (0.3,0,0) (0.35,0,0) (0.4,0,0) (0.45,0,0) (0.5,0,0) Solving for band polarization: tm. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1.00285-9, mean 4.66295, harm. mean 2.61927, 100% > 1, 45.7676% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 100% "fill" Outputting check-mu... solve_kpoint (0,0,0): tmfreqs:, k index, k1, k2, k3, kmag/2pi, tm band 1, tm band 2, tm band 3, tm band 4, tm band 5, tm band 6, tm band 7, tm band 8 Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 34 iterations. tmfreqs:, 1, 0, 0, 0, 0, 0, 0.460648, 0.542428, 0.968587, 1.01617, 1.48337, 1.48387, 1.96764 elapsed time for k point: 0 seconds. solve_kpoint (0.05,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 86 iterations. tmfreqs:, 2, 0.05, 0, 0, 0.05, 0.0231425, 0.454294, 0.548894, 0.95836, 1.02642, 1.45914, 1.50811, 1.94949 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tmfreqs:, 3, 0.1, 0, 0, 0.1, 0.0462091, 0.439084, 0.564453, 0.938267, 1.04658, 1.43467, 1.53261, 1.92548 elapsed time for k point: 0 seconds. solve_kpoint (0.15,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 4, 0.15, 0, 0, 0.15, 0.0691102, 0.420016, 0.584145, 0.91575, 1.06922, 1.41023, 1.55711, 1.9011 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 5, 0.2, 0, 0, 0.2, 0.0917239, 0.399481, 0.605653, 0.892484, 1.09267, 1.38581, 1.58161, 1.87665 elapsed time for k point: 0 seconds. solve_kpoint (0.25,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 6, 0.25, 0, 0, 0.25, 0.113863, 0.378568, 0.628026, 0.868958, 1.11645, 1.36143, 1.6061, 1.85216 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 7, 0.3, 0, 0, 0.3, 0.135213, 0.357979, 0.650805, 0.845399, 1.14036, 1.33713, 1.63058, 1.82769 elapsed time for k point: 0 seconds. solve_kpoint (0.35,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 8, 0.35, 0, 0, 0.35, 0.155194, 0.338479, 0.673671, 0.822003, 1.16428, 1.31295, 1.65502, 1.80325 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 9, 0.4, 0, 0, 0.4, 0.17268, 0.32129, 0.696193, 0.799138, 1.188, 1.28906, 1.67935, 1.77891 elapsed time for k point: 0 seconds. solve_kpoint (0.45,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 8 iterations. tmfreqs:, 10, 0.45, 0, 0, 0.45, 0.185503, 0.308628, 0.717034, 0.778102, 1.21076, 1.26621, 1.70326, 1.755 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 8 iterations. tmfreqs:, 11, 0.5, 0, 0, 0.5, 0.190416, 0.303766, 0.72859, 0.766483, 1.2253, 1.25164, 1.721, 1.73726 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.19041596916807002 at #(0.5 0.0 0.0) Band 2 range: 0.303766174808589 at #(0.5 0.0 0.0) to 0.460648275218067 at #(0.0 0.0 0.0) Band 3 range: 0.5424277685174704 at #(0.0 0.0 0.0) to 0.7285900510563766 at #(0.5 0.0 0.0) Band 4 range: 0.7664833049941241 at #(0.5 0.0 0.0) to 0.9685865863610894 at #(0.0 0.0 0.0) Band 5 range: 1.0161706920388196 at #(0.0 0.0 0.0) to 1.2253024188808164 at #(0.5 0.0 0.0) Band 6 range: 1.2516393338141774 at #(0.5 0.0 0.0) to 1.483373337947421 at #(0.0 0.0 0.0) Band 7 range: 1.483865937538404 at #(0.0 0.0 0.0) to 1.7209986559263173 at #(0.5 0.0 0.0) Band 8 range: 1.737259239629598 at #(0.5 0.0 0.0) to 1.9676376863404632 at #(0.0 0.0 0.0) Gap from band 1 (0.19041596916807002) to band 2 (0.303766174808589), 45.87385724963493% Gap from band 2 (0.460648275218067) to band 3 (0.5424277685174704), 16.305741485929598% Gap from band 3 (0.7285900510563766) to band 4 (0.7664833049941241), 5.069082902774636% Gap from band 4 (0.9685865863610894) to band 5 (1.0161706920388196), 4.794954647158872% Gap from band 5 (1.2253024188808164) to band 6 (1.2516393338141774), 2.126567159256415% Gap from band 6 (1.483373337947421) to band 7 (1.483865937538404), 0.03320255262544471% Gap from band 7 (1.7209986559263173) to band 8 (1.737259239629598), 0.9403916188076465% total elapsed time for run: 0 seconds. done. check-freqs: PASSED ************************************************************************** Test case: get-dominant-planewave. ************************************************************************** init-params: initializing eigensolver data Computing 8 bands with 1.000000e-09 tolerance. Working in 1 dimensions. Grid size is 32 x 1 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: Geometric object tree has depth 1 and 0 object nodes (vs. 0 actual objects) Initializing epsilon function... Initializing mu function... 1 k-points: (0.4,0,0) Solving for band polarization: te. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-1, mean 1, harm. mean 1, 0% > 1, 100% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 100% "fill" Outputting check-mu... solve_kpoint (0.4,0,0): tefreqs:, k index, k1, k2, k3, kmag/2pi, te band 1, te band 2, te band 3, te band 4, te band 5, te band 6, te band 7, te band 8 Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 11 iterations. tefreqs:, 1, 0.4, 0, 0, 0.4, 0.4, 0.6, 1.4, 1.6, 2.4, 2.6, 3.4, 3.6 elapsed time for k point: 0 seconds. total elapsed time for run: 0 seconds. done. check-almost-equal: PASSED ************************************************************************** Test case: Square lattice of dielectric rods in air. ************************************************************************** init-params: initializing eigensolver data Computing 8 bands with 1.000000e-09 tolerance. Working in 2 dimensions. Grid size is 32 x 32 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 11.56, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... Allocating fields... 16 k-points: (0,0,0) (0.1,0,0) (0.2,0,0) (0.3,0,0) (0.4,0,0) (0.5,0,0) (0.5,0.1,0) (0.5,0.2,0) (0.5,0.3,0) (0.5,0.4,0) (0.5,0.5,0) (0.4,0.4,0) (0.3,0.3,0) (0.2,0.2,0) (0.1,0.1,0) (0,0,0) Solving for band polarization: te. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-11.56, mean 2.327, harm. mean 1.1441, 14.5508% > 1, 12.5663% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 0% "fill" Outputting check-mu... solve_kpoint (0,0,0): tefreqs:, k index, k1, k2, k3, kmag/2pi, te band 1, te band 2, te band 3, te band 4, te band 5, te band 6, te band 7, te band 8 Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 35 iterations. tefreqs:, 1, 0, 0, 0, 0, 0, 0.561945, 0.780843, 0.780846, 0.924372, 1.00804, 1.00804, 1.09859 elapsed time for k point: 1 seconds. solve_kpoint (0.1,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 19 iterations. tefreqs:, 2, 0.1, 0, 0, 0.1, 0.0897998, 0.56021, 0.767785, 0.782438, 0.912603, 1.00889, 1.00949, 1.12115 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0,0): Solving for bands 1 to 8... iteration 40: trace = 5.671940683808276 (0.000353561% change) Finished solving for bands 1 to 8 after 92 iterations. tefreqs:, 3, 0.2, 0, 0, 0.2, 0.178853, 0.553301, 0.732522, 0.786672, 0.890463, 1.01386, 1.02078, 1.11301 elapsed time for k point: 1 seconds. solve_kpoint (0.3,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 15 iterations. tefreqs:, 4, 0.3, 0, 0, 0.3, 0.266124, 0.534865, 0.689376, 0.79205, 0.872926, 1.0209, 1.04285, 1.10753 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 20 iterations. tefreqs:, 5, 0.4, 0, 0, 0.4, 0.349588, 0.49478, 0.658508, 0.796524, 0.862684, 1.02898, 1.07044, 1.10003 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 34 iterations. tefreqs:, 6, 0.5, 0, 0, 0.5, 0.413346, 0.444623, 0.648672, 0.798265, 0.859327, 1.03311, 1.0958, 1.09742 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.1,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 19 iterations. tefreqs:, 7, 0.5, 0.1, 0, 0.509902, 0.424298, 0.448983, 0.644549, 0.802551, 0.854853, 0.98984, 1.05831, 1.11679 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.2,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 15 iterations. tefreqs:, 8, 0.5, 0.2, 0, 0.538516, 0.455354, 0.461161, 0.63323, 0.814782, 0.834251, 0.934558, 1.01274, 1.12367 elapsed time for k point: 1 seconds. solve_kpoint (0.5,0.3,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 14 iterations. tefreqs:, 9, 0.5, 0.3, 0, 0.583095, 0.478427, 0.501912, 0.61745, 0.784489, 0.83371, 0.906541, 0.967734, 1.12636 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.4,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 14 iterations. tefreqs:, 10, 0.5, 0.4, 0, 0.640312, 0.495605, 0.556867, 0.601739, 0.720584, 0.858432, 0.897772, 0.926611, 1.12762 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.5,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tefreqs:, 11, 0.5, 0.5, 0, 0.707107, 0.503606, 0.594344, 0.594357, 0.679389, 0.883665, 0.895782, 0.895787, 1.12801 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0.4,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tefreqs:, 12, 0.4, 0.4, 0, 0.565685, 0.474671, 0.549917, 0.608032, 0.745035, 0.848355, 0.895854, 0.945891, 1.12708 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0.3,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 15 iterations. tefreqs:, 13, 0.3, 0.3, 0, 0.424264, 0.373254, 0.543225, 0.646251, 0.817389, 0.830457, 0.896516, 1.0176, 1.12198 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0.2,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 20 iterations. tefreqs:, 14, 0.2, 0.2, 0, 0.282843, 0.252369, 0.551073, 0.700804, 0.797552, 0.899586, 0.90339, 1.08791, 1.10099 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0.1,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 24 iterations. tefreqs:, 15, 0.1, 0.1, 0, 0.141421, 0.12694, 0.558854, 0.755651, 0.785047, 0.909884, 0.968813, 1.04814, 1.11332 elapsed time for k point: 0 seconds. solve_kpoint (0,0,0): Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 73 iterations. tefreqs:, 16, 0, 0, 0, 0, 0, 0.561945, 0.780843, 0.780846, 0.924372, 1.00804, 1.00804, 1.09859 elapsed time for k point: 1 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.5036058015317009 at #(0.5 0.5 0.0) Band 2 range: 0.4446229135155296 at #(0.5 0.0 0.0) to 0.5943440246130545 at #(0.5 0.5 0.0) Band 3 range: 0.5943566395317801 at #(0.5 0.5 0.0) to 0.7808428124011783 at #(0.0 0.0 0.0) Band 4 range: 0.6793887413186159 at #(0.5 0.5 0.0) to 0.8173893720817362 at #(0.30000000000000004 0.30000000000000004 0.0) Band 5 range: 0.830457382271798 at #(0.30000000000000004 0.30000000000000004 0.0) to 0.9243716832604366 at #(0.0 0.0 0.0) Band 6 range: 0.8957817684697044 at #(0.5 0.5 0.0) to 1.033110413950032 at #(0.5 0.0 0.0) Band 7 range: 0.8957868745561711 at #(0.5 0.5 0.0) to 1.0958021498671253 at #(0.5 0.0 0.0) Band 8 range: 1.0974168097704617 at #(0.5 0.0 0.0) to 1.1280127657657455 at #(0.5 0.5 0.0) Gap from band 2 (0.5943440246130545) to band 3 (0.5943566395317801), 0.002122471889875288% Gap from band 4 (0.8173893720817362) to band 5 (0.830457382271798), 1.5860710536992217% Gap from band 7 (1.0958021498671253) to band 8 (1.0974168097704617), 0.1472411038798638% total elapsed time for run: 4 seconds. done. check-freqs: PASSED init-params: initializing eigensolver data Computing 8 bands with 1.000000e-09 tolerance. Working in 2 dimensions. Grid size is 32 x 32 x 1. Solving for 8 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 0.2, height 1e+20, axis (0, 0, 1) epsilon = 11.56, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... 16 k-points: (0,0,0) (0.1,0,0) (0.2,0,0) (0.3,0,0) (0.4,0,0) (0.5,0,0) (0.5,0.1,0) (0.5,0.2,0) (0.5,0.3,0) (0.5,0.4,0) (0.5,0.5,0) (0.4,0.4,0) (0.3,0.3,0) (0.2,0.2,0) (0.1,0.1,0) (0,0,0) Solving for band polarization: tm. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-11.56, mean 2.327, harm. mean 1.1441, 14.5508% > 1, 12.5663% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 0% "fill" Outputting check-mu... solve_kpoint (0,0,0): tmfreqs:, k index, k1, k2, k3, kmag/2pi, tm band 1, tm band 2, tm band 3, tm band 4, tm band 5, tm band 6, tm band 7, tm band 8 Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 21 iterations. tmfreqs:, 1, 0, 0, 0, 0, 0, 0.550336, 0.561338, 0.56134, 0.822948, 0.868842, 0.965325, 1.08938 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 21 iterations. tmfreqs:, 2, 0.1, 0, 0, 0.1, 0.0651416, 0.525004, 0.561885, 0.586307, 0.823536, 0.867346, 0.954632, 1.05688 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 3, 0.2, 0, 0, 0.2, 0.127665, 0.49365, 0.563323, 0.617312, 0.822737, 0.86353, 0.924306, 1.03882 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 4, 0.3, 0, 0, 0.3, 0.184047, 0.461592, 0.565122, 0.651397, 0.810839, 0.858984, 0.892949, 1.03589 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 5, 0.4, 0, 0, 0.4, 0.227778, 0.433361, 0.566596, 0.689045, 0.778792, 0.855425, 0.879398, 1.03874 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 6, 0.5, 0, 0, 0.5, 0.245809, 0.420657, 0.567163, 0.720092, 0.747203, 0.85409, 0.877012, 1.0408 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.1,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 8 iterations. tmfreqs:, 7, 0.5, 0.1, 0, 0.509902, 0.249299, 0.427308, 0.560221, 0.718025, 0.756122, 0.855013, 0.877106, 1.02917 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.2,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 8, 0.5, 0.2, 0, 0.538516, 0.258693, 0.445541, 0.543385, 0.711698, 0.77935, 0.858136, 0.877357, 1.00092 elapsed time for k point: 1 seconds. solve_kpoint (0.5,0.3,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 9, 0.5, 0.3, 0, 0.583095, 0.270963, 0.470322, 0.524053, 0.701396, 0.810963, 0.864394, 0.877683, 0.965108 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.4,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 10, 0.5, 0.4, 0, 0.640312, 0.281614, 0.493207, 0.508823, 0.689961, 0.846785, 0.874595, 0.878068, 0.926581 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.5,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 12 iterations. tmfreqs:, 11, 0.5, 0.5, 0, 0.707107, 0.285906, 0.502981, 0.502983, 0.684476, 0.874359, 0.883317, 0.883317, 0.892993 elapsed time for k point: 0 seconds. solve_kpoint (0.4,0.4,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 13 iterations. tmfreqs:, 12, 0.4, 0.4, 0, 0.565685, 0.276089, 0.491352, 0.508683, 0.69258, 0.839723, 0.856431, 0.907219, 0.907347 elapsed time for k point: 0 seconds. solve_kpoint (0.3,0.3,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 13, 0.3, 0.3, 0, 0.424264, 0.240239, 0.47923, 0.523499, 0.685362, 0.829266, 0.840451, 0.910791, 0.941595 elapsed time for k point: 0 seconds. solve_kpoint (0.2,0.2,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 9 iterations. tmfreqs:, 14, 0.2, 0.2, 0, 0.282843, 0.17524, 0.488533, 0.541561, 0.64751, 0.83002, 0.850443, 0.922657, 0.983925 elapsed time for k point: 0 seconds. solve_kpoint (0.1,0.1,0): Solving for bands 1 to 8... Finished solving for bands 1 to 8 after 10 iterations. tmfreqs:, 15, 0.1, 0.1, 0, 0.141421, 0.0915259, 0.516393, 0.555924, 0.601211, 0.824561, 0.86529, 0.948412, 1.03526 elapsed time for k point: 0 seconds. solve_kpoint (0,0,0): Solving for bands 2 to 8... Finished solving for bands 2 to 8 after 22 iterations. tmfreqs:, 16, 0, 0, 0, 0, 0, 0.550336, 0.561338, 0.56134, 0.822948, 0.868842, 0.965325, 1.08938 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.2859057791271611 at #(0.5 0.5 0.0) Band 2 range: 0.4206573384061857 at #(0.5 0.0 0.0) to 0.5503360755136674 at #(0.0 0.0 0.0) Band 3 range: 0.5029830978387358 at #(0.5 0.5 0.0) to 0.5671632878212806 at #(0.5 0.0 0.0) Band 4 range: 0.5613397939915359 at #(0.0 0.0 0.0) to 0.720091820469093 at #(0.5 0.0 0.0) Band 5 range: 0.7472029910634785 at #(0.5 0.0 0.0) to 0.8743593805271213 at #(0.5 0.5 0.0) Band 6 range: 0.8404509697715067 at #(0.30000000000000004 0.30000000000000004 0.0) to 0.883317372585055 at #(0.5 0.5 0.0) Band 7 range: 0.8770118718590372 at #(0.5 0.0 0.0) to 0.965325380932607 at #(0.0 0.0 0.0) Band 8 range: 0.8929933495601431 at #(0.5 0.5 0.0) to 1.089377902232541 at #(0.0 0.0 0.0) Gap from band 1 (0.2859057791271611) to band 2 (0.4206573384061857), 38.14282289442171% Gap from band 4 (0.720091820469093) to band 5 (0.7472029910634785), 3.695395142312023% total elapsed time for run: 1 seconds. done. check-freqs: PASSED **************************************************************************** Test case: square lattice of magneto-electric rods in air. **************************************************************************** init-params: initializing eigensolver data Computing 10 bands with 1.000000e-09 tolerance. Working in 2 dimensions. Grid size is 32 x 32 x 1. Solving for 10 bands at a time. Creating Maxwell data... Mesh size is 3. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: cylinder, center = (0,0,0) radius 0.11, height 1e+20, axis (0, 0, 1) Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... Allocating fields... 10 k-points: (0,0,0) (0.166667,0,0) (0.333333,0,0) (0.5,0,0) (0.5,0.166667,0) (0.5,0.333333,0) (0.5,0.5,0) (0.333333,0.333333,0) (0.166667,0.166667,0) (0,0,0) Solving for band polarization: tm. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-15, mean 1.53164, harm. mean 1.04189, 4.39453% > 1, 3.79742% "fill" Outputting check-epsilon... mu: 1-13.9136, mean 1.42353, harm. mean 1.03906, 4.39453% > 1, 3.27973% "fill" Outputting check-mu... solve_kpoint (0,0,0): tmfreqs:, k index, k1, k2, k3, kmag/2pi, tm band 1, tm band 2, tm band 3, tm band 4, tm band 5, tm band 6, tm band 7, tm band 8, tm band 9, tm band 10 Solving for bands 2 to 10... iteration 73: trace = 3.284543891383691 (0.11949% change) iteration 135: trace = 2.948968004805939 (0.0161451% change) iteration 199: trace = 2.69898322825364 (0.0253256% change) iteration 224: trace = 2.674665467966938 (0.0604414% change) iteration 282: trace = 2.452429910048918 (0.0565924% change) iteration 342: trace = 2.414033997869022 (0.00436713% change) iteration 404: trace = 2.408300336607943 (0.00700163% change) iteration 467: trace = 2.403716107831429 (0.0016442% change) iteration 521: trace = 2.400812781858495 (0.000798663% change) iteration 582: trace = 2.398066166760905 (0.000977101% change) iteration 652: trace = 2.396011000473594 (0.000459715% change) iteration 724: trace = 2.39479247032656 (0.000671889% change) iteration 763: trace = 2.394284751882646 (0.000263241% change) iteration 802: trace = 2.393821456521515 (0.000298517% change) iteration 841: trace = 2.393596245332671 (0.000100502% change) iteration 907: trace = 2.393276739907282 (0.00176964% change) iteration 959: trace = 2.393134929014358 (1.84165e-05% change) iteration 1012: trace = 2.393097865449186 (3.44397e-05% change) iteration 1069: trace = 2.393061025206864 (8.24005e-06% change) iteration 1135: trace = 2.393050981878563 (4.43877e-06% change) iteration 1188: trace = 2.393045838325383 (1.75991e-05% change) iteration 1251: trace = 2.393042463255056 (1.02795e-06% change) iteration 1311: trace = 2.393041309578154 (2.17328e-07% change) iteration 1380: trace = 2.393040388715379 (2.31234e-07% change) Finished solving for bands 2 to 10 after 1391 iterations. tmfreqs:, 1, 0, 0, 0, 0, 0, 0.265885, 0.356852, 0.36899, 0.503897, 0.506552, 0.539911, 0.635681, 0.645852, 0.660017 elapsed time for k point: 12 seconds. solve_kpoint (0.166667,0,0): Solving for bands 1 to 10... iteration 42: trace = 2.414344733390499 (0.00142812% change) iteration 96: trace = 2.412984274155959 (0.000571832% change) iteration 152: trace = 2.412702871965183 (3.53249e-05% change) iteration 209: trace = 2.412669742580899 (3.51466e-05% change) iteration 264: trace = 2.4126578641278 (2.70032e-06% change) iteration 324: trace = 2.412644932464488 (3.32172e-06% change) iteration 380: trace = 2.412640482397965 (2.87884e-06% change) iteration 433: trace = 2.412636225852249 (4.91028e-06% change) iteration 484: trace = 2.412634223661579 (1.84641e-06% change) iteration 534: trace = 2.412633040415287 (5.12986e-07% change) iteration 572: trace = 2.41263206326654 (2.21221e-06% change) Finished solving for bands 1 to 10 after 600 iterations. tmfreqs:, 2, 0.166667, 0, 0, 0.166667, 0.123306, 0.278292, 0.357366, 0.369461, 0.503797, 0.506642, 0.538217, 0.63557, 0.645595, 0.659432 elapsed time for k point: 6 seconds. solve_kpoint (0.333333,0,0): Solving for bands 1 to 10... iteration 56: trace = 2.455244994197765 (0.00191746% change) iteration 114: trace = 2.453415732188682 (0.00122397% change) iteration 164: trace = 2.452559142544338 (0.000281813% change) iteration 223: trace = 2.452235897761291 (0.000166585% change) iteration 285: trace = 2.451999658357778 (6.99726e-05% change) iteration 348: trace = 2.451902096889206 (6.4124e-05% change) iteration 406: trace = 2.45185685161993 (4.6202e-05% change) iteration 456: trace = 2.451827436393056 (2.57041e-05% change) iteration 503: trace = 2.451816227903016 (1.00689e-05% change) iteration 549: trace = 2.451806263605273 (6.55096e-06% change) iteration 592: trace = 2.451800140847527 (1.06412e-05% change) iteration 647: trace = 2.451796441509538 (1.91273e-06% change) iteration 709: trace = 2.451793184022576 (8.19139e-07% change) iteration 771: trace = 2.451791565915028 (6.02552e-07% change) Finished solving for bands 1 to 10 after 797 iterations. tmfreqs:, 3, 0.333333, 0, 0, 0.333333, 0.194655, 0.331458, 0.361679, 0.377947, 0.502693, 0.506755, 0.527227, 0.633335, 0.644215, 0.653288 elapsed time for k point: 8 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 10... iteration 51: trace = 2.455981933477009 (0.00440596% change) iteration 112: trace = 2.450757126441165 (0.00163361% change) iteration 165: trace = 2.449742916815921 (0.000376882% change) iteration 205: trace = 2.449512617940444 (0.000480813% change) iteration 228: trace = 2.449424159549771 (0.000106911% change) iteration 280: trace = 2.449366077368989 (7.03198e-06% change) iteration 334: trace = 2.449334879074217 (4.42443e-05% change) iteration 379: trace = 2.449315603452948 (7.12579e-06% change) iteration 432: trace = 2.449307403169984 (2.57347e-06% change) iteration 489: trace = 2.44930060431031 (3.36301e-06% change) iteration 549: trace = 2.449296929134784 (6.66391e-07% change) iteration 602: trace = 2.449295519391521 (3.54523e-07% change) iteration 655: trace = 2.449294132150464 (8.47062e-07% change) Finished solving for bands 1 to 10 after 699 iterations. tmfreqs:, 4, 0.5, 0, 0, 0.5, 0.208799, 0.348978, 0.365278, 0.443856, 0.477902, 0.506845, 0.509076, 0.617527, 0.63971, 0.646891 elapsed time for k point: 7 seconds. solve_kpoint (0.5,0.166667,0): Solving for bands 1 to 10... iteration 55: trace = 2.489464544406752 (0.00312308% change) iteration 109: trace = 2.487234699332965 (0.000912293% change) iteration 164: trace = 2.486209593514058 (0.000300727% change) iteration 218: trace = 2.486012788401608 (3.42504e-05% change) iteration 263: trace = 2.485938358536024 (8.06293e-05% change) iteration 312: trace = 2.48589496543296 (1.00375e-05% change) iteration 341: trace = 2.485884876844713 (2.85511e-06% change) iteration 370: trace = 2.485879110091961 (4.7741e-06% change) iteration 414: trace = 2.48587184429556 (1.76915e-06% change) iteration 473: trace = 2.485867742960644 (2.5336e-06% change) iteration 519: trace = 2.485866654055232 (2.60299e-07% change) iteration 579: trace = 2.485865728060445 (4.14613e-07% change) Finished solving for bands 1 to 10 after 598 iterations. tmfreqs:, 5, 0.5, 0.166667, 0, 0.527046, 0.210673, 0.351103, 0.364915, 0.453355, 0.500502, 0.507018, 0.510125, 0.616821, 0.642478, 0.646894 elapsed time for k point: 6 seconds. solve_kpoint (0.5,0.333333,0): Solving for bands 1 to 10... iteration 61: trace = 2.584147703201229 (0.00398781% change) iteration 121: trace = 2.579932344483964 (0.00179053% change) iteration 155: trace = 2.579249962924961 (0.000447169% change) iteration 189: trace = 2.579032773013822 (0.000163534% change) iteration 240: trace = 2.57892546489783 (1.81957e-05% change) iteration 274: trace = 2.57889938977938 (1.07981e-05% change) iteration 305: trace = 2.578889619155277 (8.78132e-06% change) iteration 339: trace = 2.578885134517451 (8.73343e-07% change) iteration 381: trace = 2.578882900952178 (1.3719e-06% change) iteration 434: trace = 2.578881958966127 (1.9612e-07% change) Finished solving for bands 1 to 10 after 444 iterations. tmfreqs:, 6, 0.5, 0.333333, 0, 0.600925, 0.214081, 0.353217, 0.365001, 0.471977, 0.507165, 0.51079, 0.563894, 0.615052, 0.646178, 0.648357 elapsed time for k point: 5 seconds. solve_kpoint (0.5,0.5,0): Solving for bands 1 to 10... iteration 60: trace = 2.660517315971257 (0.00351086% change) iteration 120: trace = 2.654222746758167 (0.00135733% change) iteration 180: trace = 2.653889537238184 (0.000171716% change) iteration 231: trace = 2.653865612669728 (2.75258e-06% change) iteration 291: trace = 2.653861461688727 (1.65938e-07% change) Finished solving for bands 1 to 10 after 308 iterations. tmfreqs:, 7, 0.5, 0.5, 0, 0.707107, 0.215649, 0.353602, 0.365269, 0.480748, 0.507335, 0.512534, 0.60792, 0.619784, 0.648908, 0.650323 elapsed time for k point: 3 seconds. solve_kpoint (0.333333,0.333333,0): Solving for bands 1 to 10... iteration 50: trace = 2.535202897662554 (0.00316097% change) iteration 94: trace = 2.532923086867295 (0.000314726% change) iteration 129: trace = 2.532442482774198 (0.000367272% change) iteration 189: trace = 2.532266376173264 (8.03707e-05% change) iteration 248: trace = 2.532232517403471 (2.72844e-06% change) iteration 309: trace = 2.532228060067403 (1.77109e-06% change) iteration 341: trace = 2.532227158000674 (1.01927e-06% change) Finished solving for bands 1 to 10 after 367 iterations. tmfreqs:, 8, 0.333333, 0.333333, 0, 0.471405, 0.210879, 0.35463, 0.361245, 0.439486, 0.507022, 0.508252, 0.543815, 0.620599, 0.646589, 0.651154 elapsed time for k point: 4 seconds. solve_kpoint (0.166667,0.166667,0): Solving for bands 1 to 10... iteration 54: trace = 2.431502621209246 (0.00594636% change) iteration 112: trace = 2.429488400673375 (0.000107952% change) iteration 173: trace = 2.429229261545035 (7.13575e-05% change) iteration 218: trace = 2.429186892539972 (9.66426e-06% change) iteration 263: trace = 2.429173188159541 (9.34178e-06% change) iteration 315: trace = 2.429165184961218 (4.78154e-06% change) iteration 366: trace = 2.429162053705261 (1.56622e-06% change) iteration 407: trace = 2.429161001525221 (1.57814e-06% change) Finished solving for bands 1 to 10 after 443 iterations. tmfreqs:, 9, 0.166667, 0.166667, 0, 0.235702, 0.162353, 0.295875, 0.356214, 0.372161, 0.504258, 0.506735, 0.537355, 0.633231, 0.645916, 0.65717 elapsed time for k point: 4 seconds. solve_kpoint (0,0,0): Solving for bands 2 to 10... iteration 74: trace = 2.393322735884219 (0.00032093% change) iteration 144: trace = 2.39306075287922 (2.30517e-05% change) iteration 217: trace = 2.393042681631062 (2.08168e-06% change) iteration 288: trace = 2.393040379176155 (6.07113e-07% change) Finished solving for bands 2 to 10 after 304 iterations. tmfreqs:, 10, 0, 0, 0, 0, 0, 0.265885, 0.356852, 0.36899, 0.503897, 0.506552, 0.539911, 0.635681, 0.645853, 0.660016 elapsed time for k point: 2 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.21564853471674153 at #(0.5 0.5 0.0) Band 2 range: 0.2658849940257783 at #(0.0 0.0 0.0) to 0.3546302418196162 at #(0.3333333333333333 0.3333333333333333 0.0) Band 3 range: 0.3562142389790025 at #(0.16666666666666666 0.16666666666666666 0.0) to 0.36527833701789714 at #(0.5 0.0 0.0) Band 4 range: 0.36899014763669635 at #(0.0 0.0 0.0) to 0.4807484735010422 at #(0.5 0.5 0.0) Band 5 range: 0.47790151572201556 at #(0.5 0.0 0.0) to 0.5073352931310028 at #(0.5 0.5 0.0) Band 6 range: 0.5065518606522352 at #(0.0 0.0 0.0) to 0.5125343710321437 at #(0.5 0.5 0.0) Band 7 range: 0.5090758763184378 at #(0.5 0.0 0.0) to 0.607920097111441 at #(0.5 0.5 0.0) Band 8 range: 0.6150518826422737 at #(0.5 0.33333333333333337 0.0) to 0.6356808127406188 at #(0.0 0.0 0.0) Band 9 range: 0.6397097834016685 at #(0.5 0.0 0.0) to 0.6489079032761239 at #(0.5 0.5 0.0) Band 10 range: 0.6468905136920686 at #(0.5 0.0 0.0) to 0.6600169180967272 at #(0.0 0.0 0.0) Gap from band 1 (0.21564853471674153) to band 2 (0.2658849940257783), 20.865196839034088% Gap from band 2 (0.3546302418196162) to band 3 (0.3562142389790025), 0.4456663031572527% Gap from band 3 (0.36527833701789714) to band 4 (0.36899014763669635), 1.0110227243500136% Gap from band 7 (0.607920097111441) to band 8 (0.6150518826422737), 1.1663039953325716% Gap from band 8 (0.6356808127406188) to band 9 (0.6397097834016685), 0.6318018453697605% total elapsed time for run: 57 seconds. done. check-freqs: PASSED ************************************************************************** Test case: 3x3 triangular lattice of rods in air, dipole defect states. ************************************************************************** ************************************************************************** Test case: fcc lattice of air spheres in dielectric. ************************************************************************** init-params: initializing eigensolver data Computing 10 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 10 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (0, 0.707107, 0.707107) (0.707107, 0, 0.707107) (0.707107, 0.707107, 0) Cell volume = 0.707107 Reciprocal lattice vectors (/ 2 pi): (-0.707107, 0.707107, 0.707107) (0.707107, -0.707107, 0.707107) (0.707107, 0.707107, -0.707107) Geometric objects: sphere, center = (0,0,0) radius 0.5 epsilon = 1, mu = 1 Geometric object tree has depth 4 and 64 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... Allocating fields... 13 k-points: (0,0.5,0.5) (0,0.5625,0.4375) (0,0.625,0.375) (0,0.5625,0.1875) (0,0.5,0) (0,0.25,0) (0,0,0) (0,0.25,0.25) (0,0.5,0.5) (0.125,0.625,0.5) (0.25,0.75,0.5) (0.3125,0.75,0.4375) (0.375,0.75,0.375) Solving for band polarization: . Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-11.56, mean 3.72661, harm. mean 1.46352, 45.5811% > 1, 25.8202% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 37.5% "fill" Outputting check-mu... solve_kpoint (0,0.5,0.5): freqs:, k index, k1, k2, k3, kmag/2pi, band 1, band 2, band 3, band 4, band 5, band 6, band 7, band 8, band 9, band 10 Solving for bands 1 to 10... iteration 14: trace = 2.422816182530671 (0.0381649% change) iteration 26: trace = 2.418028459594542 (0.00355549% change) iteration 42: trace = 2.417926671402957 (4.12491e-07% change) Finished solving for bands 1 to 10 after 49 iterations. freqs:, 1, 0, 0.5, 0.5, 0.707107, 0.370307, 0.372085, 0.38252, 0.3836, 0.493989, 0.512361, 0.522307, 0.523822, 0.595103, 0.664079 elapsed time for k point: 2 seconds. solve_kpoint (0,0.5625,0.4375): Solving for bands 1 to 10... iteration 12: trace = 2.407427410269236 (0.0321996% change) iteration 25: trace = 2.406716506193876 (4.38833e-06% change) Finished solving for bands 1 to 10 after 31 iterations. freqs:, 2, 0, 0.5625, 0.4375, 0.71807, 0.368359, 0.377487, 0.384992, 0.387438, 0.471895, 0.506546, 0.523725, 0.531091, 0.610293, 0.649714 elapsed time for k point: 1 seconds. solve_kpoint (0,0.625,0.375): Solving for bands 1 to 10... iteration 16: trace = 2.416824100020383 (1.3666e-05% change) Finished solving for bands 1 to 10 after 23 iterations. freqs:, 3, 0, 0.625, 0.375, 0.75, 0.357533, 0.38095, 0.392789, 0.401454, 0.43935, 0.495158, 0.527478, 0.540892, 0.635515, 0.644621 elapsed time for k point: 1 seconds. solve_kpoint (0,0.5625,0.1875): Solving for bands 1 to 10... iteration 15: trace = 2.409011173205657 (0.00176438% change) iteration 30: trace = 2.408973129092047 (1.24458e-07% change) Finished solving for bands 1 to 10 after 31 iterations. freqs:, 4, 0, 0.5625, 0.1875, 0.649519, 0.323742, 0.331653, 0.39821, 0.401206, 0.462868, 0.514292, 0.533282, 0.546116, 0.631046, 0.644236 elapsed time for k point: 1 seconds. solve_kpoint (0,0.5,0): Solving for bands 1 to 10... iteration 11: trace = 2.3863536100807 (0.000116095% change) Finished solving for bands 1 to 10 after 18 iterations. freqs:, 5, 0, 0.5, 0, 0.612372, 0.307529, 0.308891, 0.387977, 0.389913, 0.491606, 0.536329, 0.536613, 0.538985, 0.624681, 0.629124 elapsed time for k point: 1 seconds. solve_kpoint (0,0.25,0): Solving for bands 1 to 10... iteration 17: trace = 2.4210576593944 (3.68809e-05% change) Finished solving for bands 1 to 10 after 27 iterations. freqs:, 6, 0, 0.25, 0, 0.306186, 0.179755, 0.180421, 0.473402, 0.475968, 0.50346, 0.53577, 0.538144, 0.540345, 0.624073, 0.627676 elapsed time for k point: 1 seconds. solve_kpoint (0,0,0): Solving for bands 3 to 10... iteration 14: trace = 2.454797764127216 (0.000697847% change) iteration 29: trace = 2.454761954994568 (6.80438e-05% change) Finished solving for bands 3 to 10 after 39 iterations. freqs:, 7, 0, 0, 0, 0, 0, 0, 0.518283, 0.521124, 0.521132, 0.543942, 0.543952, 0.546505, 0.61287, 0.61397 elapsed time for k point: 1 seconds. solve_kpoint (0,0.25,0.25): Solving for bands 1 to 10... iteration 17: trace = 2.429734562412393 (0.0269389% change) iteration 33: trace = 2.428244402873749 (0.000697566% change) iteration 50: trace = 2.428132004368119 (2.18611e-06% change) Finished solving for bands 1 to 10 after 57 iterations. freqs:, 8, 0, 0.25, 0.25, 0.353553, 0.207258, 0.207841, 0.473224, 0.475268, 0.506873, 0.526874, 0.530366, 0.532023, 0.604689, 0.653955 elapsed time for k point: 2 seconds. solve_kpoint (0,0.5,0.5): Solving for bands 1 to 10... iteration 15: trace = 2.429840506213174 (0.0494696% change) iteration 28: trace = 2.418683855228275 (0.00708379% change) iteration 30: trace = 2.418447796007927 (0.00429416% change) iteration 47: trace = 2.417970501348602 (0.000224729% change) iteration 61: trace = 2.417939628661853 (3.38411e-05% change) iteration 76: trace = 2.417934651155199 (1.46663e-06% change) Finished solving for bands 1 to 10 after 83 iterations. freqs:, 9, 0, 0.5, 0.5, 0.707107, 0.370307, 0.372085, 0.38252, 0.3836, 0.493989, 0.512361, 0.522307, 0.523822, 0.595103, 0.664085 elapsed time for k point: 3 seconds. solve_kpoint (0.125,0.625,0.5): Solving for bands 1 to 10... iteration 17: trace = 2.402816909617924 (2.57746e-05% change) Finished solving for bands 1 to 10 after 27 iterations. freqs:, 10, 0.125, 0.625, 0.5, 0.728869, 0.372309, 0.377154, 0.385968, 0.392562, 0.463569, 0.502026, 0.505399, 0.549256, 0.621502, 0.639203 elapsed time for k point: 1 seconds. solve_kpoint (0.25,0.75,0.5): Solving for bands 1 to 10... iteration 17: trace = 2.423195772610677 (0.013681% change) iteration 33: trace = 2.422694284495717 (1.73988e-05% change) Finished solving for bands 1 to 10 after 45 iterations. freqs:, 11, 0.25, 0.75, 0.5, 0.790569, 0.373709, 0.386404, 0.386499, 0.410208, 0.435229, 0.490809, 0.492518, 0.566583, 0.627971, 0.654078 elapsed time for k point: 2 seconds. solve_kpoint (0.3125,0.75,0.4375): Solving for bands 1 to 10... iteration 17: trace = 2.419666612124221 (0.0026528% change) iteration 34: trace = 2.419550566011305 (2.18919e-06% change) Finished solving for bands 1 to 10 after 40 iterations. freqs:, 12, 0.3125, 0.75, 0.4375, 0.760345, 0.364261, 0.382754, 0.390011, 0.405967, 0.438258, 0.493602, 0.505736, 0.559219, 0.634896, 0.645093 elapsed time for k point: 1 seconds. solve_kpoint (0.375,0.75,0.375): Solving for bands 1 to 10... iteration 14: trace = 2.422576717217212 (0.000134415% change) Finished solving for bands 1 to 10 after 25 iterations. freqs:, 13, 0.375, 0.75, 0.375, 0.75, 0.359122, 0.380258, 0.391276, 0.402681, 0.438005, 0.49516, 0.529122, 0.541679, 0.637583, 0.645624 elapsed time for k point: 1 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.373709380455834 at #(0.25 0.75 0.5) Band 2 range: 0.0 at #(0.0 0.0 0.0) to 0.38640412149814385 at #(0.25 0.75 0.5) Band 3 range: 0.3825201383155223 at #(0.0 0.5 0.5) to 0.5182834488052885 at #(0.0 0.0 0.0) Band 4 range: 0.38360034956418215 at #(0.0 0.5 0.5) to 0.5211243460581922 at #(0.0 0.0 0.0) Band 5 range: 0.43522923853967 at #(0.25 0.75 0.5) to 0.521132474139247 at #(0.0 0.0 0.0) Band 6 range: 0.49080863941128233 at #(0.25 0.75 0.5) to 0.5439419096018507 at #(0.0 0.0 0.0) Band 7 range: 0.4925183484876648 at #(0.25 0.75 0.5) to 0.5439517867826289 at #(0.0 0.0 0.0) Band 8 range: 0.5238217224838788 at #(0.0 0.5 0.5) to 0.5665833313384158 at #(0.25 0.75 0.5) Band 9 range: 0.595103036081524 at #(0.0 0.5 0.5) to 0.6375830839471325 at #(0.375 0.75 0.375) Band 10 range: 0.6139704161698655 at #(0.0 0.0 0.0) to 0.6640850273032767 at #(0.0 0.5 0.5) Gap from band 8 (0.5665833313384158) to band 9 (0.595103036081524), 4.910052410522718% total elapsed time for run: 18 seconds. done. check-freqs: PASSED ************************************************************************** Test case: simple cubic lattice with anisotropic dielectric. ************************************************************************** init-params: initializing eigensolver data Computing 3 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 3 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: block, center = (0,0,0) size (0.313,0.313,1) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (0.313,1,0.313) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (1,0.313,0.313) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (0.313,0.313,0.313) axes (1,0,0), (0,1,0), (0,0,1) epsilon = 12, mu = 1 Geometric object tree has depth 7 and 16 object nodes (vs. 4 actual objects) Initializing epsilon function... Initializing mu function... Allocating fields... 4 k-points: (0,0,0) (0.5,0,0) (0.5,0.5,0) (0.5,0.5,0.5) Solving for band polarization: . Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-12, mean 1.51138, harm. mean 1.14094, 34.9121% > 1, 4.64894% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 100% "fill" Outputting check-mu... solve_kpoint (0,0,0): freqs:, k index, k1, k2, k3, kmag/2pi, band 1, band 2, band 3 Solving for bands 3 to 3... Finished solving for bands 3 to 3 after 41 iterations. freqs:, 1, 0, 0, 0, 0, 0, 0, 0.546635 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0,0): Solving for bands 1 to 3... Finished solving for bands 1 to 3 after 30 iterations. freqs:, 2, 0.5, 0, 0, 0.5, 0.259951, 0.259951, 0.444658 elapsed time for k point: 0 seconds. solve_kpoint (0.5,0.5,0): Solving for bands 1 to 3... Finished solving for bands 1 to 3 after 44 iterations. freqs:, 3, 0.5, 0.5, 0, 0.707107, 0.300692, 0.345674, 0.497693 elapsed time for k point: 1 seconds. solve_kpoint (0.5,0.5,0.5): Solving for bands 1 to 3... Finished solving for bands 1 to 3 after 36 iterations. freqs:, 4, 0.5, 0.5, 0.5, 0.866025, 0.362782, 0.362782, 0.502237 elapsed time for k point: 0 seconds. Band 1 range: 0.0 at #(0.0 0.0 0.0) to 0.3627824325770722 at #(0.5 0.5 0.5) Band 2 range: 0.0 at #(0.0 0.0 0.0) to 0.3627824326347104 at #(0.5 0.5 0.5) Band 3 range: 0.44465801844530395 at #(0.5 0.0 0.0) to 0.5466349524266712 at #(0.0 0.0 0.0) Gap from band 2 (0.3627824326347104) to band 3 (0.44465801844530395), 20.280278428230492% total elapsed time for run: 1 seconds. done. check-freqs: PASSED ******************************************************************************* Test case: group velocity in simple cubic lattice with anisotropic dielectric. ******************************************************************************* init-params: initializing eigensolver data Computing 3 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 3 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: block, center = (0,0,0) size (0.313,0.313,1) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (0.313,1,0.313) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (1,0.313,0.313) axes (1,0,0), (0,1,0), (0,0,1) block, center = (0,0,0) size (0.313,0.313,0.313) axes (1,0,0), (0,1,0), (0,0,1) epsilon = 12, mu = 1 Geometric object tree has depth 7 and 16 object nodes (vs. 4 actual objects) Initializing epsilon function... Initializing mu function... 1 k-points: (0.12,0.34,0.41) Solving for band polarization: . Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-12, mean 1.51138, harm. mean 1.14094, 34.9121% > 1, 4.64894% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 100% "fill" Outputting check-mu... solve_kpoint (0.12,0.34,0.41): freqs:, k index, k1, k2, k3, kmag/2pi, band 1, band 2, band 3 Solving for bands 1 to 3... Finished solving for bands 1 to 3 after 49 iterations. freqs:, 1, 0.12, 0.34, 0.41, 0.545985, 0.291848, 0.315613, 0.491642 elapsed time for k point: 1 seconds. total elapsed time for run: 1 seconds. done. check-almost-equal: PASSED check-almost-equal: PASSED ************************************************************************** Test case: symmetry transformed overlaps & inversion/mirror eigenvalues. ************************************************************************** init-params: initializing eigensolver data Computing 6 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 6 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: sphere, center = (0,0,0) radius 0.25 epsilon = 13, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... Allocating fields... 1 k-points: (0.5,0.5,0.5) Solving for band polarization: . Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-13, mean 1.78533, harm. mean 1.08703, 11.2549% > 1, 6.54442% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 33.3333% "fill" Outputting check-mu... solve_kpoint (0.5,0.5,0.5): freqs:, k index, k1, k2, k3, kmag/2pi, band 1, band 2, band 3, band 4, band 5, band 6 Solving for bands 1 to 6... Finished solving for bands 1 to 6 after 21 iterations. freqs:, 1, 0.5, 0.5, 0.5, 0.866025, 0.529861, 0.529881, 0.529888, 0.638532, 0.638544, 0.638573 elapsed time for k point: 0 seconds. total elapsed time for run: 0 seconds. done. check-almost-equal: PASSED check-almost-equal: PASSED init-params: initializing eigensolver data Computing 6 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 6 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: sphere, center = (0,0,0) radius 0.25 epsilon = 13, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... 1 k-points: (0,0,0) Solving for band polarization: zeven. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-13, mean 1.78533, harm. mean 1.08703, 11.2549% > 1, 6.54442% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 33.3333% "fill" Outputting check-mu... solve_kpoint (0,0,0): zevenfreqs:, k index, k1, k2, k3, kmag/2pi, zeven band 1, zeven band 2, zeven band 3, zeven band 4, zeven band 5, zeven band 6 Solving for bands 2 to 6... iteration 34: trace = 2.517674614600258 (2.39432e-07% change) iteration 68: trace = 2.517674409630287 (3.22847e-07% change) Finished solving for bands 2 to 6 after 78 iterations. zevenfreqs:, 1, 0, 0, 0, 0, 0, 0.552075, 0.713308, 0.713335, 0.772992, 0.773116 elapsed time for k point: 1 seconds. total elapsed time for run: 1 seconds. done. check-almost-equal: PASSED check-almost-equal: PASSED init-params: initializing eigensolver data Computing 6 bands with 1.000000e-09 tolerance. Working in 3 dimensions. Grid size is 16 x 16 x 16. Solving for 6 bands at a time. Creating Maxwell data... Mesh size is 5. Lattice vectors: (1, 0, 0) (0, 1, 0) (0, 0, 1) Cell volume = 1 Reciprocal lattice vectors (/ 2 pi): (1, -0, 0) (-0, 1, -0) (0, -0, 1) Geometric objects: sphere, center = (0,0,0) radius 0.25 epsilon = 13, mu = 1 Geometric object tree has depth 1 and 1 object nodes (vs. 1 actual objects) Initializing epsilon function... Initializing mu function... 1 k-points: (0,0,0) Solving for band polarization: zodd. Initializing fields to random numbers... elapsed time for initialization: 0 seconds. epsilon: 1-13, mean 1.78533, harm. mean 1.08703, 11.2549% > 1, 6.54442% "fill" Outputting check-epsilon... mu: 1-1, mean 1, harm. mean 1, 0% > 1, 33.3333% "fill" Outputting check-mu... solve_kpoint (0,0,0): zoddfreqs:, k index, k1, k2, k3, kmag/2pi, zodd band 1, zodd band 2, zodd band 3, zodd band 4, zodd band 5, zodd band 6 Solving for bands 2 to 6... iteration 33: trace = 2.328662645557328 (7.69489e-07% change) iteration 65: trace = 2.328660801902074 (1.01133e-05% change) Finished solving for bands 2 to 6 after 72 iterations. zoddfreqs:, 1, 0, 0, 0, 0, 0, 0.552078, 0.552093, 0.71276, 0.770964, 0.785268 elapsed time for k point: 2 seconds. total elapsed time for run: 2 seconds. done. check-almost-equal: PASSED check-almost-equal: PASSED eigensolver iterations for 87 k-points: 8-12519/10, mean = 88.4022030651341, median = 19.25 mean flops per iteration = 5151897.749640001 mean time per iteration = 0.011051890794316727 s Relative error ranged from 0.0 to 1.97981434581789e-4, with a mean of 3.6146586099975884e-6 Absolute error ranged to 1.9322135011645702e-4, with a mean of 2.717376015201954e-6 PASSED all tests. make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' Making check in examples make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/examples' make[1]: Nothing to be done for 'check'. make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/examples' Making check in utils make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/utils' make[1]: Nothing to be done for 'check'. make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/utils' make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0' make[1]: Nothing to be done for 'check-am'. make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0' >>> mpb: Entering fakeroot... libfakeroot internal error: payload not recognized! Making install in src make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make install-recursive make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' Making install in util make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' make install-am make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' make[5]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' make[5]: Nothing to be done for 'install-exec-am'. make[5]: Nothing to be done for 'install-data-am'. make[5]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/util' Making install in matrices make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrices' make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrices' make[4]: Nothing to be done for 'install-exec-am'. make[4]: Nothing to be done for 'install-data-am'. make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrices' make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrices' Making install in matrixio make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrixio' make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrixio' make[4]: Nothing to be done for 'install-exec-am'. make[4]: Nothing to be done for 'install-data-am'. make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrixio' make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/matrixio' Making install in maxwell make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/maxwell' make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/maxwell' make[4]: Nothing to be done for 'install-exec-am'. make[4]: Nothing to be done for 'install-data-am'. make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/maxwell' make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src/maxwell' make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make[4]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' /bin/mkdir -p '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/include' /bin/mkdir -p '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/lib' /bin/mkdir -p '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/include/mpb' /bin/sh ../libtool --mode=install /usr/bin/install -c libmpb.la '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/lib' /usr/bin/install -c -m 644 mpb.h '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/include' /usr/bin/install -c -m 644 matrices/eigensolver.h matrices/matrices.h matrices/scalar.h maxwell/maxwell.h util/verbosity.h '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/include/mpb' libtool: install: /usr/bin/install -c .libs/libmpb.so.6.0.0 /home/buildozer/aports/testing/mpb/pkg/mpb/usr/lib/libmpb.so.6.0.0 libtool: install: (cd /home/buildozer/aports/testing/mpb/pkg/mpb/usr/lib && { ln -s -f libmpb.so.6.0.0 libmpb.so.6 || { rm -f libmpb.so.6 && ln -s libmpb.so.6.0.0 libmpb.so.6; }; }) libtool: install: (cd /home/buildozer/aports/testing/mpb/pkg/mpb/usr/lib && { ln -s -f libmpb.so.6.0.0 libmpb.so || { rm -f libmpb.so && ln -s libmpb.so.6.0.0 libmpb.so; }; }) libtool: install: /usr/bin/install -c .libs/libmpb.lai /home/buildozer/aports/testing/mpb/pkg/mpb/usr/lib/libmpb.la libtool: install: /usr/bin/install -c .libs/libmpb.a /home/buildozer/aports/testing/mpb/pkg/mpb/usr/lib/libmpb.a libtool: install: chmod 644 /home/buildozer/aports/testing/mpb/pkg/mpb/usr/lib/libmpb.a libtool: install: ranlib /home/buildozer/aports/testing/mpb/pkg/mpb/usr/lib/libmpb.a libtool: warning: remember to run 'libtool --finish /usr/lib' make[4]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/src' Making install in tests make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/tests' make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/tests' make[2]: Nothing to be done for 'install-exec-am'. make[2]: Nothing to be done for 'install-data-am'. make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/tests' make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/tests' Making install in mpb make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' make install-am make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' make[3]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' /bin/mkdir -p '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/bin' /bin/mkdir -p '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/bin' /bin/mkdir -p '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/share/man/man1' /bin/mkdir -p '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/share/mpb' /bin/sh ../libtool --mode=install /usr/bin/install -c mpb '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/bin' /usr/bin/install -c mpb-split '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/bin' /usr/bin/install -c -m 644 mpb.scm '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/share/mpb' /usr/bin/install -c -m 644 mpb-split.1 mpb.1 '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/share/man/man1' libtool: warning: '../src/libmpb.la' has not been installed in '/usr/lib' libtool: install: /usr/bin/install -c .libs/mpb /home/buildozer/aports/testing/mpb/pkg/mpb/usr/bin/mpb make[3]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/mpb' Making install in examples make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/examples' make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/examples' make[2]: Nothing to be done for 'install-exec-am'. make[2]: Nothing to be done for 'install-data-am'. make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/examples' make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/examples' Making install in utils make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/utils' make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/utils' /bin/mkdir -p '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/bin' /bin/mkdir -p '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/share/man/man1' /bin/sh ../libtool --mode=install /usr/bin/install -c mpb-data '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/bin' /usr/bin/install -c -m 644 mpb-data.1 '/home/buildozer/aports/testing/mpb/pkg/mpb/usr/share/man/man1' libtool: warning: '../src/libmpb.la' has not been installed in '/usr/lib' libtool: install: /usr/bin/install -c .libs/mpb-data /home/buildozer/aports/testing/mpb/pkg/mpb/usr/bin/mpb-data make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/utils' make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0/utils' make[1]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0' make[2]: Entering directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0' make[2]: Nothing to be done for 'install-exec-am'. make[2]: Nothing to be done for 'install-data-am'. make[2]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0' make[1]: Leaving directory '/home/buildozer/aports/testing/mpb/src/mpb-1.12.0' >>> mpb-dev*: Running split function dev... 'usr/include' -> '/home/buildozer/aports/testing/mpb/pkg/mpb-dev/usr/include' 'usr/lib/libmpb.a' -> '/home/buildozer/aports/testing/mpb/pkg/mpb-dev/usr/lib/libmpb.a' 'usr/lib/libmpb.so' -> '/home/buildozer/aports/testing/mpb/pkg/mpb-dev/usr/lib/libmpb.so' >>> mpb-dev*: Preparing subpackage mpb-dev... >>> mpb-dev*: Stripping binaries >>> mpb-dev*: Running postcheck for mpb-dev >>> mpb-doc*: Running split function doc... 'usr/share/man' -> '/home/buildozer/aports/testing/mpb/pkg/mpb-doc/usr/share/man' >>> mpb-doc*: Preparing subpackage mpb-doc... >>> mpb-doc*: Running postcheck for mpb-doc >>> mpb*: Running postcheck for mpb >>> mpb*: Preparing package mpb... >>> mpb*: Stripping binaries >>> mpb*: Scanning shared objects >>> mpb-dev*: Scanning shared objects >>> mpb*: Tracing dependencies... openblas so:libc.musl-s390x.so.1 so:libctl.so.7 so:libctlgeom.so.7 so:libfftw3.so.3 so:libguile-3.0.so.1 so:libhdf5.so.311 >>> mpb*: Package size: 404.0 KB >>> mpb*: Compressing data... >>> mpb*: Create checksum... >>> mpb*: Create mpb-1.12.0-r0.apk >>> mpb-dev*: Tracing dependencies... mpb=1.12.0-r0 >>> mpb-dev*: Package size: 172.2 KB >>> mpb-dev*: Compressing data... >>> mpb-dev*: Create checksum... >>> mpb-dev*: Create mpb-dev-1.12.0-r0.apk >>> mpb-doc*: Tracing dependencies... >>> mpb-doc*: Package size: 5.9 KB >>> mpb-doc*: Compressing data... >>> mpb-doc*: Create checksum... >>> mpb-doc*: Create mpb-doc-1.12.0-r0.apk >>> mpb: Build complete at Tue, 30 Jun 2026 04:02:18 +0000 elapsed time 0h 2m 52s >>> mpb: Cleaning up srcdir >>> mpb: Cleaning up pkgdir >>> mpb: Uninstalling dependencies... ( 1/31) Purging .makedepends-mpb (20260630.035927) ( 2/31) Purging fftw-dev (3.3.11-r0) ( 3/31) Purging fftw-double-libs (3.3.11-r0) ( 4/31) Purging fftw-long-double-libs (3.3.11-r0) ( 5/31) Purging fftw-single-libs (3.3.11-r0) ( 6/31) Purging hdf5-dev (1.14.4.2-r1) ( 7/31) Purging zlib-dev (1.3.2-r0) ( 8/31) Purging hdf5-cpp (1.14.4.2-r1) ( 9/31) Purging hdf5-hl-cpp (1.14.4.2-r1) (10/31) Purging hdf5-hl-fortran (1.14.4.2-r1) (11/31) Purging hdf5-hl (1.14.4.2-r1) (12/31) Purging libctl-dev (4.5.1-r1) (13/31) Purging guile-dev (3.0.9-r2) (14/31) Purging guile (3.0.9-r2) (15/31) Purging gc-dev (8.2.12-r0) (16/31) Purging libgc++ (8.2.12-r0) (17/31) Purging gmp-dev (6.3.0-r4) (18/31) Purging libgmpxx (6.3.0-r4) (19/31) Purging libctl (4.5.1-r1) (20/31) Purging openblas-dev (0.3.30-r2) (21/31) Purging gfortran (15.2.0-r6) (22/31) Purging liblapack (0.3.30-r2) (23/31) Purging liblapacke (0.3.30-r2) (24/31) Purging openblas (0.3.30-r2) (25/31) Purging guile-libs (3.0.9-r2) (26/31) Purging hdf5-fortran (1.14.4.2-r1) (27/31) Purging libffi (3.5.2-r1) (28/31) Purging libgfortran (15.2.0-r6) (29/31) Purging gc (8.2.12-r0) (30/31) Purging hdf5 (1.14.4.2-r1) (31/31) Purging libsz (1.1.6-r0) Executing busybox-1.38.0-r1.trigger OK: 386.4 MiB in 107 packages >>> mpb: Updating the testing/s390x repository index... >>> mpb: Signing the index...