
 BerkeleyGW trunk revision 6943
 Epsilon code, Complex version, run on 05-Apr-2016 at 13:24:01 -0700

                                                                 ..o.          
                                                                .oxxo.         
                                                               .oxxxxo...      
                                                               oxxxxxxo.       
                                                              .oxxxxxxx.       
                                                              .ooooooxxo..     
                                                              .oooooooxo..     
                                                              .oooooxxo...     
                                                       .........oxooo......    
                                                 ............................  
                                             ................................. 
                                          .................................... 
           .          ..oo. ....  .................................oooxxxxxxxo.
     .............oxxxx@ox@@@x@x.....................o...........ooooooooooxx. 
    .o.........oox@x.oo........xxx@@............ooxxxxo..........ooooxxxxxoxo  
    .x........x@xxo...............o@xxo........oxxx@@@xoooooooooooooooxxxo...  
    .o......ox@@o..................oox@o.....ooxxx@xoooxxxxxxxoooooooooooo.... 
    o..ooooo@@xoooo....ooo...........x@o.....ooxxxxo   .oxxxxxxxxxxooooooo.... 
    . .oooo@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@....ooooox.     .oxx@@@@@xxoo........ 
      .ooooxxxxxxxooooooxooooooooooooooxo...oooooxx.        ..ox@xxxoo.........
      .ooooooxxxx@xoooooooooooooxxoooooooooooxxxxx.            .oxxooooooooxxo.
     .oooooxxxxx@@@xxxxxxxxxxxxxxxxxxxxxxxxoxxxxo.                .oxxxxxxxxo. 
   ....oooxxxxx@@@xo..oxxx@@@@@@xxxxoxxoooooooxx.                  .oxxxoo..   
  .....ooxxxx@@xo.       ........    .ooooooooxxo                              
  ..oooxxxx@@@o                       .oooooooxoo.                             
  ....oxooxxxxx.                       .ooo..oooo.                             
  .....o.ooxxxxxo.                     .oooooooxo.                             
 ......ooooxxxxxxo.                     .ooooooxoo..                           
........ooxxxxxxxo..                     .o....oxoo...                         
.......ooooxxxxxxxo.                     ........oooo.                         
.ooooooo..ooxxxxoooo.                    .........ooo...                       
..oxo...ooooxxxoooo..                    .ooo......oooo...                     
  .ooooo....o.                            .oxxxoo....ooo....                   
    .oooooo...                              ...ooooo...ooo..                   
       ...                                     .oo.......                      
                                               ....ooo...                      
                      __              __                                       
 ______              [  |            [  |                 ._____  _          _ 
|_   _ \              | |  _          | |                /  ___ \| |        | |
  | |_) | .---.  _. _.| | / |  .---.  | | .---.  _    _ / /    \_|\ \  /\  / / 
  |  __'./ /__\\[ /`\_| '' <  / /__\\ | |/ /__\\| \  | | |   _____ \ \/  \/ /  
 _| |__| | \__. | |   | |`\ \ | \___. | || \___. \ \/ / \ \.___| |  \  /\  /   
|_______/ \.__./[_]  [__|  \_] \.__./[___]\.__./  \  /   \.____./    \/  \/    
                                                  / /                          
                                                 /_/                           
--------------------------------------------------------------------------------
 Please cite the following papers when using results from BerkeleyGW:

 Jack Deslippe, Georgy Samsonidze, David A. Strubbe, Manish Jain, Marvin L.
 Cohen, and Steven G. Louie, "BerkeleyGW: A Massively Parallel Computer Package
 for the Calculation of the Quasiparticle and Optical Properties of Materials
 and Nanostructures," Comput. Phys. Commun. 183, 1269 (2012)

 Mark S. Hybertsen and Steven G. Louie, "Electron correlation in semiconductors
 and insulators: Band gaps and quasiparticle energies," Phys. Rev. B 34, 5390
 (1986)

 Michael Rohlfing and Steven G. Louie, "Electron-hole excitations and optical
 spectra from first principles," Phys. Rev. B 62, 4927 (2000)
--------------------------------------------------------------------------------

 Running MPI version (parallel)
 Running with 2 MPI task(s)

 Compilation flags:
 - Compiler: GNU
 - Para. flags: MPI
 - Math  flags: USESCALAPACK, USEFFTW3
 - Debug flags: 

 We will start a new calculation from scratch.

 Experimental Features:
  - FFTs Optimization Level:  0
  - Frequency groups:     1
  - Use HDF5 Parallel IO:     F


 We`ll perform a uniform sampling of the full BZ.

 Computing the static inverse dielectric matrix

 We are using matrix communication scheme

 We are communicating via MPI

 We are using no truncation

 Reading header of WFN
 Highest occupied band (unshifted grid) = 4
 Valence max    (unshifted grid) = 4.692353 eV
 Conduction min (unshifted grid) = 6.963707 eV
 Middle energy  (unshifted grid) = 5.828030 eV
 Fermi  energy  (unshifted grid) = 5.828030 eV

 Reading header of WFNq
 Highest occupied band (shifted grid) = 4

 Cutoff of the dielectric matrix (Ry): 10.00
 Total number of bands in the calculation: 11
 Number of fully occupied valence bands: 4
 Number of partially occ. conduction bands: 0
 Monkhorst-Pack q-grid for epsilon(q): 2 2 2

 Running with 1 valence pools
 Number of conduction bands per processor: 4
 Number of valence bands per processor: 4

 Memory available:     236.4 MB per PE
 Memory required for execution:    2.6 MB per PE
 Memory required for vcoul:    0.0 MB per PE

 Number of electrons per unit cell (from ifmax) = 8.000000
 Number of electrons per unit cell (from occupations) = 8.000000
 Plasma Frequency = 1.148512 Ry


 Started reading wavefunctions (WFN) with 44 state(s) at 13:24:01.
 [ 13:24:01 |   0% ] state  1 / 44.
 [ 13:24:01 |   5% ] state  3 / 44, remaining: 0 s.
 [ 13:24:01 |   9% ] state  5 / 44, remaining: 0 s.
 [ 13:24:01 |  18% ] state  9 / 44, remaining: 0 s.
 [ 13:24:01 |  30% ] state 14 / 44, remaining: 0 s.
 [ 13:24:01 |  39% ] state 18 / 44, remaining: 0 s.
 [ 13:24:01 |  50% ] state 23 / 44, remaining: 0 s.
 [ 13:24:01 |  59% ] state 27 / 44, remaining: 0 s.
 [ 13:24:01 |  68% ] state 31 / 44, remaining: 0 s.
 [ 13:24:01 |  80% ] state 36 / 44, remaining: 0 s.
 [ 13:24:01 |  89% ] state 40 / 44, remaining: 0 s.
 Finished reading wavefunctions (WFN) at 13:24:01.
 Elapsed time: 0 s.

 Screened Coulomb cutoff: 10.000 Ry
 Total number of bands: 11
 Number of q-points: 3

 Scissors parameters
 Valence:    es =   0.0000 eV, e0 =   0.0000 eV, edel =   0.0000
 Conduction: es =   0.0000 eV, e0 =   0.0000 eV, edel =   0.0000
 Q0 and |Q0| =  0.001000  0.000000  0.000000  0.001019
 Q  and |Q|  =  0.000000  0.000000  0.500000  0.509324
 Q  and |Q|  =  0.000000  0.500000  0.500000  0.588117

 You have a slightly shifted q0 vector and a semiconductor.
 So reading from WFNq.
 Reading header of WFNq
 Highest occupied band (shifted grid) = 4

 Started reading wavefunctions (WFNq) with 40 state(s) at 13:24:01.
 [ 13:24:01 |   0% ] state  1 / 40.
 [ 13:24:01 |   5% ] state  3 / 40, remaining: 0 s.
 [ 13:24:01 |   8% ] state  4 / 40, remaining: 0 s.
 [ 13:24:01 |  18% ] state  8 / 40, remaining: 0 s.
 [ 13:24:01 |  28% ] state 12 / 40, remaining: 0 s.
 [ 13:24:01 |  38% ] state 16 / 40, remaining: 0 s.
 [ 13:24:01 |  48% ] state 20 / 40, remaining: 0 s.
 [ 13:24:01 |  57% ] state 24 / 40, remaining: 0 s.
 [ 13:24:01 |  68% ] state 28 / 40, remaining: 0 s.
 [ 13:24:01 |  78% ] state 32 / 40, remaining: 0 s.
 [ 13:24:01 |  88% ] state 36 / 40, remaining: 0 s.
 [ 13:24:01 |  98% ] state 40 / 40, remaining: 0 s.
 Finished reading wavefunctions (WFNq) at 13:24:01.
 Elapsed time: 0 s.


 Crystal k+q WFs read from WFNq
 Number of k-points in the full BZ of WFN (nfk): 32

================================================================================
 13:24:01   Dealing with q =  0.001000  0.000000  0.000000                 1 / 3
================================================================================

 This is the special q->0 point.
 Rank of the polarizability matrix (nmtx): 169

 BLACS processor grid:   1 x   2; BLOCKSIZE =   32

 Number of k-points in the irreducible BZ(q) (nrk): 10

 Started calculation of matrix elements with 40 transition(s) at 13:24:01.
 [ 13:24:01 |   0% ] transition  1 / 40.
 [ 13:24:01 |   5% ] transition  3 / 40, remaining: 0 s.
 [ 13:24:01 |   8% ] transition  4 / 40, remaining: 0 s.
 [ 13:24:01 |  18% ] transition  8 / 40, remaining: 0 s.
 [ 13:24:01 |  28% ] transition 12 / 40, remaining: 0 s.
 [ 13:24:01 |  38% ] transition 16 / 40, remaining: 0 s.
 [ 13:24:01 |  48% ] transition 20 / 40, remaining: 0 s.
 [ 13:24:01 |  57% ] transition 24 / 40, remaining: 0 s.
 [ 13:24:01 |  68% ] transition 28 / 40, remaining: 0 s.
 [ 13:24:01 |  78% ] transition 32 / 40, remaining: 0 s.
 [ 13:24:01 |  88% ] transition 36 / 40, remaining: 0 s.
 [ 13:24:01 |  98% ] transition 40 / 40, remaining: 0 s.
 Finished calculation of matrix elements at 13:24:01.
 Elapsed time: 0 s.

 Allocating Convergence Arrays of Dim           7

 Started building polarizability matrix with 2 processor(s) at 13:24:01.
 [ 13:24:01 |   0% ] processor 1 / 2.
 [ 13:24:01 |  50% ] processor 2 / 2.
 Finished building polarizability matrix at 13:24:01.
 Elapsed time: 0 s.

 q-pt      1: Head of Epsilon         =    2.495609926128694E+001   5.961619360513543E-016
 q-pt      1: Epsilon(2,2)            =    1.417554157439670E+000  -3.284739075644060E-019
 q-pt      1: Head of Epsilon Inverse =    4.450516507940471E-002  -9.901815153402771E-019
 q-pt      1: Epsilon Inverse(2,2)    =    7.320890262414345E-001  -1.084202172485504E-018

 For q0 points, you should check the symmetry (eps(G,G') = eps*(-G,-G')) by
 using the eps0sym code. Wavefunction convergence, as well as a finite q-shift
 may cause this property of eps(G,G') to be broken.

 Writing dielectric matrix to file

 Started writing matrix with 169 column(s) at 13:24:01.
 [ 13:24:01 |   0% ] column   1 / 169.
 [ 13:24:01 |   1% ] column   3 / 169, remaining: 0 s.
 [ 13:24:01 |   9% ] column  17 / 169, remaining: 0 s.
 [ 13:24:01 |  20% ] column  34 / 169, remaining: 0 s.
 [ 13:24:01 |  30% ] column  51 / 169, remaining: 0 s.
 [ 13:24:01 |  40% ] column  68 / 169, remaining: 0 s.
 [ 13:24:01 |  50% ] column  85 / 169, remaining: 0 s.
 [ 13:24:01 |  60% ] column 102 / 169, remaining: 0 s.
 [ 13:24:01 |  70% ] column 119 / 169, remaining: 0 s.
 [ 13:24:01 |  80% ] column 136 / 169, remaining: 0 s.
 [ 13:24:01 |  90% ] column 153 / 169, remaining: 0 s.
 [ 13:24:01 |  99% ] column 169 / 169, remaining: 0 s.
 Finished writing matrix at 13:24:01.
 Elapsed time: 0 s.

 Ok

================================================================================
 13:24:01   Dealing with q =  0.000000  0.000000  0.500000                 2 / 3
================================================================================

 This is a regular non-zero q-point.
 Rank of the polarizability matrix (nmtx): 168

 BLACS processor grid:   1 x   2; BLOCKSIZE =   32

 Number of k-points in the irreducible BZ(q) (nrk): 10

 Started calculation of matrix elements with 40 transition(s) at 13:24:01.
 [ 13:24:01 |   0% ] transition  1 / 40.
 [ 13:24:01 |   5% ] transition  3 / 40, remaining: 0 s.
 [ 13:24:01 |   8% ] transition  4 / 40, remaining: 0 s.
 [ 13:24:01 |  18% ] transition  8 / 40, remaining: 0 s.
 [ 13:24:01 |  28% ] transition 12 / 40, remaining: 0 s.
 [ 13:24:01 |  38% ] transition 16 / 40, remaining: 0 s.
 [ 13:24:01 |  48% ] transition 20 / 40, remaining: 0 s.
 [ 13:24:01 |  57% ] transition 24 / 40, remaining: 0 s.
 [ 13:24:01 |  68% ] transition 28 / 40, remaining: 0 s.
 [ 13:24:01 |  78% ] transition 32 / 40, remaining: 0 s.
 [ 13:24:01 |  88% ] transition 36 / 40, remaining: 0 s.
 [ 13:24:01 |  98% ] transition 40 / 40, remaining: 0 s.
 Finished calculation of matrix elements at 13:24:01.
 Elapsed time: 0 s.

 Allocating Convergence Arrays of Dim           7

 Started building polarizability matrix with 2 processor(s) at 13:24:01.
 [ 13:24:01 |   0% ] processor 1 / 2.
 [ 13:24:01 |  50% ] processor 2 / 2.
 Finished building polarizability matrix at 13:24:01.
 Elapsed time: 0 s.

 q-pt      2: Head of Epsilon         =    2.847228384837086E+000  -1.838229419226316E-017
 q-pt      2: Epsilon(2,2)            =    2.847228412553061E+000   2.293683583364979E-017
 q-pt      2: Head of Epsilon Inverse =    3.944778658243288E-001  -2.168404344971009E-019
 q-pt      2: Epsilon Inverse(2,2)    =    3.944778624167692E-001  -3.903127820947816E-018

 Writing dielectric matrix to file

 Started writing matrix with 168 column(s) at 13:24:01.
 [ 13:24:01 |   0% ] column   1 / 168.
 [ 13:24:01 |   1% ] column   3 / 168, remaining: 0 s.
 [ 13:24:01 |  10% ] column  17 / 168, remaining: 0 s.
 [ 13:24:01 |  20% ] column  34 / 168, remaining: 0 s.
 [ 13:24:01 |  30% ] column  51 / 168, remaining: 0 s.
 [ 13:24:01 |  40% ] column  68 / 168, remaining: 0 s.
 [ 13:24:01 |  50% ] column  85 / 168, remaining: 0 s.
 [ 13:24:01 |  60% ] column 101 / 168, remaining: 0 s.
 [ 13:24:01 |  70% ] column 118 / 168, remaining: 0 s.
 [ 13:24:01 |  80% ] column 135 / 168, remaining: 0 s.
 [ 13:24:01 |  90% ] column 152 / 168, remaining: 0 s.
 Finished writing matrix at 13:24:01.
 Elapsed time: 0 s.

 Ok

================================================================================
 13:24:01   Dealing with q =  0.000000  0.500000  0.500000                 3 / 3
================================================================================

 This is a regular non-zero q-point.
 Rank of the polarizability matrix (nmtx): 150

 BLACS processor grid:   1 x   2; BLOCKSIZE =   32

 Number of k-points in the irreducible BZ(q) (nrk): 6

 Started calculation of matrix elements with 24 transition(s) at 13:24:01.
 [ 13:24:01 |   0% ] transition  1 / 24.
 [ 13:24:01 |   8% ] transition  3 / 24, remaining: 0 s.
 [ 13:24:01 |  17% ] transition  5 / 24, remaining: 0 s.
 [ 13:24:01 |  29% ] transition  8 / 24, remaining: 0 s.
 [ 13:24:01 |  38% ] transition 10 / 24, remaining: 0 s.
 [ 13:24:01 |  46% ] transition 12 / 24, remaining: 0 s.
 [ 13:24:01 |  58% ] transition 15 / 24, remaining: 0 s.
 [ 13:24:01 |  67% ] transition 17 / 24, remaining: 0 s.
 [ 13:24:01 |  79% ] transition 20 / 24, remaining: 0 s.
 [ 13:24:01 |  88% ] transition 22 / 24, remaining: 0 s.
 [ 13:24:01 |  96% ] transition 24 / 24, remaining: 0 s.
 Finished calculation of matrix elements at 13:24:01.
 Elapsed time: 0 s.

 Allocating Convergence Arrays of Dim           7

 Started building polarizability matrix with 2 processor(s) at 13:24:01.
 [ 13:24:01 |   0% ] processor 1 / 2.
 [ 13:24:01 |  50% ] processor 2 / 2.
 Finished building polarizability matrix at 13:24:01.
 Elapsed time: 0 s.

 q-pt      3: Head of Epsilon         =    2.643965054476260E+000   3.200488720974389E-018
 q-pt      3: Epsilon(2,2)            =    2.643965054476260E+000   7.878126082398496E-018
 q-pt      3: Head of Epsilon Inverse =    4.114138579427453E-001  -1.698412340886655E-019
 q-pt      3: Epsilon Inverse(2,2)    =    4.114138579427452E-001  -2.481485059464850E-018

 Writing dielectric matrix to file

 Started writing matrix with 150 column(s) at 13:24:01.
 [ 13:24:01 |   0% ] column   1 / 150.
 [ 13:24:01 |   1% ] column   3 / 150, remaining: 0 s.
 [ 13:24:01 |   9% ] column  15 / 150, remaining: 0 s.
 [ 13:24:01 |  19% ] column  30 / 150, remaining: 0 s.
 [ 13:24:01 |  29% ] column  45 / 150, remaining: 0 s.
 [ 13:24:01 |  39% ] column  60 / 150, remaining: 0 s.
 [ 13:24:01 |  49% ] column  75 / 150, remaining: 0 s.
 [ 13:24:01 |  59% ] column  90 / 150, remaining: 0 s.
 [ 13:24:01 |  69% ] column 105 / 150, remaining: 0 s.
 [ 13:24:01 |  79% ] column 120 / 150, remaining: 0 s.
 [ 13:24:01 |  89% ] column 135 / 150, remaining: 0 s.
 [ 13:24:01 |  99% ] column 150 / 150, remaining: 0 s.
 Finished writing matrix at 13:24:01.
 Elapsed time: 0 s.

 Ok


                              CPU (s)        WALL (s)           #

FULLBZ:                         0.000           0.000           1
IRRBZ:                          0.000           0.000           3
SUBGRP:                         0.000           0.000           3
RQSTAR:                         0.000           0.000          26
GMAP:                           0.001           0.002          26
GVEC:                           0.000           0.000           3
INPUT:                          0.005           0.012           1
INPUT_Q:                        0.001           0.001           1
INPUT(Q) I/O                    0.001           0.000         104
INPUT(Q) COMM                   0.001           0.001         104
GENWF:                          0.006           0.003         104
GENWF (VAL):                    0.003           0.002         104
GENWF (COND):                   0.003           0.001          26
GENWF (C-Ekin):                 0.000           0.001          26
GENWF (C-Sort):                 0.002           0.000          26
JOB SETUP:                      0.000           0.001           1
Q LOOP SETUP:                   0.000           0.000           3
INIT CUTOFF:                    0.000           0.001           3
INIT SCALAPACK:                 0.001           0.001           3
INIT ARRAYS:                    0.000           0.000           3
MTXEL:                          0.067           0.079         104
MTXEL (DENOM):                  0.000           0.000         416
MTXEL (FFT):                    0.067           0.079         104
MTXEL EXP(DENOM):               0.000           0.000           0
MTXEL EXP (FFT):                0.000           0.000           0
CHI SUM (TOTAL):                0.044           0.018           3
CHI SUM (COMM):                 0.001           0.002           6
CHI SUM (PREP):                 0.003           0.003           6
CHI SUM (ARRAY ALLOC):          0.000           0.000           0
CHI SUM (ROW):                  0.000           0.000           0
CHI SUM (COLUMN):               0.000           0.000           0
CHI SUM (zGEMM):                0.040           0.012           6
CHI SUM (HT):                   0.000           0.000           0
EPSINV (TOTAL):                 0.025           0.019           3
EPSINV (VCOUL):                 0.000           0.000           3
EPSINV (I/O)                    0.004           0.003           3
EPSINV (INVERT)                 0.018           0.013           3
EPS (I/O) COMM                  0.000           0.000           0
EPS (I/O) IO                    0.000           0.000           0
CONVERGE TESTS:                 0.000           0.001           3
OPT FFT:                        0.000           0.000           0
OPT FFT (INIT):                 0.000           0.000           0
OPT FFT (COMM_FFT):             0.000           0.000           0
OPT FFT (FFT):                  0.000           0.000           0
FFT ZERO                        0.002           0.001         520
FFT PUT                         0.003           0.002         520
FFT PLAN                        0.035           0.049         936
FFT EXEC                        0.017           0.018         936
FFT MLTPLY                      0.003           0.004         416

TOTAL:                          0.154           0.140


