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

Jump to bottom

Verlet Integration: Output standardization #855

Merged
merged 18 commits into from
Sep 4, 2021
Merged

Verlet Integration: Output standardization #855

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
18 commits
Select commit Hold shift + click to select a range
944c8d6
Standardize julia output
stormofice Sep 3, 2021
9017c00
Standardize kotlin output
stormofice Sep 3, 2021
73b472e
Standardized c output
stormofice Sep 3, 2021
3dfbf54
Standardized cpp output
stormofice Sep 3, 2021
e9d1890
Standardized lisp output
stormofice Sep 3, 2021
b1e821d
Standardized fortran output
stormofice Sep 3, 2021
5c57023
Standardized go output
stormofice Sep 3, 2021
70deb2d
Standardized java output
stormofice Sep 3, 2021
357d962
Standardize javascript output
stormofice Sep 3, 2021
70068de
Standardize nim output
stormofice Sep 3, 2021
ba11044
Standardize python output
stormofice Sep 3, 2021
886adcd
Standardize ruby output
stormofice Sep 3, 2021
fbe95d1
Standardize rust output
stormofice Sep 3, 2021
269708c
Standardize swift output
stormofice Sep 3, 2021
364dd35
Standardized haskell output
jiegillet Sep 3, 2021
03c1af5
Standardized haskell output (no quote marks)
jiegillet Sep 3, 2021
4be01a9
attempt at fix for asm
leios Sep 4, 2021
a8aadbe
Merge pull request #1 from leios/stormofice-verlet-io-overhaul
stormofice Sep 4, 2021
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
6 changes: 3 additions & 3 deletions contents/verlet_integration/code/asm-x64/verlet.s
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4,9 +4,9 @@
zero: .double 0.0
two: .double 2.0
half: .double 0.5
verlet_fmt: .string "Time for Verlet integration is: %lf\n"
stormer_fmt: .string "Time and Velocity for Stormer Verlet Integration is: %lf, %lf\n"
velocity_fmt: .string "Time and Velocity for Velocity Verlet Integration is: %lf, %lf\n"
verlet_fmt: .string "[#] Time for Verlet integration is:\n%lf\n"
stormer_fmt: .string "[#] Time for Stormer Verlet Integration is:\n%lf\n[#] Velocity for Stormer Verlet Integration is:\n%lf\n"
velocity_fmt: .string "[#] Time for Velocity Verlet Integration is:\n%lf\n[#] Velocity for Velocity Verlet Integration is:\n%lf\n"
pos: .double 5.0
acc: .double -10.0
dt: .double 0.01
Expand Down
10 changes: 5 additions & 5 deletions contents/verlet_integration/code/c++/verlet.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -64,19 +64,19 @@ int main() {
// each of these functions.

double time = verlet(5.0, -10, 0.01);
std::cout << "Time for Verlet integration is: " \
std::cout << "[#] Time for Verlet integration is:\n" \
<< time << std::endl;

timestep timestep_sv = stormer_verlet(5.0, -10, 0.01);
std::cout << "Time for Stormer Verlet integration is: " \
std::cout << "[#] Time for Stormer Verlet integration is:\n" \
<< timestep_sv.time << std::endl;
std::cout << "Velocity for Stormer Verlet integration is: " \
std::cout << "[#] Velocity for Stormer Verlet integration is:\n" \
<< timestep_sv.vel << std::endl;

timestep timestep_vv = velocity_verlet(5.0, -10, 0.01);
std::cout << "Time for velocity Verlet integration is: " \
std::cout << "[#] Time for velocity Verlet integration is:\n" \
<< timestep_vv.time << std::endl;
std::cout << "Velocity for velocity Verlet integration is: " \
std::cout << "[#] Velocity for velocity Verlet integration is:\n" \
<< timestep_vv.vel << std::endl;

return 0;
Expand Down
16 changes: 10 additions & 6 deletions contents/verlet_integration/code/c/verlet.c
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -46,16 +46,20 @@ int main() {
double time, vel;

verlet(&time, 5.0, -10, 0.01);
printf("Time for Verlet integration is: %lf\n",
time);
printf("[#] Time for Verlet integration is:\n");
printf("%lf\n", time);

stormer_verlet(&time, &vel, 5.0, -10, 0.01);
printf("Time and velocity for Stormer Verlet integration is: %lf, %lf\n",
time, vel);
printf("[#] Time for Stormer Verlet integration is:\n");
printf("%lf\n", time);
printf("[#] Velocity for Stormer Verlet integration is:\n");
printf("%lf\n", vel);

velocity_verlet(&time, &vel, 5.0, -10, 0.01);
printf("Time and velocity for velocity Verlet integration is: %lf, %lf\n",
time, vel);
printf("[#] Time for velocity Verlet integration is:\n");
printf("%lf\n", time);
printf("[#] Velocity for Stormer Verlet integration is:\n");
printf("%lf\n", vel);

return 0;
}
15 changes: 10 additions & 5 deletions contents/verlet_integration/code/clisp/verlet.lisp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -34,12 +34,17 @@
while (> p 0)
finally (return (list time vel))))

(format T "Time for Verlet integration: ~d~%" (verlet 5 -10 0.01))
(format T "[#] Time for Verlet integration:~%")
(format T "~d~%" (verlet 5 -10 0.01))

(defvar stormer-verlet-result (stormer-verlet 5 -10 0.01))
(format T "Time for Stormer Verlet integration is: ~d~%" (first stormer-verlet-result))
(format T "Velocity for Stormer Verlet integration is: ~d~%" (second stormer-verlet-result))
(format T "[#] Time for Stormer Verlet integration is:~%")
(format T "~d~%" (first stormer-verlet-result))
(format T "[#] Velocity for Stormer Verlet integration is:~%")
(format T "~d~%" (second stormer-verlet-result))

(defvar velocity-verlet-result (velocity-verlet 5 -10 0.01))
(format T "Time for velocity Verlet integration is: ~d~%" (first velocity-verlet-result))
(format T "Velocity for velocity Verlet integration is: ~d~%" (second velocity-verlet-result))
(format T "[#] Time for velocity Verlet integration is:~%")
(format T "~d~%" (first velocity-verlet-result))
(format T "[#] Velocity for velocity Verlet integration is:~%")
(format T "~d~%" (second velocity-verlet-result))
17 changes: 13 additions & 4 deletions contents/verlet_integration/code/fortran/verlet.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -91,18 +91,27 @@ SUBROUTINE velocity_verlet(pos, acc, dt, time, vel)
! Verlet
CALL verlet(pos, acc, dt, time)

WRITE(*,*) 'Time for Verlet integration: ', time
WRITE(*,*) '[#] Time for Verlet integration:'
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

At least with gfortran, this seems to add extra spaces at the start of the line, is this ok?

output looks like this:

 [#] Time for Verlet integration:
   1.0000000000000007     
 [#] Time for Stormer Verlet integration:
   1.0000000000000007     
 [#] Velocity for Stormer Verlet integration:
  -9.9999999999999805     
 [#] Time for velocity Verlet integration:
   1.0100000000000007     
 [#] Velocity for velocity Verlet integration:
  -10.099999999999980

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I was told that these spaces can be ignored

WRITE(*,*) time

! stormer Verlet
pos = 5d0
CALL stormer_verlet(pos, acc, dt, time, vel)

WRITE(*,*) 'Time for Stormer-Verlet integration: ', time
WRITE(*,*) '[#] Time for Stormer Verlet integration:'
WRITE(*,*) time
WRITE(*,*) '[#] Velocity for Stormer Verlet integration:'
WRITE(*,*) vel



! Velocity Verlet
pos = 5d0
CALL velocity_verlet(pos, acc, dt, time, vel)

WRITE(*,*) 'Time for Velocity-Verlet integration: ', time

WRITE(*,*) '[#] Time for velocity Verlet integration:'
WRITE(*,*) time
WRITE(*,*) '[#] Velocity for velocity Verlet integration:'
WRITE(*,*) vel

END PROGRAM verlet_integration
22 changes: 11 additions & 11 deletions contents/verlet_integration/code/golang/verlet.go
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -43,18 +43,18 @@ func velocityVerlet(pos, acc, dt float64) (time, vel float64) {

func main() {
time := verlet(5., -10., .01)
fmt.Println("Verlet")
fmt.Println("Time:", time)
fmt.Println()
fmt.Println("[#] Time for Verlet integration is:")
fmt.Println(time)

time, vel := stormerVerlet(5., -10., .01)
fmt.Println("Stormer-Verlet")
fmt.Println("Time:", time)
fmt.Println("Velocity:", vel)
fmt.Println()

fmt.Println("[#] Time for Stormer Verlet integration is:")
fmt.Println(time)
fmt.Println("[#] Velocity for Stormer Verlet integration is:")
fmt.Println(vel)
time, vel = velocityVerlet(5., -10., .01)
fmt.Println("Velocity Verlet")
fmt.Println("Time:", time)
fmt.Println("Velocity:", vel)
fmt.Println("[#] Time for velocity Verlet integration is:")
fmt.Println(time)
fmt.Println("[#] Velocity for velocity Verlet integration is:")
fmt.Println(vel)
}
18 changes: 14 additions & 4 deletions contents/verlet_integration/code/haskell/verlet.hs
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -48,7 +48,17 @@ main = do
dt = 0.001
freefall _ = -10
aboveGround (x, _, _, _) = x > 0
integrate m = last $ takeWhile aboveGround $ trajectory m freefall dt p0
print $ integrate verlet
print $ integrate stormerVerlet
print $ integrate velocityVerlet
timeVelocity m =
let (_, v, _, t) = last $ takeWhile aboveGround $ trajectory m freefall dt p0
in (show t, show v)

putStrLn "[#] Time for Verlet integration is:"
putStrLn $ fst $ timeVelocity verlet
putStrLn "[#] Time for Stormer Verlet integration is:"
putStrLn $ fst $ timeVelocity stormerVerlet
putStrLn "[#] Velocity for Stormer Verlet integration is:"
putStrLn $ snd $ timeVelocity stormerVerlet
putStrLn "[#] Time for velocity Verlet integration is:"
putStrLn $ fst $ timeVelocity velocityVerlet
putStrLn "[#] Velocity for velocity Verlet integration is:"
putStrLn $ snd $ timeVelocity velocityVerlet
16 changes: 11 additions & 5 deletions contents/verlet_integration/code/java/Verlet.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -65,14 +65,20 @@ static VerletValues velocity_verlet(double pos, double acc, double dt) {
public static void main(String[] args) {

double verletTime = verlet(5.0, -10, 0.01);
System.out.println("Time for Verlet integration is: " + verletTime);
System.out.println("[#] Time for Verlet integration is:");
System.out.println(verletTime);

VerletValues stormerVerlet = stormer_verlet(5.0, -10, 0.01);
System.out.println("Time for Stormer Verlet integration is: " + stormerVerlet.time);
System.out.println("Velocity for Stormer Verlet integration is: " + stormerVerlet.vel);
System.out.println("[#] Time for Stormer Verlet integration is:");
System.out.println(stormerVerlet.time);
System.out.println("[#] Velocity for Stormer Verlet integration is:");
System.out.println(stormerVerlet.vel);

VerletValues velocityVerlet = velocity_verlet(5.0, -10, 0.01);
System.out.println("Time for velocity Verlet integration is: " + velocityVerlet.time);
System.out.println("Velocity for velocity Verlet integration is: " + velocityVerlet.vel);
System.out.println("[#] Time for velocity Verlet integration is:");
System.out.println(velocityVerlet.time);
System.out.println("[#] Velocity for velocity Verlet integration is:");
System.out.println(velocityVerlet.vel);

}
}
15 changes: 10 additions & 5 deletions contents/verlet_integration/code/javascript/verlet.js
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,12 +45,17 @@ function velocityVerlet(pos, acc, dt) {
}

const time = verlet(5, -10, 0.01);
console.log(`Time for Verlet integration is: ${time}\n`);
console.log(`[#] Time for Verlet integration is:`);
console.log(`${time}`);

const stormer = stormerVerlet(5, -10, 0.01);
console.log(`Time for Stormer Verlet integration is: ${stormer.time}`);
console.log(`Velocity for Stormer Verlet integration is: ${stormer.vel}\n`);
console.log(`[#] Time for Stormer Verlet integration is:`);
console.log(`${stormer.time}`);
console.log(`[#] Velocity for Stormer Verlet integration is:`);
console.log(`${stormer.vel}`);

const velocity = velocityVerlet(5, -10, 0.01);
console.log(`Time for Velocity Verlet integration is: ${velocity.time}`);
console.log(`Velocity for Velocity Verlet integration is: ${velocity.vel}\n`);
console.log(`[#] Time for velocity Verlet integration is:`);
console.log(`${velocity.time}`);
console.log(`[#] Velocity for velocity Verlet integration is:`);
console.log(`${velocity.vel}`);
18 changes: 12 additions & 6 deletions contents/verlet_integration/code/julia/verlet.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -46,15 +46,21 @@ end

function main()
time = verlet(5.0, -10.0, 0.01);
println("Time for Verlet integration is: $(time)\n")
println("[#] Time for Verlet integration is:")
println("$(time)")

time, vel = stormer_verlet(5.0, -10.0, 0.01);
println("Time for Stormer Verlet integration is: $(time)")
println("Velocity for Stormer Verlet integration is: $(vel)\n")

println("[#] Time for Stormer Verlet integration is:")
println("$(time)")
println("[#] Velocity for Stormer Verlet integration is:")
println("$(vel)")

time, vel = velocity_verlet(5.0, -10.0, 0.01);
println("Time for velocity Verlet integration is: $(time)")
println("Velocity for velocity Verlet integration is: $(vel)\n")
println("[#] Time for velocity Verlet integration is:")
println("$(time)")
println("[#] Velocity for velocity Verlet integration is:")
println("$(vel)")

end

main()
15 changes: 10 additions & 5 deletions contents/verlet_integration/code/kotlin/verlet.kt
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -43,13 +43,18 @@ fun velocityVerlet(_pos: Double, acc: Double, dt: Double): VerletValues {

fun main(args: Array<String>) {
val verletTime = verlet(5.0, -10.0, 0.01)
println("Time for Verlet integration is: $verletTime")
println("[#] Time for Verlet integration is:")
println("$verletTime")

val stormerVerlet = stormerVerlet(5.0, -10.0, 0.01)
println("Time for Stormer Verlet integration is: $stormerVerlet.time")
println("Velocity for Stormer Verlet integration is: $stormerVerlet.vel")
println("[#] Time for Stormer Verlet integration is:")
println("${stormerVerlet.time}")
println("[#] Velocity for Stormer Verlet integration is:")
println("${stormerVerlet.vel}")

val velocityVerlet = velocityVerlet(5.0, -10.0, 0.01)
println("Time for Velocity Verlet integration is: $velocityVerlet.time")
println("Velocity for Velocity Verlet integration is: $velocityVerlet.vel")
println("[#] Time for Velocity Verlet integration is:")
println("${velocityVerlet.time}")
println("[#] Velocity for Velocity Verlet integration is:")
println("${velocityVerlet.vel}")
}
15 changes: 10 additions & 5 deletions contents/verlet_integration/code/nim/verlet.nim
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -46,12 +46,17 @@ func velocityVerlet(pos_in, acc, dt: float): (float, float) =

when isMainModule:
let timeV = verlet(5.0, -10.0, 0.01)
echo "Time for Verlet integration is: ", timeV
echo "[#] Time for Verlet integration is:"
echo timeV

let (timeSV, velSV) = stormerVerlet(5.0, -10.0, 0.01)
echo "Time for Stormer Verlet integration is: ", timeSV
echo "Velocity for Stormer Verlet integration is: ", velSV
echo "[#] Time for Stormer Verlet integration is:"
echo timeSV
echo "[#] Velocity for Stormer Verlet integration is:"
echo velSV

let (timeVV, velVV) = velocityVerlet(5.0, -10.0, 0.01)
echo "Time for velocity Verlet integration is: ", timeVV
echo "Velocity for velocity Verlet integration is: ", velVV
echo "[#] Time for velocity Verlet integration is:"
echo timeVV
echo "[#] Velocity for velocity Verlet integration is:"
echo velVV
22 changes: 11 additions & 11 deletions contents/verlet_integration/code/python/verlet.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -35,21 +35,21 @@ def velocity_verlet(pos, acc, dt):

def main():
time = verlet(5, -10, 0.01)
print("Verlet")
print("Time: {:.10f}".format(time))
print()
print("[#] Time for Verlet integration is:")
print("{:.10f}".format(time))

time, vel = stormer_verlet(5, -10, 0.01)
print("Stormer-Verlet")
print("Time: {:.10f}".format(time))
print("Velocity: {:.10f}".format(vel))
print()
print("[#] Time for Stormer Verlet integration is:")
print("{:.10f}".format(time))
print("[#] Velocity for Stormer Verlet integration is:")
print("{:.10f}".format(vel))

time, vel = velocity_verlet(5, -10, 0.01)
print("Velocity Verlet")
print("Time: {:.10f}".format(time))
print("Velocity: {:.10f}".format(vel))
print()
print("[#] Time for velocity Verlet integration is:")
print("{:.10f}".format(time))
print("[#] Velocity for velocity Verlet integration is:")
print("{:.10f}".format(vel))


if __name__ == '__main__':
main()
19 changes: 15 additions & 4 deletions contents/verlet_integration/code/ruby/verlet.rb
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -27,7 +27,7 @@ def stormer_verlet(pos, acc, dt)
vel += acc*dt
end

return time
return time, vel

end

Expand All @@ -41,10 +41,21 @@ def velocity_verlet(pos, acc, dt)
vel += acc*dt
end

return time
return time, vel

end

puts "[#] Time for Verlet integration is:"
p verlet(5.0, -10, 0.01)
p stormer_verlet(5.0, -10, 0.01)
p velocity_verlet(5.0, -10, 0.01)

time, vel = stormer_verlet(5.0, -10, 0.01)
puts "[#] Time for Stormer Verlet integration is:"
p time
puts "[#] Velocity for Stormer Verlet integration is:"
p vel

time, vel = velocity_verlet(5.0, -10, 0.01)
puts "[#] Time for velocity Verlet integration is:"
p time
puts "[#] Velocity for velocity Verlet integration is:"
p vel
21 changes: 12 additions & 9 deletions contents/verlet_integration/code/rust/verlet.rs
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -49,13 +49,16 @@ fn main() {
let (time_sv, vel_sv) = stormer_verlet(5.0, -10.0, 0.01);
let (time_vv, vel_vv) = velocity_verlet(5.0, -10.0, 0.01);

println!("Time for original Verlet integration: {}", time_v);
println!(
"Time and velocity for Stormer Verlet integration: {}, {}",
time_sv, vel_sv
);
println!(
"Time and velocity for velocity Verlet integration: {}, {}",
time_vv, vel_vv
);
println!("[#] Time for Verlet integration is:");
println!("{}", time_v);

println!("[#] Time for Stormer Verlet integration is:");
println!("{}", time_sv);
println!("[#] Velocity for Stormer Verlet integration is:");
println!("{}", vel_sv);

println!("[#] Time for velocity Verlet integration is:");
println!("{}", time_vv);
println!("[#] Velocity for velocity Verlet integration is:");
println!("{}", vel_vv);
}
Loading