← Back to team overview

fenics team mailing list archive

Re: pde.solve() command parameters

 

You seem to be using an old version of DOLFIN. Try the latest release
(0.9.0) where LinearPDE has been replaced by VariationalProblem.

-- 
Anders


On Mon, Jan 19, 2009 at 05:39:59PM -0600, Brandon Chabaud wrote:
> Hi.
> 
> I am trying to use fenics to solve a 2D elasticity problem.  I have written and
> compiled a .form file using ffc with no problem.  However, when I try to
> compile the main.cpp file that I've written, I get this error message:
> 
> main.cpp: In function 'int main(int, char**)':
> main.cpp:177: error: no matching function for call to
> 'dolfin::LinearPDE::solve(dolfin::Function&, dolfin::Function&,
> dolfin::Function&, dolfin::Function&)'
> /home/grad/chabaud/local/include/dolfin/pde/LinearPDE.h:53: note:
> candidates are: void dolfin::LinearPDE::solve(dolfin::Function&)
> /home/grad/chabaud/local/include/dolfin/pde/LinearPDE.h:56: note: void
> dolfin::LinearPDE::solve(dolfin::Function&, dolfin::Function&)
> /home/grad/chabaud/local/include/dolfin/pde/LinearPDE.h:59: note: void
> dolfin::LinearPDE::solve(dolfin::Function&, dolfin::Function&,
> dolfin::Function&)
> scons: *** [main.o] Error 1
> scons: building terminated because of errors.
> The offending line in my cpp file seems to be
>  
> pde.solve(sigma,u,p,epsilon);
>  
> The lines immediately above it are
>  
> // Set up PDE
> SolidBilinearForm a(n);
> SolidLinearForm L(f1,g,u_prev,eps_prev,u0,F,uD,n);
> LinearPDE pde(a,L,mesh);
>  
> // Solve PDE
> Function u;
> Function epsilon;
> Function sigma;
> Function p;
> 
> The solutions u, epsilon, sigma, p to the problem are defined in my .form
> file by
> 
> Element1 = VectorElement("Discontinuous Lagrange", "triangle", 0, 4)
> 
> Element2 = VectorElement("Discontinuous Lagrange", "triangle", 0, 4)
> 
> Element3 = VectorElement("Discontinuous Lagrange", "triangle", 0)
> 
> Element4 = FiniteElement("Discontinuous Lagrange", "triangle", 0)
> 
> combo = MixedElement([Element1, Element3, Element4, Element2])
> 
> (z, v, q, w) = TestFunctions(combo)
> 
> (sigma, u, p, epsilon) = TrialFunctions(combo)
> 
> Do anyone know how to fix this error? Thanks.
> 

> _______________________________________________
> FEniCS-dev mailing list
> FEniCS-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/fenics-dev

Attachment: signature.asc
Description: Digital signature


Follow ups

References