mimeticB

 Returns a m+2 by m+1 one-dimensional mimetic boundary operator

 Parameters:
                k : Order of accuracy
                m : Number of cells
 ----------------------------------------------------------------------------
 SPDX-License-Identifier: GPL-3.0-or-later
 © 2008-2024 San Diego State University Research Foundation (SDSURF).
 See LICENSE file or https://www.gnu.org/licenses/gpl-3.0.html for details.
 ----------------------------------------------------------------------------

0001 function B = mimeticB(k, m)
0002 % Returns a m+2 by m+1 one-dimensional mimetic boundary operator
0003 %
0004 % Parameters:
0005 %                k : Order of accuracy
0006 %                m : Number of cells
0007 % ----------------------------------------------------------------------------
0008 % SPDX-License-Identifier: GPL-3.0-or-later
0009 % © 2008-2024 San Diego State University Research Foundation (SDSURF).
0010 % See LICENSE file or https://www.gnu.org/licenses/gpl-3.0.html for details.
0011 % ----------------------------------------------------------------------------
0012 
0013     Q = sparse(diag(weightsQ(k, m, 1)));
0014     D = div(k, m, 1);
0015     G = grad(k, m, 1);
0016     P = sparse(diag(weightsP(k, m, 1)));
0017     
0018     B = Q*D + G'*P;
0019 end