|
Revision 495
(checked in by cbc, 12 years ago)
|
Initial import of Stark code.
|
| Line | |
|---|
| 1 |
for {i=0;i<nx+1;i++){/* waitbar((i-1)/nx) */ |
|---|
| 2 |
for {j=0;j<ny+1;j++){ |
|---|
| 3 |
if ~isnan(JXY(i,j)){ |
|---|
| 4 |
zref=FDZGRID(i,j,:); |
|---|
| 5 |
for {k=0;k<nz+1;k++){ |
|---|
| 6 |
[b2,n2,b1,n1]=basis1d2(zref,Z3D(i,j,k)); |
|---|
| 7 |
N1_3D(i,j,k)=n1; |
|---|
| 8 |
B1_3D(i,j,k)=b1; |
|---|
| 9 |
N2_3D(i,j,k)=n2; |
|---|
| 10 |
B2_3D(i,j,k)=b2; |
|---|
| 11 |
} |
|---|
| 12 |
} |
|---|
| 13 |
} |
|---|
| 14 |
} |
|---|
| 15 |
|
|---|
| 16 |
/* Set pointer to lhs */ |
|---|
| 17 |
start_of_array=(unsigned char *)mxGetPr(S); |
|---|
| 18 |
bytes_to_copy=MNP*mxGetElementSize(S); |
|---|
| 19 |
memcpy(start_of_array,sfd,bytes_to_copy); |
|---|
| 20 |
plhs[0]=S; |
|---|
| 21 |
return; |
|---|
