本文整理汇总了C++中calc_residuals函数的典型用法代码示例。如果您正苦于以下问题:C++ calc_residuals函数的具体用法?C++ calc_residuals怎么用?C++ calc_residuals使用的例子?那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了calc_residuals函数的2个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于我们的系统推荐出更棒的C++代码示例。
示例1: init_states
void model_parameters::userfunction(void)
{
f =0.0;
//*******MAIN PROGRAM********
init_states();
calc_residuals();
calc_ref_points();
calc_obj_fun();
sd_sig=sig;
ofstream ofs("par.mcmc",ios::app);
if(mceval_phase()) mcmc_stuff();
//******************************
}
开发者ID:BenjaminHCCarr,项目名称:Collapse,代码行数:14,代码来源:SPC210.cpp
示例2: while
void US_MPI_Analysis::_2dsa_worker( void )
{
bool repeat_loop = true;
MPI_Job job;
MPI_Status status;
// Use 4 here because the master will be reading 4 with the
// same instruction when reading ::READY or ::RESULTS.
int x[ 4 ];
while ( repeat_loop )
{
MPI_Send( x, // Basically don't care
4,
MPI_INT,
MPI_Job::MASTER,
MPI_Job::READY,
my_communicator ); // let master know we are ready
//if(my_rank==1)
DbgLv(1) << "w:" << my_rank << ": ready sent";
// Blocking -- Wait for instructions
MPI_Recv( &job, // get masters' response
sizeof( job ),
MPI_BYTE,
MPI_Job::MASTER,
MPI_Job::TAG0,
my_communicator,
&status ); // status not used
//if(my_rank==1)
DbgLv(1) << "w:" << my_rank << ": job_recvd length" << job.length
<< "command" << job.command;
meniscus_value = job.meniscus_value;
int offset = job.dataset_offset;
int dataset_count = job.dataset_count;
int job_length = job.length;
DbgLv(1) << "w:" << my_rank << ": offs cnt" << offset << dataset_count;
data_sets[ offset ]->run_data.meniscus = meniscus_value;
data_sets[ offset ]->simparams.meniscus = meniscus_value;
switch( job.command )
{
case MPI_Job::PROCESS: // Process solutes
{
US_SolveSim::Simulation simulation_values;
simulation_values.noisflag =
parameters[ "tinoise_option" ].toInt() > 0 ? 1 : 0;
simulation_values.noisflag +=
parameters[ "rinoise_option" ].toInt() > 0 ? 2 : 0;
simulation_values.dbg_level = dbg_level;
simulation_values.dbg_timing = dbg_timing;
//DbgLv(1) << "w:" << my_rank << ": sols size" << job.length;
//if(my_rank==1)
DbgLv(1) << "w:" << my_rank << ": sols size" << job.length;
simulation_values.solutes.resize( job.length );
MPI_Recv( simulation_values.solutes.data(), // Get solutes
job_length * solute_doubles,
MPI_DOUBLE,
MPI_Job::MASTER,
MPI_Job::TAG0,
my_communicator,
&status );
max_rss();
//*DEBUG*
//if(dbg_level>0 && my_rank==1)
//if(my_rank==1)
{
int nn = simulation_values.solutes.size() - 1;
int mm = nn/2;
DbgLv(1) << "w:" << my_rank << ": offs dscnt" << offset << dataset_count
<< "vbar s20wc bott"
<< data_sets[offset]->vbar20
<< data_sets[offset]->s20w_correction
<< data_sets[offset]->centerpiece_bottom;
DbgLv(1) << "w:" << my_rank << ": sol0 solm soln"
<< simulation_values.solutes[0].s << simulation_values.solutes[0].k
<< simulation_values.solutes[mm].s << simulation_values.solutes[mm].k
<< simulation_values.solutes[nn].s << simulation_values.solutes[nn].k;
}
//*DEBUG*
calc_residuals( offset, dataset_count, simulation_values );
// Tell master we are sending back results
int size[ 4 ] = { simulation_values.solutes.size(),
simulation_values.ti_noise.size(),
simulation_values.ri_noise.size(),
(int)max_rss() };
DbgLv(1) << "w:" << my_rank << ": result sols size" << size[0]
<< "max_rss" << size[ 3 ];
//*DEBUG*
if(dbg_level==0 && my_rank==1) {
DbgLv(1) << "w:" << my_rank << ": result sols size" << size[0]
//.........这里部分代码省略.........
开发者ID:svn2github,项目名称:UltraScan3,代码行数:101,代码来源:2dsa_worker.cpp
注:本文中的calc_residuals函数示例由纯净天空整理自Github/MSDocs等源码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。 |
请发表评论