From fd3dc913b103682e887b2613094870a3897585c5 Mon Sep 17 00:00:00 2001
From: Nicolas Marsic <nicolas.marsic@gmail.com>
Date: Fri, 22 Jun 2012 13:56:05 +0000
Subject: [PATCH] HexNodeBasis looks OK :)
---
FunctionSpace/HexNodeBasis.cpp | 57 ++++++++++++++++++++++++++++------
1 file changed, 47 insertions(+), 10 deletions(-)
diff --git a/FunctionSpace/HexNodeBasis.cpp b/FunctionSpace/HexNodeBasis.cpp
index 0e1ab6f9cc..3c866fca11 100644
--- a/FunctionSpace/HexNodeBasis.cpp
+++ b/FunctionSpace/HexNodeBasis.cpp
@@ -206,7 +206,7 @@ HexNodeBasis::~HexNodeBasis(void){
#include <cstdio>
int main(void){
- const int P = 2;
+ const int P = 3;
const double d = 0.05;
HexNodeBasis b(P);
@@ -236,7 +236,7 @@ int main(void){
printf("r = p(i, 1);\n");
printf("end\n");
printf("\n");
- /*
+
printf("d = %f;\nx = [0:d:1];\ny = x;\nz = x;\n\nlx = length(x);\nly = length(y);\nlz = length(z);\n\n", d);
for(int i = 0; i < b.getSize(); i++)
@@ -251,7 +251,7 @@ int main(void){
for(int i = 0; i < b.getSize(); i++)
printf("p%d(j, i, k) = p(%d, x(i), y(j), z(k));\n", i + 1, i + 1);
- printf("end\n");
+ printf("\nend\n");
printf("end\n");
printf("end\n");
@@ -262,15 +262,52 @@ int main(void){
printf("\n");
for(int i = b.getSize(); i > 0; i--){
printf("figure;\n");
- printf("hold on;\n");
- printf("for k = 1:lz\n");
- printf("surf(x, y, p%d(:, :, k) + z(k));\n", i);
- printf("end\n");
- printf("hold off;\n");
- }
+
+ printf("subplot(3, 2, 1);\n");
+ printf("contourf(x, y, squeeze(p%d(:, :, 1)));\n", i);
+ printf("colorbar;\n");
+ printf("title('z = 0');\n");
+ printf("ylabel('x');\n");
+ printf("xlabel('y');\n\n");
+
+ printf("subplot(3, 2, 2);\n");
+ printf("contourf(x, y, squeeze(p%d(:, :, end)));\n", i);
+ printf("colorbar;\n");
+ printf("title('z = 1');\n");
+ printf("ylabel('x');\n");
+ printf("xlabel('y');\n\n");
+
+ printf("subplot(3, 2, 3);\n");
+ printf("contourf(x, z, squeeze(p%d(:, 1, :)));\n", i);
+ printf("colorbar;\n");
+ printf("title('y = 0');\n");
+ printf("ylabel('x');\n");
+ printf("xlabel('z');\n\n");
+
+ printf("subplot(3, 2, 4);\n");
+ printf("contourf(x, z, squeeze(p%d(:, end, :)));\n", i);
+ printf("colorbar;\n");
+ printf("title('y = 1');\n");
+ printf("ylabel('x');\n");
+ printf("xlabel('z');\n\n");
+
+ printf("subplot(3, 2, 5);\n");
+ printf("contourf(y, z, squeeze(p%d(1, :, :)));\n", i);
+ printf("colorbar;\n");
+ printf("title('x = 0');\n");
+ printf("ylabel('y');\n");
+ printf("xlabel('z');\n\n");
+
+ printf("subplot(3, 2, 6);\n");
+ printf("contourf(y, z, squeeze(p%d(end, :, :)));\n", i);
+ printf("colorbar;\n");
+ printf("title('x = 1');\n");
+ printf("ylabel('y');\n");
+ printf("xlabel('z');\n\n\n");
+}
printf("\n");
- */
+
return 0;
}
*/
--
GitLab