testing api

test/test_api.jl

@@ -15,29 +15,49 @@ function test_basic()
     n2 = Node(2, [1.0, 0.0])
     n3 = Node(3, [1.0, 1.0])
     n4 = Node(4, [0.0, 1.0])
-
+    push!(model, n1, n2, n3, n4)
+    
     # create elements
     e1 = Element{Quad4}([n1, n2, n3, n4])
     e2 = Element{Seg2}([n1, n2])
 
+    # element set
+    elset = ElementSet("body", e1, e2)
+    elset2 = ElementSet("boundary", e2) 
+    push!(model, elset, elset2)
+
     # material properties
     material = Material("my material")
     material["temperature thermal conductivity"] = 6.0
     material["density"] = 36.0
+    set_material!(model, material, "Es")
 
-    # create element set body
-    body = ElementSet("body")
-    push!(body, e1, e2)
-
-    # add element set to model
-    add_element_set(model, body)
-
-    # assign material to elements in element set body
-    set_material(model, material, body)
+    # Create problem
+    field_problem = HeatProblem()
 
     # boundary conditions
+    bc = DirichletBC("body", "temperature" => 0.0)
+    # bc = DirichletBC("Es", 0.0)
+    ne = NeumannBC("boundary", "temperature flux" => ((0.0 => 0.0),
+                                                      (1.0=>600.0)))
+    
+    # LoadCase
+    case = LoadCase("Heat problem")
+    push!(case, field_problem)
+    push!(case, bc, ne)
 
+    # Get solver
+    solver = LinearSolver
+    push!(case, solver)
 
+    # add case
+    push!(model, case)
+
+    # Solve problem
+    solve!(model, 1.0, "Heat problem")
+    xi = [0.0, -1.0]
+    T = el1("temperature", xi, 1.0)
+    # T = model("temperature", [0.5, 0.0], 1)
 end
 
 end