dolfin team mailing list archive

[Marie Rognes <732543@xxxxxxxxxxxxxxxxxx>]

Public bug reported:

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())

** Affects: dolfin
     Importance: Undecided
         Status: New

-- 
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())