calculation stops

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