mersenneforum.org  

Go Back   mersenneforum.org > Math Stuff > Other Mathematical Topics

Reply
 
Thread Tools
Old 2012-09-06, 02:23   #1
ewmayer
2ω=0
 
ewmayer's Avatar
 
Sep 2002
República de California

22×3×7×139 Posts
Default Webpage on Newton-Raphson

I just uploaded a first iteration (pun!) of Iterative Methods for Solving Nonlinear Algebraic Equations, a web tutorial on the Newton-Raphson method, its higher-order-convergent analogs, and applications in areas such as solution of nonlinear differential equations and number theory. (The latter are to-be-added at a later date).

It's a long page, roughly the equivalent of a book chapter, so I suggest interested readers save to their PC and peruse the individual sections at their leisure.

I have viewed the page in Firefox and Safari; let me know if you encounter any formatting issues using your browser of choice. (Note: the page is not intended for small-screen mobile devices like smartphones, although I'll be interested to hear if it's at all readable on such.)

Comments, questions, suggestions welcome.
ewmayer is offline   Reply With Quote
Old 2012-09-06, 02:52   #2
retina
Undefined
 
retina's Avatar
 
"The unspeakable one"
Jun 2006
My evil lair

2·23·137 Posts
Default

Quote:
Originally Posted by ewmayer View Post
(Note: the page is not intended for small-screen mobile devices like smartphones, although I'll be interested to hear if it's at all readable on such.)
It is perfectly readable on one of my minions iSlave devices. Although the right edge of text wraps at about 60% of the total page width. So there ends up being a 40% wide white vertical strip on the right side for the entire length of the page except for some of the tables. Table 4(a) and 4(b) end up slightly wider than the surrounding text. The wider tables appear to be because the browser won't wrap the long numbers and will render the full number all on one line.

The line with all the asterisks (in the Hilbert space section) is rendered to almost the entire page width without a break.

Also you have an errant '-->' on the page just before "Once one has computed the coefficients of one's n-term Chebyshev ...". The aqua box just above there is rendered across the entire iSlave screen and appears to be the element that is defining the maximum page width.

ETA: The 〈 and 〉 characters render as empty square boxes on the iSlave.

Last fiddled with by retina on 2012-09-06 at 03:27 Reason: The 〈 and 〉 characters render as empty square boxes on the iSlave
retina is online now   Reply With Quote
Old 2012-09-06, 03:11   #3
WraithX
 
WraithX's Avatar
 
Mar 2006

2·35 Posts
Default

Excellent web site. It will take quite a while to go through and digest all the tasty morsels it contains.

Here are some small things I noticed on first glance:
In table 3c you have index 4 listed twice, ie 0,1,2,3,4,4, instead of 0,1,2,3,4,5.

Then, right after Inv3 you have exponent 2 listed twice, ie ε^0, ε^1, ε^2, ε^2, instead of ε^0, ε^1, ε^2, ε^3

And this one I'm not sure about, but right after you say:
"The second and third derivatives might be approximated"
you give two equations, one for f'' and one for f'''. It looks like the formula for f''' should be:
f'''(y_j) \approx \frac{f(y_{j+2})-3*f(y_{j+1})+3*f(y_{j-1})-f(y_{j-2})}{2*\delta y^3}
instead of:
f'''(y_j) \approx \frac{f(y_{j+2})-2*f(y_{j+1})+2*f(y_{j-1})-f(y_{j-2})}{2*\delta y^3}
But again, this last one is just a guess.
WraithX is offline   Reply With Quote
Old 2012-09-06, 04:01   #4
bsquared
 
bsquared's Avatar
 
"Ben"
Feb 2007

DFB16 Posts
Default

Quote:
Originally Posted by Most Excellent NR tutorial
For a discussion of this x87 implementation including assembly-code samples, optimization ideas and the resulting timing data, we refer the reader to this math-forum thread. (As it happens, the foregoing thread is what prompted the present author to finally put together some long-simmering thoughts on the topic of iterative schemes in a coherent webpage form).
I thought some the material looked familiar and then saw the above quote. Glad you took the time to write it up - nice page!
bsquared is offline   Reply With Quote
Old 2012-09-06, 18:20   #5
ewmayer
2ω=0
 
ewmayer's Avatar
 
Sep 2002
República de California

22×3×7×139 Posts
Default

Quote:
Originally Posted by retina View Post
It is perfectly readable on one of my minions iSlave devices. Although the right edge of text wraps at about 60% of the total page width. So there ends up being a 40% wide white vertical strip on the right side for the entire length of the page except for some of the tables. Table 4(a) and 4(b) end up slightly wider than the surrounding text. The wider tables appear to be because the browser won't wrap the long numbers and will render the full number all on one line.
If you or any other of the iPhone users thinks of a way to improve the formatting result - perhaps by hacking the page source to just one table and surrounding text - let me know.

Quote:
The line with all the asterisks (in the Hilbert space section) is rendered to almost the entire page width without a break.
Those are unneeded - a misguided experiment in section-highlighting - so I snipped them.

Quote:
Also you have an errant '-->' on the page just before "Once one has computed the coefficients of one's n-term Chebyshev ...". The aqua box just above there is rendered across the entire iSlave screen and appears to be the element that is defining the maximum page width.
There are actually about a half-dozen stray "-->" on the rendered page - those are thanks to my text editor (which is otherwise nice except for 2-3 such annoyances) screwing up the highlighting of about every 4th inline html-comment snip, in a fashion that causes the ensuing non-comment text to be highlighted as a comment, making it very hard to read. The only workaround I've found is to add a gratuitous html comment-closing "-->" to the comments fubared in this way. I may have to simply add an item to my pre-upload checklist, where I search for such doubled "-->-->", replace with single "-->', upload to the ftp server, then undo the edit to restore readability of my source file.

Quote:
ETA: The 〈 and 〉 characters render as empty square boxes on the iSlave.
You mean the angle-brackets, e.g. in the Hilbert space section? Not sure if I can work around that.

Thanks for the notes!
ewmayer is offline   Reply With Quote
Old 2012-09-06, 18:37   #6
ewmayer
2ω=0
 
ewmayer's Avatar
 
Sep 2002
República de California

2D9C16 Posts
Default

Quote:
Originally Posted by WraithX View Post
Excellent web site. It will take quite a while to go through and digest all the tasty morsels it contains.

Here are some small things I noticed on first glance:
In table 3c you have index 4 listed twice, ie 0,1,2,3,4,4, instead of 0,1,2,3,4,5.

Then, right after Inv3 you have exponent 2 listed twice, ie ε^0, ε^1, ε^2, ε^2, instead of ε^0, ε^1, ε^2, ε^3
Thanks - fixed the repeated-index bugs in my local version (will re-upload later, after a few more edits).

The 3rd-derivative formula is correct with the 2s ... a quick way to derive it is to use the 2nd-order-accurate central-diff approx for f'' at [j+1] and [j-1], subtract the resulting approximations f''[j+1] - f''[j-1] and divide the result by 2.delta.
ewmayer is offline   Reply With Quote
Old 2012-09-06, 19:09   #7
ewmayer
2ω=0
 
ewmayer's Avatar
 
Sep 2002
República de California

22×3×7×139 Posts
Default

Found one more bug ... formula [E1] has the k-1 and k subscripts reversed. Revised version uploaded.
ewmayer is offline   Reply With Quote
Old 2012-09-06, 21:36   #8
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3×1,993 Posts
Default

This is very nice, thank you for posting.
CRGreathouse is offline   Reply With Quote
Old 2012-09-06, 23:09   #9
retina
Undefined
 
retina's Avatar
 
"The unspeakable one"
Jun 2006
My evil lair

2·23·137 Posts
Default

Quote:
Originally Posted by ewmayer View Post
You mean the angle-brackets, e.g. in the Hilbert space section? Not sure if I can work around that.
Can you not use the ASCII equivalents? <>
retina is online now   Reply With Quote
Old 2012-09-07, 00:09   #10
ewmayer
2ω=0
 
ewmayer's Avatar
 
Sep 2002
República de California

22×3×7×139 Posts
Default

Quote:
Originally Posted by retina View Post
Can you not use the ASCII equivalents? <>
I'm not sure about the ASCII equivalent thing ... does not this note in effect say <> are not equivalents of &lang; and &rang; ?

In any case, I use <> all over the page already for less-than/greater-than ... and IMO the result of replacing angle-brackets with <> is rather fugly.

I wonder if there's a simple way to to implement the analog of a C typedef here, where &lang; and &rang; would get aliased to < and > on browsers known to have problems with the former.
ewmayer is offline   Reply With Quote
Old 2012-09-07, 02:13   #11
WraithX
 
WraithX's Avatar
 
Mar 2006

2×35 Posts
Default

I'm wondering if you could just use image files for the angle brackets? I found some images online that could be tailored for your use here:
http://www.escapecodes.info/
I would re-color the background to white and trim off the excess white-space to the left and right of the actual bracket. Then you could put in <img> tags like so:
Code:
<img src="http://www.escapecodes.info/images/9001.gif" height=12/>&fnof;, <i>g</i><img src="http://www.escapecodes.info/images/9002.gif" height=12/>
This doesn't look too bad in my testing here.

I was originally thinking that CSS could handle the C typedef / multiple browser issue, but after some searching it doesn't seem like this would catch "all" browsers.

The one issue I can think of with the image-bracket approach is that getting the height of the image to match the height of the font might be difficult across different browsers and screen resolutions. There may be an easy way around this, but this was all I had time to search for tonight.
WraithX is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
security of the webpage? Unregistered Information & Answers 4 2013-02-08 04:42
The webpage cannot be found for win7 64 benDan Software 4 2012-12-17 05:59
Gimps Results to webpage? Unregistered Information & Answers 2 2009-10-27 18:15
GIMPS webpage eratos Lounge 26 2004-01-06 20:57

All times are UTC. The time now is 20:17.


Mon Dec 6 20:17:19 UTC 2021 up 136 days, 14:46, 0 users, load averages: 1.92, 1.77, 1.88

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.