Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Remove tc thread private array #435

Merged
merged 9 commits into from
Oct 30, 2023
Merged

Remove tc thread private array #435

Show file tree
Hide file tree
Changes from 1 commit
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
5e261ae
Remove tc thread private array
marchdf Oct 26, 2023
6b104ba
Update mechs
marchdf Oct 26, 2023
b684ecf
Const some things
marchdf Oct 26, 2023
69bcaf3
update mechs
marchdf Oct 26, 2023
d782775
python updates for qss
marchdf Oct 26, 2023
321835d
update qss mechs
marchdf Oct 26, 2023
28190aa
fixes for qss mechs in python
marchdf Oct 26, 2023
18b1d9e
another python fix
marchdf Oct 26, 2023
7f8962e
update mechs
marchdf Oct 26, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
Next Next commit
update mechs
  • Loading branch information
@marchdf
marchdf committed Oct 27, 2023
commit 69bcaf32eb44fab156da04f55f6c6d87014063fe
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/Aromatic_KrNara/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45317,9 +45317,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -96938,9 +96938,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
4 changes: 2 additions & 2 deletions Support/Mechanism/Models/Aromatic_KrNara/mechanism.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -912,8 +912,8 @@ void
progressRateFR(
amrex::Real* q_f, amrex::Real* q_r, amrex::Real* sc, amrex::Real T)
{
amrex::Real invT = 1.0 / T;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real logT = log(T);
// compute the Gibbs free energy
amrex::Real g_RT[158];
gibbs(g_RT, T);
Expand Down
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/BurkeDryer/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3416,9 +3416,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -4803,9 +4803,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
4 changes: 2 additions & 2 deletions Support/Mechanism/Models/BurkeDryer/mechanism.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -82,8 +82,8 @@ void
progressRateFR(
amrex::Real* q_f, amrex::Real* q_r, amrex::Real* sc, amrex::Real T)
{
amrex::Real invT = 1.0 / T;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real logT = log(T);
// compute the Gibbs free energy
amrex::Real g_RT[13];
gibbs(g_RT, T);
Expand Down
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/C1-C2-NO/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -13509,9 +13509,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -28587,9 +28587,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
4 changes: 2 additions & 2 deletions Support/Mechanism/Models/C1-C2-NO/mechanism.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -234,8 +234,8 @@ void
progressRateFR(
amrex::Real* q_f, amrex::Real* q_r, amrex::Real* sc, amrex::Real T)
{
amrex::Real invT = 1.0 / T;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real logT = log(T);
// compute the Gibbs free energy
amrex::Real g_RT[44];
gibbs(g_RT, T);
Expand Down
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/CH4_lean/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5076,9 +5076,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -8840,9 +8840,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
4 changes: 2 additions & 2 deletions Support/Mechanism/Models/CH4_lean/mechanism.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -114,8 +114,8 @@ void
progressRateFR(
amrex::Real* q_f, amrex::Real* q_r, amrex::Real* sc, amrex::Real T)
{
amrex::Real invT = 1.0 / T;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real logT = log(T);
// compute the Gibbs free energy
amrex::Real g_RT[17];
gibbs(g_RT, T);
Expand Down
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/Davis/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3818,9 +3818,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -5790,9 +5790,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
4 changes: 2 additions & 2 deletions Support/Mechanism/Models/Davis/mechanism.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -88,8 +88,8 @@ void
progressRateFR(
amrex::Real* q_f, amrex::Real* q_r, amrex::Real* sc, amrex::Real T)
{
amrex::Real invT = 1.0 / T;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real logT = log(T);
// compute the Gibbs free energy
amrex::Real g_RT[14];
gibbs(g_RT, T);
Expand Down
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/FFCM1_Red/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6858,9 +6858,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -13374,9 +13374,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
4 changes: 2 additions & 2 deletions Support/Mechanism/Models/FFCM1_Red/mechanism.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -145,8 +145,8 @@ void
progressRateFR(
amrex::Real* q_f, amrex::Real* q_r, amrex::Real* sc, amrex::Real T)
{
amrex::Real invT = 1.0 / T;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real logT = log(T);
// compute the Gibbs free energy
amrex::Real g_RT[21];
gibbs(g_RT, T);
Expand Down
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/JL4/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2181,9 +2181,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -2457,9 +2457,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
4 changes: 2 additions & 2 deletions Support/Mechanism/Models/JL4/mechanism.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -69,8 +69,8 @@ void
progressRateFR(
amrex::Real* q_f, amrex::Real* q_r, amrex::Real* sc, amrex::Real T)
{
amrex::Real invT = 1.0 / T;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real logT = log(T);
// compute the Gibbs free energy
amrex::Real g_RT[7];
gibbs(g_RT, T);
Expand Down
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/Kolla/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -3384,9 +3384,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -4941,9 +4941,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
4 changes: 2 additions & 2 deletions Support/Mechanism/Models/Kolla/mechanism.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -83,8 +83,8 @@ void
progressRateFR(
amrex::Real* q_f, amrex::Real* q_r, amrex::Real* sc, amrex::Real T)
{
amrex::Real invT = 1.0 / T;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real logT = log(T);
// compute the Gibbs free energy
amrex::Real g_RT[12];
gibbs(g_RT, T);
Expand Down
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/LiDryer/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -2836,9 +2836,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -3961,9 +3961,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
4 changes: 2 additions & 2 deletions Support/Mechanism/Models/LiDryer/mechanism.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -80,8 +80,8 @@ void
progressRateFR(
amrex::Real* q_f, amrex::Real* q_r, amrex::Real* sc, amrex::Real T)
{
amrex::Real invT = 1.0 / T;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real logT = log(T);
// compute the Gibbs free energy
amrex::Real g_RT[9];
gibbs(g_RT, T);
Expand Down
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/LuDME/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -11167,9 +11167,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -20528,9 +20528,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
4 changes: 2 additions & 2 deletions Support/Mechanism/Models/LuDME/mechanism.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -176,8 +176,8 @@ void
progressRateFR(
amrex::Real* q_f, amrex::Real* q_r, amrex::Real* sc, amrex::Real T)
{
amrex::Real invT = 1.0 / T;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real logT = log(T);
// compute the Gibbs free energy
amrex::Real g_RT[39];
gibbs(g_RT, T);
Expand Down
12 changes: 6 additions & 6 deletions Support/Mechanism/Models/LuEthylene/mechanism.H
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -10338,9 +10338,9 @@ aJacobian_precond(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down Expand Up @@ -21681,9 +21681,9 @@ aJacobian(
val = 0.0;
}

amrex::Real invT = 1.0 / T;
amrex::Real invT2 = invT * invT;
amrex::Real logT = log(T);
const amrex::Real invT = 1.0 / T;
const amrex::Real invT2 = invT * invT;
const amrex::Real logT = log(T);

// reference concentration: P_atm / (RT) in inverse mol/m^3
amrex::Real refC = 101325 / 8.31446 / T;
Expand Down
Loading