Dear developers and others,
I experienced a strange problem when using mpb -- calculation stops at some
point.
I tried to run compiled from source and precompiled version for Ubuntu of
mpb and mpb-mpi.
But the problem does not disappear.
The end of the output file from mpb looks as follows:
iteration 9990: trace = 0.00614987443779062 (9.25949e-05% change)
iteration 9992: trace = 0.006149862923197785 (9.25559e-05% change)
iteration 9994: trace = 0.00614985218661355 (8.68906e-05% change)
iteration 9996: trace = 0.006149841476666812 (8.9275e-05% change)
iteration 9998: trace = 0.006149830919231611 (8.55977e-05% change)
The end of the output file from mpb-mpi looks as follows:
iteration 9980: trace = 0.006149329221480074 (7.39655e-05% change)
iteration 9984: trace = 0.006149309028516356 (9.35957e-05% change)
iteration 9988: trace = 0.006149288431732018 (8.2345e-05% change)
iteration 9992: trace = 0.006149265844844411 (8.87949e-05% change)
iteration 9996: trace = 0.006149244524129693 (9.49614e-05% change)
p0_14835: p4_error: : 1
p2_14845: p4_error: : 1
p1_14840: p4_error: : 1
p3_14850: p4_error: : 1
rm_l_3_14851: (26829.835938) net_send: could not write to fd=5, errno = 32
rm_l_2_14847: (26831.250000) net_send: could not write to fd=5, errno = 32
rm_l_1_14842: (26832.472656) net_send: could not write to fd=5, errno = 32
p2_14845: (26833.253906) net_send: could not write to fd=5, errno = 32
p3_14850: (26831.835938) net_send: could not write to fd=5, errno = 32
p0_14835: (26839.792969) net_send: could not write to fd=4, errno = 32
Here is my .ctl file:
-----------------------------------------------------------------------------------------
(define-param a 1) ;lattice constant (distance between neibouhring rods)
(define-param b (/ a (sqrt 3))) ;length of a side of a hexagonal unit cell
(define-param F 0.26) ; filling fraction f~(4 Pi R1^2)/(2 Sqrt[3] a^2)
(define-param fm 0.215) ; filling fraction f~(4 Pi rm1^2)/(2 Sqrt[3] a^2)
(define-param eps 14) ; dielectric
(define-param L 200)
(define-param l 26)
(define-param N (+ 1 (* 2 (round (/ (* 0.5 L) (* (sqrt 3) a)) ) )))
(define-param n (+ 1 (* 2 (round (/ (* 0.5 l) (* (sqrt 3) a)) ) )))
(define-param sx 480) ; 480=3*5*2^5
(define-param sy (* a iks))
(define-param R (* (sqrt (* F (/ (* 2 (sqrt 3)) (* 4 3.141592654)))) a)) ;
radius of rodes
(define-param rm (* (sqrt (* fm (/ (* 2 (sqrt 3)) (* 4 3.141592654)))) a))
(define-param r (* 0.6 rm)) ; radius
of small additional rods
(set! geometry-lattice (make lattice (size sx sy no-size)
(basis1 sx 0)
(basis2 0 1)))
(set! geometry
(append
(geometric-objects-duplicates (vector3 (* -1 (* (sqrt 3) a)) 0) (* -1 N)
N
(geometric-object-duplicates (vector3 0 a) 0 (/ sy a)
(make cylinder (center (* 4 (* (* -0.5 a) (sqrt 3))) (* -0.5 sy))
(radius R) (height infinity)
(material (make dielectric (epsilon eps))))
)
)
(geometric-objects-duplicates (vector3 (* -1 (* (sqrt 3) a)) 0) (* -1 N)
(- N 1)
(geometric-object-duplicates (vector3 0 a) 0 (- (/ sy a) 1)
(make cylinder (center (* 5 (* (* -0.5 a) (sqrt 3))) (+ (* -0.5 sy)
(* 0.5 a))) (radius R) (height infinity)
(material (make dielectric (epsilon eps))))
)
)
(geometric-objects-duplicates (vector3 (* -1 (* (sqrt 3) a)) 0) 0 n
(geometric-object-duplicates (vector3 0 a) 0 (/ sy a)
(make cylinder (center (* (+ 8 (* 2 N)) (* (* -0.5 a) (sqrt 3))) (* -
0.5 sy)) (radius rm) (height infinity)
(material (make dielectric (epsilon eps))))
)
)
(geometric-objects-duplicates (vector3 (* -1 (* (sqrt 3) a)) 0) 0 (- n
1)
(geometric-object-duplicates (vector3 0 a) 0 (- (/ sy a) 1)
(make cylinder (center (* (+ 9 (* 2 N)) (* (* -0.5 a) (sqrt 3))) (+
(* -0.5 sy) (* 0.5 a))) (radius rm) (height infinity)
(material (make dielectric (epsilon eps))))
)
)
(geometric-objects-duplicates (vector3 (* -1 (* (sqrt 3) a)) 0) 0 (- n
1) ; (additional small rods)
(geometric-object-duplicates (vector3 0 a) 0 (/ sy a)
(make cylinder (center (+ b (* (+ 10 (* 2 N)) (* (* -0.5 a) (sqrt
3)))) (* -0.5 sy)) (radius r) (height infinity)
(material (make dielectric (epsilon eps))))
)
)
(geometric-objects-duplicates (vector3 (* -1 (* (sqrt 3) a)) 0) 0 (- n
1) ; (additional small rods)
(geometric-object-duplicates (vector3 0 a) 0 (- (/ sy a) 1)
(make cylinder (center (+ b (* (+ 9 (* 2 N)) (* (* -0.5 a) (sqrt
3)))) (+ (* -0.5 sy) (* 0.5 a))) (radius r) (height infinity)
(material (make dielectric (epsilon eps))))
)
)
(geometric-objects-duplicates (vector3 (* (sqrt 3) a) 0) 0 n
(geometric-object-duplicates (vector3 0 a) 0 (/ sy a)
(make cylinder (center (* (+ 8 (* 2 N)) (* (* 0.5 a) (sqrt 3))) (* -
0.5 sy)) (radius rm) (height infinity)
(material (make dielectric (epsilon eps))))
)
)
(geometric-objects-duplicates (vector3 (* (sqrt 3) a) 0) 0 (- n 1)
(geometric-object-duplicates (vector3 0 a) 0 (- (/ sy a) 1)
(make cylinder (center (* (+ 9 (* 2 N)) (* (* 0.5 a) (sqrt 3))) (+
(* -0.5 sy) (* 0.5 a))) (radius rm) (height infinity)
(material (make dielectric (epsilon eps))))
)
)
(geometric-objects-duplicates (vector3 (* (sqrt 3) a) 0) 0 (- n 1) ;
(additional small rods)
(geometric-object-duplicates (vector3 0 a) 0 (/ sy a)
(make cylinder (center (+ b (* (+ 8 (* 2 N)) (* (* 0.5 a) (sqrt 3))))
(* -0.5 sy)) (radius r) (height infinity)
(material (make dielectric (epsilon eps))))
)
)
(geometric-objects-duplicates (vector3 (* (sqrt 3) a) 0) 0 (- n 1) ;
(additional small rods)
(geometric-object-duplicates (vector3 0 a) 0 (- (/ sy a) 1)
(make cylinder (center (+ b (* (+ 9 (* 2 N)) (* (* 0.5 a) (sqrt
3)))) (+ (* -0.5 sy) (* 0.5 a))) (radius r) (height infinity)
(material (make dielectric (epsilon eps))))
)
)
)
)
(set! num-bands 20)
(set! target-freq 0.486)
(set! tolerance 1e-7)
(set! resolution 8)
(set! k-points (list (vector3 0 0.5)
(vector3 0 0) ; center
)
)
(set! k-points (interpolate 19 k-points))
(run-te)
---------------------------------------------------------------------------------------------------
Thank you for the help in advance.
All the best,
Ruslan
_______________________________________________
mpb-discuss mailing list
[email protected]
http://ab-initio.mit.edu/cgi-bin/mailman/listinfo/mpb-discuss
