Skip to content

Instantly share code, notes, and snippets.

@xtalax

xtalax/PDAE_codegen.jl

Created July 15, 2022 14:51
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save xtalax/55794f65bdb67831bcce8bbc4df74e2f to your computer and use it in GitHub Desktop.
Save xtalax/55794f65bdb67831bcce8bbc4df74e2f to your computer and use it in GitHub Desktop.
Download ZIP
Raw
PDAE_codegen.jl
begin
var"##f#257" = (ModelingToolkit.ODEFunctionClosure)(function (ˍ₋arg1, ˍ₋arg2, t)
begin
begin
var"v[20](t)" = (*)(0.8414709848078965, (exp)((*)(-1, t)))
var"v[3](t)" = (+)((+)((*)(0.5, (getindex)(ˍ₋arg1, 19)), (*)(0.5, (getindex)(ˍ₋arg1, 20))), (*)(0.00014552475590851675, (exp)((*)(-1, t))))
var"v[5](t)" = (+)((+)((*)(-1.0, var"v[3](t)"), (*)(2.0, (getindex)(ˍ₋arg1, 20))), (*)(-0.0004355664317500006, (exp)((*)(-1, t))))
var"v[6](t)" = (+)((+)((*)(-1.0, (getindex)(ˍ₋arg1, 20)), (*)(-0.0005788770749661647, (exp)((*)(-1, t)))), (*)(2.0, var"v[5](t)"))
var"v[7](t)" = (+)((+)((*)(-0.0007205845507054019, (exp)((*)(-1, t))), (*)(-1.0, var"v[5](t)")), (*)(2.0, var"v[6](t)"))
var"v[8](t)" = (+)((+)((*)(-1.0, var"v[6](t)"), (*)(2.0, var"v[7](t)")), (*)(-0.0008602964080896811, (exp)((*)(-1, t))))
var"v[9](t)" = (+)((+)((*)(-1.0, var"v[7](t)"), (*)(2.0, var"v[8](t)")), (*)(-0.0009976257229936814, (exp)((*)(-1, t))))
var"v[10](t)" = (+)((+)((*)(-1.0, var"v[8](t)"), (*)(2.0, var"v[9](t)")), (*)(-0.0011321921696093797, (exp)((*)(-1, t))))
var"v[11](t)" = (+)((+)((*)(-1.0, var"v[9](t)"), (*)(-0.0012636230737369685, (exp)((*)(-1, t)))), (*)(2.0, var"v[10](t)"))
var"v[12](t)" = (+)((+)((*)(-0.0013915544448850764, (exp)((*)(-1, t))), (*)(-1.0, var"v[10](t)")), (*)(2.0, var"v[11](t)"))
var"v[13](t)" = (+)((+)((*)(-1.0, var"v[11](t)"), (*)(-0.001515631984321947, (exp)((*)(-1, t)))), (*)(2.0, var"v[12](t)"))
var"v[14](t)" = (+)((+)((*)(-0.0016355120662858347, (exp)((*)(-1, t))), (*)(-1.0, var"v[12](t)")), (*)(2.0, var"v[13](t)"))
var"v[15](t)" = (+)((+)((*)(-1.0, var"v[13](t)"), (*)(-0.0017508626896372112, (exp)((*)(-1, t)))), (*)(2.0, var"v[14](t)"))
var"v[16](t)" = (+)((+)((*)(-0.0018613643973172333, (exp)((*)(-1, t))), (*)(-1.0, var"v[14](t)")), (*)(2.0, var"v[15](t)"))
var"v[17](t)" = (+)((+)((*)(-1.0, var"v[15](t)"), (*)(2.0, var"v[16](t)")), (*)(-0.0019667111610660887, (exp)((*)(-1, t))))
var"v[18](t)" = (+)((+)((*)(-1.0, var"v[16](t)"), (*)(2.0, var"v[17](t)")), (*)(-0.002066611228951034, (exp)((*)(-1, t))))
var"v[19](t)" = (+)((+)((*)(-1.0, var"v[17](t)"), (*)(2.0, var"v[18](t)")), (*)(-0.0021607879333569546, (exp)((*)(-1, t))))
var"v[1](t)" = (+)((+)((*)(-1.0, var"v[3](t)"), (*)(2.0, (getindex)(ˍ₋arg1, 19))), (*)(-0.00014572654664285906, (exp)((*)(-1, t))))
var"u[20](t)" = (+)((+)((*)(-0.3333333333333333, (getindex)(ˍ₋arg1, 17)), (*)(1.3333333333333333, (getindex)(ˍ₋arg1, 18))), (*)(-0.029525297712557772, (exp)((*)(-1, t))))
var"u[1](t)" = (exp)((*)(-1, t))
begin
(SymbolicUtils.Code.create_array)(typeof(ˍ₋arg1), nothing, Val{1}(), Val{(20,)}(), (+)((+)((*)(361.00000000000006, var"u[1](t)"), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 2))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 1))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 1)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 3))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 2))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 2)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 3))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 4))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 3)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 5))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 4))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 4)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 5))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 6))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 5)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 7))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 6))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 6)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 7))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 8))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 7)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 9))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 8))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 8)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 10))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 9))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 9)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 10))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 11))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 10)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 11))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 12))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 11)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 12))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 13))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 12)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 14))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 13))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 13)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 14))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 15))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 14)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 16))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 15))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 15)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 17))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 16))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 16)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 18))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 17))), (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 17)), (*)(361.00000000000006, var"u[20](t)")), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 18))), (+)((+)((+)((*)(361.00000000000006, var"v[18](t)"), (*)(0.8118819450498316, (exp)((*)(-1, t)))), (*)(361.00000000000006, var"v[20](t)")), (*)(-722.0000000000001, var"v[19](t)")), (+)((+)((+)((*)(9.5, var"v[3](t)"), (*)(28.5, var"v[1](t)")), (*)(-38.0, (getindex)(ˍ₋arg1, 19))), (exp)((*)(-1, t))))
end
end
end
end, function (ˍ₋out, ˍ₋arg1, ˍ₋arg2, t)
begin
begin
var"v[20](t)" = (*)(0.8414709848078965, (exp)((*)(-1, t)))
var"v[3](t)" = (+)((+)((*)(0.5, (getindex)(ˍ₋arg1, 19)), (*)(0.5, (getindex)(ˍ₋arg1, 20))), (*)(0.00014552475590851675, (exp)((*)(-1, t))))
var"v[5](t)" = (+)((+)((*)(-1.0, var"v[3](t)"), (*)(2.0, (getindex)(ˍ₋arg1, 20))), (*)(-0.0004355664317500006, (exp)((*)(-1, t))))
var"v[6](t)" = (+)((+)((*)(-1.0, (getindex)(ˍ₋arg1, 20)), (*)(-0.0005788770749661647, (exp)((*)(-1, t)))), (*)(2.0, var"v[5](t)"))
var"v[7](t)" = (+)((+)((*)(-0.0007205845507054019, (exp)((*)(-1, t))), (*)(-1.0, var"v[5](t)")), (*)(2.0, var"v[6](t)"))
var"v[8](t)" = (+)((+)((*)(-1.0, var"v[6](t)"), (*)(2.0, var"v[7](t)")), (*)(-0.0008602964080896811, (exp)((*)(-1, t))))
var"v[9](t)" = (+)((+)((*)(-1.0, var"v[7](t)"), (*)(2.0, var"v[8](t)")), (*)(-0.0009976257229936814, (exp)((*)(-1, t))))
var"v[10](t)" = (+)((+)((*)(-1.0, var"v[8](t)"), (*)(2.0, var"v[9](t)")), (*)(-0.0011321921696093797, (exp)((*)(-1, t))))
var"v[11](t)" = (+)((+)((*)(-1.0, var"v[9](t)"), (*)(-0.0012636230737369685, (exp)((*)(-1, t)))), (*)(2.0, var"v[10](t)"))
var"v[12](t)" = (+)((+)((*)(-0.0013915544448850764, (exp)((*)(-1, t))), (*)(-1.0, var"v[10](t)")), (*)(2.0, var"v[11](t)"))
var"v[13](t)" = (+)((+)((*)(-1.0, var"v[11](t)"), (*)(-0.001515631984321947, (exp)((*)(-1, t)))), (*)(2.0, var"v[12](t)"))
var"v[14](t)" = (+)((+)((*)(-0.0016355120662858347, (exp)((*)(-1, t))), (*)(-1.0, var"v[12](t)")), (*)(2.0, var"v[13](t)"))
var"v[15](t)" = (+)((+)((*)(-1.0, var"v[13](t)"), (*)(-0.0017508626896372112, (exp)((*)(-1, t)))), (*)(2.0, var"v[14](t)"))
var"v[16](t)" = (+)((+)((*)(-0.0018613643973172333, (exp)((*)(-1, t))), (*)(-1.0, var"v[14](t)")), (*)(2.0, var"v[15](t)"))
var"v[17](t)" = (+)((+)((*)(-1.0, var"v[15](t)"), (*)(2.0, var"v[16](t)")), (*)(-0.0019667111610660887, (exp)((*)(-1, t))))
var"v[18](t)" = (+)((+)((*)(-1.0, var"v[16](t)"), (*)(2.0, var"v[17](t)")), (*)(-0.002066611228951034, (exp)((*)(-1, t))))
var"v[19](t)" = (+)((+)((*)(-1.0, var"v[17](t)"), (*)(2.0, var"v[18](t)")), (*)(-0.0021607879333569546, (exp)((*)(-1, t))))
var"v[1](t)" = (+)((+)((*)(-1.0, var"v[3](t)"), (*)(2.0, (getindex)(ˍ₋arg1, 19))), (*)(-0.00014572654664285906, (exp)((*)(-1, t))))
var"u[20](t)" = (+)((+)((*)(-0.3333333333333333, (getindex)(ˍ₋arg1, 17)), (*)(1.3333333333333333, (getindex)(ˍ₋arg1, 18))), (*)(-0.029525297712557772, (exp)((*)(-1, t))))
var"u[1](t)" = (exp)((*)(-1, t))
begin
@inbounds begin
ˍ₋out[1] = (+)((+)((*)(361.00000000000006, var"u[1](t)"), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 2))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 1)))
ˍ₋out[2] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 1)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 3))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 2)))
ˍ₋out[3] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 2)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 3))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 4)))
ˍ₋out[4] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 3)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 5))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 4)))
ˍ₋out[5] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 4)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 5))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 6)))
ˍ₋out[6] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 5)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 7))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 6)))
ˍ₋out[7] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 6)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 7))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 8)))
ˍ₋out[8] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 7)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 9))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 8)))
ˍ₋out[9] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 8)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 10))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 9)))
ˍ₋out[10] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 9)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 10))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 11)))
ˍ₋out[11] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 10)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 11))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 12)))
ˍ₋out[12] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 11)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 12))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 13)))
ˍ₋out[13] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 12)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 14))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 13)))
ˍ₋out[14] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 13)), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 14))), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 15)))
ˍ₋out[15] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 14)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 16))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 15)))
ˍ₋out[16] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 15)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 17))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 16)))
ˍ₋out[17] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 16)), (*)(361.00000000000006, (getindex)(ˍ₋arg1, 18))), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 17)))
ˍ₋out[18] = (+)((+)((*)(361.00000000000006, (getindex)(ˍ₋arg1, 17)), (*)(361.00000000000006, var"u[20](t)")), (*)(-722.0000000000001, (getindex)(ˍ₋arg1, 18)))
ˍ₋out[19] = (+)((+)((+)((*)(361.00000000000006, var"v[18](t)"), (*)(0.8118819450498316, (exp)((*)(-1, t)))), (*)(361.00000000000006, var"v[20](t)")), (*)(-722.0000000000001, var"v[19](t)"))
ˍ₋out[20] = (+)((+)((+)((*)(9.5, var"v[3](t)"), (*)(28.5, var"v[1](t)")), (*)(-38.0, (getindex)(ˍ₋arg1, 19))), (exp)((*)(-1, t)))
nothing
end
end
end
end
end)
var"##tgrad#258" = nothing
var"##jac#259" = nothing
M = [1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0]
ODEFunction{true}(var"##f#257", jac = var"##jac#259", tgrad = var"##tgrad#258", mass_matrix = M, jac_prototype = nothing, syms = [Symbol("u[2](t)"), Symbol("u[3](t)"), Symbol("u[4](t)"), Symbol("u[5](t)"), Symbol("u[6](t)"), Symbol("u[7](t)"), Symbol("u[8](t)"), Symbol("u[9](t)"), Symbol("u[10](t)"), Symbol("u[11](t)"), Symbol("u[12](t)"), Symbol("u[13](t)"), Symbol("u[14](t)"), Symbol("u[15](t)"), Symbol("u[16](t)"), Symbol("u[17](t)"), Symbol("u[18](t)"), Symbol("u[19](t)"), Symbol("v[2](t)"), Symbol("v[4](t)")], indepsym = :t)
end
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment