Changeset 1952


Ignore:
Timestamp:
Mar 30, 2011 4:50:57 PM (9 years ago)
Author:
andreasw
Message:

increased nnz per row in Hessians; scaling as used in experiments

Location:
branches/parallel/Ipopt/contrib/PetscInterface
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • branches/parallel/Ipopt/contrib/PetscInterface/Example2/IpPetscPDETempRadiation.C

    r1951 r1952  
    374374    // each processor has hole matrix, but sets only part of local constraints
    375375    Mat petsc_mat;
    376     MatCreateSeqAIJ(PETSC_COMM_SELF,NumGlobStates,NumGlobStates,8,PETSC_NULL,&petsc_mat);
     376    MatCreateSeqAIJ(PETSC_COMM_SELF,NumGlobStates,NumGlobStates,30,PETSC_NULL,&petsc_mat);
    377377    m_HessStateState = AutoPtr<PetscMatrix<Number> >(new PetscMatrix<Number>(petsc_mat));
    378378  }
     
    381381    Mat petsc_mat;  MY_DBG_PRINT("LibMeshPDEBase::calc_objective_gradient called");
    382382
    383     MatCreateSeqAIJ(PETSC_COMM_SELF,NumGlobControls,NumGlobStates,8,PETSC_NULL,&petsc_mat);
     383    MatCreateSeqAIJ(PETSC_COMM_SELF,NumGlobControls,NumGlobStates,30,PETSC_NULL,&petsc_mat);
    384384    m_HessControlState = AutoPtr<PetscMatrix<Number> >(new PetscMatrix<Number>(petsc_mat));
    385385  }
     
    387387    // each processor has hole matrix, but sets only part of local constraints
    388388    Mat petsc_mat;
    389     MatCreateSeqAIJ(PETSC_COMM_SELF,NumGlobControls,NumGlobControls,8,PETSC_NULL,&petsc_mat);
     389    MatCreateSeqAIJ(PETSC_COMM_SELF,NumGlobControls,NumGlobControls,30,PETSC_NULL,&petsc_mat);
    390390    m_HessControlControl = AutoPtr<PetscMatrix<Number> >(new PetscMatrix<Number>(petsc_mat));
    391391  }
  • branches/parallel/Ipopt/contrib/PetscInterface/IpLibMeshPDE.C

    r1940 r1952  
    139139  PG_.CreateMesh(&mesh_,PG_.GetFE_Degree());
    140140 
    141   pde_scaling_ = 1.0/pow(PG_._h,mesh_.mesh_dimension()-2);
     141  pde_scaling_ = 1e-2/pow(PG_._h,mesh_.mesh_dimension()-2);
    142142#ifdef SCALE_AUX_BOUNDS
    143   v_min_constr_scaling_ = 1.0/pow(PG_._h,mesh_.mesh_dimension()-2);
     143  v_min_constr_scaling_ = 1e-1/pow(PG_._h,mesh_.mesh_dimension()-2);
    144144#else
    145   v_min_constr_scaling_ = 1.0/pow(PG_._h,mesh_.mesh_dimension()-4);
     145  v_min_constr_scaling_ = 1e-1/pow(PG_._h,mesh_.mesh_dimension()-4);
    146146#endif
    147147  //WriteNodeFile(mesh_, "MeshGen.node");
     
    353353    Mat petsc_mat;
    354354    //MatCreateMPIAIJ(PETSC_COMM_WORLD,n_state_global,n_state_global,n_state_global,n_state_global,8,PETSC_NULL,0,PETSC_NULL,&petsc_mat);
    355     MatCreateSeqAIJ(PETSC_COMM_SELF,n_state_global,n_state_global,8,PETSC_NULL,&petsc_mat);
     355    MatCreateSeqAIJ(PETSC_COMM_SELF,n_state_global,n_state_global,30,PETSC_NULL,&petsc_mat);
    356356    hess_state_state_ = new PetscMatrix<Number>(petsc_mat);
    357357  }
     
    360360    Mat petsc_mat;
    361361    //MatCreateMPIAIJ(PETSC_COMM_WORLD,n_control_global,n_state_global,n_control_global,n_state_global,8,PETSC_NULL,0,PETSC_NULL,&petsc_mat);
    362     MatCreateSeqAIJ(PETSC_COMM_SELF,n_control_global,n_state_global,8,PETSC_NULL,&petsc_mat);
     362    MatCreateSeqAIJ(PETSC_COMM_SELF,n_control_global,n_state_global,30,PETSC_NULL,&petsc_mat);
    363363    hess_control_state_ = new PetscMatrix<Number>(petsc_mat);
    364364  }
     
    367367    Mat petsc_mat;
    368368    //MatCreateMPIAIJ(PETSC_COMM_WORLD,n_control_global,n_control_global,n_control_global,n_control_global,8,PETSC_NULL,0,PETSC_NULL,&petsc_mat);
    369     MatCreateSeqAIJ(PETSC_COMM_SELF,n_control_global,n_control_global,8,PETSC_NULL,&petsc_mat);
     369    MatCreateSeqAIJ(PETSC_COMM_SELF,n_control_global,n_control_global,30,PETSC_NULL,&petsc_mat);
    370370    hess_control_control_ = new PetscMatrix<Number>(petsc_mat);
    371371  }
Note: See TracChangeset for help on using the changeset viewer.