Re: mpb-discuss Digest, Vol 60, Issue 7 (2D Photonic Crystal Homogeneous Case)
[email protected] 22 Feb 2011 12:04:54 -0000
| Newsgroups | gmane.comp.science.photonic-bands |
|---|---|
| Message-ID | <[email protected]> |
Thanks Mahmoud Hamoud for your preciuos advice, now everything seems to be working fine.
But I have another problem.
I tried this code:
(set! num-bands 8)
(set! k-points (list (vector3 0 0 0)(vector3 0.5 0 0)(vector3 0.5 0.5 0)(vector3 0 0 0)))
(set! k-points (interpolate 4 k-points))
(set! geometry (list (make block (center 0 0 0)(material air)(size 1 1 0) ) ))
(set! geometry-lattice (make lattice (size 1 1 no-size)))
(set! resolution 32)
(run-tm)
to compute the eigenvalues of the homogeneous medium, since in that case comparisons with analytical results can be made. So, at the end of the code I tried first (run-tm), then (run-te) and finally (run), expecting the numerical results to be always the same.Here's what I got:
1) The (run-tm) and (run-te) results are the same except for the 6-th eigenvalues (computed for k(0.5,0,0) ):
tmfreqs:, 6, 0.5, 0, 0, 0.5, 0.5, 0.539878, 1.12469, 1.13339, 1.13849, 1.14212, 1.5132, 1.68302
tefreqs:, 6, 0.5, 0, 0, 0.5, 0.5, 0.5, 1.11803, 1.11803, 1.11803, 1.11803, 1.5, 1.80278
2) The (run) results are different from both (run-te) and (run-tm) results.
To be more complete I paste the output files of the three computations:
init-params: initializing eigensolver data
Computing 8 bands with 1.000000e-07 tolerance.
Working in 2 dimensions.
Grid size is 32 x 32 x 1.
Solving for 8 bands at a time.
Creating Maxwell data...
Allocating fields...
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:
block, center = (0,0,0)
size (1,1,0)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon = 1
Geometric object tree has depth 3 and 16 object nodes (vs. 1 actual objects)
Initializing dielectric 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-1, mean 1, harm. mean 1, 0% > 1, 100% "fill"
matrixio: cannot output "epsilon" (compiled without HDF)
Outputting 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 8 iterations.
tmfreqs:, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1.41421, 1.41421, 1.41421
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 2 iterations.
tmfreqs:, 2, 0.1, 0, 0, 0.1, 0.1, 0.9, 1.00499, 1.00499, 1.1, 1.34536, 1.34536, 1.48661
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 1 iterations.
tmfreqs:, 3, 0.2, 0, 0, 0.2, 0.2, 0.8, 1.0198, 1.0198, 1.2, 1.28062, 1.28062, 1.56205
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 1 iterations.
tmfreqs:, 4, 0.3, 0, 0, 0.3, 0.3, 0.7, 1.04403, 1.04403, 1.22066, 1.22066, 1.3, 1.64012
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 1 iterations.
tmfreqs:, 5, 0.4, 0, 0, 0.4, 0.4, 0.6, 1.07703, 1.07703, 1.16619, 1.16619, 1.4, 1.72047
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 2 iterations.
tmfreqs:, 6, 0.5, 0, 0, 0.5, 0.5, 0.539878, 1.12469, 1.13339, 1.13849, 1.14212, 1.5132, 1.68302
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 7 iterations.
tmfreqs:, 7, 0.5, 0.1, 0, 0.509902, 0.509902, 0.509902, 1.02956, 1.02956, 1.2083, 1.2083, 1.50333, 1.50333
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 1 iterations.
tmfreqs:, 8, 0.5, 0.2, 0, 0.538516, 0.538516, 0.538516, 0.943398, 0.943398, 1.3, 1.3, 1.51327, 1.51327
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 1 iterations.
tmfreqs:, 9, 0.5, 0.3, 0, 0.583095, 0.583095, 0.583095, 0.860233, 0.860233, 1.39284, 1.39284, 1.52971, 1.52971
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 1 iterations.
tmfreqs:, 10, 0.5, 0.4, 0, 0.640312, 0.640312, 0.640312, 0.781025, 0.781025, 1.48661, 1.48661, 1.55242, 1.55242
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 1 iterations.
tmfreqs:, 11, 0.5, 0.5, 0, 0.707107, 0.707107, 0.707107, 0.707107, 0.707107, 1.58114, 1.58114, 1.58114, 1.58114
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 2 iterations.
tmfreqs:, 12, 0.4, 0.4, 0, 0.565685, 0.565685, 0.721571, 0.722602, 0.849641, 1.45687, 1.46972, 1.52556, 1.53205
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 6 iterations.
tmfreqs:, 13, 0.3, 0.3, 0, 0.424264, 0.424264, 0.761577, 0.761577, 0.989949, 1.33417, 1.33417, 1.47648, 1.47648
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 1 iterations.
tmfreqs:, 14, 0.2, 0.2, 0, 0.282843, 0.282843, 0.824621, 0.824621, 1.13137, 1.21655, 1.21655, 1.44222, 1.44222
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 1 iterations.
tmfreqs:, 15, 0.1, 0.1, 0, 0.141421, 0.141421, 0.905539, 0.905539, 1.10454, 1.10454, 1.27279, 1.42127, 1.42127
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 1 iterations.
tmfreqs:, 16, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1.41421, 1.41421, 1.41421
elapsed time for k point: 0 seconds.
Band 1 range: 0.0 at #(0 0 0) to 0.707106781186548 at #(0.5 0.5 0)
Band 2 range: 0.509901951579547 at #(0.5 0.1 0) to 1.00000000003675 at #(0 0 0)
Band 3 range: 0.70710678118656 at #(0.5 0.5 0) to 1.124692730729 at #(0.5 0 0)
Band 4 range: 0.707106781186561 at #(0.5 0.5 0) to 1.13339282370176 at #(0.5 0 0)
Band 5 range: 1.00000000001303 at #(0 0 0) to 1.58113883008419 at #(0.5 0.5 0)
Band 6 range: 1.14212126230508 at #(0.5 0 0) to 1.58113883008419 at #(0.5 0.5 0)
Band 7 range: 1.28062484752904 at #(0.2 0 0) to 1.5811388300842 at #(0.5 0.5 0)
Band 8 range: 1.4142135623858 at #(0 0 0) to 1.72046505339032 at #(0.4 0 0)
total elapsed time for run: 0 seconds.
done.
Some deprecated features have been used. Set the environment
variable GUILE_WARN_DEPRECATED to "detailed" and rerun the
program to get more information. Set it to "no" to suppress
this message.
init-params: initializing eigensolver data
Computing 8 bands with 1.000000e-07 tolerance.
Working in 2 dimensions.
Grid size is 32 x 32 x 1.
Solving for 8 bands at a time.
Creating Maxwell data...
Allocating fields...
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:
block, center = (0,0,0)
size (1,1,0)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon = 1
Geometric object tree has depth 3 and 16 object nodes (vs. 1 actual objects)
Initializing dielectric 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: te.
Initializing fields to random numbers...
elapsed time for initialization: 0 seconds.
epsilon: 1-1, mean 1, harm. mean 1, 0% > 1, 100% "fill"
matrixio: cannot output "epsilon" (compiled without HDF)
Outputting 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 8 iterations.
tefreqs:, 1, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1.41421, 1.41421, 1.41421
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 2 iterations.
tefreqs:, 2, 0.1, 0, 0, 0.1, 0.1, 0.9, 1.00499, 1.00499, 1.1, 1.34536, 1.34536, 1.48661
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 1 iterations.
tefreqs:, 3, 0.2, 0, 0, 0.2, 0.2, 0.8, 1.0198, 1.0198, 1.2, 1.28062, 1.28062, 1.56205
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 1 iterations.
tefreqs:, 4, 0.3, 0, 0, 0.3, 0.3, 0.7, 1.04403, 1.04403, 1.22066, 1.22066, 1.3, 1.64012
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 1 iterations.
tefreqs:, 5, 0.4, 0, 0, 0.4, 0.4, 0.6, 1.07703, 1.07703, 1.16619, 1.16619, 1.4, 1.72047
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 1 iterations.
tefreqs:, 6, 0.5, 0, 0, 0.5, 0.5, 0.5, 1.11803, 1.11803, 1.11803, 1.11803, 1.5, 1.80278
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 76 iterations.
tefreqs:, 7, 0.5, 0.1, 0, 0.509902, 0.509902, 0.509902, 1.02956, 1.02956, 1.2083, 1.2083, 1.50333, 1.50333
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 1 iterations.
tefreqs:, 8, 0.5, 0.2, 0, 0.538516, 0.538516, 0.538516, 0.943398, 0.943398, 1.3, 1.3, 1.51327, 1.51327
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 1 iterations.
tefreqs:, 9, 0.5, 0.3, 0, 0.583095, 0.583095, 0.583095, 0.860233, 0.860233, 1.39284, 1.39284, 1.52971, 1.52971
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 1 iterations.
tefreqs:, 10, 0.5, 0.4, 0, 0.640312, 0.640312, 0.640312, 0.781025, 0.781025, 1.48661, 1.48661, 1.55242, 1.55242
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 1 iterations.
tefreqs:, 11, 0.5, 0.5, 0, 0.707107, 0.707107, 0.707107, 0.707107, 0.707107, 1.58114, 1.58114, 1.58114, 1.58114
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 2 iterations.
tefreqs:, 12, 0.4, 0.4, 0, 0.565685, 0.565685, 0.721901, 0.722132, 0.849091, 1.45622, 1.47371, 1.52382, 1.54419
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 6 iterations.
tefreqs:, 13, 0.3, 0.3, 0, 0.424264, 0.424264, 0.761577, 0.761577, 0.989949, 1.33417, 1.33417, 1.47648, 1.47648
elapsed time for k point: 1 seconds.
solve_kpoint (0.2,0.2,0):
Solving for bands 1 to 8...
Finished solving for bands 1 to 8 after 1 iterations.
tefreqs:, 14, 0.2, 0.2, 0, 0.282843, 0.282843, 0.824621, 0.824621, 1.13137, 1.21655, 1.21655, 1.44222, 1.44222
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 1 iterations.
tefreqs:, 15, 0.1, 0.1, 0, 0.141421, 0.141421, 0.905539, 0.905539, 1.10454, 1.10454, 1.27279, 1.42127, 1.42127
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 1 iterations.
tefreqs:, 16, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1.41421, 1.41421, 1.41421
elapsed time for k point: 0 seconds.
Band 1 range: 0.0 at #(0 0 0) to 0.707106781186547 at #(0.5 0.5 0)
Band 2 range: 0.500000000701684 at #(0.5 0 0) to 1.00000000045515 at #(0 0 0)
Band 3 range: 0.707106781186552 at #(0.5 0.5 0) to 1.11803398890731 at #(0.5 0 0)
Band 4 range: 0.707106781186553 at #(0.5 0.5 0) to 1.13137085003665 at #(0.2 0.2 0)
Band 5 range: 1.00000000000667 at #(0 0 0) to 1.58113883008419 at #(0.5 0.5 0)
Band 6 range: 1.11803398909198 at #(0.5 0 0) to 1.58113883008419 at #(0.5 0.5 0)
Band 7 range: 1.28062484762186 at #(0.2 0 0) to 1.58113883008419 at #(0.5 0.5 0)
Band 8 range: 1.41421356238112 at #(0 0 0) to 1.80277563647803 at #(0.5 0 0)
total elapsed time for run: 1 seconds.
done.
Some deprecated features have been used. Set the environment
variable GUILE_WARN_DEPRECATED to "detailed" and rerun the
program to get more information. Set it to "no" to suppress
this message.
init-params: initializing eigensolver data
Computing 8 bands with 1.000000e-07 tolerance.
Working in 2 dimensions.
Grid size is 32 x 32 x 1.
Solving for 8 bands at a time.
Creating Maxwell data...
Allocating fields...
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:
block, center = (0,0,0)
size (1,1,0)
axes (1,0,0), (0,1,0), (0,0,1)
dielectric constant epsilon = 1
Geometric object tree has depth 3 and 16 object nodes (vs. 1 actual objects)
Initializing dielectric 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: .
Initializing fields to random numbers...
elapsed time for initialization: 0 seconds.
epsilon: 1-1, mean 1, harm. mean 1, 0% > 1, 100% "fill"
matrixio: cannot output "epsilon" (compiled without HDF)
Outputting epsilon...
solve_kpoint (0,0,0):
freqs:, k index, k1, k2, k3, kmag/2pi, band 1, band 2, band 3, band 4, band 5, band 6, band 7, band 8
Solving for bands 3 to 8...
Finished solving for bands 3 to 8 after 8 iterations.
freqs:, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1
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 12 iterations.
freqs:, 2, 0.1, 0, 0, 0.1, 0.1, 0.1, 0.9, 0.9, 1.00499, 1.00499, 1.00499, 1.00499
elapsed time for k point: 1 seconds.
solve_kpoint (0.2,0,0):
Solving for bands 1 to 8...
Finished solving for bands 1 to 8 after 1 iterations.
freqs:, 3, 0.2, 0, 0, 0.2, 0.2, 0.2, 0.8, 0.8, 1.0198, 1.0198, 1.0198, 1.0198
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 1 iterations.
freqs:, 4, 0.3, 0, 0, 0.3, 0.3, 0.3, 0.7, 0.7, 1.04403, 1.04403, 1.04403, 1.04403
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 1 iterations.
freqs:, 5, 0.4, 0, 0, 0.4, 0.4, 0.4, 0.6, 0.6, 1.07703, 1.07703, 1.07703, 1.07703
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 1 iterations.
freqs:, 6, 0.5, 0, 0, 0.5, 0.5, 0.5, 0.5, 0.5, 1.11803, 1.11803, 1.11803, 1.11803
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 1 iterations.
freqs:, 7, 0.5, 0.1, 0, 0.509902, 0.509902, 0.509902, 0.509902, 0.509902, 1.02956, 1.02956, 1.2083, 1.2083
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 1 iterations.
freqs:, 8, 0.5, 0.2, 0, 0.538516, 0.538516, 0.538516, 0.538516, 0.538516, 0.943398, 0.943398, 1.3, 1.3
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 1 iterations.
freqs:, 9, 0.5, 0.3, 0, 0.583095, 0.583095, 0.583095, 0.583095, 0.583095, 0.860233, 0.860233, 1.39284, 1.39284
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 2 iterations.
freqs:, 10, 0.5, 0.4, 0, 0.640312, 0.640312, 0.640312, 0.667147, 0.70475, 0.812503, 0.8656, 1.13834, 1.50271
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 9 iterations.
freqs:, 11, 0.5, 0.5, 0, 0.707107, 0.707107, 0.707107, 0.707107, 0.707107, 0.707107, 0.707107, 0.707107, 0.707107
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 1 iterations.
freqs:, 12, 0.4, 0.4, 0, 0.565685, 0.565685, 0.565685, 0.72111, 0.72111, 0.72111, 0.72111, 0.848528, 0.848528
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 1 iterations.
freqs:, 13, 0.3, 0.3, 0, 0.424264, 0.424264, 0.424264, 0.761577, 0.761577, 0.761577, 0.761577, 0.989949, 0.989949
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 1 iterations.
freqs:, 14, 0.2, 0.2, 0, 0.282843, 0.282843, 0.282843, 0.824621, 0.824621, 0.824621, 0.824621, 1.13137, 1.13137
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 2 iterations.
freqs:, 15, 0.1, 0.1, 0, 0.141421, 0.141421, 0.141421, 0.919635, 0.922099, 0.924578, 0.925376, 1.13347, 1.14718
elapsed time for k point: 0 seconds.
solve_kpoint (0,0,0):
Solving for bands 3 to 8...
Finished solving for bands 3 to 8 after 5 iterations.
freqs:, 16, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1
elapsed time for k point: 0 seconds.
Band 1 range: 0.0 at #(0 0 0) to 0.707106781186548 at #(0.5 0.5 0)
Band 2 range: 0.0 at #(0 0 0) to 0.707106781186548 at #(0.5 0.5 0)
Band 3 range: 0.500000000005301 at #(0.5 0 0) to 1.00000000036222 at #(0 0 0)
Band 4 range: 0.500000000012381 at #(0.5 0 0) to 1.00000000042076 at #(0 0 0)
Band 5 range: 0.707106781794994 at #(0.5 0.5 0) to 1.11803398875865 at #(0.5 0 0)
Band 6 range: 0.707106781859084 at #(0.5 0.5 0) to 1.11803398876086 at #(0.5 0 0)
Band 7 range: 0.707106781868582 at #(0.5 0.5 0) to 1.39283882771848 at #(0.5 0.3 0)
Band 8 range: 0.707106782170754 at #(0.5 0.5 0) to 1.50271314409492 at #(0.5 0.4 0)
total elapsed time for run: 1 seconds.
done.
Some deprecated features have been used. Set the environment
variable GUILE_WARN_DEPRECATED to "detailed" and rerun the
program to get more information. Set it to "no" to suppress
this message.
-------------------------------------------------------------------------------
Valore legale alle tue mail
InterfreePEC - la tua Posta Elettronica Certificata
http://pec.interfree.it
-------------------------------------------------------------------------------