Skip to content
Snippets Groups Projects
Commit cfe714f1 authored by Pedro Gonnet's avatar Pedro Gonnet
Browse files

better initial guesses for h.

Former-commit-id: 1f70df3bec09555f3f0a70072b1d6c829f0f8d0b
parent 4e43037d
Branches
Tags
No related merge requests found
############################################################################### ###############################################################################
# This file is part of SWIFT. # This file is part of SWIFT.
# Coypright (c) 2012 Pedro Gonnet (pedro.gonnet@durham.ac.uk), # Coypright (c) 2013 Pedro Gonnet (pedro.gonnet@durham.ac.uk),
# Matthieu Schaller (matthieu.schaller@durham.ac.uk) # Matthieu Schaller (matthieu.schaller@durham.ac.uk)
# #
# This program is free software: you can redistribute it and/or modify # This program is free software: you can redistribute it and/or modify
...@@ -63,7 +63,7 @@ for i in range(L): ...@@ -63,7 +63,7 @@ for i in range(L):
v[index,1] = 0. v[index,1] = 0.
v[index,2] = 0. v[index,2] = 0.
m[index] = mass m[index] = mass
h[index] = 1.2 * boxSize / L h[index] = 2.251 * boxSize / L
u[index] = internalEnergy u[index] = internalEnergy
ids[index] = index ids[index] = index
......
...@@ -67,7 +67,7 @@ for i in range(L): ...@@ -67,7 +67,7 @@ for i in range(L):
v[index,1] = 0. v[index,1] = 0.
v[index,2] = 0. v[index,2] = 0.
m[index] = mass m[index] = mass
h[index] = 1.2 * boxSize / L h[index] = 2.251 * boxSize / L
u[index] = internalEnergy u[index] = internalEnergy
ids[index] = index ids[index] = index
if sqrt((x - boxSize/2.)**2 + (y - boxSize/2.)**2 + (z - boxSize/2.)**2) < 2.01 * boxSize/L: if sqrt((x - boxSize/2.)**2 + (y - boxSize/2.)**2 + (z - boxSize/2.)**2) < 2.01 * boxSize/L:
......
############################################################################### ###############################################################################
# This file is part of SWIFT. # This file is part of SWIFT.
# Coypright (c) 2012 Pedro Gonnet (pedro.gonnet@durham.ac.uk), # Coypright (c) 2013 Pedro Gonnet (pedro.gonnet@durham.ac.uk),
# Matthieu Schaller (matthieu.schaller@durham.ac.uk) # Matthieu Schaller (matthieu.schaller@durham.ac.uk)
# #
# This program is free software: you can redistribute it and/or modify # This program is free software: you can redistribute it and/or modify
...@@ -61,7 +61,7 @@ for i in range(L): ...@@ -61,7 +61,7 @@ for i in range(L):
v[index,1] = 0. v[index,1] = 0.
v[index,2] = 0. v[index,2] = 0.
m[index] = mass m[index] = mass
h[index] = 1.2 * boxSize / L h[index] = 2.251 * boxSize / L
u[index] = internalEnergy u[index] = internalEnergy
ids[index] = index ids[index] = index
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment