Dear Steven and MPB-users:
I am trying to use mpb-mpi to calculate the band structure for a special
2-dimensional photonic crystal system. The ctl file is as follows:
; supercell scheme
; * * * o o o o * * * --------layer one
;* * * o o o o o * * * --------layer two
(define-param supercell-y 21) ; the (odd) number of lateral supercell
periods,layer two
(set! geometry-lattice (make lattice
(basis1 (sqrt 2) 0)
(basis2 0 (sqrt 2))
(size 1 supercell-y no-size)))
(define-param eps 12) ; the dielectric constant of background
(define-param h 4.5) ; Size of the slab
(define-param slabw (- (/ supercell-y 2) (/ h 2)))
(define-param slaby (+ (/ h 2) (/ slabw 2)))
(set! default-material (make dielectric (epsilon eps)))
;make periodic crossing holes
(set! geometry
(list (make block (center 0 0 0) (material air) (size 0.5
0.2 infinity))
(make block (center 0 0 0) (material air) (size
0.2 0.5 infinity))
(make block (center 0.5 0.5 0) (material air)
(size 0.5 0.2 infinity))
(make block (center 0.5 0.5 0) (material air)
(size 0.2 0.5 infinity))))
(set! geometry
(append
; duplicate the crossings over the supercell:
(geometric-objects-lattice-duplicates geometry)
(list (make block (center 0 (* -1 slaby) 0) (material air) (size 1.0
slabw infinity))
(make block (center 0 slaby 0) (material air)
(size 1.0 slabw infinity)))))
(define Gamma (vector3 0 0 0))
(define X (vector3 0.5 0 0))
(define M (vector3 0.5 0.5 0))
(set! k-points (interpolate 9 (list Gamma X)))
; the bigger the supercell, the more bands you need to compute to get
; to the defect modes (the lowest band is "folded" supercell-y times):
(define-param extra-bands 40) ; number of extra bands to compute above the
gap
(set! num-bands (+ supercell-y extra-bands))
(set-param! resolution 32)
(begin-time
"total time for both TE and TM bands: "
(run-te-yeven (output-at-kpoint (vector3 0.45 0 0) output-hfield-z)))
(display-eigensolver-stats)
I used "mpirun -np 4 mpb-mpi luo.ctl >& luo.out", but it failed. The error
report is as follows:
Solving for bands 1 to 11...
[hydra:18284] *** Process received signal ***
[hydra:18284] Signal: Segmentation fault (11)
[hydra:18284] Signal code: Address not mapped (1)
[hydra:18284] Failing at address: 0x1b
[hydra:18284] [ 0] /lib64/libc.so.6 [0x7f853c700560]
[hydra:18284] [ 1]
/usr/lib64/mpi/gcc/openmpi/lib64/libmpi.so.0(MPI_Sendrecv+0x9b)
[0x7f853d7c47fb]
[hydra:18284] [ 2]
/usr/lib64/libdfftw_mpi.so.2(transpose_start_exchange_step+0x1a9)
[0x7f853f1826b9]
[hydra:18284] [ 3] /usr/lib64/libdfftw_mpi.so.2(transpose_mpi+0x33b)
[0x7f853f182b5b]
[hydra:18284] [ 4] /usr/lib64/libdfftw_mpi.so.2(fftwnd_mpi+0xd4)
[0x7f853f1836b4]
[hydra:18284] [ 5] mpb-mpi(maxwell_compute_fft+0x75) [0x428de0]
[hydra:18284] [ 6] mpb-mpi(maxwell_compute_d_from_H+0x257) [0x4291fc]
[hydra:18284] [ 7] mpb-mpi(maxwell_operator+0xe3) [0x42992c]
[hydra:18284] [ 8] mpb-mpi(eigensolver+0xf74) [0x43321e]
[hydra:18284] [ 9] mpb-mpi(solve_kpoint+0xaf3) [0x412e0f]
[hydra:18284] [10] mpb-mpi(solve_kpoint_aux+0x43) [0x4103f3]
[hydra:18284] [11] /usr/lib64/libguile.so.17 [0x7f854035b8a6]
[hydra:18284] [12] /usr/lib64/libguile.so.17 [0x7f854035a74c]
[hydra:18284] [13] /usr/lib64/libguile.so.17 [0x7f85403599cc]
[hydra:18284] [14] /usr/lib64/libguile.so.17(scm_eval_body+0xa8)
[0x7f85403619e8]
[hydra:18284] [15] /usr/lib64/libguile.so.17(scm_map+0x289) [0x7f8540362c19]
[hydra:18284] [16] /usr/lib64/libguile.so.17 [0x7f8540359f13]
[hydra:18284] [17] /usr/lib64/libguile.so.17 [0x7f85403599cc]
[hydra:18284] *** End of error message ***
[hydra:18282] *** Process received signal ***
[hydra:18282] Signal: Segmentation fault (11)
[hydra:18282] Signal code: Address not mapped (1)
[hydra:18282] Failing at address: 0x1b
[hydra:18282] [ 0] /lib64/libc.so.6 [0x7fa0920c5560]
[hydra:18282] [ 1]
/usr/lib64/mpi/gcc/openmpi/lib64/libmpi.so.0(MPI_Sendrecv+0x9b)
[0x7fa0931897fb]
[hydra:18282] [ 2]
/usr/lib64/libdfftw_mpi.so.2(transpose_start_exchange_step+0x1a9)
[0x7fa094b476b9]
[hydra:18282] [ 3] /usr/lib64/libdfftw_mpi.so.2(transpose_mpi+0x33b)
[0x7fa094b47b5b]
[hydra:18282] [ 4] /usr/lib64/libdfftw_mpi.so.2(fftwnd_mpi+0xd4)
[0x7fa094b486b4]
[hydra:18282] [ 5] mpb-mpi(maxwell_compute_fft+0x75) [0x428de0]
[hydra:18282] [ 6] mpb-mpi(maxwell_compute_d_from_H+0x257) [0x4291fc]
[hydra:18282] [ 7] mpb-mpi(maxwell_operator+0xe3) [0x42992c]
[hydra:18282] [ 8] mpb-mpi(eigensolver+0xf74) [0x43321e]
[hydra:18282] [ 9] mpb-mpi(solve_kpoint+0xaf3) [0x412e0f]
[hydra:18282] [10] mpb-mpi(solve_kpoint_aux+0x43) [0x4103f3]
[hydra:18282] [11] /usr/lib64/libguile.so.17 [0x7fa095d208a6]
[hydra:18282] [12] /usr/lib64/libguile.so.17 [0x7fa095d1f74c]
[hydra:18282] [13] /usr/lib64/libguile.so.17 [0x7fa095d1e9cc]
[hydra:18282] [14] /usr/lib64/libguile.so.17(scm_eval_body+0xa8)
[0x7fa095d269e8]
[hydra:18282] [15] /usr/lib64/libguile.so.17(scm_map+0x289) [0x7fa095d27c19]
[hydra:18282] [16] /usr/lib64/libguile.so.17 [0x7fa095d1ef13]
[hydra:18282] [17] /usr/lib64/libguile.so.17 [0x7fa095d1e9cc]
[hydra:18282] *** End of error message ***
--------------------------------------------------------------------------
mpirun noticed that process rank 2 with PID 18284 on node hydra exited on
signal 11 (Segmentation fault).
Could you please tell me what's wrong and how to fix the errors? Thanks a
lot!
Zhang, Hao
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
