From 98e42fb6a81c4642f10cfc6af962a4644fac39a5 Mon Sep 17 00:00:00 2001
From: "Pyry Runko (Puck)" <pyjorunk@student.jyu.fi>
Date: Mon, 17 Feb 2025 09:14:27 +0200
Subject: [PATCH] removed some lines, still some to remove
---
src/ipglasma.cpp | 14 ++++----------
1 file changed, 4 insertions(+), 10 deletions(-)
diff --git a/src/ipglasma.cpp b/src/ipglasma.cpp
index a683f5d..4da94cc 100644
--- a/src/ipglasma.cpp
+++ b/src/ipglasma.cpp
@@ -630,22 +630,16 @@ double IPGlasma::SumAmplitudes()
WilsonLine quark;
int thrownAway = 0;
int iWasHere = 0;
- #pragma omp parallel for reduction(+:sum,thrownAway,iWasHere) collapse(4)
+ #pragma omp parallel for reduction(+:sum)
for(int i =0; i < L; i += 1){
for(int j =0; j < L; j += 1){
-
+ int q1[2] = {i,j};
+ quark = GetWilsonLine(WilsonLineCoordinate(i, j))
+ double r1 = DistanceToOrigin(q1);
//#pragma omp parallel for reduction(+:sum) collapse(2)
for(int k =0; k < L; k +=1){
-// cout << "# Thread " << omp_get_thread_num() << " is processing (" << i << ", " << j << ", " << k <<")" << endl;
for(int l =0; l < L; l +=1){
- cout << "# Hello from inside the l loop!" << endl;
- //TODO: pointer instead of copy?
- //TODO: faster to get these two layers above?
- iWasHere += 1;
- int q1[2] = {i,j};
- //quark = GetWilsonLine(WilsonLineCoordinate(i, j));
int q2[2] = {k,l};
- double r1 = DistanceToOrigin(q1);
double r2 = DistanceToOrigin(q2);
if (r1 > 1. && r2 > 1.){
// cout << "#!!Throwing away large dipole separation!" << endl;
--
GitLab