Skip to content

Commit 8f8d7c8

Browse files
committed
correct f vector
f needs to be multiplied by dt/2
1 parent 0963026 commit 8f8d7c8

File tree

1 file changed

+1
-1
lines changed

1 file changed

+1
-1
lines changed

2013-12-03-Alternating_Direction_Implicit_Method.ipynb

+1-1
Original file line numberDiff line numberDiff line change
@@ -673,7 +673,7 @@
673673
"cell_type": "code",
674674
"collapsed": false,
675675
"input": [
676-
"f_vec = lambda U, V: numpy.multiply(dt, numpy.subtract(numpy.multiply(V, \n",
676+
"f_vec = lambda U, V: numpy.multiply(float(dt)/float(2.), numpy.subtract(numpy.multiply(V, \n",
677677
" numpy.add(k0, numpy.divide(numpy.multiply(U,U), numpy.add(1., numpy.multiply(U,U))))), U))\n",
678678
"\n",
679679
"k0 = 0.067"

0 commit comments

Comments
 (0)