PMVS三维重建

基于PMVS算法的三维重建实现，使用c++编写，简练直观，可得到较好的3D图像

/* caxpy.f -- translated by f2c （version 20061008）.
   You must link the resulting object file with libf2c:
	on Microsoft Windows system link with libf2c.lib;
	on Linux or Unix systems link with .../path/to/libf2c.a -lm
	or if you install libf2c.a in a standard place with -lf2c -lm
	-- in that order at the end of the command line as in
		cc *.o -lf2c -lm
	Source for libf2c is in /netlib/f2c/libf2c.zip e.g.

		http://www.netlib.org/f2c/libf2c.zip
*/

#include “f2c.h“
#include “blaswrap.h“

/* Subroutine */ int caxpy_（integer *n complex *ca complex *cx integer *
	incx complex *cy integer *incy）
{
    /* System generated locals */
    integer i__1 i__2 i__3 i__4;
    complex q__1 q__2;

    /* Local variables */
    integer i__ ix iy;
    extern doublereal scabs1_（complex *）;

/*     .. Scalar Arguments .. */
/*     .. */
/*     .. Array Arguments .. */
/*     .. */

/*  Purpose */
/*  ======= */

/*     CAXPY constant times a vector plus a vector. */

/*  Further Details */
/*  =============== */

/*     jack dongarra linpack 3/11/78. */
/*     modified 12/3/93 array（1） declarations changed to array（*） */

/*     .. Local Scalars .. */
/*     .. */
/*     .. External Functions .. */
/*     .. */
    /* Parameter adjustments */
    --cy;
    --cx;

    /* Function Body */
    if （*n <= 0） {
	return 0;
    }
    if （scabs1_（ca） == 0.f） {
	return 0;
    }
    if （*incx == 1 && *incy == 1） {
	goto L20;
    }

/*        code for unequal increments or equal increments */
/*          not equal to 1 */

    ix = 1;
    iy = 1;
    if （*incx < 0） {
	ix = （-（*n） + 1） * *incx + 1;
    }
    if （*incy < 0） {
	iy = （-（*n） + 1） * *incy + 1;
    }
    i__1 = *n;
    for （i__ = 1; i__ <= i__1; ++i__） {
	i__2 = iy;
	i__3 = iy;
	i__4 = ix;
	q__2.r = ca->r * cx[i__4].r - ca->i * cx[i__4].i q__2.i = ca->r * cx[
		i__4].i + ca->i * cx[i__4].r;
	q__1.r = cy[i__3].r + q__2.r q__1.i = cy[i__3].i + q__2.i;
	cy[i__2].r = q__1.r cy[i__2].i = q__1.i;
	ix += *incx;
	iy += *incy;
/* L10: */
    }
    return 0;

/*        code for both increments equal to 1 */

L20:
    i__1 = *n;
    for （i__ = 1; i__ <= i__1; ++i__） {
	i__2 = i__;
	i__3 = i__;
	i__4 = i__;
	q__2.r = ca->r * cx[i__4].r - ca->i * cx[i__4].i q__2.i = ca->r * cx[
		i__4].i + ca->i * cx[i__4].r;
	q__1.r = cy[i__3].r + q__2.r q__1.i = cy[i__3].i + q__2.i;
	cy[i__2].r = q__1.r cy[i__2].i = q__1.i;
/* L30: */
    }
    return 0;
} /* caxpy_ */

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     947776  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\OutputLinux32-gcc4.4.1\cmvs

     文件      30706  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\OutputLinux32-gcc4.4.1\genOption

     文件    1305554  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\OutputLinux32-gcc4.4.1\pmvs2

     文件     316928  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\Win32-VS2010\cmvs.exe

     文件      28160  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\Win32-VS2010\genOption.exe

     文件     446464  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\Win32-VS2010\pmvs2.exe

     文件      30208  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\Win32-VS2010\pthread.dll

     文件       2174  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\Win32-VS2010\Readme.txt

     文件     382976  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\Win64-VS2010\cmvs.exe

     文件      33280  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\Win64-VS2010\genOption.exe

     文件     607568  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\Win64-VS2010\msvcp100.dll

     文件     827728  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\Win64-VS2010\msvcr100.dll

     文件     802624  2013-09-12 07:34  PMVS_Source\binariesWin-Linux\Win64-VS2010\msvcr80.dll

     文件      68369  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\bundle.rd.out

     文件     144947  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\kermit000.rd.jpg

     文件     155824  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\kermit001.rd.jpg

     文件     152110  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\kermit002.rd.jpg

     文件     162715  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\kermit003.rd.jpg

     文件     146028  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\kermit005.rd.jpg

     文件     150450  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\kermit007.rd.jpg

     文件     156278  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\kermit008.rd.jpg

     文件     144373  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\kermit009.rd.jpg

     文件     126423  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\kermit010.rd.jpg

     文件        270  2014-01-07 11:13  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\list.rd.txt

     文件        130  2014-01-07 10:54  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\pmvs_options.txt

     文件       1774  2014-01-07 11:12  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\prep_pmvs.sh

     文件        136  2014-01-07 10:54  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\txt\00000000.txt

     文件        139  2014-01-07 10:54  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\txt\00000001.txt

     文件        139  2014-01-07 10:54  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\txt\00000002.txt

     文件        139  2014-01-07 10:54  PMVS_Source\binariesWin-Linux\Win64-VS2010\pmvs\txt\00000003.txt

............此处省略6131个文件信息