dolfin team mailing list archive

[Marie Rognes <732543@xxxxxxxxxxxxxxxxxx>]

Why does the solver automatically resize the vector?

-- 
You received this bug notification because you are a member of DOLFIN
Team, which is subscribed to DOLFIN.
https://bugs.launchpad.net/bugs/732543

Title:
  No error when attempting to solve using wrong dimensions

Status in DOLFIN:
  New

Bug description:
  Calling solve with a vector of non-matching dimensions does not
  produce an error. Something mysterious happens instead. Reproduce with
  below code.

  from dolfin import *

  mesh = UnitSquare(2, 2)

  V = FunctionSpace(mesh, "CG", 1)
  Q = FunctionSpace(mesh, "CG", 1)
  W = V * Q

  y = TrialFunction(W)
  w = TestFunction(W)

  a = inner(y, w)*dx
  L = w[0]*dx

  A = assemble(a)
  b = assemble(L)
  y = Function(W)

  info(A)
  info(b)
  u = Function(V)
  info(u.vector())
  solve(A, u.vector(), b)
  info(u.vector())