PetrKryslUCSD

# This file is machine-generated - editing it directly is not advised

julia_version = "1.8.5"
manifest_format = "2.0"
project_hash = "5fb174050415bfb2078581bbe0e257b3dbe52a7c"

[[deps.Adapt]]
deps = ["LinearAlgebra", "Requires"]
git-tree-sha1 = "cc37d689f599e8df4f464b2fa3870ff7db7492ef"
uuid = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"

PetrKryslUCSD

module mwe

using LinearAlgebra: dot, norm
using Statistics: mean
using Cthulhu
using InteractiveUtils

cmplx(r, c) = r + im * c
size3(v) = norm(v)
ssize3(v) = v[1]^2 + v[2]^2 + v[3]^2

PetrKryslUCSD

module FESpaces

using StaticArrays
using MeshCore
using MeshCore: nshapes, indextype, nrelations, retrieve, manifdim, IncRel
using MeshKeeper: Mesh, baseincrel, increl
using ..FElements: nfeatofdim, ndofsperfeat, ndofsperelem

mutable struct FEField{N, T, IT}
    dofnums::Vector{SVector{N, IT}}

PetrKryslUCSD

using PlotlyJS
using JSON

plots = scatter3d(;x = rand(6), y = rand(6), z = rand(6))
layout = Layout(width="600", height="600", autosize=true, title="Change the view, then switch to commandline",
    scene=attr(
        xaxis=attr(title="X", gridcolor="rgb(255, 255, 255)",
          zerolinecolor="rgb(255, 255, 255)",
          showbackground=true,
          backgroundcolor="rgb(25, 25, 25)"),

PetrKryslUCSD

module testinglu

using LinearAlgebra
import LinearAlgebra: lu!, generic_lufact!, BlasInt, checknonsingular

function better_generic_lufact!(A::StridedMatrix{T}, ::Val{Pivot} = Val(true);
                         check::Bool = true) where {T,Pivot}
    m, n = size(A)
    minmn = min(m,n)
    info = 0

PetrKryslUCSD

module naivelu

using BenchmarkTools
using LinearAlgebra

"""
    naivelu!(a)

Naive LU decomposition implementation.

PetrKryslUCSD

module moinbounds13
using FinEtools
using FinEtools.MatrixUtilityModule: add_mggt_ut_only!, complete_lt!
using LinearAlgebra: Transpose, mul!
using BenchmarkTools
function add_mggt_ut_only_wo!(Ke::FFltMat, gradN::FFltMat, mult::FFlt)
    @assert size(Ke, 1) == size(Ke, 2)
    Kedim = size(Ke, 1)
    nne, mdim = size(gradN)
    mx::FInt = 0

PetrKryslUCSD

Variables:
  fens1::FinEtools.FENodeSetModule.FENodeSet
  fens2::FinEtools.FENodeSetModule.FENodeSet
  tolerance::Float64
  dim::Int64
  nn1<optimized out>
  nn2<optimized out>
  xyz1<optimized out>
  id1<optimized out>
  xyz2<optimized out>

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