[arch-commits] Commit in mpfi/repos (3 files)
Felix Yan
felixonmars at archlinux.org
Tue Jul 7 17:37:16 UTC 2020
Date: Tuesday, July 7, 2020 @ 17:37:16
Author: felixonmars
Revision: 659809
archrelease: copy trunk to community-staging-x86_64
Added:
mpfi/repos/community-staging-x86_64/
mpfi/repos/community-staging-x86_64/PKGBUILD
(from rev 659808, mpfi/trunk/PKGBUILD)
mpfi/repos/community-staging-x86_64/div_ext.c
(from rev 659808, mpfi/trunk/div_ext.c)
-----------+
PKGBUILD | 29 ++++++++++++
div_ext.c | 142 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
2 files changed, 171 insertions(+)
Copied: mpfi/repos/community-staging-x86_64/PKGBUILD (from rev 659808, mpfi/trunk/PKGBUILD)
===================================================================
--- community-staging-x86_64/PKGBUILD (rev 0)
+++ community-staging-x86_64/PKGBUILD 2020-07-07 17:37:16 UTC (rev 659809)
@@ -0,0 +1,29 @@
+# Maintainer: Antonio Rojas <arojas at archlinux.org>
+# Contributor: Rémy Oudompheng <remy at archlinux.org>
+
+pkgname=mpfi
+pkgver=1.5.4
+pkgrel=3
+pkgdesc="C library for interval arithmetic"
+arch=(x86_64)
+url="http://perso.ens-lyon.fr/nathalie.revol/software.html"
+license=(GPL)
+depends=(mpfr)
+source=("https://gforge.inria.fr/frs/download.php/file/38111/mpfi-$pkgver.tgz")
+sha256sums=('3b3938595d720af17973deaf727cfc0dd41c8b16c20adc103a970f4a43ae3a56')
+
+prepare() {
+ cd $pkgname-$pkgver
+ ./autogen.sh
+}
+
+build() {
+ cd $pkgname-$pkgver
+ ./configure --prefix=/usr --enable-shared
+ make
+}
+
+package() {
+ cd $pkgname-$pkgver
+ make DESTDIR="$pkgdir" install
+}
Copied: mpfi/repos/community-staging-x86_64/div_ext.c (from rev 659808, mpfi/trunk/div_ext.c)
===================================================================
--- community-staging-x86_64/div_ext.c (rev 0)
+++ community-staging-x86_64/div_ext.c 2020-07-07 17:37:16 UTC (rev 659809)
@@ -0,0 +1,142 @@
+/* div_ext.c -- Divide two intervals:
+ when the denominator contains 0,
+ returns two intervals and a gap between them.
+
+ The returned integer is the number of returned intervals.
+
+Copyright 2018
+ AriC project, Inria Grenoble - Rhone-Alpes, France
+
+This file is part of the MPFI Library.
+
+The MPFI Library is free software; you can redistribute it and/or modify
+it under the terms of the GNU Lesser General Public License as published by
+the Free Software Foundation; either version 2.1 of the License, or (at your
+option) any later version.
+
+The MPFI Library is distributed in the hope that it will be useful, but
+WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
+or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public
+License for more details.
+
+You should have received a copy of the GNU Lesser General Public License
+along with the MPFI Library; see the file COPYING.LIB. If not, write to
+the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston,
+MA 02110-1301, USA. */
+
+#include "mpfi-impl.h"
+
+int
+mpfi_div_ext (mpfi_ptr res1, mpfi_ptr res2, mpfi_srcptr op1, mpfi_srcptr op2)
+{
+ mpfr_t tmp1, tmp2;
+ int nb_res = 0, tmp;
+
+ if ( MPFI_NAN_P (op1) || MPFI_NAN_P (op2) ) {
+ mpfr_set_nan (&(res1->left));
+ mpfr_set_nan (&(res1->right));
+ mpfr_set_nan (&(res2->left));
+ mpfr_set_nan (&(res2->right));
+ return 0;
+ }
+
+ if (MPFI_HAS_ZERO (op2)) {
+ if ( MPFR_IS_INF(&(op2->left)) && MPFR_IS_INF(&(op2->right)) ) {
+ mpfr_set_inf (&(res1->left), -1);
+ mpfr_set_inf (&(res1->right), 1);
+ mpfr_set_nan (&(res2->left));
+ mpfr_set_nan (&(res2->right));
+ nb_res = 1;
+ }
+ else if ( MPFI_HAS_ZERO_NONSTRICT (op1) ) {
+ mpfr_set_inf (&(res1->left), -1);
+ mpfr_set_inf (&(res1->right), 1);
+ mpfr_set_nan (&(res2->left));
+ mpfr_set_nan (&(res2->right));
+ nb_res = 1;
+ }
+ else if ( MPFI_IS_STRICTLY_NEG (op1) ) {
+ mpfr_init2 (tmp1, mpfi_get_prec(res1));
+ mpfr_init2 (tmp2, mpfi_get_prec(res2));
+ if ( mpfr_number_p (&(op2->left)) ) {
+ tmp = mpfr_div (&(tmp2), &(op1->right), &(op2->left), MPFI_RNDD);
+ }
+ else { /* denominator has infinite left endpoint */
+ mpfr_set_zero (&(tmp2), 1);
+ }
+
+ if ( mpfr_number_p (&(op2->right)) ) {
+ tmp = mpfr_div ( &(tmp1), &(op1->right), &(op2->right), MPFI_RNDU);
+ }
+ else { /* denominator has infinite right endpoint */
+ mpfr_set_zero( &(tmp1), -1);
+ }
+
+ mpfr_set_inf (&(res1->left), -1);
+ mpfr_set (&(res1->right), tmp1, MPFI_RNDU);
+ mpfr_set (&(res2->left), tmp2, MPFI_RNDD);
+ mpfr_set_inf (&(res2->right), 1);
+
+ mpfr_clear(tmp1);
+ mpfr_clear(tmp2);
+
+ nb_res = 2;
+ }
+ else { /* if ( MPFI_IS_STRICTLY_POS (op1) ) */
+ mpfr_init2 (tmp1, mpfi_get_prec(res1));
+ mpfr_init2 (tmp2, mpfi_get_prec(res2));
+ if ( mpfr_number_p (&(op2->left)) ) {
+ tmp = mpfr_div (&(tmp1), &(op1->left), &(op2->left), MPFI_RNDU);
+ }
+ else { /* denominator has infinite left endpoint */
+ mpfr_set_zero (&(tmp1), -1);
+ }
+
+ if ( mpfr_number_p (&(op2->right)) ) {
+ tmp = mpfr_div ( &(tmp2), &(op1->left), &(op2->right), MPFI_RNDD);
+ }
+ else { /* denominator has infinite right endpoint */
+ mpfr_set_zero( &(tmp2), 1);
+ }
+ mpfr_set_inf (&(res1->left), -1);
+ mpfr_set (&(res1->right), tmp1, MPFI_RNDU);
+ mpfr_set (&(res2->left), tmp2, MPFI_RNDD);
+ mpfr_set_inf (&(res2->right), 1);
+
+ mpfr_clear(tmp1);
+ mpfr_clear(tmp2);
+
+ nb_res = 2;
+ }
+ }
+
+ else {
+ tmp = mpfi_div(res1, op1, op2);
+ mpfr_set_nan (&(res2->left));
+ mpfr_set_nan (&(res2->right));
+ return 1;
+ }
+
+ if (nb_res >= 1) {
+ /* do not allow -0 as lower bound */
+ if (mpfr_zero_p (&(res1->left)) && mpfr_signbit (&(res1->left))) {
+ mpfr_neg (&(res1->left), &(res1->left), MPFI_RNDU);
+ }
+ /* do not allow +0 as upper bound */
+ if (mpfr_zero_p (&(res1->right)) && !mpfr_signbit (&(res1->right))) {
+ mpfr_neg (&(res1->right), &(res1->right), MPFI_RNDD);
+ }
+ }
+ if (nb_res >= 2) {
+ /* do not allow -0 as lower bound */
+ if (mpfr_zero_p (&(res2->left)) && mpfr_signbit (&(res2->left))) {
+ mpfr_neg (&(res2->left), &(res2->left), MPFI_RNDU);
+ }
+ /* do not allow +0 as upper bound */
+ if (mpfr_zero_p (&(res2->right)) && !mpfr_signbit (&(res2->right))) {
+ mpfr_neg (&(res2->right), &(res2->right), MPFI_RNDD);
+ }
+ }
+
+ return nb_res;
+}
More information about the arch-commits
mailing list