Ideally, the sequence
1 1 1
10× / + 9× / + 9× / + 9×1 + 9×10 + 9×100 + 9×1000 = 10000.
1000 100 10
However, due to the same floating point representation problem, variations arise
Similarly, the sequence
#!/usr/bin/python | |
# -*- coding: utf-8 -*- | |
## Sympy code to generate expressions for PFHub Problem 7 (MMS) | |
from sympy import Symbol, symbols, simplify | |
from sympy import Eq, sin, cos, cosh, sinh, tanh, sqrt | |
from sympy.physics.vector import divergence, gradient, ReferenceFrame, time_derivative | |
from sympy.printing import ccode, pprint | |
from sympy.abc import kappa, S, t |
time | free_energy | |
---|---|---|
0 | 326.432 | |
1000 | 72.723657 | |
2000 | 59.844391 | |
3000 | 52.347518 | |
4000 | 49.113197 | |
5000.000001 | 45.705408 | |
6000.000001 | 44.129659 | |
7000.000001 | 41.328397 | |
8000.000001 | 40.494658 |
*.aux | |
*.log | |
*.nav | |
*.out | |
*.snm | |
*.toc | |
latex
FiPy
is an object oriented, partial differential equation (PDE) solver, written in Python
, based on a standard finite volume (FV) approach. The framework has been developed in the Materials Science and Engineering Division (MSED) and Center for Theoretical and
time | free_energy | |
---|---|---|
0 | 76.4421 | |
0.05 | 179.5 | |
0.1 | 179.467 | |
0.15 | 179.436 | |
0.2 | 179.401 | |
0.25 | 179.369 | |
0.3 | 179.333 | |
0.35 | 179.3 | |
0.4 | 179.264 |
$ Database file written 1900- 1-** | |
$ From database: USER | |
ELEMENT /- ELECTRON_GAS 0.0000E+00 0.0000E+00 0.0000E+00! | |
ELEMENT VA VACUUM 0.0000E+00 0.0000E+00 0.0000E+00! | |
ELEMENT CR BCC_A2 5.1996E+01 4.0500E+03 2.3560E+01! | |
ELEMENT NB BCC_A2 9.2906E+01 5.2200E+03 3.6270E+01! | |
ELEMENT NI FCC_A1 5.8690E+01 4.7870E+03 2.9796E+01! | |
FUNCTION GLIQCR 298.15 +24339.955-11.420225*T+2.37615E-21*T**7+GHSERCR#; |
#!/bin/bash | |
#SBATCH --partition rack6i # -p, partition | |
#SBATCH --time 12:00:00 # -t, time (hh:mm:ss or dd-hh:mm:ss) | |
#SBATCH --nodes=2 # total number of machines | |
#SBATCH --ntasks=4 # 2 sockets per rack6i machine | |
#SBATCH --cpus-per-task=28 # 28 OMP threads per rack6i MPI rank | |
#SBATCH -J Alloy625 | |
#SBATCH -D /data/tnk10/phase-field/alloy625/run2/TKR4p074 | |
export OMP_NUM_THREADS=$SLURM_CPUS_PER_TASK |
time | free_energy | |
---|---|---|
0 | 319.094 | |
106.667 | 132.607 | |
213.333 | 114.85 | |
320 | 105.152 | |
426.667 | 99.4598 | |
533.333 | 95.5889 | |
640 | 92.5611 | |
746.667 | 89.9498 | |
853.333 | 87.5588 |
time | free_energy | |
---|---|---|
0 | 319.094 | |
41.6667 | 175.655 | |
83.3333 | 140.518 | |
125 | 128.394 | |
166.667 | 121.085 | |
208.333 | 115.449 | |
250 | 110.879 | |
291.667 | 107.212 | |
333.333 | 104.281 |