main.cpp

OFELI is an object oriented library of C++ classes for development of finite element codes. Its main

/*==============================================================================

                                    O  F  E  L  I

                            Object  Finite  Element  Library

  ==============================================================================

   Copyright (C) 1998 - 2004 Rachid Touzani

   This program is free software; you can redistribute it and/or modify it under
   the terms of the GNU General Public License as published by the Free 
   Software Foundation; Version 2 of the License.

   This program is distributed in the hope that it will be useful, but WITHOUT
   ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
   FOR A PARTICULAR PURPOSE. See the GNU General Public License for more 
   details.

   You should have received a copy of the GNU General Public License 
   along with this program; if not, write to the :

   Free Software Foundation
   Inc., 59 Temple Place - Suite 330
   Boston, MA  02111-1307, USA

  ==============================================================================

               An example of a Finite Element Code using OFELI

            Solution of a 2-D Poisson problem using P1 Finite elements  
  A truncated Newton Algorithm is used to solve the resulting optimization problem

  ==============================================================================*/


#include "OFELI.h"
#include "Opt.h"
#include "User.h"
using namespace OFELI;

int main(int argc, char *argv[])
{
   if (argc <= 1) {
      cout << " Usage : ex5 <par_file>\n";
      exit(1);
   }
   IPF data(argv[1]);
   Mesh ms(data.MeshFile(1));
   User ud(ms);
   int n = ms.NbDOF();

// Declare solution vector
   Vect<double> x(n),low(n),up(n);
   Vect<int> pivot(n);

// Read in boundary conditions
   Vect<double> bc(n);
   ud.SetDBC(bc);

// Define and Solve the Optimization Problem
   Opt theOpt(ms,ud);
   x = 0;
   BCAsConstraint(ms,bc,up,low);
   OptimTN<Opt>(theOpt,x,low,up,pivot,100,1.e-12,1);

// Output solution
   cout << "\nSolution :\n" << x;

#ifdef WITH_PAUSE
   system("PAUSE");
#endif
   return 0;
}