void PvaPyLogger::critical(const char* message, va_list messageArgs) const { log(LogLevelCritical, message, messageArgs); }
void WServer::Initialize() { //Start time PTBeforeLoad = boost::posix_time::microsec_clock::local_time(); /* ************************************************************************* * *********************** Connect to SQL Server ************************* * *************************************************************************/ try { log("info") << "Connecting to database backend"; Wt::Dbo::SqlConnection *SQLConnection = new Wt::Dbo::backend::MySQL("wt", "root", "", "127.0.0.1"); //Wt::Dbo::SqlConnection *SQLConnection = new Wt::Dbo::backend::Sqlite3(":memory:"); SQLConnection->setProperty("show-queries", "true"); SQLPool = new Wt::Dbo::FixedSqlConnectionPool(SQLConnection, 1); log("success") << "Successfully connected to database"; } catch(Wt::Dbo::Exception &e) { log("fatal") << "Database error connecting to database: " << e.what(); throw e; } catch(std::exception &e) { log("fatal") << "Error connecting to database: " << e.what(); throw e; } /* ************************************************************************* * *********************** Initialize Databases ************************** * *************************************************************************/ try { _DboManager = new DboDatabaseManager(this, SQLPool); _Modules = new ModulesDatabase(_DboManager); _Configurations = new ConfigurationsDatabase(_DboManager); _Languages = new LanguagesDatabase(_DboManager); _Styles = new StylesDatabase(_DboManager); _Pages = new PagesDatabase(_DboManager); _AccessPaths = new AccessPathsDatabase(_DboManager); _NavigationMenus = new NavigationMenusDatabase(_DboManager); } catch(std::exception &e) { log("fatal") << "Error initializing databases: " << e.what(); throw e; } /* ************************************************************************* * *************************** Create Tables ***************************** * *************************************************************************/ #define REINSTALLDBO 0 #if REINSTALLDBO == 1 //Drop try { _Installer = new DboInstaller(*SQLPool); _Installer->DropTables(); } catch(Wt::Dbo::Exception &e) { log("error") << "Database error dropping tables: " << e.what(); } catch(std::exception &e) { log("error") << "Error dropping tables: " << e.what(); } //Create try { _Installer->CreateTables(); } catch(Wt::Dbo::Exception &e) { log("fatal") << "Database error creating tables: " << e.what(); throw e; } catch(std::exception &e) { log("fatal") << "Error creating tables: " << e.what(); throw e; } //Insert try { _Installer->InsertRows(); } catch(Wt::Dbo::Exception &e) { log("fatal") << "Database error inserting default data: " << e.what(); throw e; } catch(std::exception &e) { log("fatal") << "Error inserting default data: " << e.what(); throw e; } #endif /* ************************************************************************* * ************************* Load DboDatabases *************************** * *************************************************************************/ // try // { log("info") << "Loading DboDatabaseManager"; _DboManager->Load(); // } // catch(Wt::Dbo::Exception &e) // { // log("fatal") << "Database error loading DboDatabaseManager: " << e.what(); // throw e; // } // catch(std::exception &e) // { // log("fatal") << "Error loading DboDatabaseManager: " << e.what(); // throw e; // } //Server localized strings setLocalizedStrings(new DboLocalizedStrings(this)); Pages()->Register404PageHandler(new Default404Page()); Pages()->RegisterPageHandler("home", ModulesDatabase::Navigation, new TestPage()); Pages()->RegisterPageHandler("sitemap", ModulesDatabase::Navigation, new AnotherPage()); /* ************************************************************************* * ********************* Create temporary XML file *********************** * *************************************************************************/ try { log("info") << "Writing XML Configuration file"; CreateWtXmlConfiguration(); log("success") << "XML Configuration file created"; } catch(std::exception &e) { log("fatal") << "Error creating XML Configuration file: " << e.what(); throw e; } //Configure authorization module ConfigureAuth(); }
void TDDSwitchTest::setUp() { log("TDD Setup is called"); log("Please write somethings"); }
void userauth_finish(Authctxt *authctxt, char *method) { int authenticated, partial; if (authctxt == NULL) fatal("%s: missing context", __func__); /* unknown method handling -- must elicit userauth failure msg */ if (authctxt->method == NULL) { authenticated = 0; partial = 0; goto done_checking; } #ifndef USE_PAM /* Special handling for root (done elsewhere for PAM) */ if (!use_privsep && authctxt->method->authenticated && authctxt->pw != NULL && authctxt->pw->pw_uid == 0 && !auth_root_allowed(method)) authctxt->method->authenticated = 0; #endif /* USE_PAM */ #ifdef _UNICOS if (authctxt->method->authenticated && cray_access_denied(authctxt->user)) { authctxt->method->authenticated = 0; fatal("Access denied for user %s.",authctxt->user); } #endif /* _UNICOS */ partial = userauth_check_partial_failure(authctxt); authenticated = authctxt->method->authenticated; #ifdef USE_PAM /* * If the userauth method failed to complete PAM work then force * partial failure. */ if (authenticated && !AUTHPAM_DONE(authctxt)) partial = 1; #endif /* USE_PAM */ /* * To properly support invalid userauth method names we set * authenticated=0, partial=0 above and know that * authctxt->method == NULL. * * No unguarded reference to authctxt->method allowed from here. * Checking authenticated != 0 is a valid guard; authctxt->method * MUST NOT be NULL if authenticated. */ done_checking: if (!authctxt->valid && authenticated) { /* * Should never happen -- if it does PAM's at fault * but we need not panic, just treat as a failure. */ authctxt->method->authenticated = 0; authenticated = 0; log("Ignoring authenticated invalid user %s", authctxt->user); auth_log(authctxt, 0, method, " ssh2"); } /* Log before sending the reply */ auth_log(authctxt, authenticated, method, " ssh2"); if (authenticated && !partial) { /* turn off userauth */ dispatch_set(SSH2_MSG_USERAUTH_REQUEST, &dispatch_protocol_ignore); packet_start(SSH2_MSG_USERAUTH_SUCCESS); packet_send(); packet_write_wait(); /* now we can break out */ authctxt->success = 1; } else { char *methods; if (authctxt->method && authctxt->method->is_initial) authctxt->init_failures++; authctxt->method = NULL; #ifdef USE_PAM /* * Keep track of last PAM error (or PERM_DENIED) for BSM * login failure auditing, which may run after the PAM * state has been cleaned up. */ authctxt->pam_retval = AUTHPAM_ERROR(authctxt, PAM_PERM_DENIED); #endif /* USE_PAM */ if (authctxt->failures++ > options.max_auth_tries) { #ifdef HAVE_BSM fatal_remove_cleanup(audit_failed_login_cleanup, authctxt); audit_sshd_login_failure(&ah, PAM_MAXTRIES); #endif /* HAVE_BSM */ packet_disconnect(AUTH_FAIL_MSG, authctxt->user); } #ifdef _UNICOS if (strcmp(method, "password") == 0) cray_login_failure(authctxt->user, IA_UDBERR); #endif /* _UNICOS */ packet_start(SSH2_MSG_USERAUTH_FAILURE); /* * If (partial) then authmethods_get() will return only * required methods, likely only "keyboard-interactive;" * (methods == NULL) implies failure, even if (partial == 1) */ methods = authmethods_get(); packet_put_cstring(methods); packet_put_char((authenticated && partial && methods) ? 1 : 0); if (methods) xfree(methods); packet_send(); packet_write_wait(); } }
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[]) { /* description: * given a set of reference stimuli and responses (forming a joint distribution), * this function will compute the log kernel density estimate over a grid of stimuli * for a set of test responses. An epanechnikov kernel is used for all dimensions * and the data is assumed to be normalized by the desired kernel width. * Note that a radial-symmetric epanechnikov kernel is used, rather than a multiplicative kernel. * * syntax: * out = kde_decode_epanechnikov( stimulus, stimulus_grid, response, test_response[, offset[, distance[, timestamp, bins]]] ); * * arguments: * stimulus = number of spikes x number of stimulus dimensions * stimulus_grid = number of grid elements x number of stimulus dimensions * response = number of spikes x number of response dimensions * test_response = number of test spikes x number of response dimensions * offset = number of grid elements * distance = optional matrix of distances. If this argument is present (and not empty), the stimulus and stimulus_grid arguments should be given as a zero-based index into the distance matrix. * timestamp = number of test spikes x 1 * bins = number of bins x 2 * * out = number of test spikes x number of stimulus grid elements (if no timestamps and bins arguments are provided) * out = number of bins x number of stimulus grid elements (if timestamps and bins arguments are provided) * */ static const double pi = 3.141592653589793238462643383279502884197; /* INPUTS */ double *stimulus; /* NxQ */ double *stimulus_grid; /* GxQ */ double *response; /* NxD or Nx0 or empty*/ double *test_response; /* MxD or Mx0 or empty*/ double *timestamp; double *bins; double *offset; /* VARIABLES */ int N, D, M, Q, G, B; int n, d, m, q, g, b; int sizeM, loopM; int next_idx, idx; double acc_a, *acc_g, *skip_g, *z; double tmp; double *pout, *pout2; mxArray *out, *out2; int idx1, idx2; mxArray *tmpmat; int *use_distance_lookup; int *NI; double **distance; double scaling_factor = 1; double v; double c1, c2; mxArray *argout=NULL, *argin[1]; /* CHECK NUMBER OF INPUTS */ if (nrhs<2) mexErrMsgTxt("This function requires at least two input arguments"); /* CHECK DIMENSIONS OF FIRST TWO ARGUMENTS */ if ( mxGetNumberOfDimensions(prhs[0])!=2 || mxGetNumberOfDimensions(prhs[1])!=2) mexErrMsgTxt("stimulus and stimulus_grid input arguments need to be matrices"); /* GET POINTERS TO FIRST TWO ARGUMENTS */ stimulus = mxGetPr( prhs[0] ); stimulus_grid = mxGetPr( prhs[1] ); /* GET ARRAY SIZES */ N = mxGetM( prhs[0] ); /* number of source (encoding) spikes */ Q = mxGetN( prhs[0] ); /* number of stimulus dimensions */ G = mxGetM( prhs[1] ); /* number of points in stimulus grid */ /* CHECK NUMBER OF DIMENSIONS IN STIMULUS GRID */ if (mxGetN(prhs[1])!=Q) mexErrMsgTxt("Incompatible size of stimulus_grid input arguments"); /* CHECK OPTIONAL INPUTS */ if (nrhs>2) { /* CHECK DIMENSIONALITY OF SPIKE RESPONSE (ENCODING)*/ if ( mxGetNumberOfDimensions(prhs[2])!=2 ) mexErrMsgTxt("Response input argument needs to be a matrix"); response = mxGetPr( prhs[2] ); D = mxGetN( prhs[2] ); /* number of response dimensions */ /* CHECK NUMBER OF SPIKES IN RESPONSE */ if ( (D>0 && mxGetM(prhs[2])!=N ) ) mexErrMsgTxt("Incompatible size of response input argument"); } else { D = 0; } if (nrhs>3) { /* CHECK DIMENSIONALITY OF TEST RESPONSE (DECODING)*/ if ( mxGetNumberOfDimensions(prhs[3])!=2 ) mexErrMsgTxt("Test_response input argument needs to be a matrix"); test_response = mxGetPr( prhs[3] ); M = mxGetM( prhs[3] ); /* number of test (decoding) spikes */ /* CHECK NUMBER OF DIMENSIONS IN TEST RESPONSE */ if ( mxGetN(prhs[3])!=D ) mexErrMsgTxt("Incompatible size of test_response input argument"); } else { M = 0; /* no test spikes specified */ } if (D==0 && M==0 ) M = 1; if (nrhs>4) { /* CHECK SIZE OF OFFSET VECTOR */ if ( !mxIsDouble( prhs[4] ) || mxGetNumberOfElements( prhs[4] )!=G ) mexErrMsgTxt("Incompatible size of offset vector"); offset = mxGetPr( prhs[4] ); } else { mexErrMsgTxt("Please provide offset vector"); } /* ALLOCATE ARRAYS FOR DISTANCE LOOK-UP TABLES */ use_distance_lookup = (int*) mxCalloc( Q, sizeof(int) ); distance = (double**) mxCalloc( Q, sizeof(double*) ); NI = (int*) mxCalloc( Q, sizeof(int) ); /* size of distance LUTs */ /* INITIALIZE ARRAY */ for (q=0;q<Q;q++) use_distance_lookup[q] = 0; if (nrhs>5) { /* CHECK CLASS AND SIZE OF DISTANCE LUTs */ if ( !mxIsCell( prhs[5] ) || mxGetNumberOfElements( prhs[5] )!=Q ) mexErrMsgTxt("Distance input argument needs to be a cell array with as many cells as stimulus dimensions"); for ( q=0 ; q<Q ; q++ ) { tmpmat = mxGetCell( prhs[5], q ); if (tmpmat==NULL || mxIsEmpty(tmpmat)) { use_distance_lookup[q] = 0; } else { /* CHECK SIZE OF DISTANCE LUT */ if ( mxGetNumberOfDimensions(tmpmat)!=2 || mxGetM( tmpmat )!=mxGetN( tmpmat ) ) mexErrMsgTxt("Distance arrays need to be square matrices"); distance[q] = (double*) mxGetPr( tmpmat ); NI[q] = mxGetM( tmpmat ); use_distance_lookup[q] = true; /* when using distance LUT, the corresponding stimulus and stimulus grid should be indices into the LUT */ /* CHECK IF STIMULUS AND STIMULUS_GRID >=0 && <NI[q] */ for ( n=0; n<N; n++ ) { if ( stimulus[n+q*N]<0 || stimulus[n+q*N]>=NI[q] ) mexErrMsgTxt("Invalid index"); } for ( g=0; g<G; g++ ) { if ( stimulus_grid[g+q*G]<0 || stimulus_grid[g+q*G]>=NI[q] ) mexErrMsgTxt("Invalid index"); } } } } B = 0; /* number of time bins */ if (nrhs>6) { /* CHECK DIMENSIONALITY OF TEST (DECODING) SPIKE TIMESTAMPS */ if ( mxGetNumberOfDimensions(prhs[6])!=2 ) mexErrMsgTxt("Timestamp input argument needs to be a matrix"); timestamp = mxGetPr( prhs[6] ); /* CHECK SIZE OF TEST (DECODING) SPIKE TIMESTAMPS */ if ( mxGetM( prhs[6] )!=M || mxGetN( prhs[6] )!=1 ) mexErrMsgTxt("Incompatible size of timestamp input argument"); } if (nrhs>7) { /* CHECK DIMENSIONALITY AND SIZE OF TIME BINS ARGUMENT */ if ( mxGetNumberOfDimensions(prhs[7])!=2 || mxGetN(prhs[7])!=2) mexErrMsgTxt("Bins input argument needs to be a Bx2 matrix"); B = mxGetM( prhs[7] ); /* number of time bins */ bins = mxGetPr( prhs[7] ); } /* COMPUTE SCALING FACTOR */ argin[0] = mxCreateDoubleScalar( 1 + 0.5*(double)(D+Q) ); /* compute volume of hypersphere */ mexCallMATLAB(1, &argout, 1, argin, "gamma" ); v = pow(pi, 0.5*(double)(D+Q))/mxGetScalar(argout); scaling_factor *= 0.5*(double)(D+Q+2)/v; /* ALLOCATE TEMPORARY ARRAYS AND OUTPUT ARRAYS */ acc_g = (double*) mxCalloc( G, sizeof(double) ); skip_g = (double*) mxCalloc( G, sizeof(double) ); if (B==0) { out = mxCreateDoubleMatrix( M, G, mxREAL ); z = mxGetPr( out ); } else { if (D==0) { z = (double*) mxCalloc( G, sizeof(double) ); } else { z = (double*) mxCalloc( M*G, sizeof(double) ); } out = mxCreateDoubleMatrix( B, G, mxREAL ); pout = mxGetPr( out ); out2 = mxCreateDoubleMatrix( B, 1, mxREAL ); pout2 = mxGetPr( out2 ); } loopM = sizeM = M; if (D==0) { loopM = 1; if (B>0) sizeM = 1; } /* COMPUTE KDE */ /* LOOP THROUGH SOURCE (ENCODING) SPIKES */ for ( n=0; n<N; n++ ) { /* LOOP THROUGH STIMULUS GRID */ for ( g=0; g<G; g++ ) { /* INITIALIZE ACCUMULATORS */ acc_g[g] = 0; skip_g[g] = 0; /* INITIALIZE INDICES */ idx1 = g; idx2 = n; /* LOOP THROUGH STIMULUS DIMENSIONS */ for ( q=0; q<Q; q++ ) { if (use_distance_lookup[q]) { tmp = distance[q][ ((int) stimulus_grid[idx1])*NI[q] + (int)stimulus[idx2] ]; } else { tmp = (stimulus_grid[idx1]-stimulus[idx2]); } tmp *= tmp; acc_g[g] += tmp; if (acc_g[g]>1) { skip_g[g] = 1; break; } /* UPDATE INDICES */ idx1 += G; idx2 += N; } } /* LOOP THROUGH TEST (DECODING) SPIKES */ for ( m=0; m<loopM; m++ ) { /* INITIALIZE ACCUMULATORS */ acc_a = 0; /* INTIALIZE INDICES */ idx1 = m; idx2 = n; /* LOOP THROUGH RESPONSE DIMENSIONS */ for ( d=0; d<D; d++ ) { tmp = (test_response[idx1]-response[idx2]); tmp *= tmp; acc_a += tmp; if (acc_a>1) { goto nextm; } /* UPDATE INDICES */ idx1 += sizeM; idx2 += N; } for ( g=0; g<G; g++ ) { if (skip_g[g] || (acc_g[g]+acc_a)>1 ) continue; z[m+g*sizeM] += (1-(acc_g[g]+acc_a)); } nextm: ; } } mxFree(acc_g); mxFree(skip_g); if (argout!=NULL) mxDestroyArray(argout); mxDestroyArray(argin[0]); /* compute log */ c1 = log(scaling_factor) - log((double)N); if (B==0 && D==0) { for (g=0;g<G;g++) z[g*sizeM]=log(z[g*sizeM] + offset[g]*(double)N/scaling_factor) + c1; } else { for (g=0; g<G; g++) { idx = g*loopM; c2 = offset[g]*(double)N/scaling_factor; for (m=0; m<loopM ; m++ ) z[m+idx] = log(z[m+idx] + c2) + c1; } } if (D==0 && B==0) { /* copy */ for (g=0; g<G; g++) { idx = g*sizeM; for (m=1; m<sizeM ; m++ ) z[m+idx] = z[idx]; } } if (B>0) { next_idx = 0; for (b=0; b<B; b++) { while (next_idx<M && (timestamp[next_idx]<bins[b])) next_idx++; idx = next_idx; while (idx<M && (timestamp[idx]<bins[b+B])) { if (D>0) { for (g=0; g<G; g++) pout[b+g*B] += z[idx+g*M]; } pout2[b]++; idx++; } if (D==0) { for (g=0; g<G; g++) { pout[b+g*B] = pout2[b] * z[g*sizeM]; } } } mxFree(z); } plhs[0] = out; plhs[1] = out2; }
int main(int argc, char **argv) { bool expert_mode = false; const char *prgname = "abrt"; abrt_init(argv); /* I18n */ setlocale(LC_ALL, ""); #if ENABLE_NLS bindtextdomain(PACKAGE, LOCALEDIR); textdomain(PACKAGE); #endif /* without this the name is set to argv[0] which confuses * desktops which uses the name to find the corresponding .desktop file * trac#180 * * env variable can be used to override the default prgname, so it's the * same as the application which is calling us (trac#303) * * note that g_set_prgname has to be called before gtk_init */ char *env_prgname = getenv("LIBREPORT_PRGNAME"); g_set_prgname(env_prgname ? env_prgname : prgname); gtk_init(&argc, &argv); /* Can't keep these strings/structs static: _() doesn't support that */ const char *program_usage_string = _( "& [-vpdx] [-e EVENT]... [-g GUI_FILE] PROBLEM_DIR\n" "\n" "GUI tool to analyze and report problem saved in specified PROBLEM_DIR" ); enum { OPT_v = 1 << 0, OPT_g = 1 << 1, OPT_p = 1 << 2, OPT_d = 1 << 3, OPT_e = 1 << 4, OPT_x = 1 << 5, }; /* Keep enum above and order of options below in sync! */ struct options program_options[] = { OPT__VERBOSE(&g_verbose), OPT_STRING('g', NULL, &g_glade_file, "FILE", _("Alternate GUI file")), OPT_BOOL( 'p', NULL, NULL, _("Add program names to log")), OPT_BOOL( 'd', "delete", NULL, _("Remove PROBLEM_DIR after reporting")), OPT_LIST( 'e', "event", &g_auto_event_list, "EVENT", _("Run only these events")), OPT_BOOL( 'x', "expert", &expert_mode, _("Expert mode")), OPT_END() }; unsigned opts = parse_opts(argc, argv, program_options, program_usage_string); argv += optind; if (!argv[0] || argv[1]) /* zero or >1 arguments */ show_usage_and_die(program_usage_string, program_options); /* Allow algorithms to add mallocated strings */ for (GList *elem = g_auto_event_list; elem; elem = g_list_next(elem)) elem->data = xstrdup((const char *)elem->data); export_abrt_envvars(opts & OPT_p); g_dump_dir_name = xstrdup(argv[0]); /* load /etc/abrt/events/foo.{conf,xml} stuff and $XDG_CACHE_HOME/abrt/events/foo.conf */ g_event_config_list = load_event_config_data(); load_event_config_data_from_user_storage(g_event_config_list); load_user_settings("report-gtk"); load_workflow_config_data(WORKFLOWS_DIR); /* list of workflows applicable to the currently processed problem */ GHashTable *possible_workflows = load_workflow_config_data_from_list( list_possible_events_glist(g_dump_dir_name, "workflow"), WORKFLOWS_DIR); /* if we have only 1 workflow, we can use the events from it as default */ if (!expert_mode && g_auto_event_list == NULL && g_hash_table_size(possible_workflows) == 1) { GHashTableIter iter; gpointer key, value; g_hash_table_iter_init(&iter, possible_workflows); if (g_hash_table_iter_next(&iter, &key, &value)) { VERB1 log("autoselected workflow: '%s'", (char *)key); g_auto_event_list = wf_get_event_names((workflow_t *)value); } } problem_data_reload_from_dump_dir(); create_assistant(expert_mode); g_custom_logger = &show_error_as_msgbox; update_gui_state_from_problem_data(UPDATE_SELECTED_EVENT); /* Enter main loop */ gtk_main(); if (opts & OPT_d) delete_dump_dir_possibly_using_abrtd(g_dump_dir_name); save_user_settings(); return 0; }
int main (void) { double y, y_expected; int e, e_expected; gsl_ieee_env_setup (); /* Test for expm1 */ y = gsl_expm1 (0.0); y_expected = 0.0; gsl_test_rel (y, y_expected, 1e-15, "gsl_expm1(0.0)"); y = gsl_expm1 (1e-10); y_expected = 1.000000000050000000002e-10; gsl_test_rel (y, y_expected, 1e-15, "gsl_expm1(1e-10)"); y = gsl_expm1 (-1e-10); y_expected = -9.999999999500000000017e-11; gsl_test_rel (y, y_expected, 1e-15, "gsl_expm1(-1e-10)"); y = gsl_expm1 (0.1); y_expected = 0.1051709180756476248117078264902; gsl_test_rel (y, y_expected, 1e-15, "gsl_expm1(0.1)"); y = gsl_expm1 (-0.1); y_expected = -0.09516258196404042683575094055356; gsl_test_rel (y, y_expected, 1e-15, "gsl_expm1(-0.1)"); y = gsl_expm1 (10.0); y_expected = 22025.465794806716516957900645284; gsl_test_rel (y, y_expected, 1e-15, "gsl_expm1(10.0)"); y = gsl_expm1 (-10.0); y_expected = -0.99995460007023751514846440848444; gsl_test_rel (y, y_expected, 1e-15, "gsl_expm1(-10.0)"); /* Test for log1p */ y = gsl_log1p (0.0); y_expected = 0.0; gsl_test_rel (y, y_expected, 1e-15, "gsl_log1p(0.0)"); y = gsl_log1p (1e-10); y_expected = 9.9999999995000000000333333333308e-11; gsl_test_rel (y, y_expected, 1e-15, "gsl_log1p(1e-10)"); y = gsl_log1p (0.1); y_expected = 0.095310179804324860043952123280765; gsl_test_rel (y, y_expected, 1e-15, "gsl_log1p(0.1)"); y = gsl_log1p (10.0); y_expected = 2.3978952727983705440619435779651; gsl_test_rel (y, y_expected, 1e-15, "gsl_log1p(10.0)"); /* Test for gsl_hypot */ y = gsl_hypot (0.0, 0.0); y_expected = 0.0; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot(0.0, 0.0)"); y = gsl_hypot (1e-10, 1e-10); y_expected = 1.414213562373095048801688e-10; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot(1e-10, 1e-10)"); y = gsl_hypot (1e-38, 1e-38); y_expected = 1.414213562373095048801688e-38; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot(1e-38, 1e-38)"); y = gsl_hypot (1e-10, -1.0); y_expected = 1.000000000000000000005; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot(1e-10, -1)"); y = gsl_hypot (-1.0, 1e-10); y_expected = 1.000000000000000000005; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot(-1, 1e-10)"); y = gsl_hypot (1e307, 1e301); y_expected = 1.000000000000499999999999e307; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot(1e307, 1e301)"); y = gsl_hypot (1e301, 1e307); y_expected = 1.000000000000499999999999e307; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot(1e301, 1e307)"); y = gsl_hypot (1e307, 1e307); y_expected = 1.414213562373095048801688e307; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot(1e307, 1e307)"); /* Test for gsl_hypot3 */ y = gsl_hypot3 (0.0, 0.0, 0.0); y_expected = 0.0; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot3(0.0, 0.0, 0.0)"); y = gsl_hypot3 (1e-10, 1e-10, 1e-10); y_expected = 1.732050807568877293527446e-10; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot3(1e-10, 1e-10, 1e-10)"); y = gsl_hypot3 (1e-38, 1e-38, 1e-38); y_expected = 1.732050807568877293527446e-38; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot3(1e-38, 1e-38, 1e-38)"); y = gsl_hypot3 (1e-10, 1e-10, -1.0); y_expected = 1.000000000000000000099; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot3(1e-10, 1e-10, -1)"); y = gsl_hypot3 (1e-10, -1.0, 1e-10); y_expected = 1.000000000000000000099; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot3(1e-10, -1, 1e-10)"); y = gsl_hypot3 (-1.0, 1e-10, 1e-10); y_expected = 1.000000000000000000099; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot3(-1, 1e-10, 1e-10)"); y = gsl_hypot3 (1e307, 1e301, 1e301); y_expected = 1.0000000000009999999999995e307; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot3(1e307, 1e301, 1e301)"); y = gsl_hypot3 (1e307, 1e307, 1e307); y_expected = 1.732050807568877293527446e307; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot3(1e307, 1e307, 1e307)"); y = gsl_hypot3 (1e307, 1e-307, 1e-307); y_expected = 1.0e307; gsl_test_rel (y, y_expected, 1e-15, "gsl_hypot3(1e307, 1e-307, 1e-307)"); /* Test for acosh */ y = gsl_acosh (1.0); y_expected = 0.0; gsl_test_rel (y, y_expected, 1e-15, "gsl_acosh(1.0)"); y = gsl_acosh (1.1); y_expected = 4.435682543851151891329110663525e-1; gsl_test_rel (y, y_expected, 1e-15, "gsl_acosh(1.1)"); y = gsl_acosh (10.0); y_expected = 2.9932228461263808979126677137742e0; gsl_test_rel (y, y_expected, 1e-15, "gsl_acosh(10.0)"); y = gsl_acosh (1e10); y_expected = 2.3718998110500402149594646668302e1; gsl_test_rel (y, y_expected, 1e-15, "gsl_acosh(1e10)"); /* Test for asinh */ y = gsl_asinh (0.0); y_expected = 0.0; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(0.0)"); y = gsl_asinh (1e-10); y_expected = 9.9999999999999999999833333333346e-11; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(1e-10)"); y = gsl_asinh (-1e-10); y_expected = -9.9999999999999999999833333333346e-11; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(1e-10)"); y = gsl_asinh (0.1); y_expected = 9.983407889920756332730312470477e-2; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(0.1)"); y = gsl_asinh (-0.1); y_expected = -9.983407889920756332730312470477e-2; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(-0.1)"); y = gsl_asinh (1.0); y_expected = 8.8137358701954302523260932497979e-1; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(1.0)"); y = gsl_asinh (-1.0); y_expected = -8.8137358701954302523260932497979e-1; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(-1.0)"); y = gsl_asinh (10.0); y_expected = 2.9982229502979697388465955375965e0; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(10)"); y = gsl_asinh (-10.0); y_expected = -2.9982229502979697388465955375965e0; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(-10)"); y = gsl_asinh (1e10); y_expected = 2.3718998110500402149599646668302e1; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(1e10)"); y = gsl_asinh (-1e10); y_expected = -2.3718998110500402149599646668302e1; gsl_test_rel (y, y_expected, 1e-15, "gsl_asinh(-1e10)"); /* Test for atanh */ y = gsl_atanh (0.0); y_expected = 0.0; gsl_test_rel (y, y_expected, 1e-15, "gsl_atanh(0.0)"); y = gsl_atanh (1e-20); y_expected = 1e-20; gsl_test_rel (y, y_expected, 1e-15, "gsl_atanh(1e-20)"); y = gsl_atanh (-1e-20); y_expected = -1e-20; gsl_test_rel (y, y_expected, 1e-15, "gsl_atanh(-1e-20)"); y = gsl_atanh (0.1); y_expected = 1.0033534773107558063572655206004e-1; gsl_test_rel (y, y_expected, 1e-15, "gsl_atanh(0.1)"); y = gsl_atanh (-0.1); y_expected = -1.0033534773107558063572655206004e-1; gsl_test_rel (y, y_expected, 1e-15, "gsl_atanh(-0.1)"); y = gsl_atanh (0.9); y_expected = 1.4722194895832202300045137159439e0; gsl_test_rel (y, y_expected, 1e-15, "gsl_atanh(0.9)"); y = gsl_atanh (-0.9); y_expected = -1.4722194895832202300045137159439e0; gsl_test_rel (y, y_expected, 1e-15, "gsl_atanh(0.9)"); /* Test for pow_int */ y = gsl_pow_2 (-3.14); y_expected = pow (-3.14, 2.0); gsl_test_rel (y, y_expected, 1e-15, "gsl_pow_2(-3.14)"); y = gsl_pow_3 (-3.14); y_expected = pow (-3.14, 3.0); gsl_test_rel (y, y_expected, 1e-15, "gsl_pow_3(-3.14)"); y = gsl_pow_4 (-3.14); y_expected = pow (-3.14, 4.0); gsl_test_rel (y, y_expected, 1e-15, "gsl_pow_4(-3.14)"); y = gsl_pow_5 (-3.14); y_expected = pow (-3.14, 5.0); gsl_test_rel (y, y_expected, 1e-15, "gsl_pow_5(-3.14)"); y = gsl_pow_6 (-3.14); y_expected = pow (-3.14, 6.0); gsl_test_rel (y, y_expected, 1e-15, "gsl_pow_6(-3.14)"); y = gsl_pow_7 (-3.14); y_expected = pow (-3.14, 7.0); gsl_test_rel (y, y_expected, 1e-15, "gsl_pow_7(-3.14)"); y = gsl_pow_8 (-3.14); y_expected = pow (-3.14, 8.0); gsl_test_rel (y, y_expected, 1e-15, "gsl_pow_8(-3.14)"); y = gsl_pow_9 (-3.14); y_expected = pow (-3.14, 9.0); gsl_test_rel (y, y_expected, 1e-15, "gsl_pow_9(-3.14)"); { int n; for (n = -9; n < 10; n++) { y = gsl_pow_int (-3.14, n); y_expected = pow (-3.14, n); gsl_test_rel (y, y_expected, 1e-15, "gsl_pow_n(-3.14,%d)", n); } } /* Test for ldexp */ y = gsl_ldexp (M_PI, -2); y_expected = M_PI_4; gsl_test_rel (y, y_expected, 1e-15, "gsl_ldexp(pi,-2)"); y = gsl_ldexp (1.0, 2); y_expected = 4.000000; gsl_test_rel (y, y_expected, 1e-15, "gsl_ldexp(1.0,2)"); y = gsl_ldexp (0.0, 2); y_expected = 0.0; gsl_test_rel (y, y_expected, 1e-15, "gsl_ldexp(0.0,2)"); y = gsl_ldexp (9.999999999999998890e-01, 1024); y_expected = GSL_DBL_MAX; gsl_test_rel (y, y_expected, 1e-15, "gsl_ldexp DBL_MAX"); y = gsl_ldexp (1e308, -2000); y_expected = 8.7098098162172166755761e-295; gsl_test_rel (y, y_expected, 1e-15, "gsl_ldexp(1e308,-2000)"); y = gsl_ldexp (GSL_DBL_MIN, 2000); y_expected = 2.554675596204441378334779940e294; gsl_test_rel (y, y_expected, 1e-15, "gsl_ldexp(DBL_MIN,2000)"); /* Test subnormals */ { int i = 0; volatile double x = GSL_DBL_MIN; y_expected = 2.554675596204441378334779940e294; x /= 2; while (x > 0) { i++ ; y = gsl_ldexp (x, 2000 + i); gsl_test_rel (y, y_expected, 1e-15, "gsl_ldexp(DBL_MIN/2**%d,%d)",i,2000+i); x /= 2; } } /* Test for frexp */ y = gsl_frexp (0.0, &e); y_expected = 0; e_expected = 0; gsl_test_rel (y, y_expected, 1e-15, "gsl_frexp(0) fraction"); gsl_test_int (e, e_expected, "gsl_frexp(0) exponent"); y = gsl_frexp (M_PI, &e); y_expected = M_PI_4; e_expected = 2; gsl_test_rel (y, y_expected, 1e-15, "gsl_frexp(pi) fraction"); gsl_test_int (e, e_expected, "gsl_frexp(pi) exponent"); y = gsl_frexp (2.0, &e); y_expected = 0.5; e_expected = 2; gsl_test_rel (y, y_expected, 1e-15, "gsl_frexp(2.0) fraction"); gsl_test_int (e, e_expected, "gsl_frexp(2.0) exponent"); y = gsl_frexp (1.0 / 4.0, &e); y_expected = 0.5; e_expected = -1; gsl_test_rel (y, y_expected, 1e-15, "gsl_frexp(0.25) fraction"); gsl_test_int (e, e_expected, "gsl_frexp(0.25) exponent"); y = gsl_frexp (1.0 / 4.0 - 4.0 * GSL_DBL_EPSILON, &e); y_expected = 0.999999999999996447; e_expected = -2; gsl_test_rel (y, y_expected, 1e-15, "gsl_frexp(0.25-eps) fraction"); gsl_test_int (e, e_expected, "gsl_frexp(0.25-eps) exponent"); y = gsl_frexp (GSL_DBL_MAX, &e); y_expected = 9.999999999999998890e-01; e_expected = 1024; gsl_test_rel (y, y_expected, 1e-15, "gsl_frexp(DBL_MAX) fraction"); gsl_test_int (e, e_expected, "gsl_frexp(DBL_MAX) exponent"); y = gsl_frexp (-GSL_DBL_MAX, &e); y_expected = -9.999999999999998890e-01; e_expected = 1024; gsl_test_rel (y, y_expected, 1e-15, "gsl_frexp(-DBL_MAX) fraction"); gsl_test_int (e, e_expected, "gsl_frexp(-DBL_MAX) exponent"); y = gsl_frexp (GSL_DBL_MIN, &e); y_expected = 0.5; e_expected = -1021; gsl_test_rel (y, y_expected, 1e-15, "gsl_frexp(DBL_MIN) fraction"); gsl_test_int (e, e_expected, "gsl_frexp(DBL_MIN) exponent"); y = gsl_frexp (-GSL_DBL_MIN, &e); y_expected = -0.5; e_expected = -1021; gsl_test_rel (y, y_expected, 1e-15, "gsl_frexp(-DBL_MIN) fraction"); gsl_test_int (e, e_expected, "gsl_frexp(-DBL_MIN) exponent"); /* Test subnormals */ { int i = 0; volatile double x = GSL_DBL_MIN; y_expected = 0.5; e_expected = -1021; x /= 2; while (x > 0) { e_expected--; i++ ; y = gsl_frexp (x, &e); gsl_test_rel (y, y_expected, 1e-15, "gsl_frexp(DBL_MIN/2**%d) fraction",i); gsl_test_int (e, e_expected, "gsl_frexp(DBL_MIN/2**%d) exponent", i); x /= 2; } } /* Test for approximate floating point comparison */ { double x, y; int i; x = M_PI; y = 22.0 / 7.0; /* test the basic function */ for (i = 0; i < 10; i++) { double tol = pow (10, -i); int res = gsl_fcmp (x, y, tol); gsl_test_int (res, -(i >= 4), "gsl_fcmp(%.5f,%.5f,%g)", x, y, tol); } for (i = 0; i < 10; i++) { double tol = pow (10, -i); int res = gsl_fcmp (y, x, tol); gsl_test_int (res, (i >= 4), "gsl_fcmp(%.5f,%.5f,%g)", y, x, tol); } } #if HAVE_IEEE_COMPARISONS /* Test for isinf, isnan, finite */ { double zero, one, inf, nan; int s; zero = 0.0; one = 1.0; inf = exp (1.0e10); nan = inf / inf; s = gsl_isinf (zero); gsl_test_int (s, 0, "gsl_isinf(0)"); s = gsl_isinf (one); gsl_test_int (s, 0, "gsl_isinf(1)"); s = gsl_isinf (inf); gsl_test_int (s, 1, "gsl_isinf(inf)"); s = gsl_isinf (-inf); gsl_test_int (s, -1, "gsl_isinf(-inf)"); s = gsl_isinf (nan); gsl_test_int (s, 0, "gsl_isinf(nan)"); s = gsl_isnan (zero); gsl_test_int (s, 0, "gsl_isnan(0)"); s = gsl_isnan (one); gsl_test_int (s, 0, "gsl_isnan(1)"); s = gsl_isnan (inf); gsl_test_int (s, 0, "gsl_isnan(inf)"); s = gsl_isnan (-inf); gsl_test_int (s, 0, "gsl_isnan(-inf)"); s = gsl_isnan (nan); gsl_test_int (s, 1, "gsl_isnan(nan)"); s = gsl_finite (zero); gsl_test_int (s, 1, "gsl_finite(0)"); s = gsl_finite (one); gsl_test_int (s, 1, "gsl_finite(1)"); s = gsl_finite (inf); gsl_test_int (s, 0, "gsl_finite(inf)"); s = gsl_finite (-inf); gsl_test_int (s, 0, "gsl_finite(-inf)"); s = gsl_finite (nan); gsl_test_int (s, 0, "gsl_finite(nan)"); } #endif { double x = gsl_fdiv (2.0, 3.0); gsl_test_rel (x, 2.0 / 3.0, 4 * GSL_DBL_EPSILON, "gsl_fdiv(2,3)"); } /* Test constants in gsl_math.h */ { double x = log(M_E); gsl_test_rel (x, 1.0, 4 * GSL_DBL_EPSILON, "ln(M_E)"); } { double x=pow(2.0,M_LOG2E); gsl_test_rel (x, exp(1.0), 4 * GSL_DBL_EPSILON, "2^M_LOG2E"); } { double x=pow(10.0,M_LOG10E); gsl_test_rel (x, exp(1.0), 4 * GSL_DBL_EPSILON, "10^M_LOG10E"); } { double x=pow(M_SQRT2, 2.0); gsl_test_rel (x, 2.0, 4 * GSL_DBL_EPSILON, "M_SQRT2^2"); } { double x=pow(M_SQRT1_2, 2.0); gsl_test_rel (x, 1.0/2.0, 4 * GSL_DBL_EPSILON, "M_SQRT1_2"); } { double x=pow(M_SQRT3, 2.0); gsl_test_rel (x, 3.0, 4 * GSL_DBL_EPSILON, "M_SQRT3^2"); } { double x = M_PI; gsl_test_rel (x, 3.1415926535897932384626433832795, 4 * GSL_DBL_EPSILON, "M_PI"); } { double x = 2 * M_PI_2; gsl_test_rel (x, M_PI, 4 * GSL_DBL_EPSILON, "2*M_PI_2"); } { double x = 4 * M_PI_4; gsl_test_rel (x, M_PI, 4 * GSL_DBL_EPSILON, "4*M_PI_4"); } { double x = pow(M_SQRTPI, 2.0); gsl_test_rel (x, M_PI, 4 * GSL_DBL_EPSILON, "M_SQRTPI^2"); } { double x = pow(M_2_SQRTPI, 2.0); gsl_test_rel (x, 4/M_PI, 4 * GSL_DBL_EPSILON, "M_SQRTPI^2"); } { double x = M_1_PI; gsl_test_rel (x, 1/M_PI, 4 * GSL_DBL_EPSILON, "M_1_SQRTPI"); } { double x = M_2_PI; gsl_test_rel (x, 2.0/M_PI, 4 * GSL_DBL_EPSILON, "M_2_PI"); } { double x = exp(M_LN10); gsl_test_rel (x, 10, 4 * GSL_DBL_EPSILON, "exp(M_LN10)"); } { double x = exp(M_LN2); gsl_test_rel (x, 2, 4 * GSL_DBL_EPSILON, "exp(M_LN2)"); } { double x = exp(M_LNPI); gsl_test_rel (x, M_PI, 4 * GSL_DBL_EPSILON, "exp(M_LNPI)"); } { double x = M_EULER; gsl_test_rel (x, 0.5772156649015328606065120900824, 4 * GSL_DBL_EPSILON, "M_EULER"); } exit (gsl_test_summary ()); }
void InPlaceReprojection::initialize() { Filter::initialize(); checkImpedance(); if (m_inferInputSRS) { m_inSRS = getPrevStage().getSpatialReference(); } #ifdef PDAL_HAVE_GDAL m_gdal_debug = boost::shared_ptr<pdal::gdal::Debug>( new pdal::gdal::Debug(isDebug(), log())); m_in_ref_ptr = ReferencePtr(OSRNewSpatialReference(0), OGRSpatialReferenceDeleter()); m_out_ref_ptr = ReferencePtr(OSRNewSpatialReference(0), OGRSpatialReferenceDeleter()); int result = OSRSetFromUserInput(m_in_ref_ptr.get(), m_inSRS.getWKT(pdal::SpatialReference::eCompoundOK).c_str()); if (result != OGRERR_NONE) { std::ostringstream msg; msg << "Could not import input spatial reference for InPlaceReprojection:: " << CPLGetLastErrorMsg() << " code: " << result << " wkt: '" << m_inSRS.getWKT() << "'"; throw std::runtime_error(msg.str()); } result = OSRSetFromUserInput(m_out_ref_ptr.get(), m_outSRS.getWKT(pdal::SpatialReference::eCompoundOK).c_str()); if (result != OGRERR_NONE) { std::ostringstream msg; msg << "Could not import output spatial reference for InPlaceReprojection:: " << CPLGetLastErrorMsg() << " code: " << result << " wkt: '" << m_outSRS.getWKT() << "'"; std::string message(msg.str()); throw std::runtime_error(message); } m_transform_ptr = TransformPtr(OCTNewCoordinateTransformation( m_in_ref_ptr.get(), m_out_ref_ptr.get()), OSRTransformDeleter()); if (!m_transform_ptr.get()) { std::ostringstream msg; msg << "Could not construct CoordinateTransformation in InPlaceReprojection:: "; std::string message(msg.str()); throw std::runtime_error(message); } #endif setSpatialReference(m_outSRS); updateBounds(); return; }
main(){ struct timeval startTime, endTime; long seconds, useconds; double mtime; double *L, *D, *U, *W, *x; double aA, bB; double ps0; double tau_j; double Dx, Dt; int i, j, k, ir; double sig=0.2; double T=0.25; double S0=1950.0, K=1975.0; double rfr=0.05; double div=0.025; double sMin = 100.0, sMax = 4100.0; int N=4000; int M=1000; Dx = (log(sMax)-log(sMin))/N; Dt = T/M; L = new double[N]; D = new double[N]; U = new double[N]; W = new double[N+1]; x = new double[N+1]; aA = sig*sig*Dt/(2.0*Dx*Dx); bB = (rfr-div-0.5*sig*sig)*Dt/(2.0*Dx); gettimeofday(&startTime, NULL); // payoff for ( i = 0; i <= N; ++i){ x[i] = log(sMin)+i*Dx; if(exp(x[i])<K) W[i] = K-exp(x[i]); else W[i] = 0.0; } for (j = M-1; j >= 0; --j) { //cout << j << endl; for (i = 1; i <=N-1; ++i){ if (i==1){ //L[i] = 0; D[i] = 1 + 2.0*aA + rfr*Dt; // - (aA-bB)* 2.0/(1.0+Dx/2.0); U[i] = -(aA+bB); // + (aA-bB)*(1.0-Dx/2)/(1.0+Dx/2.0); // Boundary condition W[1] = W[1]+(aA-bB)*(K-exp(x[0])); //W[1] = W[1]+(aA-bB)*(K*exp(-rfr*tau_j)-exp(x[0])*exp(-div*tau_j)); } else if (i==N-1){ L[i] = -(aA-bB); // + (aA+bB)*(1+Dx/2.0)/(1.0-Dx/2.0); D[i] = 1.0 + 2.0*aA + rfr*Dt; // - (aA+bB)* 2.0/(1.0-Dx/2.0); //U[i] = 0; // Boundary condition //W[N-1] = W[N-1]+0.0; } else{ L[i] = -(aA-bB); D[i] = 1 + 2.0*aA + rfr*Dt; U[i] = -(aA+bB); } } // A w_{k+1} = w_{k} tridiagSolver(L, D, U, W, N-1); } // for(i = 0; i <= N-1; ++i){ if(x[i]>log(S0)){ ir=i; break; } } ps0 = W[ir-1] + ((W[ir]-W[ir-1])/Dx)*(log(S0)-x[ir-1]); cout << "Stock Price " << S0 << endl; cout << "Put Value " << ps0 << "\n" << endl; delete [] L; delete [] D; delete [] U; delete [] W; delete [] x; gettimeofday(&endTime, NULL); seconds = endTime.tv_sec - startTime.tv_sec; useconds = endTime.tv_usec - startTime.tv_usec; mtime = ((seconds) * 1000 + useconds/1000.0); cout << "Time elapsed was: " << mtime << " (milliseconds)" << endl; return 0; }
int stream_component_open(VideoState *is, int stream_index) { AVFormatContext *pFormatCtx = is->pFormatCtx; AVCodecContext *codecCtx = NULL; AVCodec *codec = NULL; AVDictionary *optionsDict = NULL; if(stream_index < 0 || stream_index >= pFormatCtx->nb_streams) { return -1; } // Get a pointer to the codec context for the video stream codecCtx = pFormatCtx->streams[stream_index]->codec; if(codecCtx->codec_type == AVMEDIA_TYPE_AUDIO) { is->audio_callback = audio_callback; // Set audio settings from codec info AudioPlayer *player = malloc(sizeof(AudioPlayer)); is->audio_player = player; createEngine(&is->audio_player); createBufferQueueAudioPlayer(&is->audio_player, is, codecCtx->channels, codecCtx->sample_rate); //is->audio_hw_buf_size = 4096; } else if (codecCtx->codec_type == AVMEDIA_TYPE_VIDEO) { // Set video settings from codec info VideoPlayer *player = malloc(sizeof(VideoPlayer)); is->video_player = player; createVideoEngine(&is->video_player); } codec = avcodec_find_decoder(codecCtx->codec_id); if(!codec || (avcodec_open2(codecCtx, codec, &optionsDict) < 0)) { fprintf(stderr, "Unsupported codec!\n"); return -1; } switch(codecCtx->codec_type) { case AVMEDIA_TYPE_AUDIO: is->audioStream = stream_index; is->audio_st = pFormatCtx->streams[stream_index]; is->audio_buf_size = 0; is->audio_buf_index = 0; /* averaging filter for audio sync */ is->audio_diff_avg_coef = exp(log(0.01 / AUDIO_DIFF_AVG_NB)); is->audio_diff_avg_count = 0; /* Correct audio only if larger error than this */ is->audio_diff_threshold = 2.0 * SDL_AUDIO_BUFFER_SIZE / codecCtx->sample_rate; is->sws_ctx_audio = swr_alloc(); if (!is->sws_ctx_audio) { fprintf(stderr, "Could not allocate resampler context\n"); return -1; } memset(&is->audio_pkt, 0, sizeof(is->audio_pkt)); packet_queue_init(&is->audioq); break; case AVMEDIA_TYPE_VIDEO: is->videoStream = stream_index; is->video_st = pFormatCtx->streams[stream_index]; is->frame_timer = (double)av_gettime() / 1000000.0; is->frame_last_delay = 40e-3; is->video_current_pts_time = av_gettime(); packet_queue_init(&is->videoq); is->video_tid = malloc(sizeof(*(is->video_tid))); // uncomment for video pthread_create(is->video_tid, NULL, (void *) &video_thread, is); is->sws_ctx = sws_getContext ( is->video_st->codec->width, is->video_st->codec->height, is->video_st->codec->pix_fmt, is->video_st->codec->width, is->video_st->codec->height, AV_PIX_FMT_YUV420P, SWS_BILINEAR, NULL, NULL, NULL ); /*is->sws_ctx = sws_getContext ( is->video_st->codec->width, is->video_st->codec->height, is->video_st->codec->pix_fmt, //pCodecCtx->width, //pCodecCtx->height, is->video_st->codec->width, is->video_st->codec->height, TARGET_IMAGE_FORMAT, SWS_BILINEAR, NULL, NULL, NULL );*/ codecCtx->get_buffer2 = our_get_buffer; break; default: break; } return 0; }
int32_t CharsetRecog_mbcs::match_mbcs(InputText *det, const uint16_t commonChars[], int32_t commonCharsLen) const { int32_t singleByteCharCount = 0; int32_t doubleByteCharCount = 0; int32_t commonCharCount = 0; int32_t badCharCount = 0; int32_t totalCharCount = 0; int32_t confidence = 0; IteratedChar iter; while (nextChar(&iter, det)) { totalCharCount++; if (iter.error) { badCharCount++; } else { if (iter.charValue <= 0xFF) { singleByteCharCount++; } else { doubleByteCharCount++; if (commonChars != 0) { if (binarySearch(commonChars, commonCharsLen, iter.charValue) >= 0){ commonCharCount += 1; } } } } if (badCharCount >= 2 && badCharCount*5 >= doubleByteCharCount) { // Bail out early if the byte data is not matching the encoding scheme. // break detectBlock; return confidence; } } if (doubleByteCharCount <= 10 && badCharCount == 0) { // Not many multi-byte chars. if (doubleByteCharCount == 0 && totalCharCount < 10) { // There weren't any multibyte sequences, and there was a low density of non-ASCII single bytes. // We don't have enough data to have any confidence. // Statistical analysis of single byte non-ASCII charcters would probably help here. confidence = 0; } else { // ASCII or ISO file? It's probably not our encoding, // but is not incompatible with our encoding, so don't give it a zero. confidence = 10; } return confidence; } // // No match if there are too many characters that don't fit the encoding scheme. // (should we have zero tolerance for these?) // if (doubleByteCharCount < 20*badCharCount) { confidence = 0; return confidence; } if (commonChars == 0) { // We have no statistics on frequently occuring characters. // Assess confidence purely on having a reasonable number of // multi-byte characters (the more the better) confidence = 30 + doubleByteCharCount - 20*badCharCount; if (confidence > 100) { confidence = 100; } } else { // // Frequency of occurence statistics exist. // double maxVal = log((double)doubleByteCharCount / 4); /*(float)?*/ double scaleFactor = 90.0 / maxVal; confidence = (int32_t)(log((double)commonCharCount+1) * scaleFactor + 10.0); confidence = min(confidence, 100); } if (confidence < 0) { confidence = 0; } return confidence; }
binary_relation_detector binary_relation_detector_trainer:: train ( ) const { DLIB_CASSERT(num_positive_examples() > 0, "Not enough training data given."); DLIB_CASSERT(num_negative_examples() > 0, "Not enough training data given."); std::vector<sparse_vector_type> samples; std::vector<double> labels; for (unsigned long i = 0; i < pos_sentences.size(); ++i) { samples.push_back(extract_binary_relation(pos_sentences[i], pos_arg1s[i], pos_arg2s[i], tfe).feats); labels.push_back(+1); } for (unsigned long i = 0; i < neg_sentences.size(); ++i) { samples.push_back(extract_binary_relation(neg_sentences[i], neg_arg1s[i], neg_arg2s[i], tfe).feats); labels.push_back(-1); } randomize_samples(samples, labels); const int cv_folds = 6; brdt_cv_objective obj(num_threads, cv_folds, beta, samples, labels); matrix<double,2,1> params; params = 5000.0/samples.size(), 5000.0/samples.size(); // We do the parameter search in log space. params = log(params); // can't do the parameter search if we don't have enough data. So if we don't // have much data then just use the default parameters. if (pos_sentences.size() > (unsigned)cv_folds) { matrix<double,2,1> lower_params, upper_params; lower_params = 1.0/samples.size(), 1.0/samples.size(); upper_params = 100000.0/samples.size(), 100000.0/samples.size(); lower_params = log(lower_params); upper_params = log(upper_params); const double rho_begin = min(upper_params-lower_params)*0.15; const double rho_end = log(1.2/samples.size()) - log(1.0/samples.size()); find_max_bobyqa(obj, params, params.size()*2+1, lower_params, upper_params, rho_begin, rho_end, 200); } // Note that we rescale the parameters to account for the fact that the cross // validation was done on a dataset slightly smaller than the one we ultimately train // on and the C parameters of this trainer are not normalized by the number of training // samples. params = exp(params) * (cv_folds-1.0)/cv_folds; svm_c_linear_dcd_trainer<sparse_linear_kernel<sparse_vector_type> > trainer; trainer.set_c_class1(params(0)); trainer.set_c_class2(params(1)); cout << "using parameters of: " << trans(params); cout << "now doing training..." << endl; binary_relation_detector bd; bd.df = trainer.train(samples, labels); bd.relation_type = relation_name; bd.total_word_feature_extractor_fingerprint = tfe.get_fingerprint(); cout << "test on train: " << test_binary_decision_function(bd.df, samples, labels) << endl; return bd; }
void PvaPyLogger::error(const char* message, va_list messageArgs) const { log(LogLevelError, message, messageArgs); }
void PvaPyLogger::error(const std::string& message) const { log(LogLevelError, message.c_str()); }
int Logger::error(const char* message) { return log(hiPriorityBuffer(), SLOG2_ERROR, message); }
Status LegacyReplicationCoordinator::processHeartbeat(const BSONObj& cmdObj, BSONObjBuilder* resultObj) { if( cmdObj["pv"].Int() != 1 ) { return Status(ErrorCodes::BadValue, "incompatible replset protocol version"); } { string s = string(cmdObj.getStringField("replSetHeartbeat")); if (replSettings.ourSetName() != s) { log() << "replSet set names do not match, our cmdline: " << replSettings.replSet << rsLog; log() << "replSet s: " << s << rsLog; resultObj->append("mismatch", true); return Status(ErrorCodes::BadValue, "repl set names do not match"); } } resultObj->append("rs", true); if( (theReplSet == 0) || (theReplSet->startupStatus == ReplSetImpl::LOADINGCONFIG) ) { string from( cmdObj.getStringField("from") ); if( !from.empty() ) { scoped_lock lck( replSettings.discoveredSeeds_mx ); replSettings.discoveredSeeds.insert(from); } resultObj->append("hbmsg", "still initializing"); return Status::OK(); } if( theReplSet->name() != cmdObj.getStringField("replSetHeartbeat") ) { resultObj->append("mismatch", true); return Status(ErrorCodes::BadValue, "repl set names do not match (2)"); } resultObj->append("set", theReplSet->name()); MemberState currentState = theReplSet->state(); resultObj->append("state", currentState.s); if (currentState == MemberState::RS_PRIMARY) { resultObj->appendDate("electionTime", theReplSet->getElectionTime().asDate()); } resultObj->append("e", theReplSet->iAmElectable()); resultObj->append("hbmsg", theReplSet->hbmsg()); resultObj->append("time", (long long) time(0)); resultObj->appendDate("opTime", theReplSet->lastOpTimeWritten.asDate()); const Member *syncTarget = BackgroundSync::get()->getSyncTarget(); if (syncTarget) { resultObj->append("syncingTo", syncTarget->fullName()); } int v = theReplSet->config().version; resultObj->append("v", v); if( v > cmdObj["v"].Int() ) *resultObj << "config" << theReplSet->config().asBson(); Member* from = NULL; if (cmdObj.hasField("fromId")) { if (v == cmdObj["v"].Int()) { from = theReplSet->getMutableMember(cmdObj["fromId"].Int()); } } if (!from) { from = theReplSet->findByName(cmdObj.getStringField("from")); if (!from) { return Status::OK(); } } // if we thought that this node is down, let it know if (!from->hbinfo().up()) { resultObj->append("stateDisagreement", true); } // note that we got a heartbeat from this node theReplSet->mgr->send(stdx::bind(&ReplSet::msgUpdateHBRecv, theReplSet, from->hbinfo().id(), time(0))); return Status::OK(); }
int Logger::critical(const char* message) { return log(hiPriorityBuffer(), SLOG2_CRITICAL, message); }
Status LegacyReplicationCoordinator::processReplSetReconfig(OperationContext* txn, const BSONObj& newConfigObj, bool force, BSONObjBuilder* resultObj) { if( force && !theReplSet ) { replSettings.reconfig = newConfigObj.getOwned(); resultObj->append("msg", "will try this config momentarily, try running rs.conf() again in a " "few seconds"); return Status::OK(); } Status status = _checkReplEnabledForCommand(resultObj); if (!status.isOK()) { return status; } if( !force && !theReplSet->box.getState().primary() ) { return Status(ErrorCodes::NotMaster, "replSetReconfig command must be sent to the current replica set " "primary."); } try { { // just make sure we can get a write lock before doing anything else. we'll // reacquire one later. of course it could be stuck then, but this check lowers the // risk if weird things are up - we probably don't want a change to apply 30 minutes // after the initial attempt. time_t t = time(0); Lock::GlobalWrite lk(txn->lockState()); if( time(0)-t > 20 ) { return Status(ErrorCodes::ExceededTimeLimit, "took a long time to get write lock, so not initiating. " "Initiate when server less busy?"); } } scoped_ptr<ReplSetConfig> newConfig(ReplSetConfig::make(newConfigObj, force)); log() << "replSet replSetReconfig config object parses ok, " << newConfig->members.size() << " members specified" << rsLog; Status status = ReplSetConfig::legalChange(theReplSet->getConfig(), *newConfig); if (!status.isOK()) { return status; } checkMembersUpForConfigChange(*newConfig, *resultObj, false); log() << "replSet replSetReconfig [2]" << rsLog; theReplSet->haveNewConfig(*newConfig, true); ReplSet::startupStatusMsg.set("replSetReconfig'd"); } catch(const DBException& e) { log() << "replSet replSetReconfig exception: " << e.what() << rsLog; return e.toStatus(); } resetSlaveCache(); return Status::OK(); }
int piso_log2 (int a, int p) { if (p <= MAX_INT * log (2) / log (a)) return (int) p * log (a) / log (2); else return OVF; }
Status LegacyReplicationCoordinator::processReplSetInitiate(OperationContext* txn, const BSONObj& givenConfig, BSONObjBuilder* resultObj) { log() << "replSet replSetInitiate admin command received from client" << rsLog; if( !replSet ) { return Status(ErrorCodes::NoReplicationEnabled, "server is not running with --replSet"); } if( theReplSet ) { resultObj->append("info", "try querying " + rsConfigNs + " to see current configuration"); return Status(ErrorCodes::AlreadyInitialized, "already initialized"); } try { { // just make sure we can get a write lock before doing anything else. we'll // reacquire one later. of course it could be stuck then, but this check lowers the // risk if weird things are up. time_t t = time(0); Lock::GlobalWrite lk(txn->lockState()); if( time(0)-t > 10 ) { return Status(ErrorCodes::ExceededTimeLimit, "took a long time to get write lock, so not initiating. " "Initiate when server less busy?"); } /* check that we don't already have an oplog. that could cause issues. it is ok if the initiating member has *other* data than that. */ BSONObj o; if( Helpers::getFirst(txn, rsoplog, o) ) { return Status(ErrorCodes::AlreadyInitialized, rsoplog + string(" is not empty on the initiating member. " "cannot initiate.")); } } if( ReplSet::startupStatus == ReplSet::BADCONFIG ) { resultObj->append("info", ReplSet::startupStatusMsg.get()); return Status(ErrorCodes::InvalidReplicaSetConfig, "server already in BADCONFIG state (check logs); not initiating"); } if( ReplSet::startupStatus != ReplSet::EMPTYCONFIG ) { resultObj->append("startupStatus", ReplSet::startupStatus); resultObj->append("info", replSettings.replSet); return Status(ErrorCodes::InvalidReplicaSetConfig, "all members and seeds must be reachable to initiate set"); } BSONObj configObj; if (!givenConfig.isEmpty()) { configObj = givenConfig; } else { resultObj->append("info2", "no configuration explicitly specified -- making one"); log() << "replSet info initiate : no configuration specified. " "Using a default configuration for the set" << rsLog; string name; vector<HostAndPort> seeds; set<HostAndPort> seedSet; parseReplSetSeedList(replSettings.replSet, name, seeds, seedSet); // may throw... BSONObjBuilder b; b.append("_id", name); BSONObjBuilder members; members.append("0", BSON( "_id" << 0 << "host" << HostAndPort::me().toString() )); resultObj->append("me", HostAndPort::me().toString()); for( unsigned i = 0; i < seeds.size(); i++ ) { members.append(BSONObjBuilder::numStr(i+1), BSON( "_id" << i+1 << "host" << seeds[i].toString())); } b.appendArray("members", members.obj()); configObj = b.obj(); log() << "replSet created this configuration for initiation : " << configObj.toString() << rsLog; } scoped_ptr<ReplSetConfig> newConfig; try { newConfig.reset(ReplSetConfig::make(configObj)); } catch (const DBException& e) { log() << "replSet replSetInitiate exception: " << e.what() << rsLog; return Status(ErrorCodes::InvalidReplicaSetConfig, mongoutils::str::stream() << "couldn't parse cfg object " << e.what()); } if( newConfig->version > 1 ) { return Status(ErrorCodes::InvalidReplicaSetConfig, "can't initiate with a version number greater than 1"); } log() << "replSet replSetInitiate config object parses ok, " << newConfig->members.size() << " members specified" << rsLog; checkMembersUpForConfigChange(*newConfig, *resultObj, true); log() << "replSet replSetInitiate all members seem up" << rsLog; createOplog(); Lock::GlobalWrite lk(txn->lockState()); BSONObj comment = BSON( "msg" << "initiating set"); newConfig->saveConfigLocally(comment); log() << "replSet replSetInitiate config now saved locally. " "Should come online in about a minute." << rsLog; resultObj->append("info", "Config now saved locally. Should come online in about a minute."); ReplSet::startupStatus = ReplSet::SOON; ReplSet::startupStatusMsg.set("Received replSetInitiate - " "should come online shortly."); } catch(const DBException& e ) { return e.toStatus(); } return Status::OK(); }
static void input_userauth_request(int type, u_int32_t seq, void *ctxt) { Authctxt *authctxt = ctxt; Authmethod *m = NULL; char *user, *service, *method, *style = NULL; if (authctxt == NULL) fatal("input_userauth_request: no authctxt"); user = packet_get_string(NULL); service = packet_get_string(NULL); method = packet_get_string(NULL); debug("userauth-request for user %s service %s method %s", user, service, method); debug("attempt %d initial attempt %d failures %d initial failures %d", authctxt->attempt, authctxt->init_attempt, authctxt->failures, authctxt->init_failures); m = authmethod_lookup(method); if ((style = strchr(user, ':')) != NULL) *style++ = 0; authctxt->attempt++; if (m != NULL && m->is_initial) authctxt->init_attempt++; if (authctxt->attempt == 1) { /* setup auth context */ authctxt->pw = PRIVSEP(getpwnamallow(user)); /* May want to abstract SSHv2 services someday */ if (authctxt->pw && strcmp(service, "ssh-connection")==0) { /* enforced in userauth_finish() below */ authctxt->valid = 1; debug2("input_userauth_request: setting up authctxt for %s", user); } else { log("input_userauth_request: illegal user %s", user); } setproctitle("%s%s", authctxt->pw ? user : "******", use_privsep ? " [net]" : ""); authctxt->user = xstrdup(user); authctxt->service = xstrdup(service); authctxt->style = style ? xstrdup(style) : NULL; userauth_reset_methods(); if (use_privsep) mm_inform_authserv(service, style); } else { char *abandoned; /* * Check for abandoned [multi-round-trip] userauths * methods (e.g., kbdint). Userauth method abandonment * should be treated as userauth method failure and * counted against max_auth_tries. */ abandoned = authmethods_check_abandonment(authctxt, m); if (abandoned != NULL && authctxt->failures > options.max_auth_tries) { /* userauth_finish() will now packet_disconnect() */ userauth_finish(authctxt, abandoned); /* NOTREACHED */ } /* Handle user|service changes, possibly packet_disconnect() */ userauth_user_svc_change(authctxt, user, service); } authctxt->method = m; /* run userauth method, try to authenticate user */ if (m != NULL && userauth_method_can_run(m)) { debug2("input_userauth_request: try method %s", method); m->postponed = 0; m->abandoned = 0; m->authenticated = 0; if (!m->is_initial || authctxt->init_failures < options.max_init_auth_tries) m->userauth(authctxt); authmethod_count_attempt(m); if (authctxt->unwind_dispatch_loop) { /* * Method ran nested dispatch loop but was * abandoned. Cleanup and return without doing * anything else; we're just unwinding the stack. */ authctxt->unwind_dispatch_loop = 0; goto done; } if (m->postponed) goto done; /* multi-round trip userauth not finished */ if (m->abandoned) { /* multi-round trip userauth abandoned, log failure */ auth_log(authctxt, 0, method, " ssh2"); goto done; } } userauth_finish(authctxt, method); done: xfree(service); xfree(user); xfree(method); }
int Logger::debug(const char* message) { return log(lowPriorityBuffer(), SLOG2_DEBUG1, message); }
void adjustdata() /* **---------------------------------------------------------------------- ** Input: none ** Output: none ** Purpose: adjusts data after input file has been processed **---------------------------------------------------------------------- */ { int i; double ucf; /* Unit conversion factor */ Pdemand demand; /* Pointer to demand record */ /* Use 1 hr pattern & report time step if none specified */ if (Pstep <= 0) Pstep = 3600; if (Rstep == 0) Rstep = Pstep; /* Hydraulic time step cannot be greater than pattern or report time step */ if (Hstep <= 0) Hstep = 3600; if (Hstep > Pstep) Hstep = Pstep; if (Hstep > Rstep) Hstep = Rstep; /* Report start time cannot be greater than simulation duration */ if (Rstart > Dur) Rstart = 0; /* No water quality analysis for steady state run */ if (Dur == 0) Qualflag = NONE; /* If no quality timestep, then make it 1/10 of hydraulic timestep */ if (Qstep == 0) Qstep = Hstep/10; /* If no rule time step then make it 1/10 of hydraulic time step; */ /* Rule time step cannot be greater than hydraulic time step */ if (Rulestep == 0) Rulestep = Hstep/10; Rulestep = MIN(Rulestep, Hstep); /* Quality timestep cannot exceed hydraulic timestep */ Qstep = MIN(Qstep, Hstep); /* If no quality tolerance, then use default values */ if (Ctol == MISSING) { if (Qualflag == AGE) Ctol = AGETOL; else Ctol = CHEMTOL; } /* Determine unit system based on flow units */ switch (Flowflag) { case LPS: /* Liters/sec */ case LPM: /* Liters/min */ case MLD: /* megaliters/day */ case CMH: /* cubic meters/hr */ case CMD: /* cubic meters/day */ Unitsflag = SI; break; default: Unitsflag = US; } /* Revise pressure units depending on flow units */ if (Unitsflag != SI) Pressflag = PSI; else if (Pressflag == PSI) Pressflag = METERS; /* Store value of viscosity & diffusivity */ ucf = 1.0; if (Unitsflag == SI) ucf = SQR(MperFT); if (Viscos == MISSING) /* No value supplied */ Viscos = VISCOS; else if (Viscos > 1.e-3) /* Multiplier supplied */ Viscos = Viscos*VISCOS; else /* Actual value supplied */ Viscos = Viscos/ucf; if (Diffus == MISSING) Diffus = DIFFUS; else if (Diffus > 1.e-4) Diffus = Diffus*DIFFUS; else Diffus = Diffus/ucf; /* Set exponent in head loss equation and adjust flow-resistance tolerance. */ if (Formflag == HW) Hexp = 1.852; else Hexp = 2.0; /*** Updated 9/7/00 ***/ /*** No adjustment made to flow-resistance tolerance ***/ /*RQtol = RQtol/Hexp;*/ /* See if default reaction coeffs. apply */ for (i=1; i<=Nlinks; i++) { if (Link[i].Type > PIPE) continue; if (Link[i].Kb == MISSING) Link[i].Kb = Kbulk; /* Bulk coeff. */ if (Link[i].Kw == MISSING) /* Wall coeff. */ { /* Rfactor is the pipe roughness correlation factor */ if (Rfactor == 0.0) Link[i].Kw = Kwall; else if ((Link[i].Kc > 0.0) && (Link[i].Diam > 0.0)) { if (Formflag == HW) Link[i].Kw = Rfactor/Link[i].Kc; if (Formflag == DW) Link[i].Kw = Rfactor/ABS(log(Link[i].Kc/Link[i].Diam)); if (Formflag == CM) Link[i].Kw = Rfactor*Link[i].Kc; } else Link[i].Kw = 0.0; } } for (i=1; i<=Ntanks; i++) if (Tank[i].Kb == MISSING) Tank[i].Kb = Kbulk; /* Use default pattern if none assigned to a demand */ for (i=1; i<=Nnodes; i++) { for (demand = Node[i].D; demand != NULL; demand = demand->next) if (demand->Pat == 0) demand->Pat = DefPat; } /* Remove QUALITY as a reporting variable if no WQ analysis */ if (Qualflag == NONE) Field[QUALITY].Enabled = FALSE; } /* End of adjustdata */
int Logger::info(const char* message) { return log(lowPriorityBuffer(), SLOG2_INFO, message); }
long long BSONTool::processFile( const path& root ){ string fileString = root.string(); long long fileLength = file_size( root ); if ( fileLength == 0 ) { out() << "file " << fileString << " empty, skipping" << endl; return 0; } FILE* file = fopen( fileString.c_str() , "rb" ); if ( ! file ){ log() << "error opening file: " << fileString << endl; return 0; } #if !defined(__sunos__) && defined(POSIX_FADV_SEQUENTIAL) posix_fadvise(fileno(file), 0, fileLength, POSIX_FADV_SEQUENTIAL); #endif log(1) << "\t file size: " << fileLength << endl; long long read = 0; long long num = 0; long long processed = 0; const int BUF_SIZE = 1024 * 1024 * 5; boost::scoped_array<char> buf_holder(new char[BUF_SIZE]); char * buf = buf_holder.get(); ProgressMeter m( fileLength ); while ( read < fileLength ) { int readlen = fread(buf, 4, 1, file); int size = ((int*)buf)[0]; if ( size >= BUF_SIZE ){ cerr << "got an object of size: " << size << " terminating..." << endl; } uassert( 10264 , "invalid object size" , size < BUF_SIZE ); readlen = fread(buf+4, size-4, 1, file); BSONObj o( buf ); if ( _objcheck && ! o.valid() ){ cerr << "INVALID OBJECT - going try and pring out " << endl; cerr << "size: " << size << endl; BSONObjIterator i(o); while ( i.more() ){ BSONElement e = i.next(); try { e.validate(); } catch ( ... ){ cerr << "\t\t NEXT ONE IS INVALID" << endl; } cerr << "\t name : " << e.fieldName() << " " << e.type() << endl; cerr << "\t " << e << endl; } } if ( _matcher.get() == 0 || _matcher->matches( o ) ){ gotObject( o ); processed++; } read += o.objsize(); num++; m.hit( o.objsize() ); } uassert( 10265 , "counts don't match" , m.done() == fileLength ); out() << "\t " << m.hits() << " objects found" << endl; if ( _matcher.get() ) out() << "\t " << processed << " objects processed" << endl; return processed; }
int Logger::notice(const char* message) { return log(lowPriorityBuffer(), SLOG2_NOTICE, message); }
/* HTS_Engine_get_volume: get volume in db */ double HTS_Engine_get_volume(HTS_Engine * engine) { return log(engine->condition.volume) / DB; }
int Logger::warn(const char* message) { return log(lowPriorityBuffer(), SLOG2_WARNING, message); }
void TDDSwitchTest::tearDown() { log("TDD tearDown is called"); log("Please do clean up here"); }
void PvaPyLogger::critical(const std::string& message) const { log(LogLevelCritical, message.c_str()); }