638d667 | Peter Boyle | 31 August 2015, 15:33:20 UTC | Tested rms dH is ~ dt^4 numerically, so believe the ForceGradient is correct now. Paranoia makes me want to diddle with the FG step to ensure dt^2 reappears. | 31 August 2015, 15:33:20 UTC |
357c6ab | Peter Boyle | 31 August 2015, 15:32:04 UTC | Reunitarise. Complete the HMC and integrator changes. | 31 August 2015, 15:32:04 UTC |
755dca9 | Peter Boyle | 31 August 2015, 05:23:02 UTC | Added ForceGradient integrator. dH dropped so seems to work. Will only believe it is right once I have pulled a dt^4 error scaling plot out. | 31 August 2015, 05:23:02 UTC |
29fd004 | Peter Boyle | 30 August 2015, 12:39:19 UTC | Unified integrator and integrator algorithm into virtual class used as a policy for the HMC. | 30 August 2015, 12:39:19 UTC |
eed889e | Peter Boyle | 30 August 2015, 11:23:08 UTC | Update on todo list | 30 August 2015, 11:23:08 UTC |
aa52fda | Peter Boyle | 30 August 2015, 11:18:34 UTC | Global edit on HMC sector -- making GaugeField a template parameter and preparing to pass integrator, smearing, bc's as policy classes to hmc. Propose to unify "integrator" and integrator algorithm in a base/derived way to override step. Want to read through ForceGradient to ensure that abstraction covers the force gradient case. | 30 August 2015, 11:18:34 UTC |
76d7525 | Peter Boyle | 29 August 2015, 16:18:43 UTC | Started a tidy up in the HMC sector. Now comfortable with the two level integrators; to a little figure out what Guido had done & why -- but there is a neat saving of force evaluations across the nesting time boundary making use of linearity of the leapP in dt. I cleaned up the printing, reduced the volume of code, in the process sharing printing between all integrators. Placed an assert that the total integration time for all integrators must match at end of trajectory. Have now verified e-dH = 1 for nested integrators in Wilson/Wilson runs with both Omelyan and with Leapfrog so substantial confidence gained. | 29 August 2015, 16:18:43 UTC |
dc814f3 | Peter Boyle | 26 August 2015, 12:40:29 UTC | Binary IO file for generic Grid array parallel I/O. Number of IO MPI tasks can be varied by selecting which dimensions use parallel IO and which dimensions use Serial send to boss I/O. Thus can neck down from, say 1024 nodes = 4x4x8x8 to {1,8,32,64,128,256,1024} nodes doing the I/O. Interpolates nicely between ALL nodes write their data, a single boss per time-plane in processor space [old UKQCD fortran code did this], and a single node doing all I/O. Not sure I have the transfer sizes big enough and am not overly convinced fstream is guaranteed to not give buffer inconsistencies unless I set streambuf size to zero. Practically it has worked on 8 tasks, 2x1x2x2 writing /cloning NERSC configurations on my MacOS + OpenMPI and Clang environment. It is VERY easy to switch to pwrite at a later date, and also easy to send x-strips around from each node in order to gather bigger chunks at the syscall level. That would push us up to the circa 8x 18*4*8 == 4KB size write chunk, and by taking, say, x/y non parallel we get to 16MB contiguous chunks written in multi 4KB transactions per IOnode in 64^3 lattices for configuration I/O. I suspect this is fine for system performance. | 26 August 2015, 12:40:29 UTC |
612957f | Peter Boyle | 21 August 2015, 09:19:08 UTC | pull in original license. | 21 August 2015, 09:19:08 UTC |
cea8ac9 | Peter Boyle | 21 August 2015, 09:14:53 UTC | Credits to orig source where I found the macro tricks. | 21 August 2015, 09:14:53 UTC |
476da3e | Peter Boyle | 21 August 2015, 09:06:33 UTC | Separated IO reader/writers into a proper abstract base, derived relationship. Have Text/Binary/Xml versions of Reader & Writer. Any new Reader/Writer class inheriting the interface can give object serialisation to any desired format now. new file: lib/serialisation/BaseIO.h modified: lib/serialisation/BinaryIO.h modified: lib/serialisation/Serialisation.h modified: lib/serialisation/TextIO.h modified: lib/serialisation/XmlIO.h The test uses the Xml, Binary and Text formats as well as cout << Object. | 21 August 2015, 09:06:33 UTC |
35818fd | Peter Boyle | 20 August 2015, 22:04:38 UTC | Text and Binary readers | 20 August 2015, 22:04:38 UTC |
091785e | Peter Boyle | 20 August 2015, 16:19:48 UTC | Better list | 20 August 2015, 16:19:48 UTC |
77d299b | Peter Boyle | 20 August 2015, 15:30:52 UTC | Cosmetic | 20 August 2015, 15:30:52 UTC |
ab81a25 | Peter Boyle | 20 August 2015, 15:21:26 UTC | XMLReader implementation and a virtual Reader/Writer template framework. Test_serialisation has an example of *code* *free* object serialisation to both ostream and to XML using macro magic. Implementing TextReader/TextWriter, YAML, JSON etc.. should be trivial and we can use configure time options to select the default "Reader" typedef. Present done with "using XMLPolicy::Reader" to pick up the default serialisation strategy. | 20 August 2015, 15:21:26 UTC |
fdfe194 | Peter Boyle | 19 August 2015, 13:41:05 UTC | Threading bug in RNG fill fixed. | 19 August 2015, 13:41:05 UTC |
8b070ae | Peter Boyle | 19 August 2015, 10:26:01 UTC | Gparity now accepting twists through constructor | 19 August 2015, 10:26:01 UTC |
4e085dd | Peter Boyle | 19 August 2015, 09:26:07 UTC | Domain wall even-odd 2f HMC with wilson gauge and PV 2f ratio now running and giving small dH. Azusa is working hard on the rectangle term and we'll hopefully start reproducing plaquettes from RBC-UKQCD parameters soon ! My new laptop is pretty warm and is starting to groan ;) | 19 August 2015, 09:26:07 UTC |
e8d63c9 | Peter Boyle | 19 August 2015, 04:49:00 UTC | Merge branch 'master' of https://github.com/paboyle/Grid | 19 August 2015, 04:49:00 UTC |
c54c086 | Peter Boyle | 19 August 2015, 04:46:58 UTC | Even odd preconditioned one flavour ratio (no support for non-const EE schur block) | 19 August 2015, 04:46:58 UTC |
dd6bb73 | Peter Boyle | 19 August 2015, 03:58:40 UTC | Added one flavour rational ratios (unprec) | 19 August 2015, 03:58:40 UTC |
fc160ee | Peter Boyle | 19 August 2015, 03:58:40 UTC | Added one flavour rational ratios (unprec) | 19 August 2015, 03:58:40 UTC |
48db722 | Peter Boyle | 18 August 2015, 17:37:39 UTC | EvenOdd schur decomposed mpcdagmpc version of rhmc determinant. dH is also small and plaquette looks right. | 18 August 2015, 17:37:39 UTC |
570150f | Peter Boyle | 18 August 2015, 17:37:39 UTC | EvenOdd schur decomposed mpcdagmpc version of rhmc determinant. dH is also small and plaquette looks right. | 18 August 2015, 17:37:39 UTC |
9c7840c | Peter Boyle | 18 August 2015, 15:58:56 UTC | rhmc for 1+1 wilson is conserving dH~0. A good days work ;) | 18 August 2015, 15:58:56 UTC |
aef98b7 | Peter Boyle | 18 August 2015, 15:58:56 UTC | rhmc for 1+1 wilson is conserving dH~0. A good days work ;) | 18 August 2015, 15:58:56 UTC |
5c364f8 | Peter Boyle | 18 August 2015, 13:40:08 UTC | One flavour rational unprec added; untested but does compile. Moving param structs into a single header for later connection to file I/O using macromagic.h | 18 August 2015, 13:40:08 UTC |
a842a6c | Peter Boyle | 18 August 2015, 13:40:08 UTC | One flavour rational unprec added; untested but does compile. Moving param structs into a single header for later connection to file I/O using macromagic.h | 18 August 2015, 13:40:08 UTC |
2dd9ad7 | Peter Boyle | 18 August 2015, 09:43:32 UTC | Update TODO list | 18 August 2015, 09:43:32 UTC |
cd242a2 | Peter Boyle | 18 August 2015, 09:43:32 UTC | Update TODO list | 18 August 2015, 09:43:32 UTC |
9306921 | Peter Boyle | 18 August 2015, 09:37:08 UTC | Even Odd two flavour ratio added and dH == small | 18 August 2015, 09:37:08 UTC |
bdcbfe9 | Peter Boyle | 18 August 2015, 09:37:08 UTC | Even Odd two flavour ratio added and dH == small | 18 August 2015, 09:37:08 UTC |
76f3855 | Peter Boyle | 18 August 2015, 08:23:58 UTC | Merge branch 'master' of https://github.com/paboyle/Grid | 18 August 2015, 08:23:58 UTC |
8621e24 | Peter Boyle | 18 August 2015, 08:23:58 UTC | Merge branch 'master' of https://github.com/paboyle/Grid | 18 August 2015, 08:23:58 UTC |
6212807 | Peter Boyle | 18 August 2015, 08:21:29 UTC | Small dh obtained in two flavour ratio so looks ok. | 18 August 2015, 08:21:29 UTC |
7622f0c | Peter Boyle | 18 August 2015, 08:21:29 UTC | Small dh obtained in two flavour ratio so looks ok. | 18 August 2015, 08:21:29 UTC |
0bc38a6 | Peter Boyle | 18 August 2015, 08:19:42 UTC | Adding PV pseudofermion in prep for DWF HMC. Not compiled this yet, but cloned in from BFM. | 18 August 2015, 08:19:42 UTC |
25d0eae | Peter Boyle | 18 August 2015, 08:19:42 UTC | Adding PV pseudofermion in prep for DWF HMC. Not compiled this yet, but cloned in from BFM. | 18 August 2015, 08:19:42 UTC |
24382d7 | Peter Boyle | 17 August 2015, 22:14:48 UTC | Adding PV pseudofermion in prep for DWF HMC. Not compiled this yet, but cloned in from BFM. | 17 August 2015, 22:14:48 UTC |
ef6a9e6 | Peter Boyle | 17 August 2015, 22:14:48 UTC | Adding PV pseudofermion in prep for DWF HMC. Not compiled this yet, but cloned in from BFM. | 17 August 2015, 22:14:48 UTC |
353d66d | Peter Boyle | 16 August 2015, 00:41:05 UTC | Unused apparently | 16 August 2015, 00:41:05 UTC |
b8166af | Peter Boyle | 16 August 2015, 00:41:05 UTC | Unused apparently | 16 August 2015, 00:41:05 UTC |
6180487 | Peter Boyle | 15 August 2015, 23:14:10 UTC | Tidying | 15 August 2015, 23:14:10 UTC |
afeabe0 | Peter Boyle | 15 August 2015, 23:14:10 UTC | Tidying | 15 August 2015, 23:14:10 UTC |
0e088d2 | Peter Boyle | 15 August 2015, 23:13:14 UTC | Merge branch 'master' of https://github.com/paboyle/Grid | 15 August 2015, 23:13:14 UTC |
2d6b97b | Peter Boyle | 15 August 2015, 23:13:14 UTC | Merge branch 'master' of https://github.com/paboyle/Grid | 15 August 2015, 23:13:14 UTC |
53da927 | Peter Boyle | 15 August 2015, 22:59:04 UTC | Merge branch 'master' of https://github.com/paboyle/Grid | 15 August 2015, 22:59:04 UTC |
f0e32f1 | Peter Boyle | 15 August 2015, 22:59:04 UTC | Merge branch 'master' of https://github.com/paboyle/Grid | 15 August 2015, 22:59:04 UTC |
c7b50d1 | Peter Boyle | 15 August 2015, 22:56:31 UTC | Merge branch 'master' of https://github.com/paboyle/Grid | 15 August 2015, 22:56:31 UTC |
155c164 | Peter Boyle | 15 August 2015, 22:25:49 UTC | * Finished the template/policy style introduction of gparity, except the gparity force terms. So valence sector looks ok. FermionOperatorImpl.h provides the policy classes. Expect HMC will introduce a smearing policy and a fermion representation change policy template param. Will also probably need multi-precision work. * HMC is running even-odd and non-checkerboarded (checked 4^4 wilson fermion/wilson gauge). There appears to be a bug in the multi-level integrator -- <e-dH> passes with single level but not with multi-level. In any case there looks to be quite a bit to clean up. This is the "const det" style implementation that is not appropriate yet for clover since it assumes that Mee is indept of the gauge fields. Easily fixed in future. | 15 August 2015, 22:25:49 UTC |
55cfc89 | Peter Boyle | 15 August 2015, 22:25:49 UTC | * Finished the template/policy style introduction of gparity, except the gparity force terms. So valence sector looks ok. FermionOperatorImpl.h provides the policy classes. Expect HMC will introduce a smearing policy and a fermion representation change policy template param. Will also probably need multi-precision work. * HMC is running even-odd and non-checkerboarded (checked 4^4 wilson fermion/wilson gauge). There appears to be a bug in the multi-level integrator -- <e-dH> passes with single level but not with multi-level. In any case there looks to be quite a bit to clean up. This is the "const det" style implementation that is not appropriate yet for clover since it assumes that Mee is indept of the gauge fields. Easily fixed in future. | 15 August 2015, 22:25:49 UTC |
ba8c09a | Peter Boyle | 14 August 2015, 13:16:45 UTC | Reorganising the Fermion interface | 14 August 2015, 13:16:45 UTC |
f40475f | Peter Boyle | 14 August 2015, 13:16:45 UTC | Reorganising the Fermion interface | 14 August 2015, 13:16:45 UTC |
045c858 | Peter Boyle | 14 August 2015, 12:18:59 UTC | Extra test | 14 August 2015, 12:18:59 UTC |
b3b46fd | Peter Boyle | 14 August 2015, 12:18:59 UTC | Extra test | 14 August 2015, 12:18:59 UTC |
e846279 | Peter Boyle | 14 August 2015, 12:18:59 UTC | Extra test | 14 August 2015, 12:18:59 UTC |
4dc7c36 | Peter Boyle | 14 August 2015, 11:57:42 UTC | Gparity works now even if simd distributed in a Gparity twist direction. Tested by doubling lattice in t-direction. | 14 August 2015, 11:57:42 UTC |
59d66eb | Peter Boyle | 14 August 2015, 11:57:42 UTC | Gparity works now even if simd distributed in a Gparity twist direction. Tested by doubling lattice in t-direction. | 14 August 2015, 11:57:42 UTC |
cc63078 | Peter Boyle | 14 August 2015, 11:57:42 UTC | Gparity works now even if simd distributed in a Gparity twist direction. Tested by doubling lattice in t-direction. | 14 August 2015, 11:57:42 UTC |
028e206 | Peter Boyle | 13 August 2015, 23:01:04 UTC | Gparity valence test now working. Interface in FermionOperator will change a lot in future | 13 August 2015, 23:01:04 UTC |
7d3512a | Peter Boyle | 13 August 2015, 23:01:04 UTC | Gparity valence test now working. Interface in FermionOperator will change a lot in future | 13 August 2015, 23:01:04 UTC |
e6bed00 | Peter Boyle | 13 August 2015, 23:01:04 UTC | Gparity valence test now working. Interface in FermionOperator will change a lot in future | 13 August 2015, 23:01:04 UTC |
2c216a4 | Peter Boyle | 13 August 2015, 09:51:29 UTC | Gamma5 mult direct | 13 August 2015, 09:51:29 UTC |
45b0185 | Peter Boyle | 13 August 2015, 09:51:29 UTC | Gamma5 mult direct | 13 August 2015, 09:51:29 UTC |
fc9b36c | Peter Boyle | 13 August 2015, 09:51:29 UTC | Gamma5 mult direct | 13 August 2015, 09:51:29 UTC |
145b807 | Peter Boyle | 13 August 2015, 09:51:01 UTC | Gparity improvements | 13 August 2015, 09:51:01 UTC |
1c2d148 | Peter Boyle | 13 August 2015, 09:51:01 UTC | Gparity improvements | 13 August 2015, 09:51:01 UTC |
c390781 | Peter Boyle | 13 August 2015, 09:51:01 UTC | Gparity improvements | 13 August 2015, 09:51:01 UTC |
5465138 | Peter Boyle | 13 August 2015, 09:50:34 UTC | Some bug fixes for more complicated types introduced with gparity | 13 August 2015, 09:50:34 UTC |
7e9203d | Peter Boyle | 13 August 2015, 09:50:34 UTC | Some bug fixes for more complicated types introduced with gparity | 13 August 2015, 09:50:34 UTC |
8d4c433 | Peter Boyle | 13 August 2015, 09:50:34 UTC | Some bug fixes for more complicated types introduced with gparity | 13 August 2015, 09:50:34 UTC |
6ab73c5 | Peter Boyle | 12 August 2015, 08:49:33 UTC | Gparity test added; partial implementation -- this is Chris K's doubled lattice only and have to regress this with the 2 flavour implementation. | 12 August 2015, 08:49:33 UTC |
8a0be42 | Peter Boyle | 12 August 2015, 08:49:33 UTC | Gparity test added; partial implementation -- this is Chris K's doubled lattice only and have to regress this with the 2 flavour implementation. | 12 August 2015, 08:49:33 UTC |
9183380 | Peter Boyle | 12 August 2015, 08:49:33 UTC | Gparity test added; partial implementation -- this is Chris K's doubled lattice only and have to regress this with the 2 flavour implementation. | 12 August 2015, 08:49:33 UTC |
04e0e9f | Peter Boyle | 11 August 2015, 05:37:42 UTC | File list update. | 11 August 2015, 05:37:42 UTC |
c8dca58 | Peter Boyle | 11 August 2015, 05:37:42 UTC | File list update. | 11 August 2015, 05:37:42 UTC |
ded3945 | Peter Boyle | 11 August 2015, 05:37:42 UTC | File list update. | 11 August 2015, 05:37:42 UTC |
69ce87f | Peter Boyle | 11 August 2015, 05:24:53 UTC | Preconditioned conjugate residual | 11 August 2015, 05:24:53 UTC |
826fbb1 | Peter Boyle | 11 August 2015, 05:24:53 UTC | Preconditioned conjugate residual | 11 August 2015, 05:24:53 UTC |
9cd7f9e | Peter Boyle | 11 August 2015, 05:24:53 UTC | Preconditioned conjugate residual | 11 August 2015, 05:24:53 UTC |
07d672b | Peter Boyle | 11 August 2015, 05:23:38 UTC | Header | 11 August 2015, 05:23:38 UTC |
26f5ee0 | Peter Boyle | 11 August 2015, 05:23:38 UTC | Header | 11 August 2015, 05:23:38 UTC |
f165b1a | Peter Boyle | 11 August 2015, 05:23:38 UTC | Header | 11 August 2015, 05:23:38 UTC |
0a9ebac | Peter Boyle | 11 August 2015, 05:22:20 UTC | Gparity modifications in the Gparity compressor variant. | 11 August 2015, 05:22:20 UTC |
3903dfe | Peter Boyle | 11 August 2015, 05:22:20 UTC | Gparity modifications in the Gparity compressor variant. | 11 August 2015, 05:22:20 UTC |
881acaa | Peter Boyle | 11 August 2015, 05:22:20 UTC | Gparity modifications in the Gparity compressor variant. | 11 August 2015, 05:22:20 UTC |
1b3c93e | Peter Boyle | 10 August 2015, 19:47:44 UTC | Rework/global edit to enforce type templating of fermion operators. Allows multi-precision work and paves the way for alternate BC's and such like allowing for example G-parity which is important for K pipi programme. In particular, can drive an extra flavour index into the fermion fields using template types. | 10 August 2015, 19:47:44 UTC |
84a6647 | Peter Boyle | 10 August 2015, 19:47:44 UTC | Rework/global edit to enforce type templating of fermion operators. Allows multi-precision work and paves the way for alternate BC's and such like allowing for example G-parity which is important for K pipi programme. In particular, can drive an extra flavour index into the fermion fields using template types. | 10 August 2015, 19:47:44 UTC |
aeb7442 | Peter Boyle | 10 August 2015, 19:47:44 UTC | Rework/global edit to enforce type templating of fermion operators. Allows multi-precision work and paves the way for alternate BC's and such like allowing for example G-parity which is important for K pipi programme. In particular, can drive an extra flavour index into the fermion fields using template types. | 10 August 2015, 19:47:44 UTC |
2be8df9 | Peter Boyle | 07 August 2015, 07:37:15 UTC | Adding components for even odd decomposed determinant in HMC. dH not yet conserved, so something wrong in the eo force code still | 07 August 2015, 07:37:15 UTC |
a01aa15 | Peter Boyle | 07 August 2015, 07:37:15 UTC | Adding components for even odd decomposed determinant in HMC. dH not yet conserved, so something wrong in the eo force code still | 07 August 2015, 07:37:15 UTC |
ce34856 | Peter Boyle | 07 August 2015, 07:37:15 UTC | Adding components for even odd decomposed determinant in HMC. dH not yet conserved, so something wrong in the eo force code still | 07 August 2015, 07:37:15 UTC |
6ec087d | Peter Boyle | 01 August 2015, 13:48:21 UTC | Continued fraction overlap, partial fraction overlap force terms have a successful test passing. | 01 August 2015, 13:48:21 UTC |
b5a483a | Peter Boyle | 01 August 2015, 13:48:21 UTC | Continued fraction overlap, partial fraction overlap force terms have a successful test passing. | 01 August 2015, 13:48:21 UTC |
d98e836 | Peter Boyle | 01 August 2015, 13:48:21 UTC | Continued fraction overlap, partial fraction overlap force terms have a successful test passing. | 01 August 2015, 13:48:21 UTC |
742db5d | Peter Boyle | 01 August 2015, 13:30:00 UTC | Merge problem fixed | 01 August 2015, 13:30:00 UTC |
772cd81 | Peter Boyle | 01 August 2015, 13:30:00 UTC | Merge problem fixed | 01 August 2015, 13:30:00 UTC |
bb372a6 | Peter Boyle | 01 August 2015, 13:30:00 UTC | Merge problem fixed | 01 August 2015, 13:30:00 UTC |
35feb93 | Peter Boyle | 01 August 2015, 13:24:54 UTC | Merge branch 'master' of https://github.com/paboyle/Grid Conflicts: lib/Make.inc lib/qcd/hmc/HMC.h tests/Make.inc tests/Test_hmc_WilsonFermionGauge.cc | 01 August 2015, 13:24:54 UTC |
5e9bef8 | Peter Boyle | 01 August 2015, 13:24:54 UTC | Merge branch 'master' of https://github.com/paboyle/Grid Conflicts: lib/Make.inc lib/qcd/hmc/HMC.h tests/Make.inc tests/Test_hmc_WilsonFermionGauge.cc | 01 August 2015, 13:24:54 UTC |