notes 71.6 KB
Newer Older
Francois Gygi committed
1
--------------------------------------------------------------------------------
Francois Gygi committed
2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
rel1_71_1
--------------------------------------------------------------------------------
0b54b01 Merge branch 'develop'
04bf383 Fix timing msg in Bisection.C
8ad14a3 Fix bisection timer for nspin==2
4b2de68 Fix constness in BasisMapping
1a2fb01 Enable spin-polarized CP dynamics
2aad1f5 Add cstdlib header for abs()
2d6f542 Use vector swap in BasisMapping for single-node
e9d9631 Use memset to zero arrays in BasisMapping
90f6ddc Update gather/scatter in BasisMapping for MKL
4f1caa4 Fix timers in testFourierTransform
2311fdd Add fxy and fz functions in FourierTransform
984676c Cleanup FourierTransform and BasisMapping
8dc7100 Rewrite BasisMapping and FourierTransform
57f8d91 Switch names of transpose fns in BasisMapping
--------------------------------------------------------------------------------
Francois Gygi committed
19 20 21 22 23
rel1_71_0
--------------------------------------------------------------------------------
a26cbfa Update constraints before ionic step
88181a8 Change sign of force on constraints and pair force
--------------------------------------------------------------------------------
Francois Gygi committed
24 25 26 27 28 29 30
rel1_70_0
--------------------------------------------------------------------------------
c790914 Fix xsd links in upf2qso dir
2348554 Update test/Makefile
785af78 Add test/siatom test of custom occupation numbers
0af0f85 Add occ variable for custom occupation numbers
--------------------------------------------------------------------------------
Francois Gygi committed
31 32 33 34 35 36
rel1_69_1
--------------------------------------------------------------------------------
cb15c6c Fix output formatting in SpectrumCmd.C
94915e6 Fix gaussian normalization in SpectrumCmd.C
9f4d98f Remove tabs in ExchangeOperator.C
--------------------------------------------------------------------------------
Francois Gygi committed
37 38 39 40 41 42 43 44 45 46 47 48 49
rel1_69_0
--------------------------------------------------------------------------------
beb8bee Cleanup trailing spaces in util scripts
3137cbd Add scripts qbox_reduce.sh, qbox_move_subsample.py
8d434af Fix spectrum file for nspin==2
7a00811 Update spectrum command
6ba92e5 Update spectrum command help msg
80c343d Initialize eig with zero in SlaterDet
4031b74 Modify spectrum command parameters
539499d write spectrum on file, fix parallel SpectrumCmd
54058fe Add spectrum cmd to Makefile and qb.C
fc34625 Add spectrum command
--------------------------------------------------------------------------------
Francois Gygi committed
50 51 52 53
rel1_68_4
--------------------------------------------------------------------------------
d9bac45 Use syevd/heevd in JD stepper, workaround lwork size calculation
--------------------------------------------------------------------------------
Francois Gygi committed
54 55 56 57
rel1_68_3
--------------------------------------------------------------------------------
41edad7 Fix issue fail to restart MD with multiple k-points
--------------------------------------------------------------------------------
Francois Gygi committed
58 59 60 61 62 63 64
rel1_68_2
--------------------------------------------------------------------------------
260b6be Remove debug print in ExchangeOperator.C
0d66277 Fix Matrix.C calc of lwork, return values
a0ff3b7 Fix D3tensor.h
39fd798 Changed default to all atomsets in qbox_xyz.py
--------------------------------------------------------------------------------
Francois Gygi committed
65 66 67 68 69 70 71 72
rel1_68_1
--------------------------------------------------------------------------------
b2326e0 Remove unneeded loop in XCPotential.C
6ce2309 Fix white space, reorder tests in SCANFunctional.C
186b7a8 Cleanup ExchangeOperator functions params and return values
e1a0d16 Modified SCANFunctional to not crash with grad=0, or negative density
52c16a2 Set echo=true in qb.C to echo commands in server mode
--------------------------------------------------------------------------------
Francois Gygi committed
73 74 75 76 77 78 79 80
rel1_68_0
--------------------------------------------------------------------------------
05cd227 Merge branch 'scan' into develop
5c8ff65 Make Sample const in ExchangeOperator
15649b6 Make nlevels const in Bisection
8afc2a0 Add qbox_species_temp.sh utility script
4f46485 Add scripts to extract position,velocity,force
--------------------------------------------------------------------------------
Francois Gygi committed
81 82 83 84 85
rel1_67_4
--------------------------------------------------------------------------------
6e23b88 Merge branch 'develop'
c1261b7 Fix units of stress_tol test in BOSampleStepper
--------------------------------------------------------------------------------
Francois Gygi committed
86 87 88 89 90
rel1_67_3
--------------------------------------------------------------------------------
48c0d4b Merge branch 'develop'
211d894 Fix Harris-Foulkes estimate for NLCC when nspin=2
--------------------------------------------------------------------------------
Francois Gygi committed
91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106
rel1_67_2
--------------------------------------------------------------------------------
9459eed Merge branch 'develop'
1d85485 Use daxpy in XCPotential loops
136f424 Fix total energy and stress calc with NLCC
4e6d0f5 Fix non-linear core correction
c59340b Add rhocore_r and rhocore_g to ChargeDensity
d010476 Fix initialization of core correction in Species.C
593b1d1 update release string to 1.9
ddcb699 Restore 1.0 coefficient of non-local projectors to follow QE definitions
715e140 Add fast return in NonLocalPotential if no non-local species
72caf01 Add factor 0.5 in upf2qso.C for consistency of energy units
a52c966 Add <algorithm> header in upf2qso.C for g++ compilation
797585d Added ISO date stamp in upf2qso output
4f93df7 Fix upf2qso.C to work with UPF 2 pseudopotentials
--------------------------------------------------------------------------------
Francois Gygi committed
107 108 109 110
rel1_67_1
--------------------------------------------------------------------------------
00fb07e Add missing files ForceTol.h and StressTol.h
--------------------------------------------------------------------------------
Francois Gygi committed
111 112 113 114 115 116 117 118 119 120 121 122
rel1_67_0
--------------------------------------------------------------------------------
c1bad59 Restore scaling of alpha_start in CGIonicStepper.C
8ab5bfa Fix force_tol and stress_tol convergence test in BOSampleStepper.C
a7c7157 Fix test convergence test of force_tol and stress_tol
61ac1c1 Add optional VERSION macro in qb.C and Makefile
0ed5ed8 Add force_tol and stress_tol variables
ef0e7eb Adjust params in CGIonicStepper, remove rescaling of alpha_start
4a664cb Change beta_max in testCGOptimizer
66305fd Add set_nstep_max function to LineMinimizer
e800eeb Fix PlotCmd argument parsing for -vlocal with nspin=2
--------------------------------------------------------------------------------
Francois Gygi committed
123 124 125 126 127 128 129 130 131 132 133
rel1_66_3
--------------------------------------------------------------------------------
23da694 Adjust atomic radii. Update logic of non-scf loop termination.
24ec148 Fix pure virtual functions declarations in SampleStepper.h
0cae884 Fix MPI_Comm leak in SlaterDet copy ctor
ffa583a Remove old kpgen.C, replaced by util/kpgen/kpgen.cpp
2e5500e Update util/kpgen files
c67d4c3 Fix units of ext_stress
648574e Update mk file for MacBook Pro Mojave
d51cdb9 Fix precision of eigenvalue_sum output
--------------------------------------------------------------------------------
Francois Gygi committed
134 135 136 137 138
rel1_66_2
--------------------------------------------------------------------------------
2fdf056 Fix Function3d byte swapping on big endian systems
f31f37e Fix check of argc in Vext.h
--------------------------------------------------------------------------------
Francois Gygi committed
139 140 141 142
rel1_66_1
--------------------------------------------------------------------------------
2d07308 Add cstring header in FourierTransform.C
--------------------------------------------------------------------------------
Francois Gygi committed
143 144 145 146 147
rel1_66_0
--------------------------------------------------------------------------------
124fb28 Fix test causing fail on assert in FCC cell
Merge vext branch
--------------------------------------------------------------------------------
Francois Gygi committed
148 149 150 151 152 153 154 155 156 157 158 159 160 161 162
rel1_65_1
--------------------------------------------------------------------------------
a23e645 Remove static_cast and cleanup RSH and HSE functionals
de21bb8 Cleanup test programs
a312251 Add error msg in BisectionCmd if nkp > 1
cd77e32 Print np0v,np1v,np2v as XML elements in EnergyFunctional
5abca0d Compute and print np0v,np1v,np2v in StatusCmd
dc43a96 Add README and COPYING files
0827954 Rewrite and simplify UserInterface
218356e Update Makefile dependencies
df25ed5 Fix missing continuation character
58547ce Add basis_fits_in_grid function
e66b7c0 Replace zeroing loops with memset in FourierTransform.C
5af642d Add stress tests to test Makefile
--------------------------------------------------------------------------------
Francois Gygi committed
163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178
rel1_65_0
--------------------------------------------------------------------------------
fb22565 Merge branch 'develop'
939c153 Add Mac OS makefile
3ffc3f0 Fix -atomic_charge initialization when net_charge!=0
e997363 Fix UserInterface processing of empty lines
ca27b6e Add centos7 and centos7_intel makefiles
707eadd Use StrX in XML handler classes to avoid mem leaks
c4e49cf Fixed calls to align wfs in BOSampleStepper.C
0436991 Add BHandHLYP hybrid functional
f3b32ed Add qbox_move_to.py script
5298a69 Added qso_plot_species.py plot utility
fbf9477 Removed unused code from HSE and RSH functionals
75919ee Test negative density in HSE and RSH functionals
c940023 Fix text-form in SlaterDet::write
--------------------------------------------------------------------------------
Francois Gygi committed
179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197
rel1_64_0
--------------------------------------------------------------------------------
d4c5c81 Merge branch 'develop'
28c1fc2 Merge branch 'rsh-dev' into develop
38843c5 Removed unused variables
52cb12a Implemented stress for general RSH case
211d90a adjust blocking factors of square matrices
7c62abc Implement reset_rotation command
0b1e4fa fix help msg in MoveCmd.h
84416bf Reorganized stress tests as si4stress and si2stress
7e0b69a Added tests for invalid bisection parameters
5f7c1c7 Adjust strain values in test/stress files
0d593a2 Cleanup files in test/stress dir
9c5654f Add test of stress tensor
6f049b7 Fixed x_coeff in RSHFunctional.C Cleanup & comment
af533ce Fix errors in ExchangeOperator at gamma
5e24b92 Cleanup ExchangeOperator
ccd020c First RSH implementation
52e6976 Merged HSE implementation
198 199 200 201 202 203 204 205 206
fbb7421 Merge branch 'hse-dev' into develop
5b93ad8 Add cstring header in ConstraintSet.C
72e55a8 Fix force normalization in PairExtForce::energy
000b3db Fixed charge integral in ChargeDensity::update_rhor
80f2730 Fix long lines
1c14353 Merge branch 'develop' into hse-dev
e9d66bf Add zero volume test to run command
667cce0 Updated centos6.mk
f9b90ad Updated notes
Francois Gygi committed
207
--------------------------------------------------------------------------------
Francois Gygi committed
208 209 210 211 212
rel1_63_11
--------------------------------------------------------------------------------
be7d4cc rel1_63_11
4e6e1b6 adjust blocking factors of square matrices
--------------------------------------------------------------------------------
Francois Gygi committed
213 214 215 216 217 218 219 220 221 222
rel1_63_10
--------------------------------------------------------------------------------
058730d Implemented PSDA for complex wave functions
26c3495 Update printing of electronic charge & net charge
c102a32 Add testSampleReader program
8f85701 Add option to define constraint with current value
35f2374 changed output format of constraint value, force
e85f85d Changed logic of shifts: zero shifts include Gamma
25f2f96 Generate xyz file from kpoint file
718452c Changed bisection tolerance from 1.0 to 0.001
Francois Gygi committed
223
--------------------------------------------------------------------------------
Francois Gygi committed
224 225
rel1_63_9
--------------------------------------------------------------------------------
Francois Gygi committed
226
5310d81 Fix discontinuity in torsion constraint function
Francois Gygi committed
227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244
ee9efe5 Fix test for negative scf_tol value
f47d0c1 New implementation of kpgen
86be185 Added kpoint move functionality
dd1a11b Merge branch 'upf2' into develop
45aea94 Added NLCC functionality
ac8d7e6 Implemented NL potential in UPF 2 format
ff8a484 Implemented UPF 2 format for SL potentials
feb88ae Fix end of line in rseed help msg.
6309a0a Removed enthalpy_int.plt
e19b136 Add util/enthalpy.plt script
8d051ce Added -range option to cell.plt and econste_cmp.plt
f030ffc Fix file headers
--------------------------------------------------------------------------------
rel1_63_8
--------------------------------------------------------------------------------
4ee30f9 Fix discontinuity in torsion constraint function
d2a3c52 Removed $Id strings
--------------------------------------------------------------------------------
Francois Gygi committed
245 246
rel1_63_7
--------------------------------------------------------------------------------
Francois Gygi committed
247
r1920: Modified UserInterface.C to use C library, fsync()
Francois Gygi committed
248 249 250 251 252 253 254 255 256
to ensure synchronization.
r1912-1919: improvements in util/qbdriver files, twin.C
Added test programs in util/qbdriver
Added runtwin.sh script for interactive test of client-server mode
--------------------------------------------------------------------------------
rel1_63_6
--------------------------------------------------------------------------------
r1907: add sync in UserInterface.C
--------------------------------------------------------------------------------
Francois Gygi committed
257 258 259 260 261 262 263 264
rel1_63_5
--------------------------------------------------------------------------------
r1886: fix coeff of divergence correction for nspin=2 and multiple k-points
r1881,r1883: qbox_eig.py, xyz2qbox.py scripts
r1871: fix etotal_int calc with multiple k-points
r1869: removed blanks from timing name attribute
Enabled iter_cmd with modifications of positions or velocities
--------------------------------------------------------------------------------
Francois Gygi committed
265 266 267 268 269 270
rel1_63_4
--------------------------------------------------------------------------------
r1844,r1845,r1846: add qbox_maxforce.py script
r1835-r1843: B3LYP and BLYP spin-polarized implementation, SVWN updated
r1829-r1832: qbox_angle.py and qbox_torsion.py scripts
--------------------------------------------------------------------------------
Francois Gygi committed
271 272
rel1_63_3
--------------------------------------------------------------------------------
Francois Gygi committed
273
r1825: TorsionCmd.h: revert to previous sign convention for compatibility
Francois Gygi committed
274 275 276 277 278 279
with VMD and consistency with torsion constraint.
--------------------------------------------------------------------------------
r1823: Revert to use of wf (rather than dwf) in preconditioner update
--------------------------------------------------------------------------------
r1818: Added set_velocity command.
--------------------------------------------------------------------------------
Francois Gygi committed
280 281 282
r1809-r1813: modified Context to accept MPI communicators as ctor argument.
Removed pimpl idiom from implementation of Context.
--------------------------------------------------------------------------------
Francois Gygi committed
283 284 285
r1801 moved sync_positions and sync_velocities to AtomSet set_positions and
set_velocities functions
--------------------------------------------------------------------------------
Francois Gygi committed
286 287 288 289 290
rel1_63_2
--------------------------------------------------------------------------------
r1797 BOSampleStepper.C removed timer used in onpe0 section
r1796 JDWavefunctionStepper.C fixed timer name
--------------------------------------------------------------------------------
Francois Gygi committed
291 292 293 294 295 296 297 298
r1792 rel1_63_1
--------------------------------------------------------------------------------
r1792 CGIonicStepper.C: reset velocities at end of CG step
r1791 BOSampleStepper.C: modified etotal_int convergence criterion
r1788 added timers in JDWavefunctionStepper.C
--------------------------------------------------------------------------------
r1786 rel1_63_0
--------------------------------------------------------------------------------
Francois Gygi committed
299 300 301 302
r1779-r1782 Implementation of Harris-Foulkes functional for calc of etotal_int
--------------------------------------------------------------------------------
rel1_62_7 r1776
--------------------------------------------------------------------------------
Francois Gygi committed
303
r1774: Fixed output of partial_charge command.
Francois Gygi committed
304 305
Increased max iter in ConstraintSet (note: two changes in r1774. accidental)
--------------------------------------------------------------------------------
Francois Gygi committed
306
r1771 rel1_62_6
307 308 309 310
Implementation of partial_charge command
Implementation of iter_cmd and iter_cmd_period variables
JDWavefunctionStepper.C Fix: use dwf for preconditioning.
--------------------------------------------------------------------------------
Francois Gygi committed
311
rel1_62_5 r1760
Francois Gygi committed
312 313 314 315 316 317 318 319 320
r1758 GGIonicStepper.C
Reduced alpha_max to 5.0
Rescale alpha when displacement exceeds limit
r1757 testFourierTransform.C
Added MPI barriers for more accurate timing.
Added print of grid sizes.
r1756 AtomSet.C fix reset_vcm in randomize_velocities
r1754 r1755 ElectricEnthalpy.C remove omp pragmas to avoid race condition
--------------------------------------------------------------------------------
Francois Gygi committed
321 322 323 324 325 326 327 328 329 330 331 332 333 334 335
rel1_62_4 r1748
r1748 qb.C: Change URL in header to qboxcode.org
r1742-r1744 ExchangeOperator.C Removed implementation of quadratic correction
to exchange. Correction does not improve accuracy when far from asymptotic
convergence of kpoint sampling.
--------------------------------------------------------------------------------
rel1_62_3 r1735
r1733 Species.C: Fix bug in implementation of nonlocal potential f-projector
r1735 NonLocalPotential.C: Fix bug in nonlocal potential stress contribution
from f projectors
--------------------------------------------------------------------------------
rel1_62_2 r1725
ElectricEnthalpy.[Ch] Changed "polarization" to "dipole" to reflect the
correct physical definition of output. Adjusted output format.
--------------------------------------------------------------------------------
Francois Gygi committed
336 337 338 339 340 341 342
rel1_62_1 r1722
Species.C use 12 digit precision
Fix uninitialized variables
--------------------------------------------------------------------------------
rel1_62_0
Implementation of electric field
--------------------------------------------------------------------------------
343 344 345 346 347 348 349 350 351 352 353 354 355 356
rel1_61_0
Implementation of ONCV potentials
/xml/species.xsd
/xml/sample.xsd
/src/NonLocalPotential.C
/src/SpeciesHandler.C
/src/NonLocalPotential.h
/src/EnergyFunctional.C
/src/SpeciesHandler.h
/src/EnergyFunctional.h
/src/Species.C
/src/Species.h
/src/ChargeDensity.C
/src/ChargeDensity.h
Francois Gygi committed
357
--------------------------------------------------------------------------------
358 359 360 361 362 363 364 365 366
rel1_60_9
r1656 ExchangeOperator.C: changed default bisection tolerance to 1.0
--------------------------------------------------------------------------------
rel1_60_8
r1655 ExchangeOperator print exch time only if DEBUG defined
r1652 Bisection.C print nsweep only if TIMING defined
r1651 ExchangeOperator.C added omp directives in jade.C
r1649 ExchangeOperator.C, Bisection.C added BISECTION_NSWEEP debug option
--------------------------------------------------------------------------------
Francois Gygi committed
367 368 369 370 371 372 373
rel1_60_7
--------------------------------------------------------------------------------
r1645: Bisection.C: modified pair_fraction to reduce calls to overlap.
r1644: optimization of Bisection::overlap speeds up exchange on large problems
r1644: added detailed timers to ExchangeOperator.C
--------------------------------------------------------------------------------
rel1_60_6
Francois Gygi committed
374
--------------------------------------------------------------------------------
Francois Gygi committed
375 376 377 378 379 380
r1605: added -vlocal option to the plot command. Plot Vloc+VHart+Vxc.
r1570: cleanup testXMLGFPreprocessor.C
r1569: cleanup tab ExchangeOperator.C
--------------------------------------------------------------------------------
rel1_60_5
--------------------------------------------------------------------------------
Francois Gygi committed
381 382 383 384 385 386 387 388 389
r1564: modify FFTW3 use of threads
Added -pbc option to distance, angle and torsion functions
PlotCmd.C: fixed -wfs option.
r1552 CPSampleStepper.C: added stepper kinetic energy to econst
Added reset_vcm at end of randomize_v
--------------------------------------------------------------------------------
rel1_60_4 r1550
--------------------------------------------------------------------------------
added FFTW3 option and restored ESSL_FFT functionality
Francois Gygi committed
390 391
r1506: Fixed typo in AlphaPBE0.h
--------------------------------------------------------------------------------
Francois Gygi committed
392 393 394 395
rel1_60_2
CGCellStepper.C: revert to 1.60.0 use of matrix transformations
Use reduced alpha_init (0.002)
--------------------------------------------------------------------------------
Francois Gygi committed
396 397
rel1_60_1
CGCellStepper.C: fixed incorrect use of amat_inv() in calc of gradient.
Francois Gygi committed
398
Note: use was correct. later restored.
Francois Gygi committed
399 400
CGCellStepper.C: reduced initial alpha value to 0.01, max to 0.5.
--------------------------------------------------------------------------------
401 402 403 404
qbox_msd.py: print msd value only, not step index
--------------------------------------------------------------------------------
Context.C: fixed deprecated cast to const char*
--------------------------------------------------------------------------------
405 406 407 408 409 410
rel1_60_0
Added alpha_PBE0 variable
--------------------------------------------------------------------------------
rel1_59_0
Added scf_tol variable
--------------------------------------------------------------------------------
Francois Gygi committed
411 412 413
rel1_58_0
Added utility scripts
--------------------------------------------------------------------------------
Francois Gygi committed
414 415 416
rel1_57_15
--------------------------------------------------------------------------------
r1409: ExchangeOperator.C: fixed coeffs of divergent terms for empty states
Francois Gygi committed
417
r1407: XMLGFPreprocessor.C: add unistd header for close function
Francois Gygi committed
418
--------------------------------------------------------------------------------
Francois Gygi committed
419 420 421 422 423 424
rel1_57_14
--------------------------------------------------------------------------------
r1396: util/kpgen/kpgen.C: fixed odd numbered case and added test.
r1395: Species.C: fix for case of missing radial functions
r1394: BOSampleStepper.C: include weight info in eigenvalues
r1390,r1391: Fix exchange calculation with multiple k-points
Francois Gygi committed
425 426
r1388: Species.[Ch]: ensure that print outputs the same data as was read
from the pseudopotential file.
Francois Gygi committed
427 428
r1363: XMLGFPreprocessor.C: implemented sample loading using http connection.
--------------------------------------------------------------------------------
Francois Gygi committed
429 430
rel1_57_13
--------------------------------------------------------------------------------
Francois Gygi committed
431
r1358: ExchangeOperator.C: Fixed exchange contributions to stress.
Francois Gygi committed
432
KPGridConnectivity.[Ch]: Added cell_moved() member to update cell-dependent
Francois Gygi committed
433 434
quantities during cell optimization.
--------------------------------------------------------------------------------
435 436
rel1_57_12
--------------------------------------------------------------------------------
Francois Gygi committed
437
r1354: Added cell_moved() member in ExchangeOperator to resize vbasis
438
--------------------------------------------------------------------------------
Francois Gygi committed
439 440 441 442 443
rel1_57_11
--------------------------------------------------------------------------------
r1349: modified Wavefunction::operator= to allow assignment between wfs with
different cell sizes. Needed for variable cell calculations.
--------------------------------------------------------------------------------
Francois Gygi committed
444 445 446 447 448 449
rel1_57_10
--------------------------------------------------------------------------------
r1344: removed dependency on Context in multiple classes
       SlaterDet uses MPI column communicator instead of Context
--------------------------------------------------------------------------------
rel1_57_9
450
--------------------------------------------------------------------------------
Francois Gygi committed
451 452
r1326: WavefunctionHandler.[Ch]: fixed current_gfdata_pos variable
--------------------------------------------------------------------------------
Francois Gygi committed
453 454 455 456 457 458 459 460 461 462 463 464 465 466
rel1_57_8
--------------------------------------------------------------------------------
r1317: undo r1312: keep positive constraint on delta_spin
r1316: Fixed loading of gfdata into slater_det with non-zero delta_spin
r1315: Restored correct preconditioning when using stress confinement
r1314: EnergyFunctional.[Ch]: removed dependency on ChargeDensity.h from header
r1313: XCOperator.[Ch]: removed dependency on ChargeDensity.h from header
r1312: allow negative delta_spin
r1311: MDWavefunctionStepper.C: removed unused vars
r1310: NonLocalPotential.C: removed unused vars
r1309: removed deprecated conversions from string const to char* in headers
--------------------------------------------------------------------------------
rel1_57_7
--------------------------------------------------------------------------------
Francois Gygi committed
467 468 469 470
r1304: BOSampleStepper.C: removed printing of delta_ehart and delta_eigsum.
Renamed variable nscf_converged.
r1303: Wavefunction.C: fixed add_kpoint to avoid modifying occ of existing kpts
This caused incorrect occupation numbers when loading wf with multiple kpts.
Francois Gygi committed
471 472
r1301: BOSampleStepper.C: print RC1 override msg on pe0 only
--------------------------------------------------------------------------------
Francois Gygi committed
473 474 475 476
rel1_57_6
--------------------------------------------------------------------------------
r1296: SampleHandler.C: set wfv to use same ctxt as s_.wf
--------------------------------------------------------------------------------
Francois Gygi committed
477 478 479 480 481 482 483 484
rel1_57_5
--------------------------------------------------------------------------------
r1292: fixed mlwf calc for spin-polarized case
r1291: ComputeMLWFCmd.C: added spin implementation
r1290: Anderson charge mixing on both spins simultaneously
r1288: BOSampleStepper: moved enthalpy definition out of onpe0 test
--------------------------------------------------------------------------------
r1286: rel1_57_4
Francois Gygi committed
485 486 487 488
r1283: PSDAWavefunctionStepper.C: fixed extrapolate flag
r1282: rewrite of kpgen.C
--------------------------------------------------------------------------------
rel1_57_3
Francois Gygi committed
489
r1275: fix sample read/write for nspin=2. Save cmd writes full species info.
490 491 492
--------------------------------------------------------------------------------
rel1_57_2
r1272: BOSampleStepper.C: fixed nite loop test for nite==0 case
Francois Gygi committed
493
--------------------------------------------------------------------------------
Francois Gygi committed
494
rel1_57_1
Francois Gygi committed
495 496 497
r1265: BOSampleStepper.C: changed behavior for nite=1, use eigenvalue_sum for
convergence test of nite loop, changed threshold value to 0.01*delta_ehart
r1264: RunCmd.C: changed default nite value to 0
Francois Gygi committed
498
r1268: ExchangeOperator.C: Added EXCHANGE_NOQUAD debug string. Cleanup.
Francois Gygi committed
499 500 501 502 503 504
--------------------------------------------------------------------------------
rel1_57_0 r1262
Cleanup of XCOperator and ExchangeOperator to include stress at gamma point
--------------------------------------------------------------------------------
rel1_56_2 r1259
r1257: EnergyFunctional: fix nspin=2 forces, stress.
Francois Gygi committed
505 506
BOSampleStepper: fixed bug in mixing with charge_mix_ndim=0
r1212: Adaptive preconditioner.
Francois Gygi committed
507 508 509
r1211: tolerance for non-self consistent iterations
--------------------------------------------------------------------------------
rel1_55_6
Francois Gygi committed
510 511 512 513 514 515 516
Redesign of Bisection and ExchangeOperator separates bisection transformation
and load-balancing permutation.
Fixed KPGridConnectivity
Added unit_cell info in BOSampleStepper output.
SlatedDet: added occ_ and eig_ to operator=
--------------------------------------------------------------------------------
rel1_55_5
Francois Gygi committed
517
Added polar decomposition in Matrix class.
Francois Gygi committed
518 519 520 521 522 523 524
Use polar in ortho_align, align, lowdin ortho.
--------------------------------------------------------------------------------
rel1_55_4
Added OpenMP pragmas
--------------------------------------------------------------------------------
rel1_55_0
merge of hfb functionality
Francois Gygi committed
525
--------------------------------------------------------------------------------
Francois Gygi committed
526 527 528
rel1_54_6
Added USE_DFFTW macro for platforms providing the dfftw.h header only
--------------------------------------------------------------------------------
Francois Gygi committed
529 530 531 532 533 534 535 536
rel1_54_5
Added enthalpy in BOSampleStepper for use in cell optimization.
BasisMapping.C: fix bug of serial version: MLWF calc of dipole wrong.
Added randomize_r command
removed printing on stdout in rseed command
Changed number of iterations in SlaterDet::{ortho_align,align,riccati}
--------------------------------------------------------------------------------
rel1_54_4
537 538 539 540 541 542 543
Added randomize_v command to set velocities to a Maxwell-Boltzmann distribution
Modified CGIonicStepper: project forces on constraints
Modified BMDIonicStepper: rescale velocity rather than stop
Other steppers: added synchronization calls to ensure consistency of positions
and unit cell values across tasks.
--------------------------------------------------------------------------------
rel1_54_3
544 545 546 547 548 549 550 551 552 553 554 555 556 557
Rewrite of CellStepper, SDCellStepper, CGCellStepper classes.
CGOptimizer and LineMinimizer restructured.
--------------------------------------------------------------------------------
rel1_54_2
fix CG algorithm in CGCellStepper.C
--------------------------------------------------------------------------------
rel1_54_1
CGCellStepper.C: fixed CG algorithm.
--------------------------------------------------------------------------------
rel1_54_1
Basis.C: fixed bug 26: cannot build basis for kpoint close to FCC BZ boundary
fixed in r924
--------------------------------------------------------------------------------
rel1_54_0
558 559 560 561 562 563 564
LineMinimizer.h: limited the number of iterations in bracketing mode to 4
UnitCell.C: accept unit cells with orientation leading to negative volume and
define the volume as fabs(a0*(a1^a2))
CellStepper.C, CGCellStepper.C: implemented CG algorithm to optimize both
cell parameters and atomic positions. Moved implementation of common CellStepper
member functions to CellStepper.C.
--------------------------------------------------------------------------------
Francois Gygi committed
565 566
rel1_53_0
RseedCmd.h: rseed command to initialize random number generator
Francois Gygi committed
567
BOSampleStepper.C: reintroduced wf extrapolation when nite>1 (leads to
Francois Gygi committed
568 569 570 571 572 573
much improved energy conservation in liquid Si)
Matrix.[Ch]: added IntegerMatrix, permutations
ComputeMLWF.C: test if not at Gamma point
StructureFactor.C, SampleWriter.C: added <cstring> for memset
XCPotential.h: moved header dependencies to XCPotential.C
--------------------------------------------------------------------------------
574
rel1_52_4
Francois Gygi committed
575
Introduced a USE_UUID macro in qb.C, uuid_str.[Ch] (some platforms have no
576
uuid library)
Francois Gygi committed
577
Fix bug 23: Wavefunction::set_nempty: added update occ when nempty=0.
578
--------------------------------------------------------------------------------
Francois Gygi committed
579 580 581 582 583 584 585
rel1_52_3
Fixed StructuredDocumentHandler.[Ch] and WavefunctionHandler.C for compilation
with Xerces-C 3.x.x
SampleReader.C: removed debug print
--------------------------------------------------------------------------------
rel1_52_2
RunCmd.h: updated help msg
Francois Gygi committed
586 587 588 589 590 591
rel1_52_2
RunCmd.h: updated help msg
LoadCmd.C: fix bug in reading a sample with atomset only. (unit_cell not set).
Add <cmd> </cmd> tags when echoing commands executed by qbox.
Add <eigenset> tag around <eigenvalues> for each kpoint and spin.
BOSampleStepper.C: Disabled wf extrapolation if nite>1 (unstable)
Francois Gygi committed
592
--------------------------------------------------------------------------------
Francois Gygi committed
593 594 595 596 597 598 599 600 601
rel1_52_1
SlaterDet::gram use beta=0 in herk, for consistency (complex case).
(does not affect results since Matrix ctor clears data).
SlaterDet: removed constraint on Im(c(G=0)) for complex case.
Note: this fixes a bug introduced in 1.51.1 which incorrectly set Im=0 for k!=0.
This fixes bug 2010-04-04T01:14:19Z
Matrix.C: parameters alpha and beta of pzherk and zherk must be double*
Moved files uuid.[Ch] to uuid_str.[Ch] to avoid name conflicts with uuid/uuid.h
--------------------------------------------------------------------------------
Francois Gygi committed
602 603 604 605 606 607 608
rel1_52_0
Added the rescale_v command (RescaleVCmd.h) to rescale all velocities.
Implemented the Bussi-Donadio-Parrinello thermostat.
Added files sampling.C sampling.h
Modified IonicStepper.h to add ekin_stepper member
BOSampleStepper.C: <econst> includes stepper kinetic energy when using BDP.
FourierTransform.C: included some OpenMP pragmas
Francois Gygi committed
609
Known bug: 2010-04-04T01:14:19Z
Francois Gygi committed
610 611
--------------------------------------------------------------------------------
rel1_51_1
Francois Gygi committed
612
Attempt to fix bug 2010-04-04T01:14:19Z
Francois Gygi committed
613
Small imaginary part in c(G=0), initialization in SlaterDet::resize.
Francois Gygi committed
614 615 616
Incorrect fix of 2010-04-04T01:14:19Z causes incorrect results in k!=0 wfs
Later corrected in 1.52.1
Known bug: 2010-04-04T01:14:19Z
Francois Gygi committed
617
--------------------------------------------------------------------------------
Francois Gygi committed
618 619 620 621 622 623 624 625
rel1_51_0
Added external forces: atomic, pair or global.
Added capability to load samples multiple times
Added capability to redefine species
Added uuid identifier in output
SlaterDet.C, Wavefunction.C modified clear, reset, init members
Species.C: modified Species::info to print XML elements
--------------------------------------------------------------------------------
Francois Gygi committed
626 627 628
rel1_50_4
ConstraintSet.C: Fixed uninitialized ndofs_ in ctor.
--------------------------------------------------------------------------------
Francois Gygi committed
629 630 631
rel1_50_3
SampleWriter.C: Fixed size calculation for serial version
--------------------------------------------------------------------------------
632 633 634 635 636
rel1_50_2
Updated several files to enable serial version (without MPI or ScaLAPACK)
Added include<cstring> in various files to get correct memcpy declaration.
(caused trouble with gcc 4.4.1 on Ubuntu).
--------------------------------------------------------------------------------
Francois Gygi committed
637 638 639
rel1_50_1
EnergyFunctional.C:732: use sigma only if compute_stress==true
--------------------------------------------------------------------------------
Francois Gygi committed
640 641 642 643 644 645
rel1_50_0
BOSampleStepper.C, RunCmd.C: added [-atomic_density] option to the run cmd.
AndersonMixer.C: iterative LS problem with adaptive Tikhonov parameter.
Matrix.[Ch]: added getsub for arbitrary block destination
JDWavefunctionStepper.[Ch]: Jacobi-Davidson WavefunctionStepper
--------------------------------------------------------------------------------
Francois Gygi committed
646 647 648 649
rel1_49_1
EnergyFunctional.C: increased criterion for determination of neighboring
cells in esr calculation.
--------------------------------------------------------------------------------
Francois Gygi committed
650 651 652 653 654 655 656
rel1_49_0
BOSampleStepper.C: introduced separate row-weighting of Anderson LS problem
with fixed cutoff rc1=3.0 a.u.
AndersonMixer.C: modified regularization
PlotCmd.C: white space
SlaterDet.C: white space
--------------------------------------------------------------------------------
Francois Gygi committed
657
rel1_48_1
Francois Gygi committed
658 659 660 661 662
SlaterDet.C: fix bug in SlaterDet::write for -text mode
BOSampleStepper.C: modified Kerker weights. Added mixer.reset() call.
AtomCmd.h: fixed text of help message.
AndersonMixer.C: added Tikhonov regularization. Bcast results from task 0.
PlotCmd.C: cleanup include.
Francois Gygi committed
663
--------------------------------------------------------------------------------
Francois Gygi committed
664 665 666 667 668
rel1_48_0
PlotCmd.[Ch]: added plot command. Generates cube format or xyz
{LDA,PBE,BLYP}Functional: modified API for compat. with future implem. of PBE0
SlaterDet.C: added timing ifdefs in gram ortho.
--------------------------------------------------------------------------------
Francois Gygi committed
669 670 671
rel1_47_1
Added a virtual member function in Constraint that returns the number of DOFs to
be subtracted due to that constraint. That number is 1 for all constraints,
Francois Gygi committed
672 673
except for PositionConstraint, for which it is 3. This function is used
by ConstraintSet to keep track of the total number of DOFs to be removed
Francois Gygi committed
674 675
when computing the system temperature. The ConstraintSet class maintains the
total number of DOFs and modifies it only when constraints are added or deleted.
Francois Gygi committed
676 677
MDIonicStepper::ekin() modified to compute kinetic energy correctly in the
presence of constraints.
Francois Gygi committed
678 679
Emass.h: allow reset to zero
BOSampleStepper.C: fix print format of etotal_int
Francois Gygi committed
680
--------------------------------------------------------------------------------
Francois Gygi committed
681 682 683 684 685 686 687 688
rel1_47_0
Add BMDIonicStepper class.
ConstraintCmd.h ConstraintSet.C: Added PositionConstraint class to lock atoms.
Species.C: modified calculation of ndft to ensure rmax>40 a.u.
This corrects the problem that ndft would depend on the number of mesh points
given in the species file.
WavefunctionHandler.C: fix XMLString::release() call for Xerces 3.0.1
--------------------------------------------------------------------------------
Francois Gygi committed
689
rel1_46_0
Francois Gygi committed
690
AndersonMixer.C, BOSampleStepper.C: Implemented multidimensional Anderson
Francois Gygi committed
691
acceleration of charge mixing.
Francois Gygi committed
692
Added charge_mix_ndim to specify dimension of Anderson acceleration of
Francois Gygi committed
693
charge mixing. charge_mix_ndim=0 results in simple mixing.
Francois Gygi committed
694 695 696 697 698
--------------------------------------------------------------------------------
rel1_45_3
UserInterface.C: modified server mode to allow multiple commands in input file.
Use file "<inputfile>.lock" to synchronize with client.
Output file in server mode is now a valid XML fragment, with fpmd namespace.
Francois Gygi committed
699 700
--------------------------------------------------------------------------------
rel1_45_2
Francois Gygi committed
701
qb.C, UserInterface.C: implemented -server option, cleaned up destructors.
Francois Gygi committed
702 703 704 705 706 707 708
AtomSet.C: implemented destructor (fix leaks)
qb.C: added delete s at end (fix leaks)
ConstraintSet: implemented destructor (fix leaks)
Constraint.h and derived classes: implemented virtual destructors.
StrainCmd.h: fix error in usage message.
XMLGFPreprocessor.C: new macro PARALLEL_FS. Limit the number of simultaneous
fread calls if PARALLEL_FS is not defined.
Francois Gygi committed
709
--------------------------------------------------------------------------------
Francois Gygi committed
710
rel1_45_1: fixed SDCellStepper.C: make local copy of sigma to avoid
Francois Gygi committed
711 712
PGI compiler error on franklin: "must be lvalue or...".
--------------------------------------------------------------------------------
Francois Gygi committed
713 714 715 716 717 718 719 720
rel1_45_0
Added strain command for deformation of the cell.
UnitCell.C: removed member compute_deda (moved calc to SDCellStepper). Added
symmetric matrix - vector product function.
SDCellStepper.C fixed calculation of cell correction from stress tensor.
BOSampleStepper.C: reintroduced wf extrapolation when nempty>0 for MD.
RunCmd.h: fixed help msg.
--------------------------------------------------------------------------------
Francois Gygi committed
721 722 723
rel1_44_5
BOSampleStepper.C: Fix uninitialized mlwft ptr. Fixed delete statements.
--------------------------------------------------------------------------------
Francois Gygi committed
724 725
rel1_44_4
BOSampleStepper.C: reintroduced feature that keeps memory of wfv between run
Francois Gygi committed
726
commands. Allows for smooth continuation of multiple runs. Added MLWF and
Francois Gygi committed
727 728 729 730
MLWFC calculation in the niter loop.
Modified wf extrapolation to include alignment in all cases. Use lowdin
orthogonalization after extrapolation.
WfDiag.h: Added MLWF and MLWFC option for the wf_diag variable: computes MLWFs
Francois Gygi committed
731
or MLWF centers at each ionic step.
Francois Gygi committed
732 733
SaveCmd.[Ch], SampleWriter.[Ch]: added -no_wfv option: save without wf velocity.
--------------------------------------------------------------------------------
Francois Gygi committed
734 735
rel1_44_3
UserInterface.h, Base64Transcoder.C: added include <cstring>
Francois Gygi committed
736
removed trailing blanks.
Francois Gygi committed
737
--------------------------------------------------------------------------------
Francois Gygi committed
738 739 740
rel1_44_2
Added copyright headers.
--------------------------------------------------------------------------------
741 742 743 744
rel1_44_1
EnergyFunctional.C: fixed calculation of Esr and Esr stress in cases where the
unit cell is small compared with rcps (initialized in AtomSet.C).
--------------------------------------------------------------------------------
Francois Gygi committed
745 746 747 748
rel1_44_0
WavefunctionHandler.[Ch], SampleReader.[Ch]: fixed compatibility with older
versions when reading samples with smaller grid sizes.
--------------------------------------------------------------------------------
Francois Gygi committed
749 750 751 752 753 754 755 756 757 758
rel1_43_5
BOSampleStepper.C: print energy at end of wf optimization when using niter=0.
Basis.C: modified grid size to cover all grids for any k-point
--------------------------------------------------------------------------------
rel1_43_4
Basis.C: np() yields same grid size for all kpoints
XMLGFPreprocessor.C, SampleReader.C: fixed reading samples with both real and
complex wavefunctions
Wavefunction.C: removed unneeded resize in add_kpoint
--------------------------------------------------------------------------------
Francois Gygi committed
759 760 761 762 763 764
rel1_43_3
TorsionCmd.h: fixed sign of torsion angle.
SlaterDet.[Ch]: fixed const-ness of some members.
Wavefunction.C: modif allocation in add_kpoint to avoid O(nkp**2) time
x8664_gcc.mk: link statically to gfortran libs for better portability of exec.
--------------------------------------------------------------------------------
765 766 767
rel1_43_2
SlaterDet.C: rewrite of the parallel I/O section, using less memory.
--------------------------------------------------------------------------------
Francois Gygi committed
768 769 770 771 772
rel1_43_1
BOSampleStepper.C, CPSampleStepper.C: fixed initialization of ekin_ion,temp_ion.
SampleWriter.C, Wavefunction.C, SlaterDet.C: implementation of parallel I/O
using SharedFilePtr class.
--------------------------------------------------------------------------------
Francois Gygi committed
773 774 775 776 777 778
rel1_43_0
SampleWriter.C, Wavefunction.C, SlaterDet.C: implemented MPI-IO parallel I/O
using individual file pointers, no shared pointers.
SlaterDet.C: fixed bug: base64 transcoding done twice on big endian platforms
when writing in parallel.
--------------------------------------------------------------------------------
Francois Gygi committed
779 780 781 782
rel1_42_0
SlaterDet.C, Wavefunction.C: implemented dynamic resize of ecut or nempty
without loss of data. Changes of cell also preserve wavefunctions.
--------------------------------------------------------------------------------
Francois Gygi committed
783
rel1_41_0
Francois Gygi committed
784
BOSampleStepper.C, CPSampleStepper.C: remove redundant printing of unitcell
Francois Gygi committed
785 786 787 788 789
LineMinimizer.h: new LineMinimizer class used in optimization algorithms
SDIonicStepper.C: back to simple dt^2/m steepest descent (back to rel1_39_0)
SDAIonicStepper.C: steepest descent with line minimization.
CGIonicStepper.C: conjugate gradient with line minimization.
--------------------------------------------------------------------------------
Francois Gygi committed
790 791 792 793 794 795
rel1_40_0
AtomsDyn.h,CGIonicStepper.[Ch]: added CG ionic stepper.
BOSampleStepper.C: enabled niter=0 in run command.
MDIonicStepper.C: removed trailing blanks.
SDAIonicStepper: reduced theta_max, test for energy increase.
SDIonicStepper: implemented line search, removed dt, mass dependency.
Francois Gygi committed
796
Known problems: SD, SDA, CG algorithms are not working properly
Francois Gygi committed
797
--------------------------------------------------------------------------------
Francois Gygi committed
798 799 800 801 802
rel1_39_0
SampleReader.C, SpeciesReader.C: set validation scheme to Val_Auto to allow
successful parsing in the absence of a local copy of sample.xsd or species.xsd.
Added fold_in_ws command: fold all atoms in the WS cell.
MDIonicStepper.C: removed collision debug info printout.
Francois Gygi committed
803
MDIonicStepper.C: make calc of collisions on task 0, then sync in LOWE and
Francois Gygi committed
804
ANDERSEN thermostats. Ensures that all velocities are consistent on all nodes.
Francois Gygi committed
805
PSDAWavefunctionStepper.C: removed psda_enhancement_factor which caused
Francois Gygi committed
806
oscillations in some PSDA optimizations.
Francois Gygi committed
807
BOSampleStepper.C: use wf velocity only inside BOSampleStepper. wfv is not
Francois Gygi committed
808 809
saved at end of steps. wfv is not used at first step for wf extrapolation.
--------------------------------------------------------------------------------
Francois Gygi committed
810 811
rel1_38_0
MDIonicStepper.C: Modified minimum image convention in LOWE thermostat to
Francois Gygi committed
812
conserve angular momentum in clusters. Fixed LOWE and ANDERSEN thermostats to
Francois Gygi committed
813 814 815 816 817 818 819 820
have correct collision rate.
UserInterface.C: fixed endl in output to avoid double endl.
Basis.C: remove pimpl implementation.
ChargeDensity.C: use factorizable sizes for density grid.
EnergyFunctional.C: corrected upper bound of ia1 in esr calculation.
EnergyFunctional.C: implemented multiple cell in esr calculation.
AtomSet.C: changed rcps value to 1.5. Changed interface to access cell.
--------------------------------------------------------------------------------
821 822 823
rel1_37_0
Cleaned up XML formatting of output. Removed comments. Use <timing> tag.
--------------------------------------------------------------------------------
824 825 826 827 828 829
rel1_36_0
Added unit cell to the sample schema. Modified AtomSet to contain a unit cell.
Included unit cell in saved samples.
--------------------------------------------------------------------------------
rel1_35_0
Implemented parallel I/O for writing sample files. Uses MPI-IO.
Francois Gygi committed
830
This does not work on BG/L with PVFS since PVFS does not support the
831 832 833
collective calls MPI_File_write_ordered and MPI_File_write_shared.
Option -serial added to the save command to use the old serial mode.
--------------------------------------------------------------------------------
Francois Gygi committed
834 835 836 837 838 839
rel1_34_4
Base64Transcoder.C: modif print functions to use larger string buffer for I/O.
SampleReader.C, SampleWriter.C: added timer.
Corrected lda arg in various BLACS calls. Did cause errors on
AIX/pessl. (Not a bug when using blacs lib compiled from netlib distribution).
--------------------------------------------------------------------------------
Francois Gygi committed
840 841 842 843 844
rel1_34_3
Wavefunction.C: fixed incorrect XML comments
Basis.C: fixed warning message printed by all tasks when k out of BZ.
AtomSet.C, SpeciesCmd.C: modif to always print species info when a species is
defined, i.e. also when loading from a sample file.
Francois Gygi committed
845
--------------------------------------------------------------------------------
Francois Gygi committed
846 847 848 849 850 851 852 853 854 855 856 857 858
rel1_34_2
WavefunctionHandler.C: fixed for compilation on BG/L
SampleReader.C: fixed for case where loading multiple k-points with wfv
NonLocalPotential.C: fixed calc of forces with multiple k-points
EnergyFunctional.C: fixed forces and stress with multiple k-points
MDWavefunctionStepper.C: fixed CP dynamics with k-points
MLWFTransform.C: changed syntax of blas call (zcopy instead of zcopy_)
NonLocalPotential.C: fix forces with k-point, fix bug introduced in 1.34.0.
SampleReader.C: multiple k-point and wf velocity functionality.
SlaterDet.C: fixed print() and write().
Wavefunction.C: changed context allocation. Fix update_occ, entropy.
PSDWavefunctionStepper.C: fixed complex case (k!=0)
--------------------------------------------------------------------------------
Francois Gygi committed
859 860 861 862 863 864 865
rel1_34_1
Several bug fixes for kpoint functionality.
CPSampleStepper.C: create wfv if absent.
ChargeDensity.C: increased grid size to avoid aliasing with k!=0
--------------------------------------------------------------------------------
rel1_34_0 (not for distribution)
Added kpoints. This release still has bugs. Not for distribution.
Francois Gygi committed
866 867 868 869
--------------------------------------------------------------------------------
rel1_33_5
removed trailing blanks in *.[Ch] *.mk Makefile
--------------------------------------------------------------------------------
Francois Gygi committed
870 871 872
rel1_33_4
ComputeMLWFCmd.C: fixed memory leak.
blas.h: added zgemm, dnrm2.
Francois Gygi committed
873
BOSampleStepper.C: disabled Anderson acceleration of charge mixing by default,
Francois Gygi committed
874 875 876
can be enabled with debug="AND_CHMIX"
AtomSet.C: removed unused variables.
--------------------------------------------------------------------------------
Francois Gygi committed
877 878 879 880 881 882 883
rel1_33_3
BOSampleStepper.C: extrapolation: restored alignment before extrapolation
ChargeMixCoeff.h: changed default value back to 0.5
Matrix.h: added iglobal and jglobal functions
SDAIonicStepper.h: modified default theta_max and theta_nc values
Species.[Ch]: added phi function
--------------------------------------------------------------------------------
Francois Gygi committed
884 885 886
rel1_33_2
Implemented calculation of total electronic and ionic dipole in MLWF.
--------------------------------------------------------------------------------
Francois Gygi committed
887 888 889 890 891 892 893 894 895
rel1_33_1
MLWFTransform.C: fixed bug in MLWF calc when cell sizes in x and y direction
differ.
--------------------------------------------------------------------------------
rel1_33_0 candidate
Implemented compute_mlwf command
Changed spline.C and sinft.C with new implementations.
Modified Species.C to use new spline and sinft interface.
--------------------------------------------------------------------------------
Francois Gygi committed
896
rel1_32_0
897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919
Dt.h, MDIonicStepper.C: modified to allow for negative time step.
Removed Sample.h dependency in WavefunctionStepper and derived classes.
Files affected are
BOSampleStepper.C
SDWavefunctionStepper.C
PSDWavefunctionStepper.C
PSDAWavefunctionStepper.C
EnergyFunctional.C
CPSampleStepper.C
MDWavefunctionStepper.C
WavefunctionStepper.h
SDWavefunctionStepper.h
PSDWavefunctionStepper.h
PSDAWavefunctionStepper.h
MDWavefunctionStepper.h
Removed facs, sigmas from Control.h. Moved vars to ConfinementPotential.
Use hardwired values 2.0 and 0.5 in
BOSampleStepper.C, CPSampleStepper.C
Removed namespace declarations in header files.
--------------------------------------------------------------------------------
rel1_31_0
Added class SampleWriter to separate write function from SaveCmd class.
--------------------------------------------------------------------------------
Francois Gygi committed
920 921 922 923 924 925 926 927
rel1_30_1
Wavefunction.C: typos
XMLGFPreprocessor.C: fixed bug in boundary adjustment.
jacobi.C: implement special case of fewer blocks than process columns.
BOSampleStepper.C: Moved Kerker preconditioning of charge density update
before Anderson acceleration.
ChargeMixCoeff.h: Changed default value to 1.0
ExtStress.h: initialize with zero in ctor
928
MDIonicStepper.C, Thermostat.h: corrected misspelled ANDERSEN
Francois Gygi committed
929
--------------------------------------------------------------------------------
Francois Gygi committed
930 931 932 933 934 935 936
rel1_30_0
D3vector.h: modified ctor to avoid accidental use with only one argument.
AtomSet.[Ch]: added functions to get/reset center of mass velocity.
StatusCmd.h: status cmd prints center of mass velocity.
MDIonicStepper.[Ch]: new thermostats ANDERSON LOWE. Old thermo renamed
to SCALING.
--------------------------------------------------------------------------------
Francois Gygi committed
937 938 939 940 941 942 943
rel1_29_1
BOSampleStepper.C: use wfmm only when needed. ASP extrapolation implemented.
Activate NTC or ASP extrapolation using debug variable.
XMLGFPreprocessor.C: added tests to avoid div by zero if no local grid
--------------------------------------------------------------------------------
rel1_29_0
BOSampleStepper.C: major rewrite:
Francois Gygi committed
944
added a wf_stepper->preprocess before the itscf loop to
Francois Gygi committed
945
get correct extrapolation when ions move.
Francois Gygi committed
946
Modified NTC extrapolation: activated if NTC_EXTRAPOLATION is found in the
Francois Gygi committed
947 948 949 950 951
debug variable.
ReadTSC.h: added missing #define statement
Debug.h: modified to take arbitrary strings
jacobi.C: new jacobi diagonalization. activate with ifdef in Wavefunction.C
--------------------------------------------------------------------------------
Francois Gygi committed
952 953 954 955 956 957 958 959 960 961 962 963
rel1_28_2
AtomSet.h: initialize nel_ in ctor.
XMLGFPreprocessor.C: use char buffer to read file.
Include timers in SlaterDet::ortho_align with #if TIMING
qb.C: remove <> brackets in processor names for XML compatibility (BGL)
AtomSetHandler.C, SpeciesHandler.C: removed msg about xlink, replaced with href
Species.C, SaveCmd.C: removed xlink namespace declaration in sample
SlaterDet.C: added #if TIMING and timers in ortho_align
MDIonicStepper.h: corrected comments
--------------------------------------------------------------------------------
rel1_28_1e{1234567} experimental versions containing
various modifications leading to rel1_28_2
964 965 966
--------------------------------------------------------------------------------
rel1_28_0
qb.C, SaveCmd.C, qbox_xmlns.C, xmlSpecies.C: change namespace to fpmd
Francois Gygi committed
967 968 969
--------------------------------------------------------------------------------
rel1_27_1
minor updates to make single-node version work.
Francois Gygi committed
970 971 972 973 974
--------------------------------------------------------------------------------
rel1_26_1e
redefine constraint command and subcommands, output format.
Introduce constraint names.
Remove ListConstraintsCmd.h
Francois Gygi committed
975
SDAIonicStepper.h: fixed bug:
Francois Gygi committed
976 977
Use of ndofs_ incorrect when ndofs<3*na (when constraints are defined).
Use CHOLESKY_REMAP macro value for size of square context in SlaterDet.C (bgl).
Francois Gygi committed
978
--------------------------------------------------------------------------------
Francois Gygi committed
979
rel1_26_0
Francois Gygi committed
980 981 982 983 984 985
Species.C: changed spline BC for local potential.
SlaterDet.C: randomize: return if basis size is zero.
Implementation of a new logic in all IonicStepper classes.
Implementation of the move command.
Implementation of constraints.
Implementation of distance, angle and torsion measurement commands.
Francois Gygi committed
986 987 988 989 990
--------------------------------------------------------------------------------
rel1_25_0
implementation of netcharge.
SpeciesCmd.C: print species information when adding species.
bug fix: MDWavefunctionStepper.C: correct dt2bye when ecuts>0.
Francois Gygi committed
991 992 993 994 995 996
--------------------------------------------------------------------------------
rel1_24_3
FourierTransform.C: removed unused code.
NonLocalPotential.C: cleanup calculation of G=0 in update_twnl.
BOSampleStepper.C: added APC counters
qb.C: added APC counters
Francois Gygi committed
997 998
--------------------------------------------------------------------------------
rel1_24_2
Francois Gygi committed
999
NonLocalPotential.C: Fixed bug in G=0 component of l=0 non-local projector in
Francois Gygi committed
1000
semilocal form (merged bugfix from rel1_23_4)
Francois Gygi committed
1001 1002
--------------------------------------------------------------------------------
rel1_24_1
Francois Gygi committed
1003
SlaterDet.C: fixed bug introduced in 1.24.0: Cholesky decomp was skipped if
Francois Gygi committed
1004
CHOLESKY_REMAP is not defined.
Francois Gygi committed
1005 1006 1007 1008 1009 1010 1011 1012 1013 1014
--------------------------------------------------------------------------------
rel1_24_0
UnitCell.C: added space in output of lattice vectors.
Wavefunction.C: added syevx option for diagonalization (not active)
BOSampleStepper.C: added diag timer.
SlaterDet.C: Added CHOLESKY_REMAP macro to enable remapping to a square context
in SlaterDet::gram.
SlaterDet.C: rewritten communication in SlaterDet::print and SlaterDet::write
to avoid a deadlock on BGL.
Matrix.C: Added syevx and syevd functions
Francois Gygi committed
1015
Note: modif of SlaterDet.C introduced a bug: see 1.24.1
1016
--------------------------------------------------------------------------------
Francois Gygi committed
1017 1018 1019
rel1_23_4 (on branch rel1_23_3_bf)
Fixed bug in NonLocalPotential.C: G=0 component of l=0 semilocal projector.
--------------------------------------------------------------------------------
1020 1021
rel1_23_3
NonLocalPotential.C: fixed bug: double fnl_size should be int fnl_size.
Francois Gygi committed
1022 1023 1024
--------------------------------------------------------------------------------
rel1_23_2
FourierTransform.C: Fixed bug #39: When unit cells are non-orthogonal, the
Francois Gygi committed
1025
transforms along x are not done correctly. Bounds including all non-zero
Francois Gygi committed
1026 1027 1028
transforms are redefined.
Basis.C: Simplified basis set construction.
NonLocalPotential.C: removed variables used in old implementation (pre 1.23.0)
Francois Gygi committed
1029
--------------------------------------------------------------------------------
Francois Gygi committed
1030
rel1_23_1
Francois Gygi committed
1031 1032
PSDAWavefunctionStepper.C: Reset default theta for negative curvature to 0.0.
AndersonMixer: added settable default theta value for negative curvature.
Francois Gygi committed
1033 1034 1035 1036 1037
--------------------------------------------------------------------------------
rel1_23_0
Added SDA algorithm for atoms dynamics: Anderson acceleration of SD.
New implementation of NonLocalPotential::energy, recomputing anl array at
every call.
Francois Gygi committed
1038
Modified IonicStepper to make ekin and temp virtual. Implemented only in
Francois Gygi committed
1039
MDIonicStepper.
1040
--------------------------------------------------------------------------------
Francois Gygi committed
1041 1042 1043
rel1_22_0
BOSampleStepper.C: new implementation of Kerker-Anderson charge mixing.
new variables charge_mix_rcut and charge_mix_coeff introduced: determine
Francois Gygi committed
1044
the Kerker cutoff and mixing coefficient parameters.
Francois Gygi committed
1045
--------------------------------------------------------------------------------
Francois Gygi committed
1046 1047 1048 1049 1050 1051 1052 1053 1054 1055
rel1_21_0
Context.[Ch]: redesigned Context to be column-major in all cases.
ChargeDensity.C: vcontext definition modified.
SlaterDet.C: new my_col_ctxt member: store single column Context.
Use modified Context allocation for my_col_ctxt to ensure that calls to
MPI_Comm_create are done with same arguments on all tasks.
Wavefunction.C: use simpler Context constructor (default is col-major).
UserInterface.C: removed printing of newline when echo==true.
Cleaned up timers in wf steppers.
--------------------------------------------------------------------------------
Francois Gygi committed
1056 1057
rel1_20_1
NonLocalPotential.C: fixed index calculation (bug 38).
Francois Gygi committed
1058
--------------------------------------------------------------------------------
Francois Gygi committed
1059 1060 1061 1062 1063 1064 1065
rel1_20_0
Wavefunction.C: Changed allocation of sd context to column-major order to get
contiguous nodes in  context columns.
Wavefunction.C: Changed maxiter to 500 and dmu to 2 eV in update_occ.
Nrowmax.h: add update_occ when changing nrowmax.
Matrix.C:219 and similar: removed assert. Some valptr may be 0 on some blocks.
SlaterDet.C: print sdcontext dimensions in SlaterDet::info.
Francois Gygi committed
1066
ChargeDensity.C: added (void*) cast to delete statements (for xlC 7.0).
Francois Gygi committed
1067
Known bug #38: incorrect forces with nrowmax=16 with 512 tasks. Occurs when
Francois Gygi committed
1068
nrowmax <= (ntasks/nrowmax)/2.
Francois Gygi committed
1069 1070
--------------------------------------------------------------------------------
rel1_19_0 2004-10-04
Francois Gygi committed
1071
BOSampleStepper.C: Implemented simple charge mixing for SCF iterations.
Francois Gygi committed
1072
Disabled potential mixing.
Francois Gygi committed
1073 1074 1075 1076 1077 1078
SpeciesReader.C: fixed parsing of species files.
SlaterDet.C in print and write: fixed incorrect allocation of char buffer b:
caused segfaults when saving sample.
Base64Transcoder.[Ch]: added cstdio style interface for print function.
FourierTransform.[Ch]: implemented timers, activated if TIMING is defined.
testFourierTransform.C: updated to print timing info in readable form.
Francois Gygi committed
1079 1080
Known bug: Wavefunction::update_occ iteration on mu does not converge in some
systems.
Francois Gygi committed
1081
--------------------------------------------------------------------------------
Francois Gygi committed
1082
rel1_18_0 2004-09-14 (not for distribution)
Francois Gygi committed
1083
BOSampleStepper modified for scf iterations. Update of charge density and
Francois Gygi committed
1084 1085 1086 1087
calculation of VHXC is taken out of EnergyFunctional.
RunCmd modified. Use run niter nitscf nite.
Added FermiTemp variable.
Modifications for calculation of occupation numbers if nempty > 0.
Francois Gygi committed
1088
Added USE_XERCES macro. If not defined, the Xerces library is not used, and
Francois Gygi committed
1089 1090 1091 1092
only species files can be read. Loading of sample files without the Xerces lib
is not implemented yet.
SlaterDet.C modified to remove dependency on xercesc. Uses Base64Transcoder.
Makefile modified for libqb.a target for BGL.
Francois Gygi committed
1093 1094
Known bug: SpeciesReader.C incorrect if not USE_XERCES
Known bug: Base64 encoding not done in SlaterDet.C
Francois Gygi committed
1095
--------------------------------------------------------------------------------
1096
rel1_17_1 2004-08-18
Francois Gygi committed
1097
Fixed bug #27: overrun of variable tmpfion in NonLocalPotential.C is due to
Francois Gygi committed
1098 1099 1100 1101
incorrect index calculation in situations where na[is] % (ntasks/nrowmax) != 0.
This is the case for (CdSe)_110 on 128 tasks, nrowmax=32, since na[is]/4 = 27.5.
Note: bug #27 does not affect calculations with a single non-local projector
(npr[is]==1).
1102
Fixed bug #33 LDAFunctional.C: incorrect use of Perdew-Zunger caused small
1103 1104 1105 1106 1107
differences (~1.e-3 a.u.) in total energy when compared to other codes.
The adjustment of the constant gamma to make the PZ energy function continuous
was incorrect. Instead, the constants C and D must be adjusted to ensure
continuity of Ec and its derivative at rs=1.
--------------------------------------------------------------------------------
Francois Gygi committed
1108
rel1_17_0 2004-08-11
Francois Gygi committed
1109
Fixed buffer overrun in Context::string_bcast. (note: bug #27 is not fixed).
Francois Gygi committed
1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121
Fixed bug 26: inconsistent qbox xmlns namespace in species file and qbox output.
Modified qb.C: use new qbox_xmlns() function to print qbox XML namespace.
EnergyFunctional.C:909: change lda argument to 14. This was done automatically
by the BLACS lib, but corrected line conforms to the BLACS convention that
lda >= m.
ChargeDensity.C: changed the <total_electronic_charge> output to an XML comment.
SaveCmd.C: include proper xmlns info when saving a sample.
blas.h: reorganized using the ADD_ macro, removed all platform-specific ifdefs.
XCPotential.C: removed underscores in blas calls.
SlaterDet.C: removed underscores in blas calls.
FourierTransform.C: use ADD_ for zdscal.
NonLocalPotential.C: removed underscores in blas calls, removed old OSF ifdefs.
Francois Gygi committed
1122
--------------------------------------------------------------------------------
Francois Gygi committed
1123 1124 1125 1126
rel1_16_0
Implemented l=1,l=2 non-local stress.
Modified UnitCell.C to avoid infinite loops due to finite precision.
Modified mcr.mk to use library mkl_p4 instead of mkl: -lmkl_p4 before -lmkl
Francois Gygi committed
1127
This change of library accelerates the si432 calculation by 12% overall
Francois Gygi committed
1128
--------------------------------------------------------------------------------
Francois Gygi committed
1129
rel1_15_2
Francois Gygi committed
1130
Ecut.h: exit without any action if ecut is same as current value. Caused
Francois Gygi committed
1131 1132 1133 1134
zeroing of the wavefunction velocity when used in a redundant manner.
Added isodate.C/isodate.h and release.C/release.h files for separate access
to these functions from any other object.
Added a <description> tag in the sample file in SaveCmd.C including the ISO
Francois Gygi committed
1135 1136
date and the release of qbox that was used to create the file.
(Note: the tag was already allowed by the Schema specification, but the
Francois Gygi committed
1137 1138
save command did not write any description information.)
--------------------------------------------------------------------------------
Francois Gygi committed
1139 1140 1141
rel1_15_1
Fixed bug: stop on assert in NonLocalPotential for l=1 or l=2 projectors
Fixed bug: Matrix::init_size() with zero dimensions.
Francois Gygi committed
1142 1143
Known bugs:
-wf_diag variable is ignored.
Francois Gygi committed
1144
--------------------------------------------------------------------------------
Francois Gygi committed
1145 1146 1147 1148 1149
rel1_15_0
Added R option for the cell_lock variable: rescale unit cell only (keep aspect
ratio fixed).
Corrected inconsistencies between BOSampleStepper.C and CPSampleStepper.C as
regards printing of variables during a run: unit_cell, enthalpy, pv.
Francois Gygi committed
1150 1151 1152 1153
Known bugs:
-Stops on assert in NonLocalPotential::update_twnl() for l=1 or l=2 projectors
-Crashes when defining a species without using it, i.e. species with zero atoms.
 Traced to incorrect initialization of sizes in Matrix.C
1154
--------------------------------------------------------------------------------
Francois Gygi committed
1155
rel1_14_3
Francois Gygi committed
1156 1157 1158
Fixed calculation of ekin_e in MDWavefunctionStepper.C to avoid small error
when printing ekin_e at the first iteration of a CP run. Note: only the printed
ekin_e is affected, wfs are correct.
Francois Gygi committed
1159
--------------------------------------------------------------------------------
Francois Gygi committed
1160 1161 1162
rel1_14_2
Fixed Thermostat in MDIonicStepper.C
Fixed use of ionic velocities at first step in CPSampleStepper.C
Francois Gygi committed
1163 1164 1165 1166 1167 1168
Known bug:
Calculation of ekin_e in MDWavefunctionStepper.C at the first iteration uses
the wf velocity at t=0. This is second order accurate, but not consistent with
the calculation of ekin_e in subsequent steps, which introduces a small error
in the printed ekin_e at the first iteration. Ekin_e should be calculated using
a first-order expression, i.e. (wf0-wfm)/dt to be consistent.
Francois Gygi committed
1169 1170
--------------------------------------------------------------------------------
rel1_14_1a
Francois Gygi committed
1171 1172 1173 1174 1175
SlaterDet.C: Removed messages from resize about cell being enclosed by ref_cell.
Note: the determination of enclosure should be revised.
SDCellStepper.C: removed debug messages. New "S" option to preserve shape.
BOSampleStepper.C: Added output of cell dimensions if cell is not locked.
Preconditioner.C: bug fixed in calculation of diag when using confinement.
Francois Gygi committed
1176
Known bugs:
Francois Gygi committed
1177 1178
Does not use ionic velocities at first step of CP dynamics.
Thermostat is not working.
Francois Gygi committed
1179 1180
--------------------------------------------------------------------------------
rel1_14_0
Francois Gygi committed
1181
Matrix.h Reintroduced the clear() call in the resize members. Working with
Francois Gygi committed
1182 1183 1184 1185
uninitialized matrices sometimes caused NaNs to appear in the padding areas
of the danl matrix in NonLocalPotential.C.
Note that the clear() call was removed in rel1_13_0, which invalidates 1_13_0.
Matrix.C: Added assert()'s in clear() functions.
1186 1187
BOSampleStepper.C: Introduced quadratic wf extrapolation option (switch
in BOSampleStepper.C, currently disabled in 1.14.0. Added calculation
Francois Gygi committed
1188
of forces at each electronic step to monitor the convergence
1189
of the force x velocity product (disabled in 1.14.0).
Francois Gygi committed
1190
PSDWavefunctionStepper.C: Fixed bug in usage of preconditioner (affects
Francois Gygi committed
1191 1192
1_13_0 only).
PSDAWavefunctionStepper.C: Fixed bug in usage of preconditioner (see above).
Francois Gygi committed
1193
SlaterDet.[Ch]: introduced new orthogonalization and alignment functions:lowdin,
Francois Gygi committed
1194 1195 1196 1197 1198
align, ortho_align.


--------------------------------------------------------------------------------
rel1_13_0 (not for release, see Matrix clear() bug above)
1199 1200 1201 1202 1203 1204 1205 1206 1207
Implemented stress calculation and variable cell dynamics
Stress calculation works for s-only non-local potentials.
Changed names of variables tau0, taup to r0, rp (and keep tau for scaled
coordinates)
modified the IonicStepper interface to separate the calculation of rp
from the update r0 := rp.
Removed a bug in calculation of eta at first step (was not using ekin_ion)
--------------------------------------------------------------------------------
rel1_13_0a1 (not for release)
Francois Gygi committed
1208
Computes stress tensor. Limited to s-only n-l projectors. No cell eq. of
1209 1210
motion yet. Tested against GP results in cubic Si64.
--------------------------------------------------------------------------------
Francois Gygi committed
1211 1212 1213
rel1_12_1
Removed Gauss-Jacobi quadrature in NonLocalPotential introduced in 1_12_0 and
reverted to trapezoidal rule. Tests on Si64 using trapezoidal, Simpson, Gauss-
Francois Gygi committed
1214 1215 1216 1217 1218
Legendre and Gauss-Jacobi show that the trapezoidal rule is more efficient
than all other methods, likely because higher order methods are not in their
asymptotic range for reasonable numbers of points  (4-32). A possible
improvement over the present status may be to use Richardson extrapolation,
since the trapezoidal rule has O(h^2) error. However, from the example of Si64,
Francois Gygi committed
1219 1220
little improvement is obtained for small (4-8) numbers of points.
--------------------------------------------------------------------------------
Francois Gygi committed
1221
rel1_12_0
Francois Gygi committed
1222
Replaced trapezoidal integration in NonLocalPotential.C by a Gauss-Jacobi rule.
Francois Gygi committed
1223
(Note: this was later found to be less efficient. Removed in 1.12.1e)
Francois Gygi committed
1224 1225 1226
Fixed bug in PSDWavefunctionStepper.C and PSDAWavefunctionStepper.C:
use ecut if ecutprec==0.0. Caused nans when ecutprec was not set.
--------------------------------------------------------------------------------
1227
rel1_11_2 2003-12-04
Francois Gygi committed
1228
Fixed Elan allocation failure problem in SlaterDet::write. Problem was caused
1229 1230 1231
by the large number of buffers allocated by the elan lib for all the messages
sent to task 0 before a receive was posted for them. A barrier in the loop
over states limits the number of messages sent to task 0 at any given time.
Francois Gygi committed
1232 1233
Note: the problem likely doesn't show on AIX since the MPI implementation may
limit the number of buffered received messages. The barrier is added in the
1234 1235
AIX code too (SlaterDet::print).
--------------------------------------------------------------------------------
Francois Gygi committed
1236 1237 1238
rel1_11_1
Fixed uninitialized var eta_ in MDIonicStepper that caused errors on AIX.
Known problem: elan allocation failure on mcr when writing h2o64 sample
Francois Gygi committed
1239
using 140 tasks on 70 nodes. Possible bug in SlaterDet::write if USE_CSTDIO_LFS.
Francois Gygi committed
1240
When running on AIX (using SlaterDet::print instead of SlaterDet::write) no
Francois Gygi committed
1241
problem arises, and memory usage per node is low.(fixed in 1_11_2).
Francois Gygi committed
1242
--------------------------------------------------------------------------------
Francois Gygi committed
1243
rel1_11_0
Francois Gygi committed
1244
BOSampleStepper: changed to do extrapolation and ionic move before
Francois Gygi committed
1245 1246
smoothing the wf. Uses wf velocity if available. Creates it if not
Included PSDA wf stepper.
Francois Gygi committed
1247
BO dyn Works ok with PSDA.
Francois Gygi committed
1248 1249 1250 1251
Tested 1pe vs 2pe on CH4 BO dynamics: ok.
Corrected bug in SampleReader.C (reading of wf velocity, appeared in rel1_9_0)
--------------------------------------------------------------------------------
rel1_10_0 (not for release) (Note: known bugs)
1252 1253 1254
Complete rewrite of the SampleStepper hierarchy, including CP and BO classes.
Runs CP and BO/SD and BO/PSD tests ok. No extrapolation of wavefunctions yet.
Note: wavefunction extrapolation implemented yet.
Francois Gygi committed
1255
Planned work: interchange nite iterations and ionic move in BO stepper, so
1256 1257
that the wfs at the end of the iteration are the GS for the ionic positions at
the same time.
Francois Gygi committed
1258 1259
XML files: a <description> tag (optional) was added to the sample.xsd file
to include any information about the sample that a user may want to include.
Francois Gygi committed
1260
Bug in SampleReader.C: reading of density matrix for the
Francois Gygi committed
1261 1262
wavefunction velocity is incorrect. Causes out of sync MPI call and Abort
Corrected in 1_11_0 (candidate)
Francois Gygi committed
1263
--------------------------------------------------------------------------------
Francois Gygi committed
1264
rel1_9_0 (Note: known bugs)
Francois Gygi committed
1265
Large file support (LFS): LFS is not correctly supported for <iostream> in icc
Francois Gygi committed
1266
and gcc 3.2. A workaround is implemented in SaveCmd.C, Wavefunction.C,
Francois Gygi committed
1267 1268 1269 1270
SlaterDet.C. The workaround is activated by defining the USE_CSTDIO_LFS macro.
Occupation numbers moved to the SlaterDet class.
Occupation numbers are not fixed to 2.0 anymore. Load and save now use
occupation numbers.
Francois Gygi committed
1271
Bug in SampleReader.C: reading of density matrix for the
Francois Gygi committed
1272 1273
wavefunction velocity is incorrect. Causes out of sync MPI call and Abort
Corrected in 1_11_0 (candidate)
Francois Gygi committed
1274 1275
--------------------------------------------------------------------------------
rel1_8_2
Francois Gygi committed
1276 1277
Fixed problem with wfvtmp in SampleReader: cannot load sample and
then set wf_dyn = MD. Also, the order of load and set wf_dyn MD commands causes
Francois Gygi committed
1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299
a seg fault.
--------------------------------------------------------------------------------
rel1_8_1
Fixed large-file support for sample files.
Fixed resizing of wf velocity.
Defined Matrix::operator=.
Fixed bug in PSD optimization: G=0 coefficient had zero value of dt2bye.
Fixed bug in UnitCell::fold_in_ws: caused infinite loop in some cases.
--------------------------------------------------------------------------------
rel1_8_0
Added thermostat, th_time, th_temp.
Added emass variable.
--------------------------------------------------------------------------------
rel1_7_1
Fixed bug in calculation of electronic kinetic energy in SampleStepper.C which
caused small kink in ekin_e and spike in econst during restarts.
--------------------------------------------------------------------------------
rel1_7_0
Implemented parallel XML preprocessing of sample files for fast loading.
Added -serial option to the load command to force serial parsing of files.
Fixed a bug in Context::string_send: was incorrect for zero-length strings
initialization of ibuf removed, unnecessary.
Francois Gygi committed
1300
Added new dmin,dmax,imin,imax,string_send,string_recv,pmap functions in
Francois Gygi committed
1301 1302 1303 1304 1305 1306 1307 1308 1309
Context
Changed the block distribution of the slabs in FourierTransform to adopt the
same type of block distribution used in ScaLAPACK. Needed in SampleReader.C
in order to read from the gfdata matrix.
--------------------------------------------------------------------------------
rel1_6_1exp
2003-07-24:modified initialization of coefficients in SlaterDet constructor
to choose g vectors of minimum norm. Note: reverting to -O2 option using xlC
after noting that -O3 causes incorrect results in testBasis.
Francois Gygi committed
1310
Added dmin,dmax,imin,imax members to Context.C to identify process where
Francois Gygi committed
1311 1312 1313
max or min is found.
--------------------------------------------------------------------------------
rel1_6_0
Francois Gygi committed
1314
Implemented MD for wfs and for ions, print ekin_e, ekin_ion, econst in
Francois Gygi committed
1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349
SampleStepper.C. Use Stoermer's rule to compute velocities at end points.
Implemented load/save of wf velocities
--------------------------------------------------------------------------------
rel1_5_1
Fixed bug in NonLocalPotential that appeared for Contexts having multiple
columns and rows (computation of enl and fion incorrect).
--------------------------------------------------------------------------------
rel1_5_0    May 21, 2003
Parallel Load/Save XML sample implemented.
Load is slow: 6 min for 256 states on a (64)^3 grid (H2O_64, 73 Ry)
Save is ok (30 s.)
Bug found in NonLocalPotential::update_eigr: index of first local atom not
computed. Fixed. Affects calculations with multiple-column contexts.
Fixed bug in NonLocalPotential::energy: loop over fnl elements.
Added optional ecut attribute to the wavefunction element to insure
restart from xml sample in exact same conditions. If ecut attribute is not
specified in xml file, a value is inferred from the grid size.
Fixed incorrect initialization of dwf in EnergyFunctional::energy.
Implemented PSD for wf_dyn.
--------------------------------------------------------------------------------
rel1_4_0    May 16, 2003
Implemented byteswap operation when loading/saving XML base64 file on AIX.
(affects WavefunctionHandler.[Ch] and SlaterDet.[Ch].
All storage in base64 form is therefore little-endian. No byte swapping needed
on Intel platforms.
--------------------------------------------------------------------------------
rel1_3_0    May 16, 2003
New features in rel1_3_0:
XML input for sample, species, atomset, wavefunctions, including base64
encoding. New hiearchy of XML handler classes for parsing of sample files.
Parsing done with Apache xerces-c, 2.2.0.
Clean-up of FourierTransform: convention of fwd/bwd changed to be consistent
with usually accepted definitions (forward transform is Fourier analysis, i.e.
computation of Fourier coefficients, and involves scaling by 1/np012. Backward
transform is computation of real-space function from coefficients, no scaling).
Francois Gygi committed
1350
Load command includes reading of sample and appropriate resizing of all
Francois Gygi committed
1351 1352
variables.
The "locked" variable is removed from the Atom class, and will be reintroduced
Francois Gygi committed
1353
as a Constraint. This is to avoid having to include the "locked" info in the
Francois Gygi committed
1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367
XML definition of an Atom.
Builds and runs on Linux (linuc-pc_mpi_icc.mk) and on snow (aix_mpi.mk)
--------------------------------------------------------------------------------
Jan 24, 2003
Checking in release rel1_0_0:
performs SD steps with LDA and local pseudopotentials
--------------------------------------------------------------------------------
Thu Jan  9 16:01:42 PST 2003
Added UserInterface to qb. Commands Set, Print, Quit, Atom.
TO BE FIXED: SlaterDet cannot resize if more states than basis functions.
This forces the order in which ecut, cell and atoms are defined.
--------------------------------------------------------------------------------
Tue Sep 24 14:19:13 PDT 2002
SlaterDet can exist on any context.
Francois Gygi committed
1368
Developing a Wavefunction that can hold any number of SlaterDets for
Francois Gygi committed
1369 1370
arbitrary nspin and nkpoints.
--------------------------------------------------------------------------------