Alright, thanks!
Does this mean that I have to mpf_set_default_prec always before the first
mpf_... declaration?
Best,
P
On Monday, December 14, 2015 at 11:24:02 AM UTC+1, Bill Hart wrote:
>
> You are setting the default precision before initialising z. I guess this
> means z has the previous default precision, which is probably 64 bits or
> something like that.
>
> Bill.
>
> On 14 December 2015 at 11:16, <highcal...@gmail.com <javascript:>> wrote:
>
>> Dear Bill, Dear All,
>>
>> I am trying to simply assign a passed mpf_t (y in the function) to
>> another mpf_t (x). The code suggest inaccuracy after 20 digits.
>> Here is the function, ctest.c:
>>
>> #include <stdio.h>
>> #include "mpir.h"
>> #include "ctest.h"
>>
>> void ctest(mpf_t y)
>> {
>> mpf_set_default_prec(10000);
>>
>> mpf_t x; mpf_inits(x, 0);
>> gmp_printf("x, y are %.*Ff, %.*Ff \n", 70, x, 70, y);
>> mpf_set(x, y);
>> gmp_printf("x, y are %.*Ff, %.*Ff \n", 70, x, 70, y);
>>
>> mpf_clears(x, 0);
>> }
>>
>> and here is the calling routine, ctest_main.c:
>>
>> #include <stdio.h>
>> #include "mpir.h"
>> #include "ctest.h"
>>
>> int main(void) {
>> FILE *finput; char tfile[] = "t1.txt"; // contains one number:
>> 1000000000000.000 and 80 trailing zeros
>>
>> mpf_t z; mpf_inits(z, 0);
>>
>> mpf_set_default_prec(10000);
>>
>> finput = fopen(tfile, "r");
>> mpf_inp_str(z, finput, 10);
>> fclose(finput);
>>
>> ctest(z);
>>
>> mpf_clears(z, 0); return(0);
>> }
>>
>> This results in the output:
>>
>> x, y are
>> 0.0000000000000000000000000000000000000000000000000000000000000000000000,
>> 1000000000000.0000000000000000000000000000000000000000000000000000000000000000000000
>>
>>
>> x, y are
>> 1000000000000.0000000000000000001626303258728256651011179201304912567138671875000000
>> ,
>> 1000000000000.0000000000000000000000000000000000000000000000000000000000000000000000
>>
>>
>> >>
>> I am compiling (via Matlab, but this has not been an issue so far) by
>>
>> !gcc ctest_main.c ctest.c -LC:/MPIR/mpir-2.7.0/ -lmpir -o ctest_main
>>
>> Declaring y as an mpf_t in the function gave an error by the compiler. Is
>> this a bug?
>>
>> Thanks,
>> Pat
>>
>>
>> --
>> You received this message because you are subscribed to the Google Groups
>> "mpir-devel" group.
>> To unsubscribe from this group and stop receiving emails from it, send an
>> email to mpir-devel+...@googlegroups.com <javascript:>.
>> To post to this group, send email to mpir-...@googlegroups.com
>> <javascript:>.
>> Visit this group at https://groups.google.com/group/mpir-devel.
>> For more options, visit https://groups.google.com/d/optout.
>>
>
>
--
You received this message because you are subscribed to the Google Groups
"mpir-devel" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to mpir-devel+unsubscr...@googlegroups.com.
To post to this group, send email to mpir-devel@googlegroups.com.
Visit this group at https://groups.google.com/group/mpir-devel.
For more options, visit https://groups.google.com/d/optout.
