copy.c

Go to the documentation of this file.

#include <../../nrnconf.h>

/**************************************************************************
**
** 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[] = "copy.c,v 1.1 1997/12/04 17:55:17 hines Exp";
#include    <stdio.h>
#include    "matrix.h"



/* _m_copy -- copies matrix into new area */
MAT *_m_copy(in,out,i0,j0)
MAT *in,*out;
u_int   i0,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(in,out,i0)
VEC *in,*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(in,out)
PERM    *in,*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(in,i0,j0,m0,n0,out,i1,j1)
MAT *in, *out;
int i0, j0, m0, n0, i1, 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(in,i0,dim0,out,i1)
VEC *in, *out;
int i0, dim0, 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(in,i0,j0,m0,n0,out,i1)
MAT *in;
VEC *out;
int i0, j0, m0, n0, 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(in,i0,out,i1,j1,m1,n1)
VEC *in;
MAT *out;
int i0, i1, j1, m1, 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;
}