HMC algorithm with multiple time scale integration and mass preconditioning

We describe a new HMC algorithm variant we have recently introduced and extend the published results by preliminary results of a simulation with a pseudo scalar mass value of about 300 MeV. This new run confirms our expectation that simulations with such pseudo scalar mass values become feasible and affordable with our HMC variant. In addition we discuss simulations from hot and cold starts at a pseudo scalar mass value of about 300 MeV, which we performed in order to test for possible meta-stabilities.Comment: 6 pages, Talk presented at Lattice 2005 (machines and algorithms

Experiences with OpenMP in tmLQCD

An overview is given of the lessons learned from the introduction of multi-threading using OpenMP in tmLQCD. In particular, programming style, performance measurements, cache misses, scaling, thread distribution for hybrid codes, race conditions, the overlapping of communication and computation and the measurement and reduction of certain overheads are discussed. Performance measurements and sampling profiles are given for different implementations of the hopping matrix computational kernel.Comment: presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, German

On the efficient numerical solution of lattice systems with low-order couplings

We apply the Quasi Monte Carlo (QMC) and recursive numerical integration methods to evaluate the Euclidean, discretized time path-integral for the quantum mechanical anharmonic oscillator and a topological quantum mechanical rotor model. For the anharmonic oscillator both methods outperform standard Markov Chain Monte Carlo methods and show a significantly improved error scaling. For the quantum mechanical rotor we could, however, not find a successful way employing QMC. On the other hand, the recursive numerical integration method works extremely well for this model and shows an at least exponentially fast error scaling

Outflow boundary conditions for 3D simulations of non-periodic blood flow and pressure fields in deformable arteries

The simulation of blood flow and pressure in arteries requires outflow boundary conditions that incorporate models of downstream domains. We previously described a coupled multidomain method to couple analytical models of the downstream domains with 3D numerical models of the upstream vasculature. This prior work either included pure resistance boundary conditions or impedance boundary conditions based on assumed periodicity of the solution. However, flow and pressure in arteries are not necessarily periodic in time due to heart rate variability, respiration, complex transitional flow or acute physiological changes. We present herein an approach for prescribing lumped parameter outflow boundary conditions that accommodate transient phenomena. We have applied this method to compute haemodynamic quantities in different physiologically relevant cardiovascular models, including patient-specific examples, to study non-periodic flow phenomena often observed in normal subjects and in patients with acquired or congenital cardiovascular disease. The relevance of using boundary conditions that accommodate transient phenomena compared with boundary conditions that assume periodicity of the solution is discussed

Non-perturbative renormalization of moments of parton distribution functions

We compute non-perturbatively the evolution of the twist-2 operators corresponding to the average momentum of non-singlet quark densities. The calculation is based on a finite-size technique, using the Schr\"odinger Functional, in quenched QCD. We find that a careful choice of the boundary conditions, is essential, for such operators, to render possible the computation. As a by-product we apply the non-perturbatively computed renormalization constants to available data of bare matrix elements between nucleon states.Comment: Lattice2003(Matrix); 3 pages, 3 figures. Talk by A.

Complexity and Inapproximability Results for Parallel Task Scheduling and Strip Packing

We study the Parallel Task Scheduling problem $Pm|size_j|C_{\max}$ with a constant number of machines. This problem is known to be strongly NP-complete for each $m \geq 5$, while it is solvable in pseudo-polynomial time for each $m \leq 3$. We give a positive answer to the long-standing open question whether this problem is strongly $NP$-complete for $m=4$. As a second result, we improve the lower bound of $\frac{12}{11}$ for approximating pseudo-polynomial Strip Packing to $\frac{5}{4}$. Since the best known approximation algorithm for this problem has a ratio of $\frac{4}{3} + \varepsilon$, this result narrows the gap between approximation ratio and inapproximability result by a significant step. Both results are proven by a reduction from the strongly $NP$-complete problem 3-Partition

Scaling test of quenched Wilson twisted mass QCD at maximal twist

We present the results of an extended scaling test of quenched Wilson twisted mass QCD. We fix the twist angle by using two definitions of the critical mass, the first obtained by requiring the vanishing of the pseudoscalar meson mass m_PS for standard Wilson fermions and the second by requiring restoration of parity at non-zero value of the twisted mass mu and subsequently extrapolating to mu=0. Depending on the choice of the critical mass we simulate at values of beta in [5.7,6.45], for a range of pseudoscalar meson masses 250 MeV < m_PS < 1 GeV and we perform the continuum limit for the pseudoscalar meson decay constant f_PS and various hadron masses (vector meson m_V, baryon octet m_oct and baryon decuplet m_dec) at fixed value of r_0 m_PS. For both definitions of the critical mass, lattice artifacts are consistent with O(a) improvement. However, with the second definition, large O(a^2) discretization errors present at small quark mass with the first definition are strongly suppressed. The results in the continuum limit are in very good agreement with those from the Alpha and CP-PACS Collaborations.Comment: 6 pages, Talk presented at Lattice 2005, Dublin, 25-30 July 200

Hopping conductivity in heavily doped n-type GaAs layers in the quantum Hall effect regime

We investigate the magnetoresistance of epitaxially grown, heavily doped n-type GaAs layers with thickness (40-50 nm) larger than the electronic mean free path (23 nm). The temperature dependence of the dissipative resistance R_{xx} in the quantum Hall effect regime can be well described by a hopping law (R_{xx} \propto exp{-(T_0/T)^p}) with p=0.6. We discuss this result in terms of variable range hopping in a Coulomb gap together with a dependence of the electron localization length on the energy in the gap. The value of the exponent p>0.5 shows that electron-electron interactions have to be taken into account in order to explain the occurrence of the quantum Hall effect in these samples, which have a three-dimensional single electron density of states.Comment: 5 pages, 2 figures, 1 tabl