Peter Boyle
cc4ca48d13
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
bc09d7c3bd
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
9de40578d3
Two flavour pseudofermion action
2015-07-26 12:28:03 +09:00
Peter Boyle
d7e6b65a76
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
28bdc90908
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
Peter Boyle
8925845684
Merge
2015-07-21 13:56:22 +09:00
Peter Boyle
c7925e5c9b
Small pretty layout change
2015-07-21 13:53:23 +09:00
Peter Boyle
9d18773fbc
This was needed to compile on gcc
2015-07-21 13:52:59 +09:00
neo
32e6887d5f
Rearranging files in hmc
2015-07-06 16:46:43 +09:00
neo
1991852025
Added minimum norm integrator
...
Little rearrangement of HMC and integrator classes
2015-07-06 16:17:32 +09:00
neo
2718038977
HMC for Wilson Gauge action works
...
Fixed bug in momenta generation
2015-07-06 12:58:49 +09:00
neo
62d8952c0a
Debugged vector version of ProjectOnGroup
2015-07-06 02:24:58 +09:00
neo
b1f94fa292
HMC ready but untested
2015-07-04 17:47:50 +09:00
neo
30c9dc473d
More progress in the HMC construction
2015-07-04 02:43:14 +09:00
neo
9655d43017
Skeleton of HMC/Integrators
2015-07-03 16:51:41 +09:00
Peter Boyle
f41c7dffef
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
74e397b29c
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
7de5ccb879
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
dec68e5c0e
Some small steps towards a multigrid
2015-06-22 12:49:44 +01:00
Azusa Yamaguchi
945bb93e48
Variable preconditioned GCR with restarting.
...
Orthogonalisation depth and restart frequency is controllable via constructor
2015-06-21 10:58:46 +01:00
Peter Boyle
bcf1d5160f
Merge
...
Merge branch 'master' of https://github.com/paboyle/Grid
2015-06-20 22:25:31 +01:00
Peter Boyle
6ad96f7383
Patches for beginnings of an overlap multigrid
2015-06-20 22:22:56 +01:00
Azusa Yamaguchi
6cebd006d4
Merge branch 'master' of https://github.com/paboyle/Grid
2015-06-20 14:22:29 +01:00
Azusa Yamaguchi
79a9f8b9c9
Quenched works for wilson gauge
2015-06-16 14:17:11 +01:00
Azusa Yamaguchi
173b31ce05
uninitialised bug fix
2015-06-16 14:07:05 +01:00
Azusa Yamaguchi
ae0873bc77
First cut at SUN support for quenched updates
2015-06-14 01:28:54 +01:00
Azusa Yamaguchi
87cc0d4ca3
Peek poke colour/spin/complex and trace transpose support
2015-06-14 01:00:11 +01:00
Azusa Yamaguchi
66e5718610
const safety
2015-06-14 00:59:50 +01:00
Peter Boyle
6abbd35d81
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
c7152c520a
g5 and g5R5 hermitian are now differentiated
2015-06-09 22:40:58 +01:00
Peter Boyle
d8ddec86f7
Merge branch 'master' of https://github.com/paboyle/Grid
...
Conflicts:
lib/Make.inc
2015-06-09 10:27:10 +01:00
Peter Boyle
506dfd1517
Some unary ops and coarse grid support
2015-06-09 10:26:19 +01:00
neo
6b8fe04054
Experimental support for ARM
2015-06-09 15:46:21 +09:00
Peter Boyle
9e7035f5dc
Conjugate residual algorithm; some more unary functions
2015-06-08 12:04:59 +01:00
Peter Boyle
50e8b2160e
Conjugate residual added
2015-06-05 18:16:25 +01:00
Azusa Yamaguchi
ad18df92d0
Compile fix
2015-06-05 10:29:42 +01:00
Azusa Yamaguchi
1d7f9567ee
Endif terminated
2015-06-05 10:19:42 +01:00
Azusa Yamaguchi
a8b86e747b
merge to the head
2015-06-05 10:15:31 +01:00
Azusa Yamaguchi
c05fe2706c
Merge branch 'master' of https://github.com/paboyle/Grid
2015-06-05 10:04:46 +01:00
Azusa Yamaguchi
58cdcbb5e4
Adding some wilson loop support
2015-06-05 10:02:36 +01:00
Peter Boyle
b9e9777912
PartialFraction Hw with Zolo and Tanh approx converged under CG and passed EO breakdown
...
and hermiticity tests.
2015-06-04 13:28:37 +01:00
neo
b9edadc53e
Addedd Ta functionality to the tensor types
...
Merge remote-tracking branch 'upstream/master'
Conflicts:
configure
2015-06-04 18:11:32 +09:00
Peter Boyle
37aa74dfd2
CG Tests work for wilson kernel cont frac zolo and tanh
2015-06-04 06:02:00 +01:00
Peter Boyle
c327019574
Implementing the Hw kernel continued fraction 5d overlap cases
2015-06-04 00:23:16 +01:00
Peter Boyle
50bd293527
First pass at continued fraction; solver and even odd decomposition tests pass.
...
Have to make ContFrac class virtual and derive end non-abstract actions for the particular
cases.
2015-06-04 00:00:45 +01:00
Peter Boyle
4bcc319e11
Reorganise of file naming
2015-06-03 12:47:05 +01:00
Peter Boyle
8fe3d4f971
Overlap Wilson Cayley tanh & zolo
2015-06-03 11:26:54 +01:00
Peter Boyle
343d039b37
Scaled Shamir and Scaled Shamir Zolotarev aliases for special cases of Mobius.
2015-06-03 09:51:06 +01:00
Peter Boyle
5916386242
Mobius Caley form, Mobius Zolotarev operators. Pass Even Odd vs unprec test and hermiticity checks
...
in tests/Grid_any_evenodd.cc; will work on inversion tests shortly.
2015-06-03 09:36:26 +01:00
Peter Boyle
2583570e17
Domain wall fermions now invert ; have the basis set up for
...
Tanh/Zolo * (Cayley/PartFrac/ContFrac) * (Mobius/Shamir/Wilson)
Approx Representation Kernel.
All are done with space-time taking part in checkerboarding, Ls uncheckerboarded
Have only so far tested the Domain Wall limit of mobius, and at that only checked
that it
i) Inverts
ii) 5dim DW == Ls copies of 4dim D2
iii) MeeInv Mee == 1
iv) Meo+Mee+Moe+Moo == M unprec.
v) MpcDagMpc is hermitan
vi) Mdag is the adjoint of M between stochastic vectors.
That said, the RB schur solve, RB MpcDagMpc solve, Unprec solve
all converge and the true residual becomes small; so pretty good tests.
2015-06-02 16:57:12 +01:00
