/* zero.f -- translated by f2c (version 20000817).
   You must link the resulting object file with the libraries:
	-lf2c -lm   (in that order)
*/

#include "f2c.h"
#include "blaswrap.h"

/* Copyright (C) 2000, International Business Machines */
/* Corporation and others.  All Rights Reserved. */
/* Subroutine */ int rzrmat_(a, n, m)
doublereal *a;
integer *n, *m;
{
    /* System generated locals */
    integer a_dim1, a_offset, i__1, i__2;

    /* Local variables */
    static integer i__, j;

    /* Parameter adjustments */
    a_dim1 = *n;
    a_offset = 1 + a_dim1 * 1;
    a -= a_offset;

    /* Function Body */
    i__1 = *n;
    for (i__ = 1; i__ <= i__1; ++i__) {
	i__2 = *m;
	for (j = 1; j <= i__2; ++j) {
	    a[i__ + j * a_dim1] = 0.;
/* L20: */
	}
/* L10: */
    }
    return 0;
} /* rzrmat_ */

/* Subroutine */ int rzrvec_(b, n)
doublereal *b;
integer *n;
{
    /* System generated locals */
    integer i__1;

    /* Local variables */
    static integer i__;

    /* Parameter adjustments */
    --b;

    /* Function Body */
    i__1 = *n;
    for (i__ = 1; i__ <= i__1; ++i__) {
	b[i__] = 0.;
/* L10: */
    }
    return 0;
} /* rzrvec_ */

/* Subroutine */ int izrmat_(a, n, m)
integer *a, *n, *m;
{
    /* System generated locals */
    integer a_dim1, a_offset, i__1, i__2;

    /* Local variables */
    static integer i__, j;

    /* Parameter adjustments */
    a_dim1 = *n;
    a_offset = 1 + a_dim1 * 1;
    a -= a_offset;

    /* Function Body */
    i__1 = *n;
    for (i__ = 1; i__ <= i__1; ++i__) {
	i__2 = *m;
	for (j = 1; j <= i__2; ++j) {
	    a[i__ + j * a_dim1] = 0;
/* L20: */
	}
/* L10: */
    }
    return 0;
} /* izrmat_ */

/* Subroutine */ int izrvec_(b, n)
integer *b, *n;
{
    /* System generated locals */
    integer i__1;

    /* Local variables */
    static integer i__;

    /* Parameter adjustments */
    --b;

    /* Function Body */
    i__1 = *n;
    for (i__ = 1; i__ <= i__1; ++i__) {
	b[i__] = 0;
/* L10: */
    }
    return 0;
} /* izrvec_ */