Changeset 61 for trunk/ADOLC/src/adouble.cpp
 Timestamp:
 Dec 7, 2009 9:49:34 AM (10 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/ADOLC/src/adouble.cpp
r42 r61 25 25 #include <taping_p.h> 26 26 27 using namespace std; 28 27 29 /****************************************************************************/ 28 30 /* HELPFUL FUNCTIONS */ … … 1552 1554 1553 1555 ADOLC_GLOBAL_TAPE_VARS.store[locat] = 1554 ADOLC_MATH_NSP ::asinh(ADOLC_GLOBAL_TAPE_VARS.store[x.location]);1556 ADOLC_MATH_NSP_ERF::asinh(ADOLC_GLOBAL_TAPE_VARS.store[x.location]); 1555 1557 ADOLC_OPENMP_RESTORE_THREAD_NUMBER; 1556 1558 return locat; … … 1578 1580 1579 1581 ADOLC_GLOBAL_TAPE_VARS.store[locat] = 1580 ADOLC_MATH_NSP ::acosh(ADOLC_GLOBAL_TAPE_VARS.store[x.location]);1582 ADOLC_MATH_NSP_ERF::acosh(ADOLC_GLOBAL_TAPE_VARS.store[x.location]); 1581 1583 ADOLC_OPENMP_RESTORE_THREAD_NUMBER; 1582 1584 return locat; … … 1604 1606 1605 1607 ADOLC_GLOBAL_TAPE_VARS.store[locat] = 1606 ADOLC_MATH_NSP ::atanh(ADOLC_GLOBAL_TAPE_VARS.store[x.location]);1608 ADOLC_MATH_NSP_ERF::atanh(ADOLC_GLOBAL_TAPE_VARS.store[x.location]); 1607 1609 ADOLC_OPENMP_RESTORE_THREAD_NUMBER; 1608 1610 return locat; … … 1616 1618 locint locat = next_loc(); 1617 1619 1618 adouble y = exp(x*x); 1620 adouble y = 2.0 / 1621 ADOLC_MATH_NSP_ERF::sqrt(ADOLC_MATH_NSP::acos(1.0))*exp(x*x); 1619 1622 1620 1623 if (ADOLC_CURRENT_TAPE_INFOS.traceFlag) { // old: write_quad(erf_op,locat,x.location,y.location); … … 1630 1633 1631 1634 ADOLC_GLOBAL_TAPE_VARS.store[locat] = 1632 ADOLC_MATH_NSP ::erf(ADOLC_GLOBAL_TAPE_VARS.store[x.location]);1635 ADOLC_MATH_NSP_ERF::erf(ADOLC_GLOBAL_TAPE_VARS.store[x.location]); 1633 1636 ADOLC_OPENMP_RESTORE_THREAD_NUMBER; 1634 1637 return locat;
Note: See TracChangeset
for help on using the changeset viewer.