当前位置: 首页>>代码示例>>C++>>正文


C++ TPZFMatrix::VerifySymmetry方法代码示例

本文整理汇总了C++中TPZFMatrix::VerifySymmetry方法的典型用法代码示例。如果您正苦于以下问题:C++ TPZFMatrix::VerifySymmetry方法的具体用法?C++ TPZFMatrix::VerifySymmetry怎么用?C++ TPZFMatrix::VerifySymmetry使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在TPZFMatrix的用法示例。


在下文中一共展示了TPZFMatrix::VerifySymmetry方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。

示例1: ContributeHDiv

void TPZMatPoisson3d::Contribute(TPZMaterialData &data,REAL weight,TPZFMatrix<STATE> &ek,TPZFMatrix<STATE> &ef) {
	
    if(data.numberdualfunctions)
    {
        ContributeHDiv(data , weight , ek, ef);
        
        return;
    }
    
    if(fNeumann){
        
        LocalNeumanContribute(data , weight , ek, ef);
        
        return;
    }

    TPZFMatrix<REAL>  &phi = data.phi;
    TPZFMatrix<REAL> &dphi = data.dphix;
    TPZVec<REAL>  &x = data.x;
    TPZFMatrix<REAL> &axes = data.axes;
    TPZFMatrix<REAL> &jacinv = data.jacinv;
    int phr = phi.Rows();

    STATE fXfLoc = fXf;

    if(fForcingFunction) {            // phi(in, 0) = phi_in
        TPZManVector<STATE,1> res(1);
        TPZFMatrix<STATE> dres(Dimension(),1);
        fForcingFunction->Execute(x,res,dres);       // dphi(i,j) = dphi_j/dxi
        fXfLoc = res[0];
    }
    REAL delx = 0.;
    STATE ConvDirAx[3] = {0.};
    if(fC != 0.0) {
        int di,dj;
        delx = 0.;
        for(di=0; di<fDim; di++) {
            for(dj=0; dj<fDim; dj++) {
                delx = (delx<fabs(jacinv(di,dj))) ? fabs(jacinv(di,dj)) : delx;
            }
        }
        delx = 2./delx;
        
        
        switch(fDim) {
            case 1:
                ConvDirAx[0] = axes(0,0)*fConvDir[0]+axes(0,1)*fConvDir[1]+axes(0,2)*fConvDir[2];
                break;
            case 2:
                ConvDirAx[0] = axes(0,0)*fConvDir[0]+axes(0,1)*fConvDir[1]+axes(0,2)*fConvDir[2];
                ConvDirAx[1] = axes(1,0)*fConvDir[0]+axes(1,1)*fConvDir[1]+axes(1,2)*fConvDir[2];
                break;
            case 3:
                ConvDirAx[0] = axes(0,0)*fConvDir[0]+axes(0,1)*fConvDir[1]+axes(0,2)*fConvDir[2];
                ConvDirAx[1] = axes(1,0)*fConvDir[0]+axes(1,1)*fConvDir[1]+axes(1,2)*fConvDir[2];
                ConvDirAx[2] = axes(2,0)*fConvDir[0]+axes(2,1)*fConvDir[1]+axes(2,2)*fConvDir[2];
                break;
            default:
                PZError << "TPZMatPoisson3d::Contribute dimension error " << fDim << endl;
        }
    }

//    std::cout << " fxfloc " << fXfLoc << " weight " << weight << " prod " << fXfLoc*weight << std::endl;
    //Equacao de Poisson
    for( int in = 0; in < phr; in++ ) {
        int kd;
        STATE dphiic = 0;
        for(kd = 0; kd<fDim; kd++) dphiic += ConvDirAx[kd]*(STATE)dphi(kd,in);
        ef(in, 0) += - (STATE)weight * fXfLoc * ( (STATE)phi(in,0) + (STATE)(0.5*delx*fC)*fSD*dphiic );
        for( int jn = 0; jn < phr; jn++ ) {
            for(kd=0; kd<fDim; kd++) {
                ek(in,jn) += (STATE)weight * (
                                       +fK * (STATE)( dphi(kd,in) * dphi(kd,jn) ) 
                                       - (STATE)(fC* dphi(kd,in) * phi(jn)) * ConvDirAx[kd]
                                       + (STATE)(0.5 * delx * fC * dphi(kd,jn)) * fSD * dphiic * ConvDirAx[kd]
                                       );
            }
        }
    }

    if (this->IsSymetric()){    
        if ( !ek.VerifySymmetry() ) cout << __PRETTY_FUNCTION__ << "\nMATRIZ NAO SIMETRICA" << endl;
    }


}
开发者ID:labmec,项目名称:neopz,代码行数:86,代码来源:pzpoisson3d.cpp


注:本文中的TPZFMatrix::VerifySymmetry方法示例由纯净天空整理自Github/MSDocs等开源代码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。