Changeset 1954


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

setting up scaling of PDE as in AirFlow? example

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

Legend:

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

    r1952 r1954  
    342342  m_State =   AutoPtr<PetscVector<Number> >(new PetscVector<Number>(NumGlobStates, NumLocStates, m_StateDofMap->get_send_list(), GHOSTED));
    343343  m_Control = AutoPtr<PetscVector<Number> >(new PetscVector<Number>(NumGlobControls, NumLocControls));
    344   m_PDEConstrScale = 1/pow(h,m_StateMesh.mesh_dimension()-2);
     344  m_PDEConstrScale = 1e-2/pow(h,m_StateMesh.mesh_dimension()-2);
    345345  m_PDEConstr = AutoPtr<PetscVector<Number> >(new PetscVector<Number>(NumGlobPDEConstr, NumLocPDEConstr));
    346346  m_AuxConstr = AutoPtr<PetscVector<Number> >(new PetscVector<Number>(NumGlobAuxConstr, NumLocAuxConstr));
  • branches/parallel/Ipopt/contrib/PetscInterface/Example2/Problem.dat

    r1950 r1954  
    1 h=0.1
     1h=0.07
    22# MinX1, MinX2, MinX3, MaxX1, MaxX2, MaxX3
    33# MinX1, MinX2, MaxX1, MaxX2
Note: See TracChangeset for help on using the changeset viewer.