Impossible de calculer arcsinus en multiprécision

Bonjour à tous,

Je travaille avec la bibliothèque de multiprecision "boost".

Pas de problème pour calculer les sinus et cosinus avec cette bibliothèque .
Par contre, impossible de calculer les arcsinus et les arcosinus avec cette même bibliothèque.

Merci pour votre aide.

Marc75

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49

#include <boost/multiprecision/cpp_dec_float.hpp> #include <boost/math/constants/constants.hpp> #include <iostream> #include <limits> #include <vector> #include <algorithm> #include <iomanip> #include <iterator> #include <fstream> #include <string> #include <stdio.h> #include <stdlib.h> #include <math.h> #include <sys/time.h> #include <unistd.h> using namespace std;   // --------------------------------------------------------------------   int main() {   using boost::multiprecision::cpp_dec_float_50;   using boost::math::constants::pi; cpp_dec_float_50 p = pi<cpp_dec_float_50>();     int SET_PREC = 50;       cpp_dec_float_50 R = cpp_dec_float_50("1.0"); cpp_dec_float_50 C = cpp_dec_float_50("3.0");     cpp_dec_float_50 angle = cpp_dec_float_50(R/C);   cpp_dec_float_50 resu = cpp_dec_float_50(asin(angle*p/180));     std::cout << std::setprecision(SET_PREC) << std::fixed << std::showpos << resu << " " << std::endl;   }

---

Bonjour,

On ne le dira jamais assez, évitez le using namespace std.
Sinon quelle erreur obtenez vous avec ce code?

Bonjour,

Merci pour le conseil sur "using namespace std;".
J'obtiens le message (plutôt long) suivant :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21

In file included from /usr/local/include/boost/type_traits/make_unsigned.hpp:26:0, from /usr/local/include/boost/multiprecision/number.hpp:20, from /usr/local/include/boost/multiprecision/cpp_dec_float.hpp:28, from prog_cosinus_3.cpp:1: /usr/local/include/boost/math/tools/promotion.hpp: Dans l'instanciation de «*struct boost::m ath::tools::promote_args<long double, float, float, float, float, float>*»*: /usr/local/include/boost/math/special_functions/sign.hpp:124:59: requis par la substitutio n de «*template<class T> typename boost::math::tools::promote_args<T>::type boost::math::cha ngesign(const T&) [with T = long double]*» /usr/local/include/boost/lexical_cast.hpp:1031:56: requis par «*bool boost::detail::parse_ inf_nan_impl(const CharT*, const CharT*, T&, const CharT*, const CharT*, const CharT*, const CharT*, CharT, CharT) [with CharT = char; T = long double]*» /usr/local/include/boost/lexical_cast.hpp:1143:38: requis par «*bool boost::detail::parse_ inf_nan(const CharT*, const CharT*, T&) [with CharT = char; T = long double]*» /usr/local/include/boost/lexical_cast.hpp:2162:34: requis par «*bool boost::detail::lexica l_stream_limited_src<CharT, Traits, RequiresStringbuffer>::float_types_converter_internal(T& , int) [with T = long double; CharT = char; Traits = std::char_traits<char>; bool RequiresSt ringbuffer = false]*» /usr/local/include/boost/lexical_cast.hpp:2217:54: requis par «*bool boost::detail::lexica l_stream_limited_src<CharT, Traits, RequiresStringbuffer>::operator>>(long double&) [with Ch arT = char; Traits = std::char_traits<char>; bool RequiresStringbuffer = false]*» /usr/local/include/boost/lexical_cast.hpp:2350:51: [ passe outre 8 contextes d'instanciati on, utilisez -ftemplate-backtrace-limit=0 pour désactiver ] /usr/local/include/boost/multiprecision/number.hpp:1169:8: requis par «*void boost::multip recision::number<Backend, ExpressionTemplates>::do_assign_function_1(const F&, const Exp&, c onst Tag&) [with F = boost::multiprecision::detail::acos_funct<boost::multiprecision::backen ds::cpp_dec_float<50u> >; Exp = boost::multiprecision::detail::expression<boost::multiprecis ion::detail::divides, boost::multiprecision::detail::expression<boost::multiprecision::detai l::multiply_immediates, boost::multiprecision::number<boost::multiprecision::backends::cpp_d ec_float<50u> >, boost::multiprecision::number<boost::multiprecision::backends::cpp_dec_floa t<50u> >, void, void>, int, void, void>; Tag = boost::multiprecision::detail::divides; Backe nd = boost::multiprecision::backends::cpp_dec_float<50u>; boost::multiprecision::expression_ template_option ExpressionTemplates = (boost::multiprecision::expression_template_option)1u *» /usr/local/include/boost/multiprecision/number.hpp:1158:7: requis par «*void boost::multip recision::number<Backend, ExpressionTemplates>::do_assign_function(const Exp&, const mpl_::i nt_<2>&) [with Exp = boost::multiprecision::detail::expression<boost::multiprecision::detail ::function, boost::multiprecision::detail::acos_funct<boost::multiprecision::backends::cpp_d ec_float<50u> >, boost::multiprecision::detail::expression<boost::multiprecision::detail::di vides, boost::multiprecision::detail::expression<boost::multiprecision::detail::multiply_imm ediates, boost::multiprecision::number<boost::multiprecision::backends::cpp_dec_float<50u> > , boost::multiprecision::number<boost::multiprecision::backends::cpp_dec_float<50u> >, void, void>, int, void, void>, void, void>; Backend = boost::multiprecision::backends::cpp_dec_fl oat<50u>; boost::multiprecision::expression_template_option ExpressionTemplates = (boost::mu ltiprecision::expression_template_option)1u]*» /usr/local/include/boost/multiprecision/number.hpp:1063:7: requis par «*void boost::multip recision::number<Backend, ExpressionTemplates>::do_assign(const Exp&, const boost::multiprec ision::detail::function&) [with Exp = boost::multiprecision::detail::expression<boost::multi precision::detail::function, boost::multiprecision::detail::acos_funct<boost::multiprecision ::backends::cpp_dec_float<50u> >, boost::multiprecision::detail::expression<boost::multiprec ision::detail::divides, boost::multiprecision::detail::expression<boost::multiprecision::det ail::multiply_immediates, boost::multiprecision::number<boost::multiprecision::backends::cpp _dec_float<50u> >, boost::multiprecision::number<boost::multiprecision::backends::cpp_dec_fl oat<50u> >, void, void>, int, void, void>, void, void>; Backend = boost::multiprecision::bac kends::cpp_dec_float<50u>; boost::multiprecision::expression_template_option ExpressionTempl ates = (boost::multiprecision::expression_template_option)1u]*» /usr/local/include/boost/multiprecision/number.hpp:680:7: requis par «*void boost::multipr ecision::number<Backend, ExpressionTemplates>::do_assign(const boost::multiprecision::detail ::expression<tag, Arg1, Arg2, Arg3, Arg4>&, const true_&) [with tag = boost::multiprecision: :detail::function; Arg1 = boost::multiprecision::detail::acos_funct<boost::multiprecision::b ackends::cpp_dec_float<50u> >; Arg2 = boost::multiprecision::detail::expression<boost::multi precision::detail::divides, boost::multiprecision::detail::expression<boost::multiprecision: :detail::multiply_immediates, boost::multiprecision::number<boost::multiprecision::backends: :cpp_dec_float<50u> >, boost::multiprecision::number<boost::multiprecision::backends::cpp_de c_float<50u> >, void, void>, int, void, void>; Arg3 = void; Arg4 = void; Backend = boost::mu ltiprecision::backends::cpp_dec_float<50u>; boost::multiprecision::expression_template_optio n ExpressionTemplates = (boost::multiprecision::expression_template_option)1u; mpl_::true_ = mpl_::bool_<true>]*» /usr/local/include/boost/multiprecision/number.hpp:134:7: requis par «*typename boost::ena ble_if<boost::is_convertible<typename boost::multiprecision::detail::expression<tag, Arg1, A rg2, Arg3, Arg4>::result_type, boost::multiprecision::number<Backend, ExpressionTemplates> > , boost::multiprecision::number<Backend, ExpressionTemplates>&>::type boost::multiprecision: :number<Backend, ExpressionTemplates>::operator=(const boost::multiprecision::detail::expres sion<tag, Arg1, Arg2, Arg3, Arg4>&) [with tag = boost::multiprecision::detail::function; Arg 1 = boost::multiprecision::detail::acos_funct<boost::multiprecision::backends::cpp_dec_float <50u> >; Arg2 = boost::multiprecision::detail::expression<boost::multiprecision::detail::div ides, boost::multiprecision::detail::expression<boost::multiprecision::detail::multiply_imme diates, boost::multiprecision::number<boost::multiprecision::backends::cpp_dec_float<50u> >, boost::multiprecision::number<boost::multiprecision::backends::cpp_dec_float<50u> >, void, void>, int, void, void>; Arg3 = void; Arg4 = void; Backend = boost::multiprecision::backends ::cpp_dec_float<50u>; boost::multiprecision::expression_template_option ExpressionTemplates = (boost::multiprecision::expression_template_option)1u; typename boost::enable_if<boost::is _convertible<typename boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4> ::result_type, boost::multiprecision::number<Backend, ExpressionTemplates> >, boost::multipr ecision::number<Backend, ExpressionTemplates>&>::type = boost::multiprecision::number<boost: :multiprecision::backends::cpp_dec_float<50u> >&]*» /usr/local/include/boost/multiprecision/number.hpp:181:13: requis par «*boost::multiprecis ion::number<Backend, ExpressionTemplates>::number(const boost::multiprecision::detail::expre ssion<tag, Arg1, Arg2, Arg3, Arg4>&, typename boost::enable_if_c<boost::is_convertible<typen ame boost::multiprecision::detail::expression<tag, Arg1, Arg2, Arg3, Arg4>::result_type, boo st::multiprecision::number<Backend, ExpressionTemplates> >::value>::type*) [with tag = boost ::multiprecision::detail::function; Arg1 = boost::multiprecision::detail::acos_funct<boost:: multiprecision::backends::cpp_dec_float<50u> >; Arg2 = boost::multiprecision::detail::expres sion<boost::multiprecision::detail::divides, boost::multiprecision::detail::expression<boost ::multiprecision::detail::multiply_immediates, boost::multiprecision::number<boost::multipre cision::backends::cpp_dec_float<50u> >, boost::multiprecision::number<boost::multiprecision: :backends::cpp_dec_float<50u> >, void, void>, int, void, void>; Arg3 = void; Arg4 = void; Ba ckend = boost::multiprecision::backends::cpp_dec_float<50u>; boost::multiprecision::expressi on_template_option ExpressionTemplates = (boost::multiprecision::expression_template_option) 1u; typename boost::enable_if_c<boost::is_convertible<typename boost::multiprecision::detail ::expression<tag, Arg1, Arg2, Arg3, Arg4>::result_type, boost::multiprecision::number<Backen d, ExpressionTemplates> >::value>::type = void]*» prog_cosinus_3.cpp:39:60: requis depuis ici /usr/local/include/boost/math/tools/promotion.hpp:141:10: erreur*: static assertion failed: (0 == ::boost::is_same<type, long double>::value) BOOST_STATIC_ASSERT((0 == ::boost::is_same<type, long double>::value)); ^

---

Alors que le même code avec sin au lieu de asin fonctionne parfaitement ...

Bonjour,

Sur Wandbox (https://wandbox.org/), ça fonctionne avec la commande suivante g++ prog.cc -Wall -Wextra -I/opt/wandbox/boost-1.67.0/gcc-head/include -std=c++17.

Il faudrait en dire plus sur ton environnement : version de Boost, nom et version du compilateur, version de C++, etc.

Bonjour,
Mon programme a comme extension cpp comme c'est en c++.
Je compile avec les options suivantes : g++ -std=c++11 -o prog.exe prog.cpp

La version de boost utilisée est 1.54.0
La version du compilateur g++ : g++ (GCC) 6.4.0

J'ai essayé avec vos options sans succès.

Merci.

Que se trouve en prog_cosinus_3.cpp [ligne] 39 ([caractère] 60) ?
L'erreur semble indiquer que tu utilises des long double alors qu'il attend autre chose, en fait tout sauf ça : BOOST_STATIC_ASSERT((0 == ::boost::is_same<type, long double>::value));

Sur ma debian boost 1.62, g++ 6.3 ça marche aussi...
-peut être votre boost est un peu ancien, et que les versions plus récentes acceptent les long double..-

Le programme prog_cosinus_3.cpp est en fait celui affiché.

La ligne 39 correspond à :

cpp_dec_float_50 resu = cpp_dec_float_50(asin(angle*p/180));

Le caractère 60 correspond à p (qui correspond à pi qui est définit en multiprécision) :

using boost::math::constants::pi;
cpp_dec_float_50 p = pi<cpp_dec_float_50>();

Je vais upgrader les versions de g++ et de boost au niveau indiqué par kajtp. Merci pour l'info.

C'est l'upgrade de la bibliothèque boost en 1.62.0 qui a corrigé le problème.