Skip to content

Commit

Permalink
Removed the B Matrix changes...?
Browse files Browse the repository at this point in the history
  • Loading branch information
@sap8b
sap8b committed Jun 19, 2021
1 parent 1406adf commit c18270e
Show file tree
Hide file tree
Showing 2 changed files with 48 additions and 44 deletions.
88 changes: 44 additions & 44 deletions DiffusionSimulators_2D_MatrixD.cs
View file
Original file line number Diff line number Diff line change
Expand Up @@ -930,28 +930,28 @@ public void TwoD_ADI(int n_time_steps, int output_interval)

RVector f12s = f12.GetColVector(j); //

switch (BCs[0].TypeBC)
{
case BoundaryConditions.Dirichlet:
B[j][0, 0] = 1.0;
B[j][0, 1] = 0.0;
break;
default:
B[j][0, 0] = 1.0;
B[j][0, 1] = 0.0;
break;
}
switch (BCs[3].TypeBC)
{
case BoundaryConditions.Dirichlet:
B[j][ncols - 1, ncols - 1] = 1.0;
B[j][ncols - 1, ncols - 2] = 0.0;
break;
default:
B[j][ncols - 1, ncols - 1] = 1.0;
B[j][ncols - 1, ncols - 2] = 0.0;
break;
}
//switch (BCs[0].TypeBC)
//{
// case BoundaryConditions.Dirichlet:
// B[j][0, 0] = 1.0;
// B[j][0, 1] = 0.0;
// break;
// default:
// B[j][0, 0] = 1.0;
// B[j][0, 1] = 0.0;
// break;
//}
//switch (BCs[3].TypeBC)
//{
// case BoundaryConditions.Dirichlet:
// B[j][ncols - 1, ncols - 1] = 1.0;
// B[j][ncols - 1, ncols - 2] = 0.0;
// break;
// default:
// B[j][ncols - 1, ncols - 1] = 1.0;
// B[j][ncols - 1, ncols - 2] = 0.0;
// break;
//}
RVector u12 = TridiagonalMatrix.Thomas_Algorithm(B[j], v1 + f12s);
C_Im1.ReplaceCol(u12, j);
}
Expand Down Expand Up @@ -1030,28 +1030,28 @@ public void TwoD_ADI(int n_time_steps, int output_interval)
b[0] = CL[k]; //nu *
b[ncols - 1] = CR[k]; //nu *

switch (BCs[0].TypeBC)
{
case BoundaryConditions.Dirichlet:
B[k][0, 0] = 1.0;
B[k][0, 1] = 0.0;
break;
default:
B[k][0, 0] = 1.0;
B[k][0, 1] = 0.0;
break;
}
switch (BCs[3].TypeBC)
{
case BoundaryConditions.Dirichlet:
B[k][ncols - 1, ncols - 1] = 1.0;
B[k][ncols - 1, ncols - 2] = 0.0;
break;
default:
B[k][ncols - 1, ncols - 1] = 1.0;
B[k][ncols - 1, ncols - 2] = 0.0;
break;
}
//switch (BCs[0].TypeBC)
//{
// case BoundaryConditions.Dirichlet:
// B[k][0, 0] = 1.0;
// B[k][0, 1] = 0.0;
// break;
// default:
// B[k][0, 0] = 1.0;
// B[k][0, 1] = 0.0;
// break;
//}
//switch (BCs[3].TypeBC)
//{
// case BoundaryConditions.Dirichlet:
// B[k][ncols - 1, ncols - 1] = 1.0;
// B[k][ncols - 1, ncols - 2] = 0.0;
// break;
// default:
// B[k][ncols - 1, ncols - 1] = 1.0;
// B[k][ncols - 1, ncols - 2] = 0.0;
// break;
//}
RVector u1 = TridiagonalMatrix.Thomas_Algorithm(B[k], b);
u1[0] = CL[k]; //nu *
u1[ncols - 1] = CR[k]; //nu *
Expand Down
4 changes: 4 additions & 0 deletions log.txt
View file
Original file line number Diff line number Diff line change
Expand Up @@ -37,3 +37,7 @@
{"message":"Invalid file link:(~/2D_Region.png).","source":"BuildCommand.BuildCore.Build Document.LinkPhaseHandlerWithIncremental.ConceptualDocumentProcessor.Save","file":"README.md","line":"0","date_time":"2021-06-18T18:39:35.936992Z","message_severity":"warning","code":"InvalidFileLink"}
{"message":"Invalid file link:(~/2D_Region.png).","source":"BuildCommand.BuildCore.Build Document.LinkPhaseHandlerWithIncremental.ConceptualDocumentProcessor.Save","file":"README.md","line":"0","date_time":"2021-06-18T18:54:21.9024564Z","message_severity":"warning","code":"InvalidFileLink"}
{"message":"Invalid file link:(~/2D_Region.png).","source":"BuildCommand.BuildCore.Build Document.LinkPhaseHandlerWithIncremental.ConceptualDocumentProcessor.Save","file":"README.md","line":"0","date_time":"2021-06-18T18:54:40.3941345Z","message_severity":"warning","code":"InvalidFileLink"}
{"message":"Invalid file link:(~/2D_Region.png).","source":"BuildCommand.BuildCore.Build Document.LinkPhaseHandlerWithIncremental.ConceptualDocumentProcessor.Save","file":"README.md","line":"0","date_time":"2021-06-19T12:31:22.1021961Z","message_severity":"warning","code":"InvalidFileLink"}
{"message":"Invalid file link:(~/2D_Region.png).","source":"BuildCommand.BuildCore.Build Document.LinkPhaseHandlerWithIncremental.ConceptualDocumentProcessor.Save","file":"README.md","line":"0","date_time":"2021-06-19T12:31:32.8836769Z","message_severity":"warning","code":"InvalidFileLink"}
{"message":"Invalid file link:(~/2D_Region.png).","source":"BuildCommand.BuildCore.Build Document.LinkPhaseHandlerWithIncremental.ConceptualDocumentProcessor.Save","file":"README.md","line":"0","date_time":"2021-06-19T13:29:39.9687672Z","message_severity":"warning","code":"InvalidFileLink"}
{"message":"Invalid file link:(~/2D_Region.png).","source":"BuildCommand.BuildCore.Build Document.LinkPhaseHandlerWithIncremental.ConceptualDocumentProcessor.Save","file":"README.md","line":"0","date_time":"2021-06-19T13:29:51.4077891Z","message_severity":"warning","code":"InvalidFileLink"}

0 comments on commit c18270e

Please sign in to comment.