Geomi
Problem.hpp
1 #ifndef DEF_PROBLEM
2 #define DEF_PROBLEM
3 
4 /*
5 template <typename M, typename Q, typename Alg>
6 class Problem
7 {
8 private:
9  // The actual physical configurations
10  DiscSyst<M,Q>& m_ds;
11 
12  ProblemInterface& m_interface;
13 
14 public:
15 
16  Problem (DiscSyst<M,Q>& ds, ProblemInterface& interface) :
17  m_ds(ds),
18  m_interface(interface)
19  { }
20 
21  ~Problem ()
22  { }
23 
24  bool
25  computeA (Alg& A, const Q& x)
26  {
27  return m_interface.computeA(A,y);
28  }
29 
30  bool
31  computeJacobianA (std::vector<Alg>& JA, const Q& x)
32  {
33  return m_interface.computeJacobianA(JA,y);
34  }
35 
36 };
37 */
38 
39 #endif
40