17 #if defined(ADDEBUG_PRINT)
18 extern int addebug_count;
40 double * xd=
x.get_u_dot();
106 #if defined(CHECK_COUNT)
107 if (
ncount >= ncount_check)
111 constexpr
size_t sizeofdf1b2_header =
sizeof(
df1b2_header);
112 size_t total_bytes = sizeofdf1b2_header + sizeofdf1b2_header;
114 #if defined(SAFE_ALL)
116 int slen=strlen(ids);
121 #if defined(SAFE_ALL)
152 cerr <<
"illegal value for df1b2variable::pass = "
177 #if defined(SAFE_ALL)
181 constexpr
size_t sizeofdf1b2_header =
sizeof(
df1b2_header);
184 list.
bptr += sizeofdf1b2_header;
186 list.
bptr += sizeofdf1b2_header;
199 constexpr
size_t sizeofdouble =
sizeof(double);
200 size_t total_bytes = 2 * nvar * sizeofdouble;
202 #if defined(SAFE_ALL)
204 int slen=strlen(ids);
209 void * tmpptr2=list2.
bptr;
211 #if defined(SAFE_ALL)
215 memcpy(list2,pz->get_u_bar(), nvar * sizeofdouble);
216 memcpy(list2,pz->get_u_dot_bar(), nvar * sizeofdouble);
222 #if defined(PRINT_DERS)
228 double* pxu_bar = px->
u_bar;
229 double* pzu_bar = pz->u_bar;
230 for (
size_t i=0;i<nvar;i++)
232 *pxu_bar += *pzu_bar;
236 #if defined(ADDEBUG_PRINT)
238 if (addebug_count == 49)
240 cout <<
"trap" <<
endl;
242 cout << px->
u_bar[i] <<
" " << pz->u_bar[i] <<
" " << addebug_count <<
endl;
246 double* pzu_dot_bar = pz->u_dot_bar;
247 for (
size_t i=0;i<nvar;i++)
249 *pxu_dot_bar += *pzu_dot_bar;
254 size_t size = nvar * sizeofdouble;
255 memset(pz->u_bar, 0, size);
256 memset(pz->u_dot_bar, 0, size);
257 #if defined(PRINT_DERS)
277 constexpr
size_t sizeofdf1b2_header =
sizeof(
df1b2_header);
278 size_t total_bytes = sizeofdf1b2_header + sizeofdf1b2_header;
279 #if defined(SAFE_ALL)
281 int slen=strlen(ids);
298 int num_bytes2=*nlist2.
bptr;
303 #if defined(SAFE_ALL)
312 list.
bptr += sizeofdf1b2_header;
323 double * z_bar_tilde=pz->get_u_bar_tilde();
324 double * z_dot_bar_tilde=pz->get_u_dot_bar_tilde();
327 #if defined(PRINT_DERS)
333 constexpr
size_t sizeofdouble =
sizeof(double);
334 memset(z_bar_tilde, 0, nvar * sizeofdouble);
335 memset(z_dot_bar_tilde, 0, nvar * sizeofdouble);
337 double* pz_bar_tilde = z_bar_tilde;
338 double* px_bar_tilde = x_bar_tilde;
339 for (
size_t i=0;i<nvar;i++)
341 *pz_bar_tilde += *px_bar_tilde;
346 double* pz_dot_bar_tilde = z_dot_bar_tilde;
347 double* px_dot_bar_tilde = x_dot_bar_tilde;
348 for (
size_t i=0;i<nvar;i++)
350 *pz_dot_bar_tilde += *px_dot_bar_tilde;
355 #if defined(PRINT_DERS)
381 #if defined(SAFE_ALL)
384 constexpr
size_t sizeofdf1b2_header =
sizeof(
df1b2_header);
387 list.
bptr+=sizeofdf1b2_header;
391 #if defined(PRINT_DERS)
398 #
if defined(ADDEBUG_PRINT)
400 if (addebug_count == 49)
402 cout <<
"trap" <<
endl;
404 cout << *(px->
u_tilde) <<
" " << *(pz->u_tilde) <<
" " << addebug_count
408 double* pz_u_dot_tilde = pz->u_dot_tilde;
409 for (
size_t i=0;i<nvar;i++)
411 *px_u_dot_tilde += *pz_u_dot_tilde;
415 #if defined(PRINT_DERS)
420 constexpr
size_t sizeofdouble =
sizeof(double);
421 memset(pz->u_dot_tilde, 0, nvar * sizeofdouble);
void ad_read_pass1_eq(void)
Description not yet available.
void check_buffer_size(const size_t)
Description not yet available.
df1b2_gradlist * f1b2gradlist
void memcpy(test_smartlist &dest, void *source, const size_t nsize)
memcpy for test_smartlist
void read_pass1_eq_2(void)
Description not yet available.
df1b2variable & operator=(const df3_one_variable &v)
#define ADUNCONST(type, obj)
Creates a shallow copy of obj that is not CONST.
int adptr_diff(void *x, void *y)
void read_pass1_eq_3(void)
Description not yet available.
void ncount_checker(int ncount, int ncount_check)
void restoreposition(void)
Description not yet available.
Description not yet available.
Description not yet available.
prnstream & endl(prnstream &)
void print_derivatives(const adstring &s, double f, double df, double d2f, double d3f, int bflag)
Description not yet available.
static int no_derivatives
int write_pass1_eq(const df1b2variable *px, df1b2variable *pz)
Description not yet available.
Description not yet available.
void initialize(void)
Initialize.
Description not yet available.
void checkidentiferstring(const char *ids, test_smartlist &list)
Description not yet available.
int write_save_pass2_tilde_values(const df1b2variable *px)
Description not yet available.
void read_pass1_eq_1(void)
Description not yet available.