Michael Marshall
2b1fcd78c3
Fixes post review with Peter: a) Correct bug in isRegularShape - detect 3d matrix where 1st slice is 2x2 and second slice is 2x1; b) Synchronisation of EigenResizeCounter done by checking we're the OMP primary thread; c) Move definition of EigenResizeCounter to new file, BaseIO.cc
2021-05-31 22:24:54 +01:00
Michael Marshall
0a4e0b49a0
BaseIO: Added "EigenResizeCounter" to keep track of any allocations/deallocations to Eigen tensors during readback. On read, if the tensor is resized, EigenResizeCounter += delta memory (in bytes)
2021-05-31 12:49:56 +01:00
Michael Marshall
76af169f05
Add global namespace to Writer<T> and Reader<T> inside GRID_SERIALIZABLE_CLASS_MEMBERS (so that "using Grid" not necessary).
...
Fix issue with output of Grid::iMatrix so that M<3>{{148,149,150,} {151,152,153,} {154155156}} becomes M<3>{{148,149,150} {151,152,153} {154,155,156}}
2021-05-31 08:43:02 +01:00
Michael Marshall
7b89232251
Extended HDF5 serialisation of std::vector<T> where T now also includes Grid scalar/vector/matrix
...
Changed VectorUtils element traits to is_flattenable, because: a) contract changed on what it does; and b) no other Grid dependencies on element. Needs review.
Initial tests work ... needs proper regression testing.
2021-05-30 20:27:53 +01:00
Peter Boyle
b5aeae526f
Make Cshift fields static to avoid repeated reallocaate overhead
2021-05-28 16:33:08 +02:00
Michael Marshall
ef0ddd5d04
std::vector serialisation in hdf5 uses a different format if the vector is ragged. When reading back std::vector we need to check which format we're reading (since we don't know a priori) and this involves looking for attributes that may not exist. The c++ API: a) throws; and b) prints voluminous logging. Switched to non-throwing, non-logging, C version of the API after code review.
2021-05-24 18:43:55 +01:00
Michael Marshall
9b73dacf50
First row might still be ragged if multi dimensional. attrExists() doesn't throw, but easier to wrap in try ... catch than to explain in comment.
2021-05-22 04:34:32 +01:00
Michael Marshall
244b4aa07f
Serialise std::vector of numeric types as multidimensional object if size is regular ... or individually if ragged
2021-05-21 20:08:56 +01:00
Christopher Kelly
86f08c6b9a
Added a check that the initial EOFA action agrees with |eta|^2, thus checking the quality of the rational approximation in the heatbath
2021-05-18 13:57:44 -04:00
Christopher Kelly
9f0271039f
Completed implementation of Meofa method of ExactOneFlavourRatio pseudofermion action
...
Added tests to tests/forces/Test_mobius_force_eofa.cc testing that the EOFA heatbath results in Phi = M^{-1/2} eta
2021-05-18 12:27:51 -04:00
Christopher Kelly
24df770f74
Added tests/IO/Test_field_array_io.cc testing/demonstrating parallel IO of an array of 5D fermion fields
2021-05-13 12:32:45 -04:00
Christopher Kelly
45b6c7effc
Added a test code forces/Test_gpdwf_force_1f_2f that compares the action and force for DWF, EOFA and DSDR actions between the 1f and 2f implementations of G-parity BCs
...
Broke up ExactOneFlavourRatio refresh into a virtual routine that generates eta and one that uses it as with the ratio and RHMC actions
Added accessors to the pseudofermion field to TwoFlavourEvenOddRatio and ExactOneFlavourRatio
2021-05-12 16:34:07 -04:00
Quadro
1c70d8c4d9
Warning remove
2021-05-05 19:56:04 -04:00
Quadro
f0e9a5299f
Happy on GCC I hope
2021-05-05 19:55:34 -04:00
Quadro
f1b8ba45e7
Warning on GCC suppress unrelated to my code so why doesn't it shut up about its ABI fix
2021-05-05 19:54:21 -04:00
Peter Boyle
fe998ab578
Merge branch 'feature/gparity_HMC' of https://github.com/paboyle/Grid into feature/gparity_HMC
2021-05-05 17:36:51 -04:00
Peter Boyle
c2ee2b5fd1
Random chhanges
2021-05-05 17:36:38 -04:00
Peter Boyle
3b734ee397
two point function example
2021-05-05 17:36:19 -04:00
Peter Boyle
8637a9512a
Freeze Gaussian implementation
2021-05-05 17:34:54 -04:00
Peter Boyle
7f6e2ee03e
Drop normal_distribution, standardise
2021-05-05 17:34:17 -04:00
u61464
8cfc7342cd
staggered hand unroll read coalesce
2021-05-05 14:17:18 -07:00
Peter Boyle
7b02acb2bd
Merge branch 'feature/gparity_HMC' of https://github.com/paboyle/Grid into feature/gparity_HMC
2021-05-04 13:45:11 -04:00
Peter Boyle
86948c6ea0
CRC for finger print fields - aids debug / version diff
2021-05-04 13:44:38 -04:00
Peter Boyle
53d226924a
CRC added
2021-05-04 13:44:07 -04:00
Christopher Kelly
80176b1b39
RHMC now outputs some initial norms to the logs
...
Fixed DWF+I Gparity binaries not correctly assigning twist directions (thanks Peter!)
2021-05-04 13:12:23 -04:00
u61464
15ae317858
Merge branch 'develop' of https://github.com/paboyle/Grid into develop
2021-05-04 08:40:38 -07:00
u61464
834f536b5f
Fastest option on SyCL is now std::complex
2021-05-04 08:40:18 -07:00
Christopher Kelly
29ddafd0fc
Added variant of G-parity DWF+I ensemble gen code using double prec RHMC
2021-04-30 13:12:24 -04:00
Peter Boyle
c332d9f08b
Merge pull request #356 from felixerben/bugfix/stoutSmearing
...
Jamie's fix
2021-04-27 14:10:49 -04:00
cf2923d5dd
Jamie's fix
2021-04-27 16:53:37 +01:00
Peter Boyle
0e4413ddde
Merge pull request #355 from felixerben/bugfix/stoutSmearing
...
bugfix 3D stout smearing
2021-04-27 08:01:55 -04:00
Peter Boyle
0f08364e4f
Mom filter refresh sRNG
2021-04-26 23:18:11 +02:00
Peter Boyle
a198d59381
Merge branch 'feature/gparity_HMC' of https://github.com/paboyle/Grid into feature/gparity_HMC
2021-04-26 21:05:52 +02:00
009ccd581e
bugfix 3D stout smearing
2021-04-26 10:36:33 +01:00
Peter Boyle
8cd4263974
Tests compile
2021-04-25 22:20:37 -04:00
Peter Boyle
d45c868656
Change interface
2021-04-25 10:53:34 -04:00
Peter Boyle
955a8113de
Expose label only to reduce number of parameters
2021-04-25 10:36:38 -04:00
Peter Boyle
dbe210dd53
Open the ens_id
2021-04-25 10:25:59 -04:00
Peter Boyle
3a4f5f2324
Merge develop, strengthen force tests
2021-04-22 18:54:00 -04:00
Peter Boyle
824d84473f
Merge branch 'develop' into feature/gparity_HMC
2021-04-22 16:32:41 -04:00
Peter Boyle
38964a4076
Switch twist direction
2021-04-22 15:57:37 -04:00
Peter Boyle
0d9aa87228
Reduce momentum to the GP plane
2021-04-22 15:56:59 -04:00
Peter Boyle
0e959d9b94
Update plaquette analysis
2021-04-22 15:55:47 -04:00
Peter Boyle
752f70cd48
Merge branch 'develop' into feature/gparity_HMC
2021-04-22 01:58:11 +02:00
54c6b1376d
Quick fix of conserved current implementation in CayleyFermion5D. Now function treats current insertion with appropriate periodic boundary conditions in the mu=3 direction.
2021-04-21 16:56:46 +01:00
Peter Boyle
86e11743ca
set twists
2021-04-20 10:19:11 -04:00
f3f11b586f
Tadpole sign now in front of forward hopping term to be consistent with previous implementation and analytic form.
2021-04-17 12:44:27 +01:00
8083e3f7e8
Sign factor for tadpole implementation corrected.
2021-04-15 11:14:31 +01:00
Christopher Kelly
e0e42873c1
Const correctness for Lattice::Replicate
...
Adapted GeneralEvenOddRationalRatio and Test_rhmc_EOWilsonRatio_doubleVsMixedPrec to recent changes that require passing in serial RNG
For GeneralEvenOddRationalRatio and TwoFlavourEvenOddRatio, broke refresh into two stages, the first of which generates the random field and the second that computes the pseudofermion field.
This allows derived classes to override the generation of the random field, for example in testing.
Test_dwf_gpforce now uses Gparity in x-direction and APBC in time as opposed to G-parity in time
Added Test_action_dwf_gparity2fvs1f that compares the DWF fermion action with the 2f and the 1f (doubled-lattice) implementations of Gparity
2021-04-14 16:41:27 -04:00
Peter Boyle
980e721f6e
Update MetaData.h
2021-04-13 09:33:01 -04:00
