On Sat, Sep 11, 2010 at 2:14 PM, Brian Gough <[email protected]> wrote: > At Wed, 8 Sep 2010 08:01:23 +0530, > Awhan Patnaik wrote: >> What can I do to get as much accuracy as the Octave method? > > Did you compute the norm from the printed output data or the data in > memory? Note that it is only printed to 7 digits, so you would lose > ~11 digits of precision in that case. >
That is precisely where I was going wrong. The loss of significant digits while transferring data to Octave through a data file resulted in the loss of accuracy. When the norm is computed in the C++ code all is right. Many thanks to Francesco Abbate for helping me in this. _______________________________________________ Help-gsl mailing list [email protected] http://lists.gnu.org/mailman/listinfo/help-gsl
