Homec4science
Diffusion COSOlver (new_coulomb) Browse

cosolver/matlab/debug_speed_functions.nb
41e7a23b16aenew_coulomb

Last Change
Raw File

debug_speed_functions.nb
View Options

(* Content-type: application/vnd.wolfram.mathematica *)
(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)
(* CreatedBy='Mathematica 12.1' *)
(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[ 158, 7]
NotebookDataLength[ 11394, 318]
NotebookOptionsPosition[ 9724, 280]
NotebookOutlinePosition[ 10118, 296]
CellTagsIndexPosition[ 10075, 293]
WindowFrame->Normal*)
(* Beginning of Notebook Content *)
Notebook[{
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"\[IndentingNewLine]",
RowBox[{"ClearAll", "\[IndentingNewLine]",
RowBox[{
RowBox[{"k", "=", "1"}], ";"}], "\[IndentingNewLine]",
RowBox[{"N", "[",
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{
RowBox[{"Erf", "[", " ", "x", "]"}], " ",
SuperscriptBox["x",
RowBox[{
RowBox[{"2", " ", "k"}], " ", "+", "1", " "}]],
RowBox[{"Exp", "[",
RowBox[{"-",
SuperscriptBox["x", "2"]}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "0", ",", "Infinity"}], "}"}]}], "]"}], "]"}],
"\[IndentingNewLine]",
RowBox[{"N", "[",
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{
FractionBox["2",
RowBox[{"Sqrt", "[", "\[Pi]", "]"}]], " ",
RowBox[{"Exp", "[",
RowBox[{"-",
SuperscriptBox["x", "2"]}], "]"}], " ",
SuperscriptBox["x",
RowBox[{"2", " ", "k"}]],
RowBox[{"Exp", "[",
RowBox[{"-",
SuperscriptBox["x", "2"]}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "0", ",", "Infinity"}], "}"}]}], "]"}], "]"}],
"\[IndentingNewLine]"}]}]], "Input",
CellChangeTimes->CompressedData["
1:eJwdyU0og3EAx/Gn0UqhzRprpDZlHBzszQ5adtmBYknMTmzJlidbTmNelsNy
sYsTUgqNtrLkMi8rNSsTq6XN27LUZtMykbUo/H/P4dvn8BUZrf1jLIqihCTo
WpBeTRvyGqEnHoU9ZYkEPFSnnmET2/oCS5qVd1i9o2PNEC1VD+XQ6b4QwsW0
Xwxrf05DUCZpPIeu3NEj5G84ktC/Z8xBOU/3CuMqM9dBjJx118BAgVMPld7t
BlhBh+WwuXOgA7Ycr/VCFSfMuJ7uG2E+vWyCRcvbNRze+mT82AzEoJZeuoG2
QvIXcjJTf5A1TrNnieYJrQTavielMHOnl0HqyTg0R9yvuzTAlH13FAY9ChPs
4ptC0KW2R+B9pTMKlbzWGNTpvdx54kneLYCioq8Nlg78jLeCbDv88okVzA/m
s6tERWowB/8BInPSjg==
"],
CellLabel->"In[6]:=",ExpressionUUID->"ba9b04bc-1a20-4a44-a090-9a16874068e6"],
Cell[BoxData["ClearAll"], "Output",
CellChangeTimes->{{3.805540766349869*^9, 3.805540781140292*^9},
3.805597515004232*^9},
CellLabel->"Out[6]=",ExpressionUUID->"be376b5e-a26a-47cb-8e5f-dae035dbab87"],
Cell[BoxData["0.44194173824159216`"], "Output",
CellChangeTimes->{{3.805540766349869*^9, 3.805540781140292*^9},
3.8055975151453943`*^9},
CellLabel->"Out[8]=",ExpressionUUID->"b7464d1e-93a5-4576-83aa-0d9332f5a09e"],
Cell[BoxData["0.17677669529663687`"], "Output",
CellChangeTimes->{{3.805540766349869*^9, 3.805540781140292*^9},
3.805597515252719*^9},
CellLabel->"Out[9]=",ExpressionUUID->"0b65ed30-0896-4a7f-9819-fcac4d127532"]
}, Open ]],
Cell[CellGroupData[{
Cell["Test speed integrals of collision frequency", "Chapter",
CellChangeTimes->{{3.8055407219448757`*^9,
3.8055407470245028`*^9}},ExpressionUUID->"ee21ce83-6aa4-4257-8b1c-\
8197fd57ee12"],
Cell[CellGroupData[{
Cell[BoxData[{
RowBox[{
RowBox[{
RowBox[{"\[Nu]D", "[", "s_", "]"}], " ", ":=", " ",
RowBox[{
FractionBox["1",
RowBox[{
SuperscriptBox["s", "3"], " "}]],
RowBox[{"(",
RowBox[{
RowBox[{
RowBox[{"(",
RowBox[{"1", " ", "-", " ",
FractionBox["1",
RowBox[{"2", " ",
SuperscriptBox["s", "2"]}]]}], ")"}], " ",
RowBox[{"Erf", "[", "s", "]"}]}], " ", "+", " ",
RowBox[{
FractionBox["1",
RowBox[{"2", " ", "s", " "}]], " ",
RowBox[{"D", "[",
RowBox[{
RowBox[{"Erf", "[", "s", "]"}], ",", "s"}], "]"}]}]}], ")"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"\[Nu]SD", "[", "s_", "]"}], " ", ":=", " ",
RowBox[{
FractionBox[
RowBox[{"4", " "}], "s"], " ",
FractionBox["1",
RowBox[{"2", " ",
SuperscriptBox["s", "2"]}]],
RowBox[{"(", " ",
RowBox[{
RowBox[{"Erf", "[", " ", "s", "]"}], " ", "-", " ",
RowBox[{"s", " ",
RowBox[{"D", "[",
RowBox[{
RowBox[{"Erf", "[", "s", "]"}], ",", "s"}], "]"}]}]}], ")"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"\[Nu]Para", "[", "s_", "]"}], " ", ":=", " ",
RowBox[{
FractionBox["2",
SuperscriptBox["s", "3"]], " ",
FractionBox["1",
RowBox[{"2", " ",
SuperscriptBox["s", "2"]}]],
RowBox[{"(",
RowBox[{
RowBox[{"Erf", "[", "s", "]"}], " ", "-", " ",
RowBox[{"s", " ",
RowBox[{"D", "[",
RowBox[{
RowBox[{"Erf", "[", "s", "]"}], ",", "s"}], "]"}]}]}], ")"}]}]}],
";"}], "\[IndentingNewLine]",
RowBox[{
RowBox[{
RowBox[{"k", " ", "=", " ", "0"}], ";"}], "\[IndentingNewLine]",
RowBox[{"(*",
RowBox[{"N", "[",
RowBox[{
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{
RowBox[{"\[Nu]D", "[", "x", "]"}],
SuperscriptBox["x",
RowBox[{"2", " ", "k"}]],
RowBox[{"Exp", "[",
RowBox[{"-",
SuperscriptBox["x", "2"]}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "0", ",", "Infinity"}], "}"}]}], "]"}], ",", "60"}],
"]"}], "*)"}]}], "\[IndentingNewLine]",
RowBox[{
RowBox[{"(*",
RowBox[{"N", "[",
RowBox[{
RowBox[{"Integrate", "[",
RowBox[{
RowBox[{
RowBox[{"\[Nu]SD", "[", "x", "]"}],
SuperscriptBox["x",
RowBox[{"2", " ", "k"}]],
RowBox[{"Exp", "[",
RowBox[{"-",
SuperscriptBox["x", "2"]}], "]"}]}], ",",
RowBox[{"{",
RowBox[{"x", ",", "0", ",", "Infinity"}], "}"}]}], "]"}], ",", "60"}],
"]"}], "*)"}]}]}], "Input",
CellChangeTimes->{{3.805540802707239*^9, 3.805540871127615*^9}, {
3.805540929672986*^9, 3.805540931624302*^9}, {3.805541686686306*^9,
3.805541708289136*^9}, {3.805541787769302*^9, 3.805541866997121*^9}, {
3.8055419033207693`*^9, 3.805541935699361*^9}, {3.80554196655737*^9,
3.805541970782301*^9}, {3.805544189652068*^9, 3.805544209799753*^9}, {
3.805544280933981*^9, 3.805544293475382*^9}, {3.805544399495862*^9,
3.80554440747244*^9}, 3.805544488086038*^9, {3.805544525517054*^9,
3.805544528718832*^9}, {3.805544595990273*^9, 3.80554459746484*^9}, {
3.805544783651104*^9, 3.805544808321331*^9}, {3.805544984396348*^9,
3.805544984485919*^9}, {3.8055466746876993`*^9, 3.805546678633816*^9}, {
3.8055467424437838`*^9, 3.805546760474421*^9}, {3.805546846041645*^9,
3.805546898942337*^9}, {3.8055470187748823`*^9, 3.8055470229592943`*^9}, {
3.805547138133713*^9, 3.805547165585322*^9}, {3.805547651003269*^9,
3.805547656254814*^9}, {3.805547756181013*^9, 3.805547821869564*^9}, {
3.805548479069907*^9, 3.805548479152157*^9}, {3.805548532597632*^9,
3.8055485449708776`*^9}, {3.8055488057500467`*^9, 3.805548809844821*^9}, {
3.8055999444962063`*^9, 3.805599948169345*^9}, {3.8055999875364027`*^9,
3.8055999890733547`*^9}, {3.805600026058917*^9, 3.8056000513153667`*^9}, {
3.805600261712356*^9, 3.805600292812903*^9}, {3.805600702358691*^9,
3.805600704935557*^9}},ExpressionUUID->"212634bb-b305-45e1-9a8a-\
4e904931f545"],
Cell[BoxData["1.\
0656799507071040471381587984598115390830230942306253248467682586747105883193`\
60."], "Output",
CellChangeTimes->{{3.805546754072839*^9, 3.8055467738392572`*^9}, {
3.805546860270019*^9, 3.805546912258359*^9}, {3.805547019653969*^9,
3.8055470363599358`*^9}, {3.805547154444459*^9, 3.805547179292563*^9},
3.805547693734625*^9, {3.805547757534422*^9, 3.805547760463541*^9},
3.8055478241202393`*^9, 3.8055484936057997`*^9, {3.805548551302459*^9,
3.8055485518428183`*^9}, {3.805548807231986*^9, 3.805548824291758*^9},
3.8055974769325533`*^9, 3.805597529412094*^9, {3.805600010656227*^9,
3.805600060257059*^9}},
CellLabel->
"Out[359]=",ExpressionUUID->"f4397ada-13a9-4771-a47c-9ff75547f728"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{"NumberForm", "[",
RowBox[{"1.0656799507141321`", ",", "16"}], "]"}]], "Input",
NumberMarks->False,
CellLabel->
"In[354]:=",ExpressionUUID->"b7fede7f-bcb5-492a-bfa8-7ee7fbc84622"],
Cell[BoxData[
TagBox[
InterpretationBox["\<\"1.065679950714132\"\>",
1.0656799507141321`,
AutoDelete->True],
NumberForm[#, 16]& ]], "Output",
CellChangeTimes->{3.805600046970869*^9},
CellLabel->
"Out[354]//NumberForm=",ExpressionUUID->"92c5178e-e8a9-4682-9f12-\
9c62b51c6952"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[""], "Input",
CellChangeTimes->{{3.805547078294024*^9, 3.8055471358988657`*^9}},
NumberMarks->False,ExpressionUUID->"bda3ed77-9b1d-411f-adbb-24c845176628"],
Cell[BoxData[
TagBox[
InterpretationBox["\<\"0.6149535993032899\"\>",
0.6149535993032899,
AutoDelete->True],
NumberForm[#, 19]& ]], "Output",
CellChangeTimes->{{3.8055470726856813`*^9, 3.805547111094249*^9}},
CellLabel->
"Out[1704]//NumberForm=",ExpressionUUID->"e9de6471-262b-40a9-8fb2-\
b89d61ca9676"]
}, Open ]],
Cell[CellGroupData[{
Cell[BoxData[
RowBox[{
RowBox[{"ResourceFunction", "[", "\"\<NthDigit\>\"", "]"}], "[",
RowBox[{"0.6149535993032899`", ",", "32"}], "]"}]], "Input",
NumberMarks->False,
CellLabel->
"In[1705]:=",ExpressionUUID->"356670b0-73a2-4495-b405-41411783c585"],
Cell[BoxData["3"], "Output",
CellChangeTimes->{3.805547126374133*^9},
CellLabel->
"Out[1705]=",ExpressionUUID->"fe5727bc-cb78-45a4-9e09-57a0caf04a10"]
}, Open ]]
}, Open ]]
},
WindowSize->{808, 755},
WindowMargins->{{53, Automatic}, {Automatic, 16}},
FrontEndVersion->"12.1 for Mac OS X x86 (64-bit) (March 18, 2020)",
StyleDefinitions->"Default.nb",
ExpressionUUID->"6f44b032-b0b5-4eeb-b0d0-af73b2ffebad"
]
(* End of Notebook Content *)
(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 1747, 46, 160, "Input",ExpressionUUID->"ba9b04bc-1a20-4a44-a090-9a16874068e6"],
Cell[2330, 70, 205, 3, 34, "Output",ExpressionUUID->"be376b5e-a26a-47cb-8e5f-dae035dbab87"],
Cell[2538, 75, 219, 3, 34, "Output",ExpressionUUID->"b7464d1e-93a5-4576-83aa-0d9332f5a09e"],
Cell[2760, 80, 217, 3, 34, "Output",ExpressionUUID->"0b65ed30-0896-4a7f-9819-fcac4d127532"]
}, Open ]],
Cell[CellGroupData[{
Cell[3014, 88, 192, 3, 69, "Chapter",ExpressionUUID->"ee21ce83-6aa4-4257-8b1c-8197fd57ee12"],
Cell[CellGroupData[{
Cell[3231, 95, 4207, 112, 204, "Input",ExpressionUUID->"212634bb-b305-45e1-9a8a-4e904931f545"],
Cell[7441, 209, 734, 12, 34, "Output",ExpressionUUID->"f4397ada-13a9-4771-a47c-9ff75547f728"]
}, Open ]],
Cell[CellGroupData[{
Cell[8212, 226, 209, 5, 30, "Input",ExpressionUUID->"b7fede7f-bcb5-492a-bfa8-7ee7fbc84622"],
Cell[8424, 233, 291, 9, 45, "Output",ExpressionUUID->"92c5178e-e8a9-4682-9f12-9c62b51c6952"]
}, Open ]],
Cell[CellGroupData[{
Cell[8752, 247, 170, 2, 30, "Input",ExpressionUUID->"bda3ed77-9b1d-411f-adbb-24c845176628"],
Cell[8925, 251, 318, 9, 45, "Output",ExpressionUUID->"e9de6471-262b-40a9-8fb2-b89d61ca9676"]
}, Open ]],
Cell[CellGroupData[{
Cell[9280, 265, 259, 6, 30, "Input",ExpressionUUID->"356670b0-73a2-4495-b405-41411783c585"],
Cell[9542, 273, 154, 3, 34, "Output",ExpressionUUID->"fe5727bc-cb78-45a4-9e09-57a0caf04a10"]
}, Open ]]
}, Open ]]
}
]
*)
c4science · Help