Prev

Next

Index-> contents reference index search external

Up-> CppAD multi_thread thread_test.cpp multi_newton.cpp multi_newton_common

multi_thread-> parallel_ad thread_test.cpp

thread_test.cpp-> a11c_openmp.cpp a11c_bthread.cpp a11c_pthread.cpp simple_ad_openmp.cpp simple_ad_bthread.cpp simple_ad_pthread.cpp team_example.cpp harmonic.cpp multi_atomic.cpp multi_newton.cpp team_thread.hpp

multi_newton.cpp-> multi_newton_common multi_newton_setup multi_newton_worker multi_newton_takedown multi_newton_run multi_newton_time

multi_newton_common

Headings-> Purpose Source

# include <cppad/cppad.hpp>
# include <cppad/utility/time_test.hpp>
# include <cmath>
# include <cstring>
# include "multi_newton.hpp"
# include "team_thread.hpp"
# define USE_THREAD_ALLOC_FOR_WORK_ALL 1

namespace {
     using CppAD::thread_alloc; // fast multi-threadeding memory allocator
     using CppAD::vector;       // uses thread_alloc

     // number of threads, set by multi_newton_time.
     size_t num_threads_ = 0;

     // function we are finding zeros of, set by multi_newton_time
     void (*fun_)(double x, double& f, double& df) = 0;

     // convergence criteria, set by multi_newton_setup
     double epsilon_ = 0.;

     // maximum number of iterations, set by  multi_newton_setup
     size_t max_itr_ = 0;

     // length for all sub-intervals
     double sub_length_ = 0.;

     // structure with information for one thread
     typedef struct {
          // number of sub intervals (worker input)
          size_t num_sub;
          // beginning of interval (worker input)
          double xlow;
          // end of interval (worker input)
          double xup;
          // vector of zero candidates (worker output)
          // after call to multi_newton_setup:    x.size() == 0
          // after call to multi_newton_work:     x.size() is number of zeros
          // after call to multi_newton_takedown: x.size() == 0
          vector<double> x;
          // false if an error occurs, true otherwise (worker output)
          bool   ok;
     } work_one_t;
     // vector with information for all threads
     // after call to multi_newton_setup:    work_all.size() == num_threads
     // after call to multi_newton_takedown: work_all.size() == 0
     // (use pointers instead of values to avoid false sharing)
     vector<work_one_t*> work_all_;
}