Peter Boyle
48db72259e
EvenOdd schur decomposed mpcdagmpc version of rhmc determinant.
...
dH is also small and plaquette looks right.
2015-08-18 18:37:39 +01:00
Peter Boyle
9c7840c3a7
rhmc for 1+1 wilson is conserving dH~0.
...
A good days work ;)
2015-08-18 16:58:56 +01:00
Peter Boyle
5c364f8082
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
2015-08-18 14:40:08 +01:00
Peter Boyle
2dd9ad7b0f
Update TODO list
2015-08-18 10:43:32 +01:00
Peter Boyle
bdcbfe9310
Even Odd two flavour ratio added and dH == small
2015-08-18 10:37:08 +01:00
Peter Boyle
76f3855629
Merge branch 'master' of https://github.com/paboyle/Grid
2015-08-18 09:23:58 +01:00
Peter Boyle
6212807a77
Small dh obtained in two flavour ratio so looks ok.
2015-08-18 09:21:29 +01:00
Peter Boyle
0bc38a69ce
Adding PV pseudofermion in prep for DWF HMC.
...
Not compiled this yet, but cloned in from BFM.
2015-08-18 09:19:42 +01:00
Peter Boyle
24382d77bb
Adding PV pseudofermion in prep for DWF HMC.
...
Not compiled this yet, but cloned in from BFM.
2015-08-17 23:14:48 +01:00
Peter Boyle
353d66def1
Unused apparently
2015-08-16 01:41:05 +01:00
Peter Boyle
afeabe0d23
Tidying
2015-08-16 00:14:10 +01:00
Peter Boyle
2d6b97be06
Merge branch 'master' of https://github.com/paboyle/Grid
2015-08-16 00:13:14 +01:00
Peter Boyle
53da927c3c
Merge branch 'master' of https://github.com/paboyle/Grid
2015-08-15 23:59:04 +01:00
Peter Boyle
c7b50d18e7
Merge branch 'master' of https://github.com/paboyle/Grid
2015-08-15 23:56:31 +01:00
Peter Boyle
155c164b0c
* 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.
2015-08-15 23:25:49 +01:00
Peter Boyle
f40475f382
Reorganising the Fermion interface
2015-08-14 14:16:45 +01:00
Peter Boyle
045c85823b
Extra test
2015-08-14 13:18:59 +01:00
Peter Boyle
cc63078de5
Gparity works now even if simd distributed in a Gparity twist direction.
...
Tested by doubling lattice in t-direction.
2015-08-14 12:57:42 +01:00
Peter Boyle
e6bed000c3
Gparity valence test now working.
...
Interface in FermionOperator will change a lot in future
2015-08-14 00:01:04 +01:00
Peter Boyle
fc9b36c769
Gamma5 mult direct
2015-08-13 10:51:29 +01:00
Peter Boyle
c39078162e
Gparity improvements
2015-08-13 10:51:01 +01:00
Peter Boyle
7e9203d8e0
Some bug fixes for more complicated types introduced with gparity
2015-08-13 10:50:34 +01:00
Peter Boyle
6ab73c5512
Gparity test added; partial implementation -- this is Chris K's doubled lattice only
...
and have to regress this with the 2 flavour implementation.
2015-08-12 09:49:33 +01:00
Peter Boyle
c8dca58e6d
File list update.
2015-08-11 06:37:42 +01:00
Peter Boyle
826fbb18c4
Preconditioned conjugate residual
2015-08-11 06:24:53 +01:00
Peter Boyle
07d672baeb
Header
2015-08-11 06:23:38 +01:00
Peter Boyle
3903dfe6a5
Gparity modifications in the Gparity compressor variant.
2015-08-11 06:22:20 +01:00
Peter Boyle
1b3c93e22a
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.
2015-08-10 20:47:44 +01:00
Peter Boyle
2be8df93ad
Adding components for even odd decomposed determinant in HMC.
...
dH not yet conserved, so something wrong in the eo force code still
2015-08-07 08:37:15 +01:00
Peter Boyle
b5a483ae60
Continued fraction overlap, partial fraction overlap force terms have a successful
...
test passing.
2015-08-01 22:48:21 +09:00
Peter Boyle
bb372a6a8a
Merge problem fixed
2015-08-01 22:30:00 +09:00
Peter Boyle
5e9bef8a1b
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
2015-08-01 22:24:54 +09:00
Peter Boyle
bb7d8535aa
Test for DWF force term passes
2015-08-01 22:07:05 +09:00
Peter Boyle
848104b1a9
Changes making force term test for DWF pass.
2015-08-01 22:06:07 +09:00
Peter Boyle
8627e237c8
Jackson smoothed chebyshev and (untested) completion of force terms
...
for Cayley, Partial and Cont fraction dwf and overlap.
have even odd and unprec forces.
2015-08-01 05:58:35 +09:00
Peter Boyle
68d9463be5
Bug in two flav pseudofermion corrected to reimport gauge field upon rejection.
...
exp(-DeltaH) = 1 now, and plaquette is sensible. Will reproduce an old Wilson Gauge
Wilson Fermion SCRI plaquette with precision in mass matching shortly.
2015-07-29 21:02:07 +09:00
Peter Boyle
0b603225d1
Two flavour HMC for Wilson/Wilson is conserving energy.
...
Still to check plaq and <e(-dH)>, but nevertheless this is
progress
2015-07-29 17:53:39 +09:00
Peter Boyle
f4c74e34d8
Committing incomplete work for parameter file I/O.
...
MacroMagic.h is central. Guido and I plan to move
over to generating virtual (XML, JSON, YAML, text, binary) encoding
from macro based system.
2015-07-27 18:32:28 +09:00
Peter Boyle
51031dd46c
Files renamed
2015-07-27 18:30:19 +09:00
Peter Boyle
63015324c1
Two flavour pseudofermion action
2015-07-26 12:28:03 +09:00
Peter Boyle
36b8f35eed
Elemental force term for Wilson dslash added and tests thereof passing.
...
Now need to construct pseudofermion two flavour, ratio, one flavour, ratio
action fragments.
2015-07-26 10:54:38 +09:00
Peter Boyle
a4953610bb
Merge branch 'master' of https://github.com/paboyle/Grid
2015-07-24 01:33:19 +09:00
Peter Boyle
5e370db6c5
Sizable improvement in multigrid for unsquared.
...
6000 matmuls CG unprec
2000 matmuls CG prec (4000 eo muls)
1050 matmuls PGCR on 16^3 x 32 x 8 m=.01
Substantial effort on timing and logging infrastructure
2015-07-24 01:31:13 +09:00
paboyle
e0ed320827
Bug work around
2015-07-21 22:49:36 -07:00
paboyle
ee37130b5e
Removed troublesome macros
2015-07-21 22:41:01 -07:00
Peter Boyle
5b475d5e08
5x speed up now
2015-07-22 00:30:05 +09:00
Peter Boyle
b382660425
INSTALL
2015-07-21 13:58:57 +09:00
Peter Boyle
987801c86d
Merge
2015-07-21 13:56:22 +09:00
Peter Boyle
f8be0aeed1
No changes shown on git diff
2015-07-21 13:54:09 +09:00
Peter Boyle
9651ab661b
Small pretty layout change
2015-07-21 13:53:23 +09:00
Peter Boyle
fb65953d82
This was needed to compile on gcc
2015-07-21 13:52:59 +09:00
Peter Boyle
98dfc9254f
This file is being developed and will remain hacky until the new algorithm
...
is complete
2015-07-21 13:52:23 +09:00
Peter Boyle
44cf212720
Printing change
2015-07-21 13:51:56 +09:00
Peter Boyle
cbec16fedd
More info
2015-07-21 13:48:57 +09:00
Peter Boyle
64703207d4
Tweaks to subspace set up to put in g5 r5 hermiticity
2015-07-21 12:13:03 +09:00
Peter Boyle
021478af3b
verbose
2015-07-21 12:12:29 +09:00
Peter Boyle
2da20f1443
This file drives me crazy
2015-07-11 23:06:31 +09:00
Peter Boyle
2c9ceaef94
No compile fix
2015-07-02 02:03:09 +01:00
Peter Boyle
55e313bd08
Cleaning up the recursion for traceIndex<n> after the changes the enable G++ to
...
compile it again.
2015-07-01 23:43:57 +01:00
Peter Boyle
ef0ec1d0b2
Merge branch 'master' of https://github.com/paboyle/Grid
2015-07-01 22:51:04 +01:00
paboyle
1e29d9f778
Some useful XC30 commands
2015-07-01 22:50:13 +01:00
paboyle
46cf661ecd
More xc30 config commansd
2015-07-01 22:48:58 +01:00
paboyle
cb9da7b371
Temporarily disable gmp dependency simply because Cray XC30's I'm benchmarking
...
have a downlevel gmp version that chokes on ::max_align_t where gmp had a
bug as far as I recall.
2015-07-01 22:47:33 +01:00
paboyle
71e6733d4d
Modified memory bw test to display word size
2015-07-01 22:46:53 +01:00
Peter Boyle
98b84d230a
Change the SIMD command correctly with precision = double vs. single and
...
connect the "Real" default precisoin to a configure flag.
Have RealF, RealD and Real types, where Real is compile target dependent single/double,
RealF is single and RealD is double etc..
2015-07-01 22:45:15 +01:00
paboyle
c4777879e6
Remove dependency on wrong file
2015-07-01 13:04:02 +01:00
Peter Boyle
e618c609fe
Merge branch 'master' of https://github.com/paboyle/Grid
2015-06-30 15:17:46 +01:00
Peter Boyle
e164ed6f12
Big commit fixing nocompiles in defective C++11 compilers (gcc, icpc). stared getting to
...
near the bleeding edge I guess
2015-06-30 15:17:27 +01:00
Peter Boyle
95ecf81d42
big commit fixing nocompiles in defective C++11 compilers (gcc, icpc). stared getting to
...
near the bleeding edge I guess
2015-06-30 15:03:11 +01:00
Peter Boyle
f36dbfffe5
VPGCR updates
2015-06-30 15:02:27 +01:00
Peter Boyle
8581c05ab3
big commit fixing nocompiles in defective C++11 compilers (gcc, icpc). stared getting to
...
near the bleeding edge I guess
2015-06-30 15:01:44 +01:00
Peter Boyle
59cd42c164
big commit fixing nocompiles in defective C++11 compilers (gcc, icpc). stared getting to
...
near the bleeding edge I guess
2015-06-30 15:01:26 +01:00
Peter Boyle
a4369e1db6
big commit fixing nocompiles in defective C++11 compilers (gcc, icpc). stared getting to
...
near the bleeding edge I guess
2015-06-30 15:00:19 +01:00
Peter Boyle
d6c79bbadb
Update Benchmark_comms.cc
2015-06-25 10:59:53 +01:00
Peter Boyle
5f8f0bc792
Some small steps towards a multigrid
2015-06-22 12:49:44 +01:00
Peter Boyle
eace9051e8
Merge
...
Merge branch 'master' of https://github.com/paboyle/Grid
2015-06-20 22:25:31 +01:00
Peter Boyle
f1916a7515
Will start this as a two level algorithm
2015-06-20 22:24:21 +01:00
Peter Boyle
5ccbac7db1
HDCG but this is not complete and placeholder for later completion
2015-06-20 22:23:57 +01:00
Peter Boyle
aba5c8595a
Patches for beginnings of an overlap multigrid
2015-06-20 22:22:56 +01:00
Peter Boyle
2534199ee7
Merge branch 'master' of https://github.com/paboyle/Grid
2015-06-19 17:24:05 +01:00
Peter Boyle
8ebff0b7c0
Merge branch 'master' of https://github.com/paboyle/Grid
2015-06-14 01:29:41 +01:00
Peter Boyle
0f6ed6b633
Merge branch 'master' of https://github.com/paboyle/Grid
2015-06-14 01:27:07 +01:00
Peter Boyle
f0d190024a
Updates to ldop tests
2015-06-10 12:26:25 +01:00
Peter Boyle
0ad72b38ab
file list
2015-06-10 12:25:15 +01:00
Peter Boyle
e79a57b423
Merge branch 'master' of https://github.com/paboyle/Grid
...
Not sure what changed in master
2015-06-09 22:51:10 +01:00
Peter Boyle
5a53086e41
Solver converges
2015-06-09 22:51:02 +01:00
Peter Boyle
c4f204f440
Solver converges
2015-06-09 22:50:45 +01:00
Peter Boyle
2022a7205c
Remove extra layers of checks now it works
2015-06-09 22:43:41 +01:00
Peter Boyle
4357f5acdc
Some extra tests
2015-06-09 22:43:10 +01:00
Peter Boyle
4963f7356a
5d OpDir direction interface refers to the 5d dims, not 4d to present a
...
sensible and consistent external interface.
2015-06-09 22:41:59 +01:00
Peter Boyle
3a52bc4ce1
G5R5 update
2015-06-09 22:41:27 +01:00
Peter Boyle
708d4f7533
g5 and g5R5 hermitian are now differentiated
2015-06-09 22:40:58 +01:00
Peter Boyle
61770d4472
Files update
2015-06-09 22:40:12 +01:00
Peter Boyle
a57ca0bbe7
Prettier reporting
2015-06-09 22:39:37 +01:00
Peter Boyle
1f2498918d
Got this sorted with the promote working in a test
2015-06-09 22:39:13 +01:00
Peter Boyle
25bfff370a
Starting to use
2015-06-09 22:38:13 +01:00
Peter Boyle
6403e3021f
Debugged finally. A silly mistake in permute cost me a day of debug.
2015-06-09 22:37:21 +01:00
Peter Boyle
4d11fc0330
silly change
2015-06-09 22:36:48 +01:00
Peter Boyle
64e2f36d40
Merge pull request #14 from mspraggs/master
...
Removed std::string calls from NerscIO map indexing
2015-06-09 10:31:52 +01:00
Peter Boyle
e8b43944e7
Merge branch 'master' of https://github.com/paboyle/Grid
...
Conflicts:
lib/Make.inc
2015-06-09 10:27:10 +01:00