Hallo zusammen,

ich spiele gerade damit herum ein Chat Interface in HTML zu erzeugen.
Das klappt soweit auch ganz gut und ich habe in meinem 'Main' ein <p> Element.


<section class="main">
<p id="output"></p>
</section>


Nun soll dieses ueber ein Python-Script immer anders befuellt werden ...

Ich habe eine Lösung gefunden.


class Right(SubDomain):
def inside(self, x, on_boundary):
return on_boundary and near(x[0], 1)
rn = Expression("0.", degree=1)

boundaries = FacetFunction("size_t", mesh)
boundaries.set_all(0)
rightbound = Right()
rightbound.mark(boundaries, 1)
ds = Measure('ds ...

Hallo zusammen,

Zunächst weiß ich leider nicht, wie man im Forum mathematische Formeln erzeugt,
deswegen leider hier als Text.

Integral _ (Neumann-Rand) u^2 * ds

Also das Integral von u^2 auf dem Neumannrand ausgewertet.
Dies kann man als L2 Norm auffassen und somit

|| u^2 ||^2 _ Neumannrand ...

from dolfin import *
import numpy

mesh = UnitSquareMesh(50,50)
F = FunctionSpace(mesh, "Lagrange", 1)
x = TrialFunction(F)
y = TestFunction(F)

m = inner(grad(x),grad(y))*dx
M = assemble(m)

n = F.dim()
d = mesh.geometry().dim()
dof = F.tabulate_dof_coordinates().reshape(n,d)

o = numpy.array ...

Ich glaube, dass es so wie du sagtest funktioniert!
Gibt es eine Möglichkeit dieses nun auch auf einen Vektor anzuwenden?
In dem Array stecken alle Indizes die 0 gesetzt werden sollen.
Der Befehl v.zero(o) geht leider nicht.

from dolfin import *
import numpy

mesh = UnitSquareMesh(50,50)
F = FunctionSpace(mesh, "Lagrange", 1)
x = TrialFunction(F)
y = TestFunction(F)

m = inner(grad(x),grad(y))*dx
M = assemble(m)

n = F.dim()
d = mesh.geometry().dim()
dof = F.tabulate_dof_coordinates().reshape(n,d)

o = numpy.array ...

Hallo zusammen,

Mit Hilfe von FENICS möchte ich eine Zeile einer Matrix auf 0 setzen.
Hierzu soll die zero() Funktion verwendet werden:
https://fenicsproject.org/documentation/dolfin/1.0.1/python/programmers-reference/cpp/GenericMatrix.html

/edit: Alternativ kann das Ganze natürlich auch mit ...