PetrKryslUCSD

using JuAFEM
# using UnicodePlots

t0 = time()
grid = generate_grid(Quadrilateral, (1000,1000))
addnodeset!(grid, "boundary", x -> abs(x[1]) ≈ 1 ||  abs(x[2]) ≈ 1);

dim = 2
ip = Lagrange{dim, RefCube, 1}()
qr = QuadratureRule{dim, RefCube}(2)

PetrKryslUCSD

using FinEtools
using FinEtools.MeshExportModule

println("""Heat conduction example from JuAFEM.""")
t0 = time()

A = 2.0
thermal_conductivity = eye(2,2); # conductivity matrix
function getsource!(forceout::FFltVec, XYZ::FFltMat, tangents::FFltMat, fe_label::FInt)
  forceout[1] = 1.0; #heat source

PetrKryslUCSD

module Driver

using FESetModule
using FEMMBaseModule

function run1()
J = ones(2, 2)
loc = zeros(1,3)
conn = reshape([1,2,3], 1, 3)
N = zeros(3,1)

PetrKryslUCSD

module FEMMBaseModule

typealias FInt Int
typealias FFlt Float64
typealias FFltVec Vector{FFlt}
typealias FIntVec Vector{FInt}
typealias FFltMat Matrix{FFlt}
typealias FIntMat Matrix{FInt}

abstract FESet