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

voidcmdpkt_beacontimerinterrupt_819xusb(	struct net_device *dev){	struct r8192_priv *priv = ieee80211_priv(dev);	u16 tx_rate;	{		//		// 070117, rcnjko: 87B have to S/W beacon for DTM encryption_cmn.		//		if(priv->ieee80211->current_network.mode == IEEE_A  ||			priv->ieee80211->current_network.mode == IEEE_N_5G ||			(priv->ieee80211->current_network.mode == IEEE_N_24G  && (!priv->ieee80211->pHTInfo->bCurSuppCCK)))		{			tx_rate = 60;			DMESG("send beacon frame  tx rate is 6Mbpm/n");		}		else		{			tx_rate =10;			DMESG("send beacon frame  tx rate is 1Mbpm/n");		}		rtl819xusb_beacon_tx(dev,tx_rate); // HW Beacon	}}

static void cmpk_handle_interrupt_status(struct net_device *dev, u8 *pmsg){	struct cmpk_intr_sta rx_intr_status;	struct r8192_priv *priv = rtllib_priv(dev);	DMESG("---> cmpk_Handle_Interrupt_Status()/n");	rx_intr_status.length = pmsg[1];	if (rx_intr_status.length != (sizeof(struct cmpk_intr_sta) - 2)) {		DMESG("cmpk_Handle_Interrupt_Status: wrong length!/n");		return;	}	if (priv->rtllib->iw_mode == IW_MODE_ADHOC) {		rx_intr_status.interrupt_status = *((u32 *)(pmsg + 4));		DMESG("interrupt status = 0x%x/n",		      rx_intr_status.interrupt_status);		if (rx_intr_status.interrupt_status & ISR_TxBcnOk) {			priv->rtllib->bibsscoordinator = true;			priv->stats.txbeaconokint++;		} else if (rx_intr_status.interrupt_status & ISR_TxBcnErr) {			priv->rtllib->bibsscoordinator = false;			priv->stats.txbeaconerr++;		}		if (rx_intr_status.interrupt_status & ISR_BcnTimerIntr)			cmdpkt_beacontimerinterrupt_819xusb(dev);	}	DMESG("<---- cmpk_handle_interrupt_status()/n");}

/*----------------------------------------------------------------------------- * Function:    cmpk_handle_interrupt_status() * * Overview:    The function is responsible for extract the message from *				firmware. It will contain dedicated info in *				ws-07-0063-v06-rtl819x-command-packet-specification-070315.doc. * 				Please refer to chapter "Interrupt Status Element". * * Input:       struct net_device *dev, *			u8*	pmsg		-	Message Pointer of the command packet. * * Output:      NONE * * Return:      NONE * * Revised History: *  When			Who			Remark *  05/12/2008	amy		Add this for rtl8192 porting from windows code. * *---------------------------------------------------------------------------*/static	voidcmpk_handle_interrupt_status(	struct net_device *dev,	u8*	pmsg){	cmpk_intr_sta_t		rx_intr_status;	/* */	struct r8192_priv *priv = ieee80211_priv(dev);	DMESG("---> cmpk_Handle_Interrupt_Status()/n");	/* 0. Display received message. */	//cmpk_Display_Message(CMPK_RX_BEACON_STATE_SIZE, pMsg);	/* 1. Extract TX feedback info from RFD to temp structure buffer. */	/* It seems that FW use big endian(MIPS) and DRV use little endian in	   windows OS. So we have to read the content byte by byte or transfer	   endian type before copy the message copy. */	//rx_bcn_state.Element_ID 	= pMsg[0];	//rx_bcn_state.Length 		= pMsg[1];	rx_intr_status.length = pmsg[1];	if (rx_intr_status.length != (sizeof(cmpk_intr_sta_t) - 2))	{		DMESG("cmpk_Handle_Interrupt_Status: wrong length!/n");		return;	}	// Statistics of beacon for ad-hoc mode.	if(	priv->ieee80211->iw_mode == IW_MODE_ADHOC)	{		//2 maybe need endian transform?		rx_intr_status.interrupt_status = *((u32 *)(pmsg + 4));		//rx_intr_status.InterruptStatus = N2H4BYTE(*((UINT32 *)(pMsg + 4)));		DMESG("interrupt status = 0x%x/n", rx_intr_status.interrupt_status);		if (rx_intr_status.interrupt_status & ISR_TxBcnOk)		{			priv->ieee80211->bibsscoordinator = true;			priv->stats.txbeaconokint++;		}		else if (rx_intr_status.interrupt_status & ISR_TxBcnErr)		{			priv->ieee80211->bibsscoordinator = false;			priv->stats.txbeaconerr++;		}		if (rx_intr_status.interrupt_status & ISR_BcnTimerIntr)		{			cmdpkt_beacontimerinterrupt_819xusb(dev);		}	}	 // Other informations in interrupt status we need?	DMESG("<---- cmpk_handle_interrupt_status()/n");}	/* cmpk_handle_interrupt_status */

static int r8180_wx_set_retry(struct net_device *dev,                              struct iw_request_info *info,                              union iwreq_data *wrqu, char *extra){    struct r8180_priv *priv = ieee80211_priv(dev);    int err = 0;    if (priv->ieee80211->bHwRadioOff)        return 0;    down(&priv->wx_sem);    if (wrqu->retry.flags & IW_RETRY_LIFETIME ||            wrqu->retry.disabled)	{        err = -EINVAL;        goto exit;    }    if (!(wrqu->retry.flags & IW_RETRY_LIMIT))	{        err = -EINVAL;        goto exit;    }    if (wrqu->retry.value > R8180_MAX_RETRY)	{        err = -EINVAL;        goto exit;    }    if (wrqu->retry.flags & IW_RETRY_MAX) {        priv->retry_rts = wrqu->retry.value;        DMESG("Setting retry for RTS/CTS data to %d", wrqu->retry.value);    }	else {        priv->retry_data = wrqu->retry.value;        DMESG("Setting retry for non RTS/CTS data to %d", wrqu->retry.value);    }    /* FIXME !     * We might try to write directly the TX config register     * or to restart just the (R)TX process.     * I'm unsure if whole reset is really needed     */    rtl8180_commit(dev);    /*    if(priv->up){    	rtl8180_rtx_disable(dev);    	rtl8180_rx_enable(dev);    	rtl8180_tx_enable(dev);    }    */exit:    up(&priv->wx_sem);    return err;}

static void cmdpkt_beacontimerinterrupt_819xusb(struct net_device *dev){	struct r8192_priv *priv = rtllib_priv(dev);	if ((priv->rtllib->current_network.mode == IEEE_A)  ||	    (priv->rtllib->current_network.mode == IEEE_N_5G) ||	    ((priv->rtllib->current_network.mode == IEEE_N_24G)  &&	    (!priv->rtllib->pHTInfo->bCurSuppCCK)))		DMESG("send beacon frame  tx rate is 6Mbpm/n");	else		DMESG("send beacon frame  tx rate is 1Mbpm/n");}

static int r8192_wx_set_crcmon(struct net_device *dev,			       struct iw_request_info *info,			       union iwreq_data *wrqu, char *extra){	struct r8192_priv *priv = ieee80211_priv(dev);	int *parms = (int *)extra;	int enable = (parms[0] > 0);	short prev = priv->crcmon;	down(&priv->wx_sem);	if(enable)		priv->crcmon=1;	else		priv->crcmon=0;	DMESG("bad CRC in monitor mode are %s",	      priv->crcmon ? "accepted" : "rejected");	if(prev != priv->crcmon && priv->up){		//rtl8180_down(dev);		//rtl8180_up(dev);	}	up(&priv->wx_sem);	return 0;}

int r8180_wx_set_key(struct net_device *dev, struct iw_request_info *info, 		     union iwreq_data *wrqu, char *key){	struct r8180_priv *priv = ieee80211_priv(dev);	struct iw_point *erq = &(wrqu->encoding);		if (erq->flags & IW_ENCODE_DISABLED) {	}		/*	i = erq->flags & IW_ENCODE_INDEX;	if (i < 1 || i > 4)*/			if (erq->length > 0) {		//int len = erq->length <= 5 ? 5 : 13;				u32* tkey= (u32*) key;		priv->key0[0] = tkey[0];		priv->key0[1] = tkey[1];		priv->key0[2] = tkey[2];		priv->key0[3] = tkey[3] &0xff;		DMESG("Setting wep key to %x %x %x %x", 		      tkey[0],tkey[1],tkey[2],tkey[3]);		rtl8180_set_hw_wep(dev);	}	return 0;}

int SQLAccessFactory::MakeSQLExecuteFactoryRequest(std::string &dran                     , std::string &command                     , SQLExecuteFactoryRequest **req){  SQLExecuteFactoryRequest*  seReq      = soap_new__wsdair__SQLExecuteFactoryRequest(&soap,-1);  *req = seReq;  seReq->wsdair__SQLExpression = soap_new_wsdair__SQLExpressionType(&soap, -1) ;  seReq->wsdair__SQLExpression->Expression = command ;  seReq->wsdai__DataResourceAbstractName = dran;  DMESG("DRAN : " << seReq->wsdai__DataResourceAbstractName << endl);  DMESG("EXPR : " << seReq->wsdair__SQLExpression->Expression << endl);  return SOAP_OK;}

    // Destructor    FileIO::~FileIO()  {       // Free all HDF5 resources       /* End access to the dataset and release resources used by it. */       // if(asciiOutput == true) asciiOutputFile.close();     // close some extra stuff     check( H5Tclose(complex_tid),   DMESG("H5Tclose"));     check( H5Tclose(timing_tid),   DMESG("H5Tclose"));     check( H5Tclose(species_tid),  DMESG("H5Tclose"));     check( H5Tclose(s256_tid),     DMESG("H5Tclose"));              // close file     check( H5Fclose(file) , DMESG("Unable to close file ..."));     delete cfl_table;}

Array2c Analysis::getSpectrum(unsigned int dir) {       Array2c spectrum;       if     (dir == SPEC_XZ)  spectrum.reference(spectrumXZ);       else if(dir == SPEC_YZ)  spectrum.reference(spectrumXZ);       else if(dir == SPEC_XY)  spectrum.reference(spectrumXY);       else   check(-1, DMESG("No such direction for get spectrum"));           return spectrum;   };

static int r8192_wx_set_retry(struct net_device *dev,				struct iw_request_info *info,				union iwreq_data *wrqu, char *extra){	struct r8192_priv *priv = rtllib_priv(dev);	int err = 0;	if (priv->bHwRadioOff == true)		return 0;	down(&priv->wx_sem);	if (wrqu->retry.flags & IW_RETRY_LIFETIME ||	    wrqu->retry.disabled) {		err = -EINVAL;		goto exit;	}	if (!(wrqu->retry.flags & IW_RETRY_LIMIT)) {		err = -EINVAL;		goto exit;	}	if (wrqu->retry.value > R8192_MAX_RETRY) {		err = -EINVAL;		goto exit;	}	if (wrqu->retry.flags & IW_RETRY_MAX) {		priv->retry_rts = wrqu->retry.value;		DMESG("Setting retry for RTS/CTS data to %d",		      wrqu->retry.value);	} else {		priv->retry_data = wrqu->retry.value;		DMESG("Setting retry for non RTS/CTS data to %d",		      wrqu->retry.value);	}	rtl8192_commit(dev);exit:	up(&priv->wx_sem);	return err;}

/*----------------------------------------------------------------------------- * Function:    cmpk_handle_interrupt_status() * * Overview:    The function is responsible for extract the message from *		firmware. It will contain dedicated info in *		ws-07-0063-v06-rtl819x-command-packet-specification-070315.doc. *		Please refer to chapter "Interrupt Status Element". * * Input:       struct net_device *dev *              u8 *pmsg		- Message Pointer of the command packet. * * Output:      NONE * * Return:      NONE * * Revised History: *  When		Who	Remark *  05/12/2008		amy	Add this for rtl8192 porting from windows code. * *--------------------------------------------------------------------------- */static void cmpk_handle_interrupt_status(struct net_device *dev, u8 *pmsg){	cmpk_intr_sta_t		rx_intr_status;	/* */	struct r8192_priv *priv = ieee80211_priv(dev);	DMESG("---> cmpk_Handle_Interrupt_Status()/n");	/* 1. Extract TX feedback info from RFD to temp structure buffer. */	/* It seems that FW use big endian(MIPS) and DRV use little endian in	 * windows OS. So we have to read the content byte by byte or transfer	 * endian type before copy the message copy.	 */	rx_intr_status.length = pmsg[1];	if (rx_intr_status.length != (sizeof(cmpk_intr_sta_t) - 2)) {		DMESG("cmpk_Handle_Interrupt_Status: wrong length!/n");		return;	}	/* Statistics of beacon for ad-hoc mode. */	if (priv->ieee80211->iw_mode == IW_MODE_ADHOC) {		/* 2 maybe need endian transform? */		rx_intr_status.interrupt_status = *((u32 *)(pmsg + 4));		DMESG("interrupt status = 0x%x/n",		      rx_intr_status.interrupt_status);		if (rx_intr_status.interrupt_status & ISR_TxBcnOk) {			priv->ieee80211->bibsscoordinator = true;			priv->stats.txbeaconokint++;		} else if (rx_intr_status.interrupt_status & ISR_TxBcnErr) {			priv->ieee80211->bibsscoordinator = false;			priv->stats.txbeaconerr++;		}		if (rx_intr_status.interrupt_status & ISR_BcnTimerIntr)			cmdpkt_beacontimerinterrupt_819xusb(dev);	}	/* Other informations in interrupt status we need? */	DMESG("<---- cmpk_handle_interrupt_status()/n");}

int S_SQLAccessFactory::commitDirectoryToIndex2(EntryProps &prop                           , const std::string &indexTable                           , std::string &table                           , Statement &statement){   std::ostringstream query;    query << "UPDATE " << indexTable         << " SET /"table_name/" = '" << table  << "', "        << " /"directory/" = '" << dran << "',"        << " /"flags/" = '" << prop.flags << "'"        << " WHERE /"id/"=" << prop.id        << ";";  DMESG("SQL: >" << query.str() << "<" << std::endl);    if(statement.exec(query.str())){    DMESG("16 Directory exists/n");    return -1;  }   return prop.id;}

static int r8180_wx_set_enc(struct net_device *dev, 			    struct iw_request_info *info, 			    union iwreq_data *wrqu, char *key){	struct r8180_priv *priv = ieee80211_priv(dev);	int ret;		down(&priv->wx_sem);			DMESG("Setting SW wep key");		ret = ieee80211_wx_set_encode(priv->ieee80211,info,wrqu,key);				up(&priv->wx_sem);	return ret;}

static int r8180_wx_set_beaconinterval(struct net_device *dev, struct iw_request_info *aa,			  union iwreq_data *wrqu, char *b){	int *parms = (int *)b;	int bi = parms[0];		struct r8180_priv *priv = ieee80211_priv(dev);		down(&priv->wx_sem);	DMESG("setting beacon interval to %x",bi);		priv->ieee80211->beacon_interval=bi;	rtl8180_commit(dev);	up(&priv->wx_sem);			return 0;	}

void SQLAccessFactory::PrintEndpointReference(EndpointReference* er,    std::string &result){  if ( er == NULL ) return ;  if (er->Address != NULL )  {    DMESG("EndPoint : " << er->Address->__item << std::endl);    std::cout << er->Address->__item << std::endl;  }  if (er->wsa__ReferenceParameters != NULL )  {    std::string rp_dran;    std::string str = er->wsa__ReferenceParameters->__item;    dran_deserialize(str, rp_dran);    std::cout << "DRAN : " << rp_dran << std::endl;    result = rp_dran;  }}

int r8180_wx_set_key(struct net_device *dev, struct iw_request_info *info,		     union iwreq_data *wrqu, char *key){	struct r8180_priv *priv = ieee80211_priv(dev);	struct iw_point *erq = &(wrqu->encoding);	if (priv->ieee80211->bHwRadioOff)		return 0;	if (erq->length > 0) {		u32* tkey = (u32*) key;		priv->key0[0] = tkey[0];		priv->key0[1] = tkey[1];		priv->key0[2] = tkey[2];		priv->key0[3] = tkey[3] & 0xff;		DMESG("Setting wep key to %x %x %x %x",		      tkey[0], tkey[1], tkey[2], tkey[3]);		rtl8180_set_hw_wep(dev);	}	return 0;}

static int r8192_wx_set_crcmon(struct net_device *dev,			       struct iw_request_info *info,			       union iwreq_data *wrqu, char *extra){	struct r8192_priv *priv = ieee80211_priv(dev);	int *parms = (int *)extra;	int enable = (parms[0] > 0);	mutex_lock(&priv->wx_mutex);	if (enable)		priv->crcmon = 1;	else		priv->crcmon = 0;	DMESG("bad CRC in monitor mode are %s",	      priv->crcmon ? "accepted" : "rejected");	mutex_unlock(&priv->wx_mutex);	return 0;}

void interp(uchar code){	uchar c = 0;	uchar d = 0;	DMESG("{");			if (code == 0xe0){		e0=1;		DMESG("SETE0,");		return;	}			if (code == 0xf0) {		DMESG("SETRELEASE,");		release=1;			return;	}			if (e0){			c = xlat(ktab_eo, code);		DMESGF("xlateo %x->%x", code,c);	}	else {		c = xlat(ktab, code);		DMESGF("xlat %x->%x", code,c);	}	if (c==K_SHIFT) {		if (release){			DMESG("UNSHIFT,");			shift=0;  //handle letting go of shift key		}		else {			DMESG("SHIFT,");			shift=1;  //handle pressing shift key		}			release=0;		e0=0;		DMESG("unRELEASE0,unE0");		return;	}	DMESG("nE0");	e0=0;	if (release) {		DMESG("ISREL,unREL");		release=0;		return;  //release codes don't generate characters	}	release=0;		if (c==0)		return;	//generate capitals	if (shift) {		DMESG("isshift,");		if ((c>='a') && (c<='z'))	{			DMESGF("alpha %c,", c);			c=c-'a'+'A';		} else {			d = xlat(ktab_shift, c); //shift the key			DMESGF("shifted %c->%c", c, d);			if(d)				c=d;		}	}			DMESGF("out %c}", c);	nextchar=c;		}

    int FileIO::create(Setup *setup) {             hid_t file_plist = H5Pcreate(H5P_FILE_ACCESS);#ifdef GKC_PARALLEL_MPI   //       pass some information onto the underlying MPI_File_open call           MPI_Info file_info;          check(MPI_Info_create(&file_info), DMESG("File info"));          /*           H5Pset_sieve_buf_size(file_plist, 262144);           H5Pset_alignment(file_plist, 524288, 262144);                          MPI_Info_set(file_info, (char *) "access_style"        , (char *) "write_once");          MPI_Info_set(file_info, (char *) "collective_buffering", (char *) "true");          MPI_Info_set(file_info, (char *) "cb_block_size"       , (char *) "1048576");          MPI_Info_set(file_info, (char *) "cb_buffer_size"      , (char *) "4194304");           * */          check( H5Pset_fapl_mpio(file_plist, parallel->Comm[DIR_ALL], file_info), DMESG("Set MPI Property"));#endif        file = check(H5Fcreate(outputFileName.c_str(), (overwriteFile ? H5F_ACC_TRUNC : H5F_ACC_EXCL),                        H5P_DEFAULT, file_plist ), DMESG("H5FCreate : HDF5 File (File already exists ? use -f to overwrite) : " + outputFileName));        check( H5Pclose(file_plist),   DMESG("H5Pclose"));#ifdef GKC_PARALLEL_MPI        MPI_Info_free(&file_info);#endif                 //////////////////////////////////////////////////////////////// Info Group ////////////////////////////////////////////////////////          hid_t infoGroup = check(H5Gcreate(file, "/Info",H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT), DMESG("Error creating group file for Phasespace : H5Gcreate"));         check(H5LTset_attribute_string(infoGroup, ".", "Output", outputFileName.c_str()), DMESG("H5LTset_attribute"));         check(H5LTset_attribute_string(infoGroup, ".", "Input",  inputFileName.c_str()), DMESG("H5LTset_attribute"));                           check(H5LTset_attribute_string(infoGroup, ".", "Version", PACKAGE_VERSION), DMESG("H5LTset_attribute"));         // Some Simulation specific stuff         //check(H5LTset_attribute_string(infoGroup, ".", "Solver", ((setup->Solver & VL_LIN) ? "Linear" : "Non-Linear")), DMESG("H5LTset_attribute"));         //heck(H5LTset_attribute_string(infoGroup, ".", "Type",   ((setup->VlasovType   & VLASOV_LOCAL ) ? "Local"  : "Global"    )), DMESG("H5LTset_attribute"));         //heck(H5LTset_attribute_string(infoGroup, ".", "FFTSolverS",   ((setup->VlasovType   & VLASOV_LOCAL ) ? "Local"  : "Global"    )), DMESG("H5LTset_attribute"));         //check(H5LTset_attribute_string(infoGroup, ".", "Initial Condition", setup->PerturbationMethod.c_str()), DMESG("H5LTset_attribute"));         check(H5LTset_attribute_string(infoGroup, ".", "Info", info.c_str()), DMESG("H5LTset_attribute"));                  check(H5LTset_attribute_string(infoGroup, ".", "Config", setup->configFileString.c_str()), DMESG("H5LTset_attribute"));         H5Gclose(infoGroup);                           /// Wrote setup constants, ugly here ////         hid_t constantsGroup = check(H5Gcreate(file, "/Constants",H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT), DMESG("Error creating group file for Phasespace : H5Gcreate"));         //         if (!setup->parser_constants.empty()) {                        std::vector<std::string> const_vec = Setup::split(setup->parser_constants, ",");            for(int s = 0; s < const_vec.size(); s++) {                 std::vector<std::string> key_value = Setup::split(const_vec[s],"=");                double value = Setup::string_to_double(key_value[1]);                int dim[] = { 1 };   //           check(H5LTmake_dataset_double(constantsGroup, Setup::trimLower(key_value[0], false).c_str(), 1, dim, &value ), DMESG("Write Constants Attributes"));                check(H5LTset_attribute_double(constantsGroup, ".", Setup::trimLower(key_value[0], false).c_str(), &value, 1), DMESG("H5LTset_attribute"));                //check(H5LTset_attribute_double(constantsGroup, ".", Setup::trimLower(key_value[0], false).c_str(), &(Setup::string_to_double(key_value[1])), 1), DMESG("H5LTset_attribute"));            };                  }                   H5Gclose(constantsGroup);                  // ********************* setup Table for CFL   *****************88         cfl_table = new CFLTable();                  cfl_offset[0] =  HOFFSET( CFLTable, timeStep );         cfl_offset[1] =  HOFFSET( CFLTable, time );         cfl_offset[2] =  HOFFSET( CFLTable, Fx );         cfl_offset[3] =  HOFFSET( CFLTable, Fy );         cfl_offset[4] =  HOFFSET( CFLTable, Fz  );         cfl_offset[5] =  HOFFSET( CFLTable, Fv );         cfl_offset[6] =  HOFFSET( CFLTable, total );                   for(int i = 1; i < 7; i++)  cfl_sizes[i] = sizeof(double); cfl_sizes[0] = sizeof(int);         hid_t   cfl_type[7]; for(int i = 1; i < 7; i++)  cfl_type [i] = H5T_NATIVE_DOUBLE; cfl_type[0] = H5T_NATIVE_INT;         const char *cfl_names[7];         cfl_names[0] = "timeStep";         cfl_names[1] = "time";         cfl_names[2] = "Fx"; cfl_names[3] = "Fy"; cfl_names[4] = "Fz"; cfl_names[5] = "Fv"; cfl_names[6] = "Total";          check(H5TBmake_table("cflTable", file, "cfl", (hsize_t) 7, (hsize_t) 0, sizeof(CFLTable), (const char**) cfl_names,                               cfl_offset, cfl_type, 32, NULL, 0, cfl_table ), DMESG("H5Tmake_table : cfl"));                  return HELIOS_SUCCESS;    }

void Grid::initData(FileIO *fileIO) {  hid_t gridGroup = fileIO->newGroup("/Grid");    // set lengths  check(H5LTset_attribute_double(gridGroup, ".", "Lx", &Lx, 1), DMESG("HDF-5 Error"));  check(H5LTset_attribute_double(gridGroup, ".", "Ly", &Ly, 1), DMESG("HDF-5 Error"));  check(H5LTset_attribute_double(gridGroup, ".", "Lz", &Lz, 1), DMESG("HDF-5 Error"));  check(H5LTset_attribute_double(gridGroup, ".", "Lv", &Lv, 1), DMESG("HDF-5 Error"));  check(H5LTset_attribute_double(gridGroup, ".", "Lm", &Lm, 1), DMESG("HDF-5 Error"));           // set grid point number  check(H5LTset_attribute_int(gridGroup, ".", "Nx", &Nx , 1), DMESG("HDF-5 Error"));  check(H5LTset_attribute_int(gridGroup, ".", "Nky",&Nky, 1), DMESG("HDF-5 Error"));  check(H5LTset_attribute_int(gridGroup, ".", "Nz", &Nz , 1), DMESG("HDF-5 Error"));  check(H5LTset_attribute_int(gridGroup, ".", "Nv", &Nv , 1), DMESG("HDF-5 Error"));  check(H5LTset_attribute_int(gridGroup, ".", "Nm", &Nm , 1), DMESG("HDF-5 Error"));  check(H5LTset_attribute_int(gridGroup, ".", "Ns", &Ns , 1), DMESG("HDF-5 Error"));           // set grids   check(H5LTset_attribute_double(gridGroup, ".", "X", &X[NxGlD], Nx), DMESG("HDF-5 Error"));  check(H5LTset_attribute_double(gridGroup, ".", "Z", &Z[NzGlD], Nz), DMESG("HDF-5 Error"));  check(H5LTset_attribute_double(gridGroup, ".", "V", &V[NvGlD], Nv), DMESG("HDF-5 Error"));  check(H5LTset_attribute_double(gridGroup, ".", "M", &M[NmGlD], Nm), DMESG("HDF-5 Error"));    check(H5LTset_attribute_string(gridGroup, ".", "MuIntegrationType", muIntegrationType.c_str()), DMESG("HDF-5 Error"));     H5Gclose(gridGroup);}

int S_SQLAccessFactory::handleQuery(){  // get what command is requested  MDInterpreter interpreter;  MDInterpreter::CommandToken cmdToken = interpreter.parseCommand (query);  // check whether the command is supported  int ret = checkInvalidCommand(soap ,cmdToken, _ACCESS_FACTORY);  if ( ret != SOAP_OK ) return ret;  ret = checkDBConnection(soap, helper);  if (ret != SOAP_OK ) return ret;  ret = checkTooManyIndirectServiceForUser();  if ( ret != SOAP_OK ) return ret;    Statement statement(*(helper->dbConn->dbServer->dbConn), debug);  // initialze prop  EntryProps prop;  prop.directoryTable = "SQLRespone";  prop.indexTable = "masterindex";  prop.flags = 0;  std::string table;  int error = 0 ;  do {    if( statement.beginTransaction()) { error = 1; break; }    // parse requested query    std::string translatedQuery;    helper->dbConn->dbServer->sqlParse(query, cmdToken, translatedQuery);    int isPlain = helper->dbConn->dbServer->asp->isPlain ;    DBResult dbResult;    std::vector<std::string> results;    helper->dbConn->getResult(dbResult,results);    if (dbResult.errorCode != 0) {   // error      return returnErrorDBResult(soap, dbResult);    }    std::stringstream query_s ;    // add an directory entry into masterindex    prop.id = helper->dbConn->dbServer->addDirectoryToIndex("SQLResponse"                , prop, statement);    if ( prop.id < 0 ) { error = 1; break; }    switch(cmdToken) {      case MDInterpreter::cSELECT:      {        table = "responses" + StringOf(prop.id) + "_"                 + StringOf(helper->dbConn->dbServer->asp->numColumn) ;        query_s << "CREATE VIEW " << table << " AS " << translatedQuery << ";";        break;      }      case MDInterpreter::cUPDATE:      case MDInterpreter::cDELETE:      case MDInterpreter::cINSERT:      {        query_s <<  translatedQuery << ";";        table = "responseu" + StringOf(prop.id) + "_"  ;        break;      }      default:        ERRLOG("Unsupported command token: " << cmdToken);      }    DMESG("Translated query : " <<  query_s.str() << std::endl);    if( statement.exec(query_s.str())) { error = 1; break; }    switch(cmdToken) {      case MDInterpreter::cUPDATE:      case MDInterpreter::cDELETE:      case MDInterpreter::cINSERT:      {        table += StringOf(statement.numRows() ) ;        break;      }      default:        ERRLOG("Unsupported command token: " << cmdToken);          }    dran = "/SQLResponse/" + table ;    prop.flags = EPT_DIR ;    if( isPlain ) prop.flags += EP_PLAIN;    prop.id = commitDirectoryToIndex2(prop, "masterindex", table, statement);  } while ( false ) ;  if (error)  {    return generateSoapError(soap, "Invalid Expression Fault" , NULL            , SOAP_TYPE_wsdai__InvalidExpressionFaultType );  }  if ( helper->dbConn->dbServer->finalizeDirectoryToIndex(prop.id, statement) < 0) {    return generateSoapError(soap, "Invalid Expression Fault" , NULL                , SOAP_TYPE_wsdai__InvalidExpressionFaultType );//.........这里部分代码省略.........

Grid:: Grid (Setup *setup, Parallel *parallel, FileIO *fileIO) {    muIntegrationType = setup->get("Grid.MuIntegrationType", "Gauss-Legendre");      // Set Initial conditions   // Global Domain Grid Number  Nx  = setup->get("Grid.Nx", 32);  Nky = setup->get("Grid.Nky", 9);  Nz  = setup->get("Grid.Nz", 8 );  Nv  = setup->get("Grid.Nv", 16);  Nm  = setup->get("Grid.Nm", 1 );  Ns  = setup->get("Grid.Ns", 1 );    // Global Domain Grid Lengths  Lx  = setup->get("Grid.Lx", 16.   );  Ly  = setup->get("Grid.Ly", 32.   );  Lz  = setup->get("Grid.Lz", 2. * M_PI);  Lv  = setup->get("Grid.Lv", 4.  );  Lm  = setup->get("Grid.Lm", 7.  );    int NmBoundary  = setup->get("Grid.NmBoundary", 0 );    const double x0  = setup->get("Grid.x0", Lx/2.);  //const double z0  = setup->get("Grid.z0", 0.   );      int Ny = 2 * Nky - 2;        check(((Nm == 1) && (Lm != 1.)) ? -1 : 0, DMESG("For Nm=1, set Lm = 1 !"));    // Calculate grid distances values (equidistant for X, Y, Z, V)  dx = (Nx > 1) ? Lx/((double) (Nx-1)) : Lx;  dy = (Ny > 1) ? Ly/((double) (Ny-1)) : Ly;  dz = (Nz > 1) ? Lz/((double) (Nz  )) : Lz;  dv = (Nv > 1) ? 2.* Lv/ ( (double) (Nv-1)) : Lv;        // Set Global Global lower/upper Domain/Boundary index  NxGlD = 3; NxGuD=Nx+2; NxGlB=1; NxGuB=Nx+4;   NyGlD = 3; NyGuD=Ny+2; NyGlB=1; NyGuB=Ny+4;   NzGlD = 3; NzGuD=Nz+2; NzGlB=1; NzGuB=Nz+4;   NvGlD = 3; NvGuD=Nv+2; NvGlB=1; NvGuB=Nv+4;   NmGlD = 1; NmGuD=Nm;   NmGlB=1; NmGuB=Nm;   NsGlD = 1; NsGuD=Ns;   NsGlB=1; NsGuB=Ns;          NkyGlD = 0  ; NkyGuD=Nky-1; NkyGlB=0; NkyGuB=Nky-1;   NkyLlD = 0  ; NkyLuD=Nky-1; NkyLlB=0; NkyLuB=Nky-1;  NkyLD  = Nky; NkyGD = Nky;     // Number of Ghost/Halo cells  NxGC = 2; NyGC = 2; NzGC = 2; NvGC = 2; NmGC = NmBoundary == 1 ? 2 : 0;   //////////////////   Grid decomposition  ////////////////////////  // Set local (L) lower (l) and upper (u) bounds  // Currently we only have equal partition decomposition  // Set local decomposition indices for X  NxLlD = (NxGuD - NxGlD + 1)/parallel->decomposition[DIR_X] * parallel->Coord[DIR_X] + NxGC + 1;  NxLuD = (NxGuD - NxGlD + 1)/parallel->decomposition[DIR_X] * (parallel->Coord[DIR_X]+1) + NxGC;  NxLlB = NxLlD - NxGC; NxLuB = NxLuD + NxGC;  // Set local decomposition indices for Y (wtf?)  NyLlD = (NyGuD - NyGlD + 1)/parallel->decomposition[DIR_Y] *  parallel->Coord[DIR_Y] + NyGC + 1;  NyLuD = (NyGuD - NyGlD + 1)/parallel->decomposition[DIR_Y] * (parallel->Coord[DIR_Y]+ 1) + NyGC;  NyLlB = NyLlD - NyGC; NyLuB = NyLuD + NyGC;  // Set local decomposition indices for Z  NzLlD = (NzGuD - NzGlD + 1)/parallel->decomposition[DIR_Z] *  parallel->Coord[DIR_Z] + NzGC + 1;  NzLuD = (NzGuD - NzGlD + 1)/parallel->decomposition[DIR_Z] * (parallel->Coord[DIR_Z]+1) + NzGC;  NzLlB =  NzLlD - NzGC; NzLuB = NzLuD + NzGC;      // Set local decomposition indices for V  NvLlD = (NvGuD - NvGlD + 1)/parallel->decomposition[DIR_V] *  parallel->Coord[DIR_V] + NvGC + 1;  NvLuD = (NvGuD - NvGlD + 1)/parallel->decomposition[DIR_V] * (parallel->Coord[DIR_V] + 1) + NvGC;  NvLlB = NvLlD - NvGC; NvLuB = NvLuD + NvGC;      // Set local decomposition indices for M (no boundary)  NmLlD = (NmGuD - NmGlD + 1)/parallel->decomposition[DIR_M] *  parallel->Coord[DIR_M] + 1;  NmLuD = (NmGuD - NmGlD + 1)/parallel->decomposition[DIR_M] * (parallel->Coord[DIR_M] + 1) ;  NmLlB = NmLlD - NmGC; NmLuB = NmLuD + NmGC;   // Set local decomposition indices for S (no boundary)  NsLlD = (NsGuD - NsGlD + 1)/parallel->decomposition[DIR_S] *  parallel->Coord[DIR_S] + 1;  NsLuD = (NsGuD - NsGlD + 1)/parallel->decomposition[DIR_S] * (parallel->Coord[DIR_S] + 1) ;  NsLlB = NsLlD; NsLuB = NsLuD;   // Set number of local domain points  NxLD = NxLuD - NxLlD + 1; NyLD = NyLuD - NyLlD + 1;  NzLD = NzLuD - NzLlD + 1; NvLD = NvLuD - NvLlD + 1;  NmLD = NmLuD - NmLlD + 1; NsLD = NsLuD - NsLlD + 1;  // Set number of local boundary points   NxLB = NxLuB - NxLlB + 1; NyLB = NyLuB - NyLlB + 1;  NzLB = NzLuB - NzLlB + 1; NvLB = NvLuB - NvLlB + 1;  NmLB = NmLuB - NmLlB + 1; NsLB = NsLuB - NsLlB + 1;      // Set number of local/global Domain points   NxGB = Nx + 2*NxGC  ; NyGB = Ny + 2*NyGC;  NzGB = Nz + 2*NxGC  ; NvGB = Nv + 2*NxGC;  NmGB = Nm           ; NsGB = Ns;//.........这里部分代码省略.........

void Analysis::initDataOutput(Setup *setup, FileIO *fileIO) {             analysisGroup = fileIO->newGroup(fileIO->getFileID(), "/Analysis");             hsize_t offset0[] = { 0, 0, 0, 0, 0, 0, 0 };         //###################################### Analysis - Heat fluxes ################################          // Heat Flux ky and Particle FluxKy ( per species)      hid_t fluxGroup = fileIO->newGroup(analysisGroup, "Flux");          hsize_t FSky_dim[]       = { plasma->nfields, Nky, Ns  , 1 };      hsize_t FSky_maxdim[]    = { plasma->nfields, Nky, Ns  , H5S_UNLIMITED} ;     hsize_t FSky_chunkdim[]  = { plasma->nfields, Nky, NsLD, 1 };     hsize_t FSky_chunkBdim[] = { plasma->nfields, Nky, NsLD, 1 };     hsize_t FSky_offset[]    = { 0, 0  , NsLlD-1, 0  };//     std::cout << "NsLD : " << NsLD << std::endl;//     std::cout << parallel->Coord(DIR_XYZVM) << std::endl << std::flush;//     parallel->barrier();     //check(-1, DMESG("STOP"));     FA_heatKy      = new FileAttr("Heat"   , fluxGroup, 4, FSky_dim, FSky_maxdim, FSky_chunkdim, offset0,  FSky_chunkBdim, FSky_offset, parallel->Coord(DIR_XYZVM == 0));     FA_particleKy  = new FileAttr("Density", fluxGroup, 4, FSky_dim, FSky_maxdim, FSky_chunkdim, offset0,  FSky_chunkBdim, FSky_offset, parallel->Coord(DIR_XYZVM == 0));         H5Gclose(fluxGroup);          //###################################### Moments - Heat fluxes ################################     hid_t momentGroup = check(H5Gcreate(fileIO->getFileID(), "/Moments",H5P_DEFAULT, H5P_DEFAULT, H5P_DEFAULT), DMESG("Error creating group file for Phi : H5Gcreate"));     hsize_t moment_dim[]       =  { grid->NzGD , NkyLD      , grid->NxGD , NsLD , 1};     hsize_t moment_maxdim[]    =  { grid->NzGD , NkyLD      , grid->NxGD , Ns   , H5S_UNLIMITED};     hsize_t moment_chunkBdim[] =  { NzLD       , NkyLD      , NxLD       , NsLD , 1          };     hsize_t moment_chunkdim[]  =  { NzLD       , NkyLD      , NxLD       , NsLD , 1};     hsize_t moment_moffset[]   =  { 0, 0, 0, 0, 0 };     hsize_t moment_offset[]    =  { NzLlD-3, 0     , NxLlD-3, 0     ,  0  };          bool momWrite = (parallel->Coord(DIR_VM) == 0);               FA_Mom_Tp        = new FileAttr("Temperature_v", momentGroup, 5, moment_dim , moment_maxdim   , moment_chunkdim   , moment_moffset    ,  moment_chunkBdim  , moment_offset, momWrite, fileIO->complex_tid);     FA_Mom_HeatFlux  = new FileAttr("HeatFlux"     , momentGroup, 5, moment_dim , moment_maxdim   , moment_chunkdim   , moment_moffset    ,  moment_chunkBdim  , moment_offset, momWrite, fileIO->complex_tid);     FA_Mom_Density   = new FileAttr("Density"      , momentGroup, 5, moment_dim , moment_maxdim   , moment_chunkdim   , moment_moffset    ,  moment_chunkBdim  , moment_offset, momWrite, fileIO->complex_tid);     FA_Mom_Time  = fileIO->newTiming(momentGroup);             H5Gclose(momentGroup);           dataOutputMoments   = Timing(setup->get("DataOutput.Moments.Step", -1)       , setup->get("DataOutput.Moments.Time", -1.));     //###################################### Power Spectrum  ################################################     // X-scalarValue      hid_t growGroup = fileIO->newGroup(analysisGroup, "PowerSpectrum");     hsize_t grow_x_dim[]       = { plasma->nfields, Nx/2+1, 1 };      hsize_t grow_x_maxdim[]    = { plasma->nfields, Nx/2+1, H5S_UNLIMITED} ;     hsize_t grow_x_chunkdim[]  = { plasma->nfields, Nx/2+1, 1 };     hsize_t grow_x_chunkBdim[] = { plasma->nfields, Nx/2+1, 1 };     FA_grow_x  = new FileAttr("X", growGroup, 3, grow_x_dim, grow_x_maxdim, grow_x_chunkdim, offset0,  grow_x_chunkBdim, offset0, parallel->myRank == 0);     // Y-scalarValue     hsize_t grow_y_dim[]       = { plasma->nfields, Nky, 1 };     hsize_t grow_y_maxdim[]    = { plasma->nfields, Nky, H5S_UNLIMITED };     hsize_t grow_y_chunkdim[]  = { plasma->nfields, Nky, 1 };     hsize_t grow_y_chunkBdim[] = { plasma->nfields, Nky, 1 };        FA_grow_y  = new FileAttr("Y", growGroup, 3, grow_y_dim, grow_y_maxdim, grow_y_chunkdim, offset0,  grow_y_chunkBdim, offset0, parallel->myRank == 0);     FA_grow_t  = fileIO->newTiming(growGroup);     H5Gclose(growGroup);               hid_t freqGroup = fileIO->newGroup(analysisGroup, "PhaseShift");     FA_freq_x  = new FileAttr("X", freqGroup, 3, grow_x_dim, grow_x_maxdim, grow_x_chunkdim, offset0,  grow_x_chunkBdim, offset0, parallel->myRank == 0);     FA_freq_y  = new FileAttr("Y", growGroup, 3, grow_y_dim, grow_y_maxdim, grow_y_chunkdim, offset0,  grow_y_chunkBdim, offset0, parallel->myRank == 0);     FA_freq_t  = fileIO->newTiming(freqGroup);     H5Gclose(freqGroup);      //////////////////////////////////////////////////////////////// Setup Table for scalar data ////////////////////////////////////////////////////////                    ScalarValues scalarValues;           size_t SV_offset[] = { HOFFSET( ScalarValues, timestep ) , HOFFSET( ScalarValues, time     ), HOFFSET( ScalarValues, phiEnergy ),                            HOFFSET( ScalarValues, ApEnergy ),  HOFFSET( ScalarValues, BpEnergy ), HOFFSET( ScalarValues, particle_number    ),                            HOFFSET( ScalarValues, kinetic_energy ), HOFFSET( ScalarValues, entropy ), HOFFSET( ScalarValues, heat_flux ),                            HOFFSET( ScalarValues, particle_flux ) };      size_t SV_sizes[] = { sizeof(scalarValues.timestep), sizeof(scalarValues.time    ), sizeof(scalarValues.phiEnergy),                            sizeof(scalarValues.ApEnergy), sizeof(scalarValues.BpEnergy), Ns * sizeof(scalarValues.particle_number[0]), Ns * sizeof(scalarValues.kinetic_energy[0]),                            Ns * sizeof(scalarValues.entropy[0]), Ns * sizeof(scalarValues.heat_flux[0]), Ns * sizeof(scalarValues.particle_flux[0])};      hid_t SV_types[] = { H5T_NATIVE_INT, H5T_NATIVE_DOUBLE, H5T_NATIVE_DOUBLE, H5T_NATIVE_DOUBLE, H5T_NATIVE_DOUBLE,                            fileIO->species_tid, fileIO->species_tid, fileIO->species_tid, fileIO->species_tid, fileIO->species_tid } ;        const char *SV_names[] = { "Timestep", "Time", "phiEnergy", "ApEnergy", "BpEnergy", "ParticleNumber", "KineticEnergy", "Entropy", "HeatFlux", "ParticleFlux" };      SVTable = new TableAttr(analysisGroup, "scalarValues", 10, SV_names, SV_offset, SV_types, SV_sizes, &scalarValues);       dataOutputStatistics  = Timing(setup->get("DataOutput.Statistics.Step", -1), setup->get("DataOutput.Statistics.Time", -1.));//.........这里部分代码省略.........