File: FatBoundary.cpp

package info (click to toggle)
freefem3d 1.0pre10-3.4
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 25,016 kB
  • ctags: 8,675
  • sloc: cpp: 57,204; sh: 8,788; yacc: 2,975; makefile: 1,149; ansic: 508; perl: 110
file content (67 lines) | stat: -rw-r--r-- 1,989 bytes parent folder | download | duplicates (6)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
//  This file is part of ff3d - http://www.freefem.org/ff3d
//  Copyright (C) 2001, 2002, 2003 Stphane Del Pino

//  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; either version 2, or (at your option)
//  any later version.

//  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.  

//  $Id: FatBoundary.cpp,v 1.7 2005/09/25 14:16:05 delpinux Exp $


#include <FatBoundary.hpp>
#include <PDESystem.hpp>

#include <Structured3DMesh.hpp>
#include <FEMDiscretization.hpp>

#include <MatrixManagement.hpp>

void FatBoundary::Discretize (ConstReferenceCounting<Problem> Pb)
{
  throw ErrorHandler(__FILE__,__LINE__,
		     "FBM is not implemented yet ",
		     ErrorHandler::unexpected);

//   //! Sets the problem.
//   __problem = Pb;

//   MemoryManager MM;
//   bool performAssembling
//     = MM.ReserveMatrix(__A,
// 		       problem().numberOfUnknown(),
// 		       mesh().numberOfVertices());

//   MM.ReserveVector(__b, 
// 		   problem().numberOfUnknown(),
// 		   mesh().numberOfVertices());

//   ffout(2) << "Discretizing the PDE Problem ... ";

//   DegreeOfFreedomSet __degreeOfFreedomSet(0,0,0,0);
//   FEMDiscretization<Structured3DMesh>
//     FEM(problem(), mesh(), *__A, *__b, __degreeOfFreedomSet);

//   if (performAssembling) {
//     FEM.assembleMatrix();
//   } else {
//     ffout(2) << "keeping previous operator discretization\n";
//   }
//   FEM.assembleSecondMember();
}

void FatBoundary::Compute ( Solution& u )
{
  ;
}