Saving unfinished changes

Author: fverdugo

src/FESpaces/FEOperators.jl

@@ -178,6 +178,15 @@ function LinearFEOperator(
 
 end
 
+function LinearFEOperator(
+  testfesp::FESpaceLike,
+  trialfesp::FESpaceLike,
+  terms::Vararg{<:AffineFETerm})
+
+  assem = SparseMatrixAssembler(testfesp,trialfesp)
+  LinearFEOperator(testfesp,trialfesp,assem,terms...)
+end
+
 function LinearFEOperator(
   biform::Function,
   liform::Function,
@@ -264,6 +273,14 @@ function NonLinearFEOperator(
   NonLinearFEOperator(testfesp,trialfesp,assem,terms)
 end
 
+function NonLinearFEOperator(
+  testfesp::FESpaceLike,
+  trialfesp::FESpaceLike,
+  terms::Vararg{<:FETerm})
+  assem = SparseMatrixAssembler(testfesp,trialfesp)
+  NonLinearFEOperator(testfesp,trialfesp,assem,terms...)
+end
+
 function NonLinearFEOperator(
   res::Function,
   jac::Function,

test/FESpacesTests/FEOperatorsTests.jl

@@ -84,6 +84,10 @@ eh1 = sqrt(sum( integrate(h1(e),trian,quad) ))
 @test el2 < 1.e-8
 @test eh1 < 1.e-8
 
+# Further tests
+op = LinearFEOperator(V,U,t_Ω,t_ΓN,t_ΓR,t_ΓS)
+op = LinearFEOperator(V,U,t_Ω)
+
 # For triangles
 
 model = simplexify(model)

test/MultiFieldTests/MultiFEOperatorsTests.jl

@@ -99,7 +99,7 @@ e2h1 = sqrt(sum( integrate(h1(e2),trian,quad) ))
 
 # Further tests
 
-# This is only to stress the single term API
 op = LinearFEOperator(a,b,V,U,assem,trian,quad)
+op = LinearFEOperator(V,U,t_Ω,t_Γ)
 
 end