自学教程：C++ w2函数代码示例

void SolveShort(Geometry *geometry, Stencil *st_uw, Stencil *st_dw,		int k, int kend, int l, int m, double I_uw,		double *chi, double *S, double *I, double *Psi){  int    local;  double chi_uw, chi_dw, S_uw, S_dw, dS_dw, dS_uw,         dtau_uw, dtau_dw, w[3], c1, c2;   local = k*geometry->Nplane + m*geometry->Nx + l;  /* --- The upwind quantities --                      -------------- */  chi_uw  = Interpolate_3D(chi, geometry, st_uw, l, m);  dtau_uw = 0.5 * (chi_uw + chi[local]) * st_uw->ds;  S_uw    = Interpolate_3D(S, geometry, st_uw, l, m);  dS_uw   = (S_uw - S[local]) / dtau_uw;  if (k == kend) {    w2(dtau_uw, w);    /* --- Piecewise linear integration in last layer -- ------------ */    c1 = (S_uw - S[local]) / dtau_uw;    I[local]   = I_uw*(1.0 - w[0]) + w[0]*S[local] + w[1]*c1;    if (Psi) Psi[local] = w[0] - w[1]/dtau_uw;  } else {    w3(dtau_uw, w);    /* --- The downwind quantities --                  -------------- */    chi_dw  = Interpolate_3D(chi, geometry, st_dw, l, m);    dtau_dw = 0.5 * (chi[local] + chi_dw) * st_dw->ds;    S_dw    = Interpolate_3D(S, geometry, st_dw, l, m);    dS_dw   = (S[local] - S_dw) / dtau_dw;    /* --- Piecewise quadratic integration --          -------------- */    c1 = (dS_uw*dtau_dw + dS_dw*dtau_uw);    c2 = (dS_uw - dS_dw);    I[local] = I_uw*(1.0 - w[0]) + w[0]*S[local] +       (w[1]*c1 + w[2]*c2) / (dtau_uw + dtau_dw);    /* --- Try piecewise linear if quadratic gives negative           monochromatic intensity --                  -------------- */     if (I[local] < 0.0) {      c1 = dS_uw;      I[local] = (1.0 - w[0])*I_uw + w[0]*S[local] + w[1]*c1;      if (Psi) Psi[local] = w[0] - w[1]/dtau_uw;    } else {      if (Psi) {	c1 = dtau_uw - dtau_dw;	Psi[local] = w[0] + (w[1]*c1 - w[2]) / (dtau_uw * dtau_dw);      }    }  }}

static int fit2_read_regr( PIA *pi, int cont, int regr ){	int  a, b, r;	if (cont) {	  if (regr != 6) return 0xff;	  r = 7;	} else r = regr + 0x10;	w2(0xc); w0(r); w2(4); w2(5);	         w0(0); a = r1();	         w0(1); b = r1();	w2(4);	return j44(a,b);}

	void testVecDot() {		glam::vec4 v1(1, 2, 3.5, 100);		glam::vec4 w1(0.5, 80, 20, 0.1);		float u1 = glam::dot(v1, w1);		TS_ASSERT_EQUALS(u1, 240.5f);		glam::Vector<double, 7> v2(1, 2, 3, 4, 5, 6, 0);		glam::Vector<double, 7> w2(0, 6, 5, 4, 3, 2, 1);		TS_ASSERT_EQUALS(glam::dot(v2, w2), 70);	}

int main(){    sem_init(&f,0,10);    sem_init(&s,0,10);    w1();    w2(); //   pthread_mutex_init(&f,NULL); //   pthread_mutex_init(&s,NULL);}

  shared_ptr<OsmMap> createTestMap()  {    shared_ptr<OsmMap> map(new OsmMap());    _map = map;    shared_ptr<Node> n1 = createNode(0.0, 0.0);    n1->setTag("building", "yes");    n1->setTag("name", "n1");    shared_ptr<Way> w1(new Way(Status::Unknown1, map->createNextWayId(), 13.0));    w1->setTag("area", "yes");    w1->setTag("building", "yes");    w1->setTag("name", "w1");    w1->addNode(createNode(0.1, 0.0)->getId());    w1->addNode(createNode(0.2, 0.0)->getId());    w1->addNode(createNode(0.2, 0.1)->getId());    w1->addNode(w1->getNodeId(0));    map->addWay(w1);    shared_ptr<Way> w2(new Way(Status::Unknown1, map->createNextWayId(), 13.0));    w2->setTag("highway", "track");    w2->setTag("name", "w2");    w2->addNode(createNode(0.3, 0.0)->getId());    w2->addNode(createNode(0.3, 0.1)->getId());    map->addWay(w2);    shared_ptr<Way> w3(new Way(Status::Unknown1, map->createNextWayId(), 13.0));    w3->setTag("highway", "road");    w3->setTag("name", "w3");    w3->addNode(createNode(0.4, 0.0)->getId());    w3->addNode(createNode(0.4, 0.1)->getId());    map->addWay(w3);    shared_ptr<Way> w4(new Way(Status::Unknown1, map->createNextWayId(), 13.0));    w4->addNode(createNode(0.5, 0.0)->getId());    w4->addNode(createNode(0.7, 0.0)->getId());    w4->addNode(createNode(0.6, 0.1)->getId());    w4->addNode(w4->getNodeId(0));    map->addWay(w4);    shared_ptr<Way> w5(new Way(Status::Unknown1, map->createNextWayId(), 13.0));    w5->addNode(createNode(0.55, 0.01)->getId());    w5->addNode(createNode(0.65, 0.01)->getId());    w5->addNode(createNode(0.6, 0.05)->getId());    w5->addNode(w5->getNodeId(0));    map->addWay(w5);    shared_ptr<Relation> r1(new Relation(Status::Unknown1, 1, 15.0, "multipolygon"));    r1->setTag("building", "yes");    r1->setTag("name", "r1");    r1->addElement("outer", w4->getElementId());    r1->addElement("inner", w5->getElementId());    map->addRelation(r1);    return map;  }

void test1() {    worker<int> w1(boost::bind<int>(&calculatefib, 7));    worker<int> w2(boost::bind<int>(&calculatefib, 8));        worker<int> w3(boost::bind<int>(&fib_job::calc, boost::shared_ptr<fib_job>(new fib_job(11))));        int sum = w1.get() + w2.get() + w3.get();        std::cout << "Done, sum=" << sum << std::endl;}

/*  reference counting having cyclic dependencies (widget implementation)//*/void testRCIPtr(){    RCWidget w1(10);    RCWidget w2(w1);    w2.doThis();    std::cout << w1.showThat() << '/n';       // prints 10    std::cout << w2.showThat() << '/n';       // prints -1}

  void test() {    Wrapper w1;    clang_analyzer_eval(w1.value == 42); // expected-warning{{TRUE}}    Wrapper w2(50);    clang_analyzer_eval(w2.value == 50); // expected-warning{{TRUE}}    Wrapper w3(false);    clang_analyzer_eval(w3.value == 42); // expected-warning{{TRUE}}  }

static void  ktti_write_regr( PIA *pi, int cont, int regr, int val){	int r;	r = regr + cont_map[cont];	w0(r); w2(0xb); w2(0xa); w2(3); w2(6); 	w0(val); w2(3); w0(0); w2(6); w2(0xb);}

static void on20_connect ( PIA *pi){	pi->saved_r0 = r0();        pi->saved_r2 = r2();	w2(4);w0(0);w2(0xc);w2(4);w2(6);w2(4);w2(6);w2(4); 	if (pi->mode) { op(2); vl(8); op(2); vl(9); }	       else   { op(2); vl(0); op(2); vl(8); }}

static void epat_write_regr( PIA *pi, int cont, int regr, int val){	int r;	r = regr + cont_map[cont];	switch (pi->mode) {	case 0:	case 1:	case 2:	w0(0x60+r); w2(1); w0(val); w2(4);		break;	case 3:	case 4:	case 5: w3(0x40+r); w4(val);		break;	}}

static void epia_write_regr( PIA *pi, int cont, int regr, int val){       int  r;	regr += cont_map[cont];        switch (pi->mode)  {        case 0:        case 1:        case 2: r = regr^0x19;                w0(r); w2(1); w0(val); w2(3); w2(4);                break;	case 3:	case 4:        case 5: r = regr^0x40;                w3(r); w4(val); w2(4);                break;        }}

void test_fraction(){    Fraction w1(1, 2);    assert(w1.Numerator() == 1);    assert(w1.Denominator() == 2);    Fraction w2(12, -16);    assert(w2.Numerator() == -3);    assert(w2.Denominator() == 4);    std::cerr << "test_fraction ok." << std::endl;}

	bool matrix_eq(const M& w, const N& w2, float prec)	{ 		boost::test_tools::percent_tolerance_t <float> pt(prec);		boost::test_tools::close_at_tolerance<typename M::value_type> cmp(pt);		for(int i=0;i<w.shape()[0];i++){			for(int j=0;j<w.shape()[1];j++){				if( !cmp(w(i,j), w2(i,j)))					return false;			}		}		return true;	}

static void fit2_read_block( PIA *pi, char * buf, int count ){	int  k, a, b, c, d;	w2(0xc); w0(0x10);	for (k=0;k<count/4;k++) {		w2(4); w2(5);		w0(0); a = r1(); w0(1); b = r1();		w0(3); c = r1(); w0(2); d = r1();		buf[4*k+0] = j44(a,b);		buf[4*k+1] = j44(d,c);                w2(4); w2(5);                       a = r1(); w0(3); b = r1();                w0(1); c = r1(); w0(0); d = r1();                buf[4*k+2] = j44(d,c);                buf[4*k+3] = j44(a,b);	}	w2(4);}

static int comm_read_regr( PIA *pi, int cont, int regr ){       int     l, h, r;        r = regr + cont_map[cont];        switch (pi->mode)  {        case 0: w0(r); P1; w0(0);        	w2(6); l = r1(); w0(0x80); h = r1(); w2(4);                return j44(l,h);        case 1: w0(r+0x20); P1;         	w0(0); w2(0x26); h = r0(); w2(4);                return h;	case 2:	case 3:        case 4: w3(r+0x20); (void)r1();        	w2(0x24); h = r4(); w2(4);                return h;        }        return -1;}       

void Grid::setSize(unsigned int width, unsigned int height, unsigned int rows, unsigned int columns, unsigned int xSep, unsigned int ySep){    // Delete previous grid    mcs.clear();    // Int/Int always returns int    int w = (int)width / columns;    int h = (int)height / rows;    for(size_t i = 0; i < rows; i++)    {        for(size_t j = 0; j < columns; j++)        {            // Added counter clockwise            Coordinate w2(w * j + xSep, h * i + ySep);            Coordinate x(w * j + xSep, h * (i + 1) - ySep);            Coordinate y(w * (j + 1) - xSep, h * (i + 1) - ySep);            Coordinate z(w * (j + 1) - xSep, h * i + ySep);            mcs.push_back(w2);            mcs.push_back(x);            mcs.push_back(x);            mcs.push_back(y);            if (xSep > 0)            {                mcs.push_back(y);                mcs.push_back(z);            }            if (ySep > 0)            {                mcs.push_back(z);                mcs.push_back(w2);            }        }    }    if (xSep == 0)    {        mcs.push_back(Coordinate(0,0));        mcs.push_back(Coordinate(width,0));    }    if (ySep == 0)    {        mcs.push_back(Coordinate(width,0));        mcs.push_back(Coordinate(width,height));    }}

static void dstr_write_regr(  PIA *pi, int cont, int regr, int val ){       int  r;        r = regr + cont_map[cont];	w0(0x81); P1; 	if (pi->mode >= 2) { w0(0x11); } else { w0(1); }	P2; w0(r); P1;	        switch (pi->mode)  {        case 0:        case 1: w0(val); w2(5); w2(7); w2(5); w2(4);		break;	case 2:	case 3:        case 4: w4(val);                 break;        }}

static void ktti_write_block( PIA *pi, char * buf, int count ){	int k;	for (k=0;k<count/2;k++) {		w0(0x10); w2(0xb); w2(0xa); w2(3); w2(6);		w0(buf[2*k]); w2(3);		w0(buf[2*k+1]); w2(6);		w2(0xb);	}}

static void aten_write_block( PIA *pi, char * buf, int count ){	int k;	w0(0x88); w2(0xe); w2(6);	for (k=0;k<count/2;k++) {		w0(buf[2*k+1]); w2(0xe); w2(6);		w0(buf[2*k]); w2(7); w2(6);	}	w2(0xc);}

static void  fit3_write_regr( PIA *pi, int cont, int regr, int val){	if (cont == 1) return;	switch (pi->mode) {	case 0:	case 1: w2(0xc); w0(regr); w2(0x8); w2(0xc);		w0(val); w2(0xd);		w0(0);   w2(0xc);		break;	case 2: w2(0xc); w0(regr); w2(0x8); w2(0xc);		w4(val); w4(0);		w2(0xc);		break;	}}

int _tmain(int argc, _TCHAR* argv[]){	Widget w1;		// 
C++ w_assert1函数代码示例
C++ w0函数代码示例