This is a static copy of a profile reportHome
integrateHGsubmatrices_GLQ (255 calls, 1.960 sec)
Generated 18-Mar-2011 23:31:44 using cpu time.
M-function in file /Users/Robert/Documents/Work/Cardiff/postdoc/myrepo/NURBS/isoBEM/integrateHGsubmatrices_GLQ.m
Copy to new window for comparing multiple runs
Parents (calling functions)
No parent Lines where the most time was spent
Line Number | Code | Calls | Total Time | % Time | Time Plot |
6 | [gpt gwt]=lgwt(ngp,-1,1); | 255 | 0.705 s | 36.0% |  |
31 | [Ttemp(i,j) Utemp(i,j)]=DBIEke... | 6120 | 0.653 s | 33.3% |  |
22 | [N(lclBasis) dN(lclBasis)]=NUR... | 4590 | 0.238 s | 12.2% |  |
25 | [jacob_xi, normals, r, dr, drd... | 1530 | 0.156 s | 7.9% |  |
36 | Gsubmatrix=Gsubmatrix + [N(1)*... | 1530 | 0.031 s | 1.6% |  |
All other lines | | | 0.176 s | 9.0% |  |
Totals | | | 1.960 s | 100% | |
Children (called functions)
M-Lint results
No M-Lint messages.Coverage results
[ Show coverage for parent directory ]
Total lines in function | 39 |
Non-code lines (comments, blank lines) | 14 |
Code lines (lines that can run) | 25 |
Code lines that did run | 25 |
Code lines that did not run | 0 |
Coverage (did run/can run) | 100.00 % |
Function listing
time calls line
1 function [ Hsubmatrix, Gsubmatrix ] = integrateHGsubmatrices_GLQ( ngp, elcoords, bsFnConn, collocCoords, range)
2 % integrate the H and G submatrices using Gauss Legendre quadrature
3
0.02 255 4 global controlPts p knotVec
5
0.71 255 6 [gpt gwt]=lgwt(ngp,-1,1);
7
255 8 Hsubmatrix=zeros(2,6);
255 9 Gsubmatrix=zeros(2,6);
10
255 11 numBasisFns=length(bsFnConn);
255 12 N=zeros(1,numBasisFns); dN=zeros(1,numBasisFns);
13
255 14 jacob_param=(range(2)-range(1)) / 2; % jacobian from parent to parameter space
15
255 16 for pt=1:ngp % integrate using Gaussian quadrature
17
0.02 1530 18 xi_param=convertToParamSpace( gpt(pt), range); % get the gauss point in parameter space
19
1530 20 for lclBasis=1:numBasisFns
4590 21 i=bsFnConn(lclBasis);
0.24 4590 22 [N(lclBasis) dN(lclBasis)]=NURBSbasis(i, p, xi_param, knotVec, controlPts(:,3)' );
0.01 4590 23 end
24
0.16 1530 25 [jacob_xi, normals, r, dr, drdn] = getKernelParameters( elcoords, collocCoords, N, dN );
1530 26 jacob=jacob_xi*jacob_param; % the final jacobian we use
27
1530 28 Ttemp=zeros(2,2); Utemp=zeros(2,2);
0.01 1530 29 for i=1:2
0.02 3060 30 for j=1:2
0.65 6120 31 [Ttemp(i,j) Utemp(i,j)]=DBIEkernels(i,j,r,dr,drdn,normals);
0.03 6120 32 end
0.03 3060 33 end
34
1530 35 Hsubmatrix=Hsubmatrix + [N(1)*Ttemp N(2)*Ttemp N(3)*Ttemp]*jacob*gwt(pt);
0.03 1530 36 Gsubmatrix=Gsubmatrix + [N(1)*Utemp N(2)*Utemp N(3)*Utemp]*jacob*gwt(pt);
37
0.02 1530 38 end
255 39 end
Other subfunctions in this file are not included in this listing.