Commit Graph

906 Commits

Author SHA1 Message Date
Bård Skaflestad
9f99ac3265 Disable diagnostic output. 2010-08-09 11:04:26 +00:00
Bård Skaflestad
1e1788dcba Uhhh, not quite (pointy hat: bska). 2010-08-09 09:32:56 +00:00
Bård Skaflestad
cafb4b5f93 Allocate correct amount of space... 2010-08-09 09:32:12 +00:00
Bård Skaflestad
a39f74ad40 Use fixed box geometry
(tests G.faces.areas ~= 1 and G.cells.volumes ~= 1).
2010-08-09 09:19:09 +00:00
Bård Skaflestad
e07d447558 Add Ubuntu's default install directory for libsuitesparse-dev package. 2010-08-09 09:13:56 +00:00
Bård Skaflestad
b558bf929f Update documentation and example to reflect reality.
Compile in optimized mode.
2010-08-09 08:56:02 +00:00
Bård Skaflestad
32a3f8462e Update inner product call (mex_ip_simple) to reflect
nconn ~= DIFF(G.cells.facePos).
2010-08-09 08:54:55 +00:00
Bård Skaflestad
803511db6e Update documentation to reflect reality. 2010-08-09 08:52:39 +00:00
Bård Skaflestad
b196437776 Add support for number of connections different from number of faces.
This is a precursor to supporting wells as faces.
2010-08-09 08:51:01 +00:00
Jostein R. Natvig
cd829252ca Additiona and modifications to mimetic assembly routines in C. Add
code to call libmwumfpack from mex code.
2010-08-09 08:02:56 +00:00
Jostein R. Natvig
7e3fe07012 Add system assembly. 2010-08-06 10:53:07 +00:00
Jostein R. Natvig
f7e93692f9 free->mxFree. createNumericMatrix->createDoubleMatrix. 2010-08-05 12:16:15 +00:00
Jostein R. Natvig
10582ea32a Convert to using MRST API. 2010-08-05 12:14:39 +00:00
Jostein R. Natvig
97565f3bb7 Return pointers from getXX(...). Transpose centroids and normals.
And coordinates.
2010-08-05 12:13:59 +00:00
Jostein R. Natvig
53de3c03d1 Some cleanup. 'Transpose' getFaceCellNeighbors. 2010-08-05 07:33:59 +00:00
Bård Skaflestad
0d0f8eb69c The geometry routines are presently in 'geometry.c', not
'mimetic_geometry.c'.  Reflect this in the build command.
2010-08-04 20:35:32 +00:00
Bård Skaflestad
02f6a33282 Assign cell volumes to '.cells' sub-structure. Having volumes in the
'.faces' sub-structure is slightly counter intuitive.
2010-08-04 20:33:57 +00:00
Bård Skaflestad
2b943b8f3a Don't spam command window with verbose build messages. While useful
in debugging, this has no place in production.
2010-08-04 19:46:25 +00:00
Bård Skaflestad
f75b20430e Expand SVN keywords. 2010-08-04 16:28:45 +00:00
Bård Skaflestad
98fcce29f6 Add MEX wrapper for computing pressure and fluxes given interface
pressures and output from 'mex_schur_comp_symm'.
2010-08-04 16:27:13 +00:00
Jostein R. Natvig
de74b2565c Add geometry processing in C. 2010-08-04 14:11:37 +00:00
Jostein R. Natvig
a118c2f8ef A first attempt at a C API to MRST grids. 2010-08-04 14:10:27 +00:00
Bård Skaflestad
a2af73ea08 Correct typo in comment. 2010-08-04 12:29:49 +00:00
Bård Skaflestad
235f534376 Appease MLINT. 2010-08-04 11:51:58 +00:00
Bård Skaflestad
53489b4e72 Add MEX wrapper (and test script) for computing the hybrid system
component matrices using compiled C code.
2010-08-04 11:44:48 +00:00
Bård Skaflestad
2f30e6f0c5 Initial implementation of hybrid system infrastructure. Actual
assembly and system solve not currently implemented.  This is untested.
2010-08-03 17:01:33 +00:00
Bård Skaflestad
2cb51aac48 Don't require build(er|system) to #define COMPILING_FOR_MATLAB in
order to compile this file.
2010-08-03 14:27:13 +00:00
Bård Skaflestad
287d230f8a Ignore anything named *.mexa64 as these files are compiler output. 2010-07-09 15:36:29 +00:00
Bård Skaflestad
4ea66c06b1 Remove macro definition which is no longer needed. 2010-07-09 10:39:26 +00:00
Bård Skaflestad
eccc63d566 Assert copyright. Forgotten in previous revision. 2010-07-09 10:16:55 +00:00
Bård Skaflestad
945253df6b Remove statements disabled in earlier revision. 2010-07-09 09:34:29 +00:00
Bård Skaflestad
81d00cc324 Add documentation and building quirk (-DCOMPILING_FOR_MATLAB=1). 2010-07-09 09:32:21 +00:00
Bård Skaflestad
5024034657 Use mxMalloc() and mxFree() over their ANSI C counterparts in order
to register memory usage with the M memory manager (e.g., for Ctrl-C
  or error purposes).
2010-07-09 09:31:26 +00:00
Bård Skaflestad
5cc657c89e Add first fully functioning M-callable C implementation of the
'ip_simple' mimetic inner product.  Coincides with the results from
  'computeMimeticIP' to (roughly) the order of round-off on a
  non-trivial test case.

  The calling interface is

      BI = mex_ip_simple(G, rock)

  and some latitude has been extended towards multiple data types in
  the various G fields (e.g., G.cells.facePos may be an 'int32').

  Additional clean-up and optimisation is likely.
2010-07-08 23:28:33 +00:00
Bård Skaflestad
fc5bf7f555 Provide call interface to MEX routine for easier M testing. 2010-07-04 22:24:37 +00:00
Bård Skaflestad
9eca2f1e92 Implement cell loop (tentatively named mim_ip_simple_all()) and a
MEX gateway for easy testing from M.  Builds with fairly strict
  warnings, but is not tested yet.
2010-07-04 21:44:04 +00:00
Bård Skaflestad
49242b7e2f The LAPACK and BLAS operators do not modify their (char*)
parameters.  Declare these operators as taking (const char*)s.
2010-07-03 11:07:56 +00:00
Bård Skaflestad
05d4a5fdef Correct a few embarrasing matrix dimension errors (in the "linear
pressure" part).
2010-07-02 22:19:38 +00:00
Bård Skaflestad
56a7db04d3 Split mim_ip_simple() up into two co-operating procedures,
- mim_ip_span_nullspace()
        Constructs orthogonal basis for mimetic IP null space.
    - mim_ip_linpress_exact()
        Adds term to impose exactness for linear pressure fields.

  Re-implement mim_ip_simple() in terms of these procedures.
2010-06-29 21:18:23 +00:00
Bård Skaflestad
01edb43f47 Allow builder to -Define expansion of MAT_SIZE_T . 2010-06-29 19:11:26 +00:00
Bård Skaflestad
c5379fb191 Appease the -pedantic -ansi gods. 2010-06-29 15:32:03 +00:00
Bård Skaflestad
08d53ece42 Don't fall prey to refactoring errors. 2010-06-28 22:54:31 +00:00
Bård Skaflestad
5441bc3efe For kicks and giggles, add an (untested) C+BLAS/LAPACK
implementation of the 'ip_simple' mimetic inner product.

  Suggested by: jrn.
2010-06-28 22:47:55 +00:00
Bård Skaflestad
ee9a8fdc59 Implement G.type for OPM processing. 2010-03-18 16:31:10 +00:00
Bård Skaflestad
b3f61acd81 Assert Copyright for 2010 in files missed in OPM -r762. 2010-03-18 15:57:02 +00:00
Jostein R. Natvig
89128edfac Add hack to reverse sign in ZCORN if ZCORN decrease along pillar.
The sign of z-coordinates in the processed grid is restored to ensure
right-handedness in relation between face geometry and face-cell topology.
2010-03-03 06:49:51 +00:00
Jostein R. Natvig
d86ba9ab8d Add sanity chech of input. 2010-01-15 13:20:25 +00:00
Atgeirr Flø Rasmussen
97718768fd Changed to ensure a convex combination is taken when points
are between pillar extremes. Original code kept in #if 0 block.
2010-01-06 09:22:42 +00:00
Jostein R. Natvig
f795a753ed Catch up with change in grid structure. 2010-01-05 15:43:59 +00:00
Jostein R. Natvig
d5c900de11 *Change some fields from double to int.
*Add copyright.

Lacking functionality to split disconnected grid.
2009-12-16 13:49:14 +00:00