----------------------------------------------------------------------------------
@MSGID: 1@dont-email.me> e40fd4aa
@REPLY: zW7d.89812@fx43.iad>
e3b47898
@REPLYADDR David Jones
@REPLYTO 2:5075/128 David Jones
@CHRS: CP866 2
@RFC: 1 0
@RFC-Message-ID: 1@dont-email.me>
@RFC-References:
<aa2d38fc-1df2-4703-b1b3-67adf45e2793n@googlegroups.com> <8a0d90c2-1912-41f9-8cf1-76c8f12698adn@googlegroups.com>
<c8ac7613-a6bc-4804-a668-40dbc9de519en@googlegroups.com> zW7d.89812@fx43.iad>
@TZUTC: -0000
@PID: XanaNews/1.21-f3fb89f (x86; Portable
ISpell)
@TID: FIDOGATE-5.12-ge4e8b94
Ron Shepard wrote:
> On 9/3/23 3:34 PM, Farzad Tatar wrote:
>
> I think you are on the right track with examination of the condition
> number of the matrix. As you already note, the solution to this
> problem may involve reformulations using different basis functions,
> or maybe the same type of basis functions but different orders of the
> function, or different grid spacing, and so on.
>
> Here are a couple of comments about your posted code.
>
> > dist_ = SQRT((xj_ - xi_)**2.0_dp + (yj_ - yi_)**2.0_dp)
>
>
> You should almost never use a floating point exponent with an integer
> value in fortran. You should use an integer exponent here. The
> floating point exponent expression is much more expensive (it
> requires a LOG() and and EXP() evaluation internally), and it is
> potentially less accurate.
>
.. and in the next statement the result is effectively squared, so the
SQRT calculation is unnecessary.
--- XanaNews/1.21-f3fb89f (x86; Portable ISpell)
* Origin: A noiseless patient Spider (2:5075/128)
SEEN-BY: 5001/100 5005/49 5015/255 5019/40 5020/715
848 1042 4441 12000
SEEN-BY: 5030/49 1081 5058/104 5075/128
@PATH: 5075/128 5020/1042 4441