Revision afb150384c671c72502d79b81811b1923ee6084e authored by Edzer J. Pebesma on 12 December 2008, 17:59:42 UTC, committed by cran-robot on 12 December 2008, 17:59:42 UTC
1 parent eb020f8
copy.c
/**************************************************************************
**
** Copyright (C) 1993 David E. Steward & Zbigniew Leyk, all rights reserved.
**
** Meschach Library
**
** This Meschach Library is provided "as is" without any express
** or implied warranty of any kind with respect to this software.
** In particular the authors shall not be liable for any direct,
** indirect, special, incidental or consequential damages arising
** in any way from use of the software.
**
** Everyone is granted permission to copy, modify and redistribute this
** Meschach Library, provided:
** 1. All copies contain this copyright notice.
** 2. All modified copies shall carry a notice stating who
** made the last modification and the date of such modification.
** 3. No charge is made for this software or works derived from it.
** This clause shall not be construed as constraining other software
** distributed on the same medium as this software, nor is a
** distribution fee considered a charge.
**
***************************************************************************/
static char rcsid[] = "$Id: copy.c,v 1.1.1.1 2003-06-23 18:31:34 cees Exp $";
#include <stdio.h>
#include "matrix.h"
/* _m_copy -- copies matrix into new area */
MAT *_m_copy(MAT *in,MAT *out,u_int i0,u_int j0)
{
u_int i /* ,j */;
if ( in==MNULL )
error(E_NULL,"_m_copy");
if ( in==out )
return (out);
if ( out==MNULL || out->m < in->m || out->n < in->n )
out = m_resize(out,in->m,in->n);
for ( i=i0; i < in->m; i++ )
MEM_COPY(&(in->me[i][j0]),&(out->me[i][j0]),
(in->n - j0)*sizeof(Real));
/* for ( j=j0; j < in->n; j++ )
out->me[i][j] = in->me[i][j]; */
return (out);
}
/* _v_copy -- copies vector into new area */
VEC *_v_copy(VEC *in,VEC *out,u_int i0)
{
/* u_int i,j; */
if ( in==VNULL )
error(E_NULL,"_v_copy");
if ( in==out )
return (out);
if ( out==VNULL || out->dim < in->dim )
out = v_resize(out,in->dim);
MEM_COPY(&(in->ve[i0]),&(out->ve[i0]),(in->dim - i0)*sizeof(Real));
/* for ( i=i0; i < in->dim; i++ )
out->ve[i] = in->ve[i]; */
return (out);
}
/* px_copy -- copies permutation 'in' to 'out' */
PERM *px_copy(PERM *in,PERM *out)
{
/* int i; */
if ( in == PNULL )
error(E_NULL,"px_copy");
if ( in == out )
return out;
if ( out == PNULL || out->size != in->size )
out = px_resize(out,in->size);
MEM_COPY(in->pe,out->pe,in->size*sizeof(u_int));
/* for ( i = 0; i < in->size; i++ )
out->pe[i] = in->pe[i]; */
return out;
}
/*
The .._move() routines are for moving blocks of memory around
within Meschach data structures and for re-arranging matrices,
vectors etc.
*/
/* m_move -- copies selected pieces of a matrix
-- moves the m0 x n0 submatrix with top-left cor-ordinates (i0,j0)
to the corresponding submatrix of out with top-left co-ordinates
(i1,j1)
-- out is resized (& created) if necessary */
MAT *m_move(MAT *in,int i0,int j0,int m0,int n0,MAT *out,int i1,int j1)
{
int i;
if ( ! in )
error(E_NULL,"m_move");
if ( i0 < 0 || j0 < 0 || i1 < 0 || j1 < 0 || m0 < 0 || n0 < 0 ||
i0+m0 > in->m || j0+n0 > in->n )
error(E_BOUNDS,"m_move");
if ( ! out )
out = m_resize(out,i1+m0,j1+n0);
else if ( i1+m0 > out->m || j1+n0 > out->n )
out = m_resize(out,max(out->m,i1+m0),max(out->n,j1+n0));
for ( i = 0; i < m0; i++ )
MEM_COPY(&(in->me[i0+i][j0]),&(out->me[i1+i][j1]),
n0*sizeof(Real));
return out;
}
/* v_move -- copies selected pieces of a vector
-- moves the length dim0 subvector with initial index i0
to the corresponding subvector of out with initial index i1
-- out is resized if necessary */
VEC *v_move(VEC *in,int i0,int dim0,VEC *out,int i1)
{
if ( ! in )
error(E_NULL,"v_move");
if ( i0 < 0 || dim0 < 0 || i1 < 0 ||
i0+dim0 > in->dim )
error(E_BOUNDS,"v_move");
if ( (! out) || i1+dim0 > out->dim )
out = v_resize(out,i1+dim0);
MEM_COPY(&(in->ve[i0]),&(out->ve[i1]),dim0*sizeof(Real));
return out;
}
/* mv_move -- copies selected piece of matrix to a vector
-- moves the m0 x n0 submatrix with top-left co-ordinate (i0,j0) to
the subvector with initial index i1 (and length m0*n0)
-- rows are copied contiguously
-- out is resized if necessary */
VEC *mv_move(MAT *in,int i0,int j0,int m0,int n0,VEC *out,int i1)
{
int dim1, i;
if ( ! in )
error(E_NULL,"mv_move");
if ( i0 < 0 || j0 < 0 || m0 < 0 || n0 < 0 || i1 < 0 ||
i0+m0 > in->m || j0+n0 > in->n )
error(E_BOUNDS,"mv_move");
dim1 = m0*n0;
if ( (! out) || i1+dim1 > out->dim )
out = v_resize(out,i1+dim1);
for ( i = 0; i < m0; i++ )
MEM_COPY(&(in->me[i0+i][j0]),&(out->ve[i1+i*n0]),n0*sizeof(Real));
return out;
}
/* vm_move -- copies selected piece of vector to a matrix
-- moves the subvector with initial index i0 and length m1*n1 to
the m1 x n1 submatrix with top-left co-ordinate (i1,j1)
-- copying is done by rows
-- out is resized if necessary */
MAT *vm_move(VEC *in,int i0,MAT *out,int i1,int j1,int m1,int n1)
{
int dim0, i;
if ( ! in )
error(E_NULL,"vm_move");
if ( i0 < 0 || i1 < 0 || j1 < 0 || m1 < 0 || n1 < 0 ||
i0+m1*n1 > in->dim )
error(E_BOUNDS,"vm_move");
if ( ! out )
out = m_resize(out,i1+m1,j1+n1);
else
out = m_resize(out,max(i1+m1,out->m),max(j1+n1,out->n));
dim0 = m1*n1;
for ( i = 0; i < m1; i++ )
MEM_COPY(&(in->ve[i0+i*n1]),&(out->me[i1+i][j1]),n1*sizeof(Real));
return out;
}
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...